diff --git a/include/boost/math/special_functions/detail/bessel_jy.hpp b/include/boost/math/special_functions/detail/bessel_jy.hpp index cbc19900b..d3414496b 100644 --- a/include/boost/math/special_functions/detail/bessel_jy.hpp +++ b/include/boost/math/special_functions/detail/bessel_jy.hpp @@ -556,7 +556,7 @@ namespace boost { namespace math { if((cp != 0) && (tools::max_value() * fabs(Yv_scale) < fabs(cp * Yv))) *Y = org_kind & need_y ? T(-sign(cp) * sign(Yv) * sign(Yv_scale) * policies::raise_overflow_error(function, 0, pol)) : T(0); else - *Y = sp * Jv + (cp == 0 ? T(0) : T((cp * Yv) / Yv_scale)); + *Y = (sp != 0 ? sp * Jv : T(0)) + (cp == 0 ? T(0) : T((cp * Yv) / Yv_scale)); } else {