From 72706c2a8182d46e77146502964f4ed0385ed02e Mon Sep 17 00:00:00 2001
From: Oliver Sander <oliver.sander@tu-dresden.de>
Date: Wed, 9 May 2018 06:00:49 +0200
Subject: [PATCH] Don't hardwire quad weights to 'double'
---
dune/gfe/harmonicenergystiffness.hh | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/dune/gfe/harmonicenergystiffness.hh b/dune/gfe/harmonicenergystiffness.hh
index 645b519a..fa3cce66 100644
--- a/dune/gfe/harmonicenergystiffness.hh
+++ b/dune/gfe/harmonicenergystiffness.hh
@@ -49,11 +49,11 @@ energy(const typename Basis::LocalView& localView,
// Local position of the quadrature point
const Dune::FieldVector<double,gridDim>& quadPos = quad[pt].position();
- const double integrationElement = element.geometry().integrationElement(quadPos);
+ const auto integrationElement = element.geometry().integrationElement(quadPos);
const auto jacobianInverseTransposed = element.geometry().jacobianInverseTransposed(quadPos);
- double weight = quad[pt].weight() * integrationElement;
+ auto weight = quad[pt].weight() * integrationElement;
// The derivative of the local function defined on the reference element
auto referenceDerivative = localInterpolationRule.evaluateDerivative(quadPos);
--
GitLab