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