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