Topology vroom 1

docs/src/reference/grid.md

@@ -46,6 +46,8 @@ get_cell_coordinates!
 Ferrite.ExclusiveTopology
 Ferrite.getneighborhood
 Ferrite.faceskeleton
+Ferrite.vertex_star_stencils
+Ferrite.getstencil
 ```
 
 ### Grid Sets Utility

src/Dofs/DofHandler.jl

@@ -620,6 +620,17 @@ function sortedge(edge::Tuple{Int,Int})
     a < b ? (return (edge, PathOrientationInfo(true))) : (return ((b, a), PathOrientationInfo(false)))
 end
 
+"""
+sortedge_fast(edge::Tuple{Int,Int})
+
+Returns the unique representation of an edge.
+Here the unique representation is the sorted node index tuple.
+"""
+function sortedge_fast(edge::Tuple{Int,Int})
+    a, b = edge
+    a < b ? (return edge) : (return (b, a))
+end
+
 """
     sortface(face::Tuple{Int})
     sortface(face::Tuple{Int,Int})
@@ -633,6 +644,20 @@ so the unique representation is always a tuple length 3.
 """
 sortface(face::Tuple{Int,Int}) = sortedge(face) # Face in 2D is the same as edge in 3D.
 
+
+"""
+    sortface_fast(face::Tuple{Int})
+    sortface_fast(face::Tuple{Int,Int})
+    sortface_fast(face::Tuple{Int,Int,Int})
+    sortface_fast(face::Tuple{Int,Int,Int,Int})
+
+Returns the unique representation of a face.
+Here the unique representation is the sorted node index tuple.
+Note that in 3D we only need indices to uniquely identify a face,
+so the unique representation is always a tuple length 3.
+"""
+sortface_fast(face::Tuple{Int,Int}) = sortedge_fast(face) # Face in 2D is the same as edge in 3D.
+
 """
     !!!NOTE TODO implement me.
 
@@ -669,6 +694,16 @@ function sortface(face::Tuple{Int,Int,Int})
     return (a, b, c), SurfaceOrientationInfo() # TODO fill struct
 end
 
+
+function sortface_fast(face::Tuple{Int,Int,Int})
+    a, b, c = face
+    b, c = minmax(b, c)
+    a, c = minmax(a, c)
+    a, b = minmax(a, b)
+    return (a, b, c)
+end
+
+
 function sortface(face::Tuple{Int,Int,Int,Int})
     a, b, c, d = face
     c, d = minmax(c, d)
@@ -680,7 +715,21 @@ function sortface(face::Tuple{Int,Int,Int,Int})
     return (a, b, c), SurfaceOrientationInfo() # TODO fill struct
 end
 
+
+function sortface_fast(face::Tuple{Int,Int,Int,Int})
+    a, b, c, d = face
+    c, d = minmax(c, d)
+    b, d = minmax(b, d)
+    a, d = minmax(a, d)
+    b, c = minmax(b, c)
+    a, c = minmax(a, c)
+    a, b = minmax(a, b)
+    return (a, b, c)
+end
+
+
 sortface(face::Tuple{Int}) = face, nothing
+sortface_fast(face::Tuple{Int}) = face
 
 """
     find_field(dh::DofHandler, field_name::Symbol)::NTuple{2,Int}

src/Ferrite.jl

@@ -89,6 +89,8 @@ include("FEValues/face_integrals.jl")
 
 # Grid
 include("Grid/grid.jl")
+include("Grid/topology.jl")
+include("Grid/utils.jl")
 include("Grid/grid_generators.jl")
 include("Grid/coloring.jl")