diff --git a/dune/gfe/harmonicenergystiffness.hh b/dune/gfe/harmonicenergystiffness.hh
index bb1e2abdd050f0c855da25b3f68b4deec14c1927..79d0da07f36b4f68277bbb5300add4f813146914 100644
--- a/dune/gfe/harmonicenergystiffness.hh
+++ b/dune/gfe/harmonicenergystiffness.hh
@@ -63,7 +63,7 @@ energy(const Entity& element,
                                                                                                       localSolution);
 
     int quadOrder = (element.type().isSimplex()) ? (localFiniteElement.localBasis().order()-1) * 2
-                                                 : (localFiniteElement.localBasis().order()-1) * 2 * gridDim;
+                                                 : localFiniteElement.localBasis().order() * 2 * gridDim;