|
From: | Andreas Schneider |
Subject: | [Bug-gsl] Wrong algorithm in gsl_rng_r250 |
Date: | Sun, 25 Apr 2004 13:23:39 +0200 |
User-agent: | Mozilla/5.0 (X11; U; Linux i686; en; Stable) Gecko/20020827 Beonex/0.9-pre-1 |
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.
Andreas Schneider
[Prev in Thread] | Current Thread | [Next in Thread] |