Ticket #7779: bug_uniform_on_sphere.cpp

#include <boost/circular_buffer.hpp>
#include <boost/random/variate_generator.hpp>
#include <boost/random/uniform_on_sphere.hpp>
#include <boost/random.hpp>
#include <iostream>
#include <fstream>
using namespace std;

int main(int argc, const char *argv[])
{
    typedef boost::mt19937                                              GeneratorType;
    typedef boost::uniform_on_sphere<float>                             DistributionType;
    typedef boost::variate_generator<GeneratorType, DistributionType >  VariateType;
    typedef boost::circular_buffer<DistributionType::result_type>       BufferType;
    const int           BUFSIZE = 10;
    const int           DIM = 2;
    GeneratorType       gen;
    DistributionType    dist(DIM);
    VariateType         variate(gen,dist);

    gen.seed(11);
    BufferType buf(BUFSIZE);
    long n(0);
    DistributionType::result_type tmp;
    while (1){
        cout << "n: " << n << "\r" << flush;
        tmp = variate();
        buf.push_back(tmp);
        if (isnan(tmp[0])) {
            cout << "n: " << n << "       " << endl;
            for (size_t i=0; i<DIM; i++)
                cout << tmp[i] << " , ";
            cout << endl;
            ofstream fout("debug.out");
            for (int i=0; i<BUFSIZE; i++)
                fout << buf[i][0] << " " << buf[i][1] << endl;
            fout.close();
            ofstream gen_out("gen.out");
            gen_out << gen;
            gen_out.close();
            exit(1);
        }
        n++;
    }

    return 0;
}