diff --git a/dune/gfe/harmonicenergystiffness.hh b/dune/gfe/harmonicenergystiffness.hh
index 8a8ea8f7bf9b71895fa9182ad21bb92476bbaa89..7a85534b4c836c99a5bc076774bc5824fc6c6027 100644
--- a/dune/gfe/harmonicenergystiffness.hh
+++ b/dune/gfe/harmonicenergystiffness.hh
@@ -55,7 +55,7 @@ energy(const Entity& element,
     
     LocalGeodesicFEFunction<gridDim, double, TargetSpace> localGeodesicFEFunction(localSolution);
 
-    int quadOrder = 1;//gridDim;
+    int quadOrder = gridDim;
 
     const Dune::QuadratureRule<double, gridDim>& quad 
         = Dune::QuadratureRules<double, gridDim>::rule(element.type(), quadOrder);
@@ -121,7 +121,7 @@ assembleEmbeddedGradient(const Entity& element,
     LocalGeodesicFEFunction<gridDim, double, TargetSpace> localGeodesicFEFunction(localSolution);
 
     // I am not sure about the correct quadrature order
-    int quadOrder = 1;//gridDim;
+    int quadOrder = gridDim;
 
     // numerical quadrature loop
     const Dune::QuadratureRule<double, gridDim>& quad