Fix when search for root in interval containing zero.

copasi · Mar 13, 2019 · 11ba5f0 · 11ba5f0
1 parent 9fb9c65
commit 11ba5f0
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/copasi/utilities/CBrent.cpp b/copasi/utilities/CBrent.cpp
@@ -68,7 +68,7 @@ bool CBrent::findRoot(double a,            // Left border
   while ((fabs(b - a) > vTol || fabs(fb - fa) > fTol) &&
          fb != 0.0 &&
          fs != 0.0 &&
-         fabs(b - a) > fabs(b + a) * 50.0 * std::numeric_limits< double >::epsilon())
+         fabs(b - a) > (fabs(b) + fabs(a)) * 50.0 * std::numeric_limits< double >::epsilon())
     {
       if (fa != fc && fb != fc)
         {