pyop3: RelabeledMesh

firedrake/mesh.py

@@ -5013,7 +5013,7 @@ def RelabeledMesh(mesh, indicator_functions, subdomain_ids, **kwargs):
         # Clear label stratum; this is a copy, so safe to change.
         plex1.clearLabelStratum(dmlabel_name, subid)
         dmlabel = plex1.getLabel(dmlabel_name)
-        section = f.topological.function_space().dm.getSection()
+        section = f.topological.function_space().local_section
         dmcommon.mark_points_with_function_array(plex, section, height, f.dat.data_ro_with_halos.real.astype(IntType), dmlabel, subid)
     distribution_parameters_noop = {"partition": False,
                                     "overlap_type": (DistributedMeshOverlapType.NONE, 0)}

tests/regression/test_mark_entities.py

@@ -1,8 +1,8 @@
 import pytest
 from firedrake import *
-pytest.skip(allow_module_level=True, reason="pyop3 TODO")
 
 
+@pytest.mark.skip(reason="pyop3 TODO: Hitting minor issues in layout tabulation")
 @pytest.mark.parallel(nprocs=8)
 def test_mark_entities_mark_points_with_function_array():
     # Mark cells in f0.
@@ -48,6 +48,7 @@ def test_mark_entities_mark_points_with_function_array():
     assert abs(v - (4 * .5 + 4 * .5 * sqrt(2))) < 1.e-10
 
 
+@pytest.mark.skip(reason="pyop3 TODO: Hitting minor issues in layout tabulation")
 @pytest.mark.parallel(nprocs=7)
 def test_mark_entities_overlapping_facet_subdomains():
     my_facet_label = 777