diff --git a/dune/gfe/localgeodesicfefunction.hh b/dune/gfe/localgeodesicfefunction.hh
index a8519b7892ef7955ccfcca6ff0cf9be6afa2c69d..a6178e45f08560abc5742987c327d7c01e635c59 100644
--- a/dune/gfe/localgeodesicfefunction.hh
+++ b/dune/gfe/localgeodesicfefunction.hh
@@ -191,11 +191,11 @@ evaluate(const Dune::FieldVector<ctype, dim>& local) const
 
     solver.setup(&assembler,
                  coefficients_[0],   // initial iterate
-                 1e-8,    // tolerance
-                 20,      // maxTrustRegionSteps
+                 1e-14,    // tolerance
+                 100,      // maxTrustRegionSteps
                  2,       // initial trust region radius
-                 20,      // inner iterations
-                 1e-8     // inner tolerance
+                 100,      // inner iterations
+                 1e-14     // inner tolerance
                  );
 
     solver.solve();