Glassvision/JuliaFEM.jl
1
0
Fork 0
mirror of https://github.com/JuliaFEM/JuliaFEM.jl.git synced 2026-09-23 11:02:36 +00:00
Code Issues Packages Projects Releases Wiki Activity

refactor(topology): Update Tetrahedron reference_coordinates to SVector of Vec

Browse Source 
- Changed Tet4 and Tet10 return types in docstrings: NTuple → SVector{N, Vec{3,Float64}}
- Changed implementations: tuple of tuples → SVector of Vec{3,Float64}
- Updated reference_coordinates for Tet4 (4 nodes) and Tet10 (10 nodes)
This commit is contained in:
Jukka Aho
2025-11-21 00:31:28 +02:00
parent fa02393718
commit 342990ce55
1 changed files with 18 additions and 18 deletions
Download Patch File Download Diff File Expand all files Collapse all files
src/topology/tetrahedra.jl
+18 -18
View File
@@ -161,7 +161,7 @@ dim(::Tetrahedron{N}) where {N} = 3
# ============================================================================
"""
reference_coordinates(::Tetrahedron{4}) -> NTuple{4, NTuple{3, Float64}}
reference_coordinates(::Tetrahedron{4}) -> SVector{4, Vec{3,Float64}}
Return reference coordinates for linear tetrahedron (Tet4) - 4 corner nodes only.
@@ -188,16 +188,16 @@ Tuple of 4 coordinate triples: ((ξ₁, η₁, ζ₁), (ξ₂, η₂, ζ₂), (
- Node 4: (0.0, 0.0, 1.0) - Along ζ-axis
"""
function reference_coordinates(::Tetrahedron{4})
return (
(0.0, 0.0, 0.0), # N1: Corner at origin
(1.0, 0.0, 0.0), # N2: Corner along ξ
(0.0, 1.0, 0.0), # N3: Corner along η
(0.0, 0.0, 1.0) # N4: Corner along ζ
return SVector(
Vec{3,Float64}((0.0, 0.0, 0.0)), # N1: Corner at origin
Vec{3,Float64}((1.0, 0.0, 0.0)), # N2: Corner along ξ
Vec{3,Float64}((0.0, 1.0, 0.0)), # N3: Corner along η
Vec{3,Float64}((0.0, 0.0, 1.0)) # N4: Corner along ζ
)
end
"""
reference_coordinates(::Tetrahedron{10}) -> NTuple{10, NTuple{3, Float64}}
reference_coordinates(::Tetrahedron{10}) -> SVector{10, Vec{3,Float64}}
Return reference coordinates for quadratic tetrahedron (Tet10) - 10 nodes total.
@@ -211,17 +211,17 @@ Return reference coordinates for quadratic tetrahedron (Tet10) - 10 nodes total.
- Node 10: Edge midpoint between N3-N4 (0.0, 0.5, 0.5)
"""
function reference_coordinates(::Tetrahedron{10})
return (
(0.0, 0.0, 0.0), # N1: Corner
(1.0, 0.0, 0.0), # N2: Corner
(0.0, 1.0, 0.0), # N3: Corner
(0.0, 0.0, 1.0), # N4: Corner
(0.5, 0.0, 0.0), # N5: Midpoint edge 1-2
(0.5, 0.5, 0.0), # N6: Midpoint edge 2-3
(0.0, 0.5, 0.0), # N7: Midpoint edge 3-1
(0.0, 0.0, 0.5), # N8: Midpoint edge 1-4
(0.5, 0.0, 0.5), # N9: Midpoint edge 2-4
(0.0, 0.5, 0.5) # N10: Midpoint edge 3-4
return SVector(
Vec{3,Float64}((0.0, 0.0, 0.0)), # N1: Corner
Vec{3,Float64}((1.0, 0.0, 0.0)), # N2: Corner
Vec{3,Float64}((0.0, 1.0, 0.0)), # N3: Corner
Vec{3,Float64}((0.0, 0.0, 1.0)), # N4: Corner
Vec{3,Float64}((0.5, 0.0, 0.0)), # N5: Midpoint edge 1-2
Vec{3,Float64}((0.5, 0.5, 0.0)), # N6: Midpoint edge 2-3
Vec{3,Float64}((0.0, 0.5, 0.0)), # N7: Midpoint edge 3-1
Vec{3,Float64}((0.0, 0.0, 0.5)), # N8: Midpoint edge 1-4
Vec{3,Float64}((0.5, 0.0, 0.5)), # N9: Midpoint edge 2-4
Vec{3,Float64}((0.0, 0.5, 0.5)) # N10: Midpoint edge 3-4
)
end