bug-gsl
[Top][All Lists]
Advanced

[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




reply via email to

[Prev in Thread] Current Thread [Next in Thread]