Fix for mixed topology when mixing cells with different facet types i…

…n 3D (#3572)
FEniCS · Dec 20, 2024 · 4d70f2f · 4d70f2f
1 parent 9ee46f0
commit 4d70f2f
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/cpp/dolfinx/fem/dofmapbuilder.cpp b/cpp/dolfinx/fem/dofmapbuilder.cpp
@@ -279,6 +279,10 @@ build_basic_dofmaps(
 
         const std::vector<std::vector<int>>& e_dofs_d = entity_dofs[d];
 
+        // Skip over undefined topology, e.g. quad facets of tetrahedra
+        if (d < D and !topology.connectivity({D, i}, {d, et}))
+          continue;
+
         // Iterate over each entity of current dimension d and type et
         std::span<const std::int32_t> c_to_e
             = d < D ? topology.connectivity({D, i}, {d, et})->links(c)