diff --git a/dune/gfe/harmonicenergystiffness.hh b/dune/gfe/harmonicenergystiffness.hh
index b39d882f2b8bc6d4382185039925a83891061b6d..645b519a5ad3ce17fc4cea5648d9aff54627cfc8 100644
--- a/dune/gfe/harmonicenergystiffness.hh
+++ b/dune/gfe/harmonicenergystiffness.hh
@@ -38,7 +38,7 @@ energy(const typename Basis::LocalView& localView,
     LocalInterpolationRule localInterpolationRule(localFiniteElement,localSolution);
 
     int quadOrder = (localFiniteElement.type().isSimplex()) ? (localFiniteElement.localBasis().order()-1) * 2
-                                                 : localFiniteElement.localBasis().order() * 2 * gridDim;
+                                                 : (localFiniteElement.localBasis().order() * gridDim - 1) * 2;
 
     const auto element = localView.element();