[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Bug-gsl] Problem in gsl_ran_gaussian_ziggurat
From: |
Gareth Haslip |
Subject: |
[Bug-gsl] Problem in gsl_ran_gaussian_ziggurat |
Date: |
Sat, 3 Jul 2010 19:19:45 +0100 |
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
- [Bug-gsl] Problem in gsl_ran_gaussian_ziggurat,
Gareth Haslip <=