Non-glued boundary spaces: gluing functions on the same volumetric patch

rafavzqz · Sep 5, 2024 · 8918059 · 8918059
1 parent 6769fce
commit 8918059
Showing 1 changed file with 13 additions and 1 deletion.
diff --git a/geopdes/inst/multipatch/@sp_multipatch/sp_multipatch.m b/geopdes/inst/multipatch/@sp_multipatch/sp_multipatch.m
@@ -141,7 +141,19 @@
       side_number  = msh.boundary.side_numbers(iptc);
       sp_bnd{iptc} = spaces{patch_number}.boundary(side_number);
     end
-    sp.boundary = sp_multipatch (sp_bnd, msh.boundary, boundary_interfaces);
+    if (gluing)
+      sp.boundary = sp_multipatch (sp_bnd, msh.boundary, boundary_interfaces);
+    else
+      glued_bnd_interfaces = true (numel(boundary_interfaces), 1);
+      vol_patch_numbers = msh.boundary.patch_numbers;
+      bnd_patches = [[boundary_interfaces.patch1]; [boundary_interfaces.patch2]];
+      for ii = 1:numel(boundary_interfaces)
+        if (vol_patch_numbers(bnd_patches(1,ii)) ~= vol_patch_numbers(bnd_patches(2,ii)))
+          glued_bnd_interfaces(ii) = false;
+        end
+      end
+      sp.boundary = sp_multipatch (sp_bnd, msh.boundary, boundary_interfaces(glued_bnd_interfaces));
+    end
 
     dofs = zeros (sp.boundary.ndof, 1);
     boundary_orientation = []; %boundary_orient = zeros (sp.boundary.ndof, 1);