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); |