boost/math
2
0
Fork 0
mirror of https://github.com/boostorg/math.git synced 2026-01-29 07:42:11 +00:00
Code Issues Projects Releases Wiki Activity

Merge pull request #1 from jzmaddock/patch-1

Browse Source 
change log(2) to constants::ln_two<T>()
This commit is contained in:
Anton Bikineev
2013-12-05 07:32:39 -08:00
parent 1edb61ffa3 fe28fea2de
commit 8e805f02dd
1 changed files with 3 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

6
boost/math/special_functions/detail/bessel_jy_derivatives_series.hpp
View File

@@ -55,7 +55,7 @@ inline T bessel_j_derivative_small_z_series(T v, T x, const Policy& pol)
}
else
{
prefix = (v - 1) * log(x / 2) - log(2) - boost::math::lgamma(v+1, pol);
prefix = (v - 1) * log(x / 2) - constants::ln_two<T>() - boost::math::lgamma(v+1, pol);
prefix = exp(prefix);
}
if (0 == prefix)
@@ -165,7 +165,7 @@ inline T bessel_y_derivative_small_z_series(T v, T x, const Policy& pol)
else
{
gam = boost::math::lgamma(v, pol);
p = (v + 1) * p + log(2);
p = (v + 1) * p + constants::ln_two<T>();
prefix = gam - log(boost::math::constants::pi<T>()) - p;
if (boost::math::tools::log_max_value<T>() < prefix)
{
@@ -197,7 +197,7 @@ inline T bessel_y_derivative_small_z_series(T v, T x, const Policy& pol)
else
{
int sgn;
prefix = boost::math::lgamma(-v, &sgn, pol) + (v - 1) * log(x / 2) - log(2);
prefix = boost::math::lgamma(-v, &sgn, pol) + (v - 1) * log(x / 2) - constants::ln_two<T>();
prefix = exp(prefix) * sgn / boost::math::constants::pi<T>();
}
bessel_y_derivative_small_z_series_term_b<T, Policy> s2(v, x);