diff --git a/MathLib/Nonlinear/NewtonRaphson-impl.h b/MathLib/Nonlinear/NewtonRaphson-impl.h index a0eee7e0af81139033bdd42ccf19c8dbf78b894e..82190b5bc2459508b4703ba0bf5295e55e8078bd 100644 --- a/MathLib/Nonlinear/NewtonRaphson-impl.h +++ b/MathLib/Nonlinear/NewtonRaphson-impl.h @@ -30,9 +30,9 @@ NewtonRaphson::NewtonRaphson() template<class F_RESIDUAL, class F_DX, class T_VALUE> bool NewtonRaphson::solve(F_RESIDUAL &f_residual, F_DX &f_dx, const T_VALUE &x0, T_VALUE &x_new) { - const bool checkAbsResidualError = (_r_abs_tol<std::numeric_limits<double>::max()); - const bool checkRelResidualError = (_r_rel_tol<std::numeric_limits<double>::max()); - const bool checkRelDxError = (_dx_rel_tol>.0); + const bool checkAbsResidualError = (_r_abs_tol < std::numeric_limits<double>::max()); + const bool checkRelResidualError = (_r_rel_tol < std::numeric_limits<double>::max()); + const bool checkRelDxError = (_dx_rel_tol > .0); const bool needXNorm = (checkRelResidualError || checkRelDxError); INFO("------------------------------------------------------------------");