diff --git a/python/dune/perftool/pdelab/geometry.py b/python/dune/perftool/pdelab/geometry.py index 5b2049b6cafb821a85bf99d1dc5fa68d1e9f1290..9d8a1c28d1bfbd50af29e9c3f6c850b60a949de1 100644 --- a/python/dune/perftool/pdelab/geometry.py +++ b/python/dune/perftool/pdelab/geometry.py @@ -294,6 +294,10 @@ def define_constant_jacobian_inveser_transposed(name, restriction): pos = name_localcenter() dim = name_dimension() + if restriction: + geo_in = name_in_cell_geometry(restriction) + pos = "{}.global({})".format(geo_in, pos) + globalarg(name, dtype=np.float64, shape=(dim, dim), managed=False) return 'auto {} = {}.jacobianInverseTransposed({});'.format(name,