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

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

Browse Source 
- Changed Triangle{3,6,7,10} return types in docstrings: NTuple → SVector{N, Vec{2,Float64}}
- Changed all implementations: tuple of tuples → SVector of Vec{2,Float64}
- Updated reference_coordinates for Tri3, Tri6, Tri7, Tri10 (4 functions)
This commit is contained in:
Jukka Aho
2025-11-21 00:30:52 +02:00
parent a1d97bb1ef
commit 0926a7cdbe
1 changed files with 34 additions and 34 deletions
Download Patch File Download Diff File Expand all files Collapse all files
src/topology/triangles.jl
+34 -34
View File
@@ -184,7 +184,7 @@ dim(::Triangle{N}) where {N} = 2
# ============================================================================
"""
reference_coordinates(::Triangle{3}) -> NTuple{3, NTuple{2, Float64}}
reference_coordinates(::Triangle{3}) -> SVector{3, Vec{2,Float64}}
Return reference coordinates for linear triangle (Tri3) - 3 corner nodes only.
@@ -210,15 +210,15 @@ Tuple of 3 coordinate pairs: ((ξ₁, η₁), (ξ₂, η₂), (ξ₃, η₃))
- Node 3: (0.0, 1.0) - Along η-axis
"""
function reference_coordinates(::Triangle{3})
return (
(0.0, 0.0), # N1: Corner at origin
(1.0, 0.0), # N2: Corner along ξ
(0.0, 1.0) # N3: Corner along η
return SVector(
Vec{2,Float64}((0.0, 0.0)), # N1: Corner at origin
Vec{2,Float64}((1.0, 0.0)), # N2: Corner along ξ
Vec{2,Float64}((0.0, 1.0)) # N3: Corner along η
)
end
"""
reference_coordinates(::Triangle{6}) -> NTuple{6, NTuple{2, Float64}}
reference_coordinates(::Triangle{6}) -> SVector{6, Vec{2,Float64}}
Return reference coordinates for quadratic triangle (Tri6) - 6 nodes total.
@@ -230,18 +230,18 @@ Return reference coordinates for quadratic triangle (Tri6) - 6 nodes total.
- Node 6: Edge midpoint between N3-N1 (0.0, 0.5)
"""
function reference_coordinates(::Triangle{6})
return (
(0.0, 0.0), # N1: Corner
(1.0, 0.0), # N2: Corner
(0.0, 1.0), # N3: Corner
(0.5, 0.0), # N4: Midpoint of edge 1-2
(0.5, 0.5), # N5: Midpoint of edge 2-3
(0.0, 0.5) # N6: Midpoint of edge 3-1
return SVector(
Vec{2,Float64}((0.0, 0.0)), # N1: Corner
Vec{2,Float64}((1.0, 0.0)), # N2: Corner
Vec{2,Float64}((0.0, 1.0)), # N3: Corner
Vec{2,Float64}((0.5, 0.0)), # N4: Midpoint of edge 1-2
Vec{2,Float64}((0.5, 0.5)), # N5: Midpoint of edge 2-3
Vec{2,Float64}((0.0, 0.5)) # N6: Midpoint of edge 3-1
)
end
"""
reference_coordinates(::Triangle{7}) -> NTuple{7, NTuple{2, Float64}}
reference_coordinates(::Triangle{7}) -> SVector{7, Vec{2,Float64}}
Return reference coordinates for quadratic triangle with centroid (Tri7).
@@ -252,19 +252,19 @@ Return reference coordinates for quadratic triangle with centroid (Tri7).
- Node 7: Face centroid (1/3, 1/3)
"""
function reference_coordinates(::Triangle{7})
return (
(0.0, 0.0), # N1: Corner
(1.0, 0.0), # N2: Corner
(0.0, 1.0), # N3: Corner
(0.5, 0.0), # N4: Midpoint edge 1-2
(0.5, 0.5), # N5: Midpoint edge 2-3
(0.0, 0.5), # N6: Midpoint edge 3-1
(1.0 / 3.0, 1.0 / 3.0) # N7: Face centroid
return SVector(
Vec{2,Float64}((0.0, 0.0)), # N1: Corner
Vec{2,Float64}((1.0, 0.0)), # N2: Corner
Vec{2,Float64}((0.0, 1.0)), # N3: Corner
Vec{2,Float64}((0.5, 0.0)), # N4: Midpoint edge 1-2
Vec{2,Float64}((0.5, 0.5)), # N5: Midpoint edge 2-3
Vec{2,Float64}((0.0, 0.5)), # N6: Midpoint edge 3-1
Vec{2,Float64}((1.0 / 3.0, 1.0 / 3.0)) # N7: Face centroid
)
end
"""
reference_coordinates(::Triangle{10}) -> NTuple{10, NTuple{2, Float64}}
reference_coordinates(::Triangle{10}) -> SVector{10, Vec{2,Float64}}
Return reference coordinates for cubic triangle (Tri10) - 10 nodes total.
@@ -278,17 +278,17 @@ Return reference coordinates for cubic triangle (Tri10) - 10 nodes total.
- Node 10: Face centroid (1/3, 1/3)
"""
function reference_coordinates(::Triangle{10})
return (
(0.0, 0.0), # N1: Corner
(1.0, 0.0), # N2: Corner
(0.0, 1.0), # N3: Corner
(1.0 / 3.0, 0.0), # N4: Edge 1-2, 1/3
(2.0 / 3.0, 0.0), # N5: Edge 1-2, 2/3
(2.0 / 3.0, 1.0 / 3.0), # N6: Edge 2-3, 1/3
(1.0 / 3.0, 2.0 / 3.0), # N7: Edge 2-3, 2/3
(0.0, 2.0 / 3.0), # N8: Edge 3-1, 1/3
(0.0, 1.0 / 3.0), # N9: Edge 3-1, 2/3
(1.0 / 3.0, 1.0 / 3.0) # N10: Face centroid
return SVector(
Vec{2,Float64}((0.0, 0.0)), # N1: Corner
Vec{2,Float64}((1.0, 0.0)), # N2: Corner
Vec{2,Float64}((0.0, 1.0)), # N3: Corner
Vec{2,Float64}((1.0 / 3.0, 0.0)), # N4: Edge 1-2, 1/3
Vec{2,Float64}((2.0 / 3.0, 0.0)), # N5: Edge 1-2, 2/3
Vec{2,Float64}((2.0 / 3.0, 1.0 / 3.0)), # N6: Edge 2-3, 1/3
Vec{2,Float64}((1.0 / 3.0, 2.0 / 3.0)), # N7: Edge 2-3, 2/3
Vec{2,Float64}((0.0, 2.0 / 3.0)), # N8: Edge 3-1, 1/3
Vec{2,Float64}((0.0, 1.0 / 3.0)), # N9: Edge 3-1, 2/3
Vec{2,Float64}((1.0 / 3.0, 1.0 / 3.0)) # N10: Face centroid
)
end