fix for subgrid in case there is empty boundary data in the parent gr…

…id (or probably also when the subgrid lies completely in the interior)
WIAS-PDELib · Jan 17, 2024 · d49f4c1 · d49f4c1
1 parent fd3587a
commit d49f4c1
Showing 1 changed file with 7 additions and 2 deletions.
diff --git a/src/subgrid.jl b/src/subgrid.jl
@@ -194,8 +194,13 @@ function subgrid(parent,
 
         subgrid[BFaceRegions]=sub_bfaceregions
         subgrid[BFaceGeometries]=sub_bfacetypes
-        subgrid[BFaceNodes]=tryfix(sub_bfacenodes)
-        subgrid[NumBFaceRegions]=maximum(sub_bfaceregions)
+        if length(sub_bfaceregions) > 1
+            subgrid[BFaceNodes]=tryfix(sub_bfacenodes)
+            subgrid[NumBFaceRegions]=maximum(sub_bfaceregions)
+        else
+            subgrid[BFaceNodes]=zeros(Ti, 2, 0)
+            subgrid[NumBFaceRegions]=0
+        end
     end
     subgrid[CoordinateSystem]=parent[CoordinateSystem]