diff options
Diffstat (limited to 'boost/math/special_functions/ellint_d.hpp')
| -rw-r--r-- | boost/math/special_functions/ellint_d.hpp | 5 |
1 files changed, 5 insertions, 0 deletions
diff --git a/boost/math/special_functions/ellint_d.hpp b/boost/math/special_functions/ellint_d.hpp index bc5a4b2a56..fa5c53db18 100644 --- a/boost/math/special_functions/ellint_d.hpp +++ b/boost/math/special_functions/ellint_d.hpp @@ -82,8 +82,12 @@ T ellint_d_imp(T phi, T k, const Policy& pol) s = -1; rphi = constants::half_pi<T>() - rphi; } + BOOST_MATH_INSTRUMENT_VARIABLE(rphi); + BOOST_MATH_INSTRUMENT_VARIABLE(m); T sinp = sin(rphi); T cosp = cos(rphi); + BOOST_MATH_INSTRUMENT_VARIABLE(sinp); + BOOST_MATH_INSTRUMENT_VARIABLE(cosp); T c = 1 / (sinp * sinp); T cm1 = cosp * cosp / (sinp * sinp); // c - 1 T k2 = k * k; @@ -99,6 +103,7 @@ T ellint_d_imp(T phi, T k, const Policy& pol) { // http://dlmf.nist.gov/19.25#E10 result = s * ellint_rd_imp(cm1, T(c - k2), c, pol) / 3; + BOOST_MATH_INSTRUMENT_VARIABLE(result); } if(m != 0) result += m * ellint_d_imp(k, pol); |