diff -up boost_1_54_0/boost/math/distributions/inverse_gaussian.hpp~ boost_1_54_0/boost/math/distributions/inverse_gaussian.hpp
|
old
|
new
|
namespace detail
|
| 285 | 285 | |
| 286 | 286 | // Define the distribution, using gamma_nooverflow: |
| 287 | 287 | typedef gamma_distribution<RealType, no_overthrow_policy> gamma_nooverflow; |
| | 288 | gamma_nooverflow g(static_cast<RealType>(0.5), static_cast<RealType>(1.)); |
| 288 | 289 | |
| 289 | | gamma_distribution<RealType, no_overthrow_policy> g(static_cast<RealType>(0.5), static_cast<RealType>(1.)); |
| 290 | | |
| 291 | | // gamma_nooverflow g(static_cast<RealType>(0.5), static_cast<RealType>(1.)); |
| 292 | 290 | // R qgamma(0.2, 0.5, 1) 0.0320923 |
| 293 | 291 | RealType qg = quantile(complement(g, p)); |
| 294 | 292 | //RealType qg1 = qgamma(1.- p, 0.5, 1.0, true, false); |