BUG: Corrected initialization of tolerances

jheylmun · jheylmun · commit c039b6071d59 · 2025-09-30T09:04:42.000-05:00
diff --git a/src/numerics/rootSolvers/rootSolver/rootSolver.C b/src/numerics/rootSolvers/rootSolver/rootSolver.C
@@ -42,7 +42,7 @@ void Foam::rootSolver::initialise(const scalarList& x) const
     tolerances_ = absTolerances_;
     forAll(x, i)
     {
-        tolerances_[i] = max(relTolerances_[i]*mag(x[i]), tolerances_[i]);
+        tolerances_[i] = max(relTolerances_[i]*mag(x[i]), absTolerances_[i]);
     }
 }
 

Original file line number	Diff line number	Diff line change
`@@ -42,7 +42,7 @@ void Foam::rootSolver::initialise(const scalarList& x) const`
`42`	`42`	`tolerances_ = absTolerances_;`
`43`	`43`	`forAll(x, i)`
`44`	`44`	`{`
`45`		`- tolerances_[i] = max(relTolerances_[i]*mag(x[i]), tolerances_[i]);`
	`45`	`+ tolerances_[i] = max(relTolerances_[i]*mag(x[i]), absTolerances_[i]);`
`46`	`46`	`}`
`47`	`47`	`}`
`48`	`48`