equil init: reduce the tolerance of the regula falsi solvers to 1e-10

This commit is contained in:
Andreas Lauser 2018-01-02 12:43:56 +01:00
parent 2f00e7804a
commit 794fc3f588

View File

@ -675,7 +675,7 @@ double satFromPc(const MaterialLawManager& materialLawManager,
assert(f0 > 0 && f1 < 0);
const int maxIter = 60;
const double tol = 1e-6;
const double tol = 1e-10;
// regula falsi with the "Pegasus" method to avoid stagnation
for (int iter = 0; iter < maxIter; ++ iter) {
@ -784,7 +784,7 @@ double satFromSumOfPcs(const MaterialLawManager& materialLawManager,
assert(f0 > 0.0 && f1 < 0.0);
const int maxIter = 60;
const double tol = 1e-6;
const double tol = 1e-10;
// regula falsi with the "Pegasus" method to avoid stagnation
for (int iter = 0; iter < maxIter; ++ iter) {