import new upstream nautilus stable release 14.2.8

[ceph.git] / ceph / src / boost / boost / math / special_functions / detail / bessel_jy.hpp

diff --git a/ceph/src/boost/boost/math/special_functions/detail/bessel_jy.hpp b/ceph/src/boost/boost/math/special_functions/detail/bessel_jy.hpp
index b67d989b689d7d6da64c0d6d78fb80ccf22e8603..315b932df1a520ec21ef2cdc2311b24313275a4f 100644 (file)
--- a/ceph/src/boost/boost/math/special_functions/detail/bessel_jy.hpp
+++ b/ceph/src/boost/boost/math/special_functions/detail/bessel_jy.hpp
@@ -561,11 +561,11 @@ namespace boost { namespace math {
          if (reflect)
          {
             if((sp != 0) && (tools::max_value<T>() * fabs(Yv_scale) < fabs(sp * Yv)))
-               *J = org_kind & need_j ? T(-sign(sp) * sign(Yv) * sign(Yv_scale) * policies::raise_overflow_error<T>(function, 0, pol)) : T(0);
+               *J = org_kind & need_j ? T(-sign(sp) * sign(Yv) * (Yv_scale != 0 ? sign(Yv_scale) : 1) * policies::raise_overflow_error<T>(function, 0, pol)) : T(0);
             else
                *J = cp * Jv - (sp == 0 ? T(0) : T((sp * Yv) / Yv_scale));     // reflection formula
             if((cp != 0) && (tools::max_value<T>() * fabs(Yv_scale) < fabs(cp * Yv)))
-               *Y = org_kind & need_y ? T(-sign(cp) * sign(Yv) * sign(Yv_scale) * policies::raise_overflow_error<T>(function, 0, pol)) : T(0);
+               *Y = org_kind & need_y ? T(-sign(cp) * sign(Yv) * (Yv_scale != 0 ? sign(Yv_scale) : 1) * policies::raise_overflow_error<T>(function, 0, pol)) : T(0);
             else
                *Y = (sp != 0 ? sp * Jv : T(0)) + (cp == 0 ? T(0) : T((cp * Yv) / Yv_scale));
          }