diff --git a/src/harmonicmaps.cc b/src/harmonicmaps.cc
index 4b1c1159d2af3e4a811d6efaccbd8b8575ace0b5..6a2335f790d64449bccf9a71a1f01fc09317ffc6 100644
--- a/src/harmonicmaps.cc
+++ b/src/harmonicmaps.cc
@@ -252,7 +252,7 @@ int main (int argc, char *argv[]) try
       GFE::EmbeddedGlobalGFEFunction<FEBasis, TargetSpace> numericalSolution(feBasis, x);
 
       // QuadratureRule for the integral of the L^2 error
-      QuadratureRuleKey quadKey(dim,3);
+      QuadratureRuleKey quadKey(dim,6);
 
       // Compute the embedded L^2 error
       double l2Error = DiscretizationError<GridType::LeafGridView>::computeL2Error(&numericalSolution,