diff --git a/src/localgeodesicfefunction.hh b/src/localgeodesicfefunction.hh
index 994de48d6e6c017af05328390be9d0800010396f..ea1672991eb05bb1ac1c234e03e6566cf15692c4 100644
--- a/src/localgeodesicfefunction.hh
+++ b/src/localgeodesicfefunction.hh
@@ -68,7 +68,7 @@ evaluate(const Dune::FieldVector<ctype, dim>& local)
 
     solver.setup(&assembler,
                  coefficients_[0],   // initial iterate
-                 1e-5,    // tolerance
+                 1e-8,    // tolerance
                  20,      // maxTrustRegionSteps
                  2,       // initial trust region radius
                  20,      // inner iterations