DirichletBC: support string subdomain_id

Author: pbrubeck

firedrake/functionspacedata.py

@@ -478,6 +478,8 @@ def __str__(self):
 
     @PETSc.Log.EventDecorator()
     def boundary_nodes(self, V, sub_domain):
+        region_names = V.mesh().topology_dm.getAttr("face_region_names") or {}
+        sub_domain = region_names.get(sub_domain, sub_domain)
         if sub_domain in ["bottom", "top"]:
             if not V.extruded:
                 raise ValueError("Invalid subdomain '%s' for non-extruded mesh",

firedrake/mg/mesh.py

@@ -129,6 +129,7 @@ def MeshHierarchy(mesh, refinement_levels,
     dm_cell_type, = mesh.dm_cell_types
     tdim = mesh.topology_dm.getDimension()
     cdm = dmcommon.submesh_create(mesh.topology_dm, tdim, "celltype", dm_cell_type, True)
+    cdm.setAttr("face_region_names", mesh.topology_dm.getAttr("face_region_names"))
     cdm.removeLabel("pyop2_core")
     cdm.removeLabel("pyop2_owned")
     cdm.removeLabel("pyop2_ghost")
@@ -142,6 +143,7 @@ def MeshHierarchy(mesh, refinement_levels,
         if i % refinements_per_level == 0:
             before(cdm, i)
         rdm = cdm.refine()
+        rdm.setAttr("face_region_names", cdm.getAttr("face_region_names"))
         if i % refinements_per_level == 0:
             after(rdm, i)
         dms.append(rdm)

firedrake/mg/netgen.py

@@ -155,6 +155,7 @@ def uniformRefinementRoutine(ngmesh, cdm):
     logger.info(f"\t\t\t[{time.time()}]Refining the plex")
     cdm.setRefinementUniform(True)
     rdm = cdm.refine()
+    rdm.setAttr("face_region_names", cdm.getAttr("face_region_names"))
     rdm.removeLabel("pyop2_core")
     rdm.removeLabel("pyop2_owned")
     rdm.removeLabel("pyop2_ghost")
@@ -267,6 +268,7 @@ def NetgenHierarchy(mesh, levs, flags, distribution_parameters=None):
     # Make a plex (cdm) without overlap.
     dm_cell_type, = mesh.dm_cell_types
     cdm = dmcommon.submesh_create(mesh.topology_dm, tdim, "celltype", dm_cell_type, True)
+    cdm.setAttr("face_region_names", mesh.topology_dm.getAttr("face_region_names"))
     cdm.removeLabel("pyop2_core")
     cdm.removeLabel("pyop2_owned")
     cdm.removeLabel("pyop2_ghost")