boost/math
2
0
Fork 0
mirror of https://github.com/boostorg/math.git synced 2026-01-28 19:32:08 +00:00
Code Issues Projects Releases Wiki Activity

Repaired some typos in the multiprecision code in gamma_imp_bernoulli().

Browse Source
This commit is contained in:
Christopher Kormanyos
2014-01-02 14:29:48 +01:00
parent 6db6405d54
commit 7dff191552
1 changed files with 3 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

7
include/boost/math/special_functions/gamma.hpp
View File

@@ -236,7 +236,6 @@ T gamma_imp_bernoulli(T x, const Policy& pol)
T gamma_value = exp(log_gamma_value);
// Rescale the result using downward recursion if necessary.
// Rescale the result using downward recursion if necessary.
for(int k = 0; k < n_recur; ++k)
{
@@ -251,8 +250,8 @@ T gamma_imp_bernoulli(T x, const Policy& pol)
}
else
{
const T sin_pi_xx = (is_near_a_negative_integer ? sinpx(xx)
: sin(boost::math::constants::pi<T>() * xx));
const T sin_pi_xx = (is_near_a_negative_integer ? sinpx(xx)
: sin(boost::math::constants::pi<T>() * xx));
return -boost::math::constants::pi<T>() / (xx * gamma_value * sin_pi_xx);
}
@@ -532,7 +531,7 @@ T gamma_imp(T z, const Policy& pol, const lanczos::undefined_lanczos&)
template <class T, class Policy>
T lgamma_imp(T z, const Policy& pol, const lanczos::undefined_lanczos&, int*)
{
return lgamma_imp_bernoulli(z,pol);
return lgamma_imp_bernoulli(z, pol);
}
//