[Bug-gsl] Problem in gsl_ran_gaussian_ziggurat

[Gareth Haslip]

Hi,

I think the line:

j = gsl_rng_uniform_int (r, 16777216);  /* sample from 2^24 */

should be j = gsl_rng_uniform_int (r, 16777216-1);  /* sample from 2^24 */

as it crashes on 2^24 (exceeds upped limit for gsl_rng_uniform_int.

Thanks,

Gareth