Fix to make sure adjacency graph is limited to faces

sandialabs · Sep 16, 2024 · a70776b · a70776b
1 parent 1671ba2
commit a70776b
Show file tree

Hide file tree

Showing 2 changed files with 15 additions and 3 deletions.
diff --git a/packages/seacas/libraries/ioss/src/text_mesh/Iotm_TextMeshAdjacencyGraph.h b/packages/seacas/libraries/ioss/src/text_mesh/Iotm_TextMeshAdjacencyGraph.h
@@ -364,7 +364,7 @@ namespace Iotm {
         std::vector<EntityId> sideNodes =
             get_sorted_side_nodes(adjacency.elementIndex, adjacency.side);
 
-        for (int otherSide = 1; otherSide <= get_element_topology(neighborElementIndex).num_sides();
+        for (int otherSide = 1; otherSide <= get_element_topology(neighborElementIndex).num_face_sides();
              ++otherSide) {
           std::vector<EntityId> otherSideNodes =
               get_sorted_side_nodes(neighborElementIndex, otherSide);
@@ -738,7 +738,7 @@ namespace Iotm {
         initialize_side_connectivity_graph(elementIndices);
 
         for (size_t elementIndex : elementIndices) {
-          int numSides = get_element_topology(elementIndex).num_sides();
+          int numSides = get_element_topology(elementIndex).num_face_sides();
           for (int side = 1; side <= numSides; ++side) {
             if (m_indexGraph[elementIndex].sideReference[side - 1] == 0) {
               CurrentAdjacency adjacency(elementIndex, side);
@@ -752,7 +752,7 @@ namespace Iotm {
       {
         for (size_t elementIndex : elementIndices) {
           m_indexGraph[elementIndex] =
-              FaceConnections(get_element_topology(elementIndex).num_sides());
+              FaceConnections(get_element_topology(elementIndex).num_face_sides());
         }
       }
 

diff --git a/packages/seacas/libraries/ioss/src/text_mesh/Iotm_TextMeshTopologyMapping.h b/packages/seacas/libraries/ioss/src/text_mesh/Iotm_TextMeshTopologyMapping.h
@@ -80,6 +80,18 @@ namespace Iotm {
 
     bool operator!=(const TopologyMapEntry &rhs) const { return !(*this == rhs); }
 
+    int num_face_sides() const
+    {
+      if (topology->is_shell()) {
+        // Only interested in face boundaries, not edges
+        if (topology->parametric_dimension() == 2) {
+          return topology->number_faces();
+        }
+      }
+
+      return topology->number_boundaries();
+    }
+
     int num_sides() const
     {
       return sideTopologies.size();