diff --git a/dune/gfe/harmonicenergystiffness.hh b/dune/gfe/harmonicenergystiffness.hh
index 645b519a5ad3ce17fc4cea5648d9aff54627cfc8..fa3cce665a0564d7e6d26d1d4159c3355e9139a5 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);