diff --git a/python/dune/codegen/pdelab/geometry.py b/python/dune/codegen/pdelab/geometry.py
index 9f943595a065f9296ecd48447093bee3a8bdb6f5..d37a315f2926ffe2f914589925eb77eee9352d5e 100644
--- a/python/dune/codegen/pdelab/geometry.py
+++ b/python/dune/codegen/pdelab/geometry.py
@@ -272,8 +272,13 @@ class EquidistantGeometryMixin(AxiparallelGeometryMixin):
         rft = lop_template_range_field()
         return "{} = {}.gridView().template begin<0>()->geometry().integrationElement(Dune::FieldVector<{}, {}>());".format(name, gfs, rft, world_dimension())
 
-    @class_member(classtag="operator")
     def define_jacobian_inverse_transposed(self, name, restriction):
+        dim = world_dimension()
+        globalarg(name, shape=(dim, dim), managed=False)
+        self._define_jacobian_inverse_transposed(name, restriction)
+
+    @class_member(classtag="operator")
+    def _define_jacobian_inverse_transposed(self, name, restriction):
         dim = world_dimension()
         self._define_jacobian_inverse_transposed_eval(name)
         from dune.codegen.pdelab.localoperator import lop_template_ansatz_gfs