do not use Inf for Newton Bounds

boostorg · Jul 25, 2023 · e7ada7f · e7ada7f
1 parent 6bf5125
commit e7ada7f
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/include/boost/math/distributions/skew_normal.hpp b/include/boost/math/distributions/skew_normal.hpp
@@ -676,8 +676,8 @@ namespace boost{ namespace math{
 
     // refine the result by numerically searching the root of (p-cdf)
 
-    const RealType search_min = range(dist).first;
-    const RealType search_max = range(dist).second;
+    const RealType search_min = support(dist).first;
+    const RealType search_max = support(dist).second;
 
     const int get_digits = policies::digits<RealType, Policy>();// get digits from policy,
     std::uintmax_t m = policies::get_max_root_iterations<Policy>(); // and max iterations.