[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Bug-gsl] Problem in gsl_ran_gaussian_ziggurat
From: |
Brian Gough |
Subject: |
Re: [Bug-gsl] Problem in gsl_ran_gaussian_ziggurat |
Date: |
Mon, 05 Jul 2010 16:34:13 +0100 |
User-agent: |
Wanderlust/2.15.6 (Almost Unreal) Emacs/23.2 Mule/6.0 (HANACHIRUSATO) |
At Sat, 3 Jul 2010 19:19:45 +0100,
Gareth Haslip wrote:
> 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.
>
Can you send an example program which shows the problem. Thank you.
--
Brian Gough