diff --git a/test/localgeodesicfefunctiontest.cc b/test/localgeodesicfefunctiontest.cc
index b832101e136c01545c57e9ca5d581ba005707804..1016ae59861065073271b52015dcbbfbc7ca9f71 100644
--- a/test/localgeodesicfefunctiontest.cc
+++ b/test/localgeodesicfefunctiontest.cc
@@ -139,7 +139,7 @@ void testDerivativeOfValueWRTCoefficients(const std::vector<TargetSpace>& corner
     LocalGeodesicFEFunction<domainDim,double,TargetSpace> f(corners);
 
     // A quadrature rule as a set of test points
-    int quadOrder = 1;
+    int quadOrder = 3;
     
     const Dune::QuadratureRule<double, domainDim>& quad 
         = Dune::QuadratureRules<double, domainDim>::rule(GeometryType(GeometryType::simplex,domainDim), quadOrder);