id,summary,reporter,owner,description,type,status,milestone,component,version,severity,resolution,keywords,cc
13620,Issue with computing quantiles for skew-normal,Ben Bales <bbbales2@…>,John Maddock,"The 0.075 skew-normal quantile breaks (returns an inf) with certain, seemingly reasonable parameters. Other surrounding values are okay. The complemented version works for some reason.

Here's some test code
{{{
#include <iostream>
#include <boost/math/distributions.hpp> 

int main()
{
    boost::math::skew_normal_distribution<> dist(2.0, 1.0, -2.5);

    std::vector<double> qs = { 0.05, 0.075, 0.1 };
    for(auto q : qs) {
        std::cout << ""q : "" << q << std::endl;
        std::cout << ""w/o complement: "" << quantile(dist, q) << std::endl;
        std::cout << ""w complement: "" << quantile(complement(dist, 1 - q)) << std::endl << std::endl;
    }
}
}}}

Dunno how easy this is to fix. I can work around it where I'm using it, but it seems like something that could be fixed.

Thanks!",Bugs,closed,To Be Determined,math,Boost 1.66.0,Problem,obsolete,,