This is the mail archive of the
gsl-discuss@sources.redhat.com
mailing list for the GSL project.
Re: bug report on random number generators
- From: Brian Gough <bjg at network-theory dot co dot uk>
- To: Heiko Bauke <Heiko dot Bauke at Physik dot Uni-Magdeburg dot DE>
- Cc: gsl-discuss at sources dot redhat dot com
- Date: Fri, 21 Feb 2003 20:33:53 +0000
- Subject: Re: bug report on random number generators
- References: <1045775039.3e5542bf72f2d@webmail.uni-magdeburg.de>
Heiko Bauke writes:
> I have found some bugs in the initialization routine for generators
> gsl_rng_transputer, gsl_rng_randu, gsl_rng_borosh13 and some other
> serious bugs.
>
Thanks for the bug report and patches.
Sorry for not replying to your previous message. I forwarded it to
the original author of the Knuth routines, Carlo Perassi, and he is
working on fixing them.
Ideally these bugs should have been caught before they were included
by the test of computing the 10,000th value analytically but I must
have let them slip through without checking that.
I will correct the period in the comments from 2^32 to 2^30. Thanks
for pointing that out.
I don't use (s|1) since it produces identical sequences for even and
odd seeds. Although this allows seeds with sub-maximal period this is
the intention -- the 'randu' type generators are provided as a way to
produce poor-quality sequences.
--
Brian