UpsilonNumworks/Upsilon
1
0
Fork 0
mirror of https://github.com/UpsilonNumworks/Upsilon.git synced 2026-01-18 16:27:34 +01:00
Code Issues Packages Projects Releases Wiki Activity

[poincare/solver] Fix typo in BrentRoot algorithm

Browse Source
This commit is contained in:
Ruben Dashyan
2020-02-26 16:47:01 +01:00
committed by EmilieNumworks
parent 74f7705413
commit bb24c21320
1 changed files with 2 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
poincare/src/solver.cpp
View File

@@ -138,7 +138,7 @@ double Solver::BrentRoot(double ax, double bx, double precision, ValueAtAbscissa
return b;
}
}
if (std::fabs(e) >= tol1 && std::fabs(fa) > std::fabs(b)) {
if (std::fabs(e) >= tol1 && std::fabs(fa) > std::fabs(fb)) {
double s = fb/fa;
double p = 2.0*xm*s;
double q = 1.0-s;
@@ -168,7 +168,7 @@ double Solver::BrentRoot(double ax, double bx, double precision, ValueAtAbscissa
if (std::fabs(d) > tol1) {
b += d;
} else {
b += xm > 0.0 ? tol1 : tol1;
b += xm > 0.0 ? tol1 : -tol1;
}
fb = evaluation(b, context, complexFormat, angleUnit, context1, context2, context3);
}