From 75bc27d4881e9c76cd7dfb48d0ba1bab752baf86 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Wed, 1 Jun 2011 10:08:18 +0000
Subject: [PATCH] Solve the minimization problems for the evaluation of gfe
functions really down to machine precision. This appears to make the overall
algorithm much more stable
[[Imported from SVN: r7355]]
---
dune/gfe/localgeodesicfefunction.hh | 8 ++++----
1 file changed, 4 insertions(+), 4 deletions(-)
diff --git a/dune/gfe/localgeodesicfefunction.hh b/dune/gfe/localgeodesicfefunction.hh
index a8519b78..a6178e45 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();
--
GitLab