[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [Bug-gsl] Wrong algorithm in gsl_rng_r250

From: Brian Gough
Subject: Re: [Bug-gsl] Wrong algorithm in gsl_rng_r250
Date: Wed, 5 May 2004 17:57:35 +0100

Andreas Schneider writes:
 > The documentation claims that gsl_rng_r250 does
 >    x_n = x_{n-103} ^ x_{n-250},
 > but actually the implementation does
 >    x_n = x_{n-147} ^ x_{n-250}.
 > This error seems to be very common. It must have been introduced
 > long time ago.
 > I fed some data from both variants into Marsaglia's diehard and
 > found that the documented version is better in most tests. Thus the
 > error is probably in the implementation and the documentation is
 > right.

Thanks for the bug report.  I'll make a note of it in the BUGS file.

I have checked the original paper and there is a discrepancy there
between the source code and the paper itself.

Brian Gough

reply via email to

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