diff options
Diffstat (limited to 'boost/math/distributions/inverse_gamma.hpp')
-rw-r--r-- | boost/math/distributions/inverse_gamma.hpp | 9 |
1 files changed, 6 insertions, 3 deletions
diff --git a/boost/math/distributions/inverse_gamma.hpp b/boost/math/distributions/inverse_gamma.hpp index 88083e084f..fa5d357ac7 100644 --- a/boost/math/distributions/inverse_gamma.hpp +++ b/boost/math/distributions/inverse_gamma.hpp @@ -91,13 +91,13 @@ public: typedef RealType value_type; typedef Policy policy_type; - inverse_gamma_distribution(RealType shape = 1, RealType scale = 1) - : m_shape(shape), m_scale(scale) + inverse_gamma_distribution(RealType l_shape = 1, RealType l_scale = 1) + : m_shape(l_shape), m_scale(l_scale) { RealType result; detail::check_inverse_gamma( "boost::math::inverse_gamma_distribution<%1%>::inverse_gamma_distribution", - scale, shape, &result, Policy()); + l_scale, l_shape, &result, Policy()); } RealType shape()const @@ -259,6 +259,9 @@ inline RealType cdf(const complemented2_type<inverse_gamma_distribution<RealType if(false == detail::check_inverse_gamma_x(function, c.param, &result, Policy())) return result; + if(c.param == 0) + return 1; // Avoid division by zero + //result = 1. - gamma_q(shape, c.param / scale, Policy()); result = gamma_p(shape, scale/c.param, Policy()); return result; |