From 5a1a0bd24d56b6630d47e69a1b435bd11aa1e7f7 Mon Sep 17 00:00:00 2001
From: Dominic Kempf <dominic.kempf@iwr.uni-heidelberg.de>
Date: Tue, 29 Jan 2019 13:55:30 +0100
Subject: [PATCH] [fixup] Fix equidistant mixin after previous change to
axiparallel one
---
python/dune/codegen/pdelab/geometry.py | 7 ++++++-
1 file changed, 6 insertions(+), 1 deletion(-)
diff --git a/python/dune/codegen/pdelab/geometry.py b/python/dune/codegen/pdelab/geometry.py
index 9f943595..d37a315f 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
--
GitLab