From 147bfc92c6a1337400785f88e57e28e23b62e9bf Mon Sep 17 00:00:00 2001
From: DennisGlaeser <dennis.glaeser@iws.uni-stuttgart.de>
Date: Wed, 29 Nov 2017 15:39:05 +0100
Subject: [PATCH] [evalGradients] fix remaining errors
---
dumux/discretization/evalgradients.hh | 8 ++++----
1 file changed, 4 insertions(+), 4 deletions(-)
diff --git a/dumux/discretization/evalgradients.hh b/dumux/discretization/evalgradients.hh
index 0384c97c0b..ebce370cde 100644
--- a/dumux/discretization/evalgradients.hh
+++ b/dumux/discretization/evalgradients.hh
@@ -61,15 +61,15 @@ evalGradients(const Element& element,
// evaluate gradients using the local finite element basis
const auto& localBasis = fvGridGeometry.feCache().get(geometry.type()).localBasis();
- // the inverse transposed of the jacobian matrix
- const auto jacInvT = g.jacobianInverseTransposed(local);
-
// evaluate the shape function gradients at the scv center
- using ShapeJacobian = typename decltype(localBasis)::Traits::JacobianType;
+ using ShapeJacobian = typename std::decay_type< decltype(localBasis) >::Traits::JacobianType;
const auto localPos = geometry.local(globalPos);
std::vector< ShapeJacobian > shapeJacobian;
localBasis.evaluateJacobian(localPos, shapeJacobian);
+ // the inverse transposed of the jacobian matrix
+ const auto jacInvT = geometry.jacobianInverseTransposed(localPos);
+
// interpolate the gradients
Dune::FieldVector<GlobalPosition, PrimaryVariables::dimension> result( PrimaryVariables(0.0) );
for (int i = 0; i < element.subEntities(Element::Geometry::mydimension); ++i)
--
GitLab