Merge pull request #6387 from jyoo1042/abcd

Change the error condition for bounds on TOMS748
sxs-collaboration · Dec 3, 2024 · 881ac83 · 881ac83
2 parents 98eb769 + 38eb651
commit 881ac83
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/NumericalAlgorithms/RootFinding/TOMS748.hpp b/src/NumericalAlgorithms/RootFinding/TOMS748.hpp
@@ -316,7 +316,7 @@ std::pair<T, T> toms748_solve(F f, const T& ax, const T& bx, const T& fax,
                               size_t& max_iter) {
   // Main entry point and logic for Toms Algorithm 748
   // root finder.
-  if (UNLIKELY(simd::any(ax >= bx))) {
+  if (UNLIKELY(simd::any(ax > bx))) {
     throw std::domain_error("Lower bound is larger than upper bound");
   }