[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Bug-gsl] gsl_rng_uniform_pos generates boundary value
From: |
Giorgos |
Subject: |
[Bug-gsl] gsl_rng_uniform_pos generates boundary value |
Date: |
Thu, 18 Oct 2012 23:54:21 +0100 |
User-agent: |
Mozilla/5.0 (X11; Linux x86_64; rv:16.0) Gecko/20121011 Thunderbird/16.0.1 |
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
- [Bug-gsl] gsl_rng_uniform_pos generates boundary value,
Giorgos <=