Hello,According to the GSL documentation, the gsl_rng_uniform_pos call should not return a boundary value (0 or 1). However, I noticed that for a given state of the RNG it generates the value 1.0. I have saved the RNG state in a file which I attach here, but I am not sure whether it is portable... Here is the program to generate the uniform
#include <stdio.h>
#include <gsl/gsl_randist.h>
int main()
{
FILE * seed = fopen("seed", "r");
gsl_rng * RNG = gsl_rng_alloc(gsl_rng_mt19937);
gsl_rng_fread(seed, RNG);
printf("%f\n", gsl_rng_uniform(RNG));
return 0;
}
I use GSL-1.15, GCC-4.7. If you would like me to give any more
information, please ask :)
seed
Description: Binary data
