fem: Fix CoordinateMapping changes for interpolation

No integral_type available in that case, so check that integration dimension matches the cell dimension.
firedrakeproject · May 4, 2021 · f6bdf1e · f6bdf1e
1 parent 94aa4d5
commit f6bdf1e
Showing 1 changed file with 3 additions and 2 deletions.
diff --git a/tsfc/fem.py b/tsfc/fem.py
@@ -136,9 +136,10 @@ def preprocess(self, expr, context):
         :arg context: The translation context.
         :returns: A new UFL expression
         """
-        ifacet = self.interface.integral_type.startswith("interior_facet")
+        is_facet = (self.interface.integration_dim !=
+                    self.interface.fiat_cell.get_dimension())
         return preprocess_expression(expr, complex_mode=context.complex_mode,
-                                     do_apply_restrictions=ifacet)
+                                     do_apply_restrictions=is_facet)
 
     @property
     def config(self):