Periodic boundaries aren't tagged

movement/monge_ampere.py

@@ -329,19 +329,6 @@ def l2_projector(self):
             for dirichlet_bc in self._l2_projector_bcs(boundary_tag)
         ]
 
-        # Check for an attempt to fix periodic boundary segments
-        subdomains = set()
-        for bc in bcs:
-            if not hasattr(bc, "sub_domain"):
-                subdomain = bc._F.sub_domain  # EquationBC doesn't have sub_domain attr
-            else:
-                subdomain = bc.sub_domain
-            if not isinstance(subdomain, Iterable):
-                subdomain = [subdomain]
-            subdomains = subdomains.union(subdomain)
-        if set(self.fixed_boundary_segments).difference(subdomains):
-            raise ValueError("Cannot fix boundary nodes for periodic segments.")
-
         # Create solver
         problem = firedrake.LinearVariationalProblem(a, L, self._grad_phi, bcs=bcs)
         self._l2_projector = firedrake.LinearVariationalSolver(

test/test_monge_ampere.py

@@ -109,16 +109,6 @@ def test_invalid_plane_valuerror(self):
         msg = "Could not determine a plane for the provided points."
         self.assertEqual(str(cm.exception), msg)
 
-    def test_periodic_fixed_boundaries_valueerror(self):
-        mesh = self.mesh(n=3, periodic=True)
-        mover = MongeAmpereMover_Relaxation(
-            mesh, ring_monitor, fixed_boundary_segments=["on_boundary"]
-        )
-        with self.assertRaises(ValueError) as cm:
-            mover.move()
-        msg = "Cannot fix boundary nodes for periodic segments."
-        self.assertEqual(str(cm.exception), msg)
-
     def test_curved_notimplementederror(self):
         coords = Function(VectorFunctionSpace(UnitTriangleMesh(), "CG", 2))
         coords.interpolate(coords.function_space().mesh().coordinates)