From f6bdf1eb36a9c769c0ebe28dc15956473eec965f Mon Sep 17 00:00:00 2001 From: Lawrence Mitchell Date: Tue, 4 May 2021 09:53:42 +0100 Subject: [PATCH] fem: Fix CoordinateMapping changes for interpolation No integral_type available in that case, so check that integration dimension matches the cell dimension. --- tsfc/fem.py | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) diff --git a/tsfc/fem.py b/tsfc/fem.py index 7dd79f6d..63bcb911 100644 --- 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):