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