From octave-maintainers-request at bevo dot che dot wisc dot edu Fri Jan 23 07:42:30 2004
Subject: Re: randn benchmarks
From: Paul Kienzle <pkienzle at users dot sourceforge dot net>
To: octave-maintainers at bevo dot che dot wisc dot edu
Date: Fri, 23 Jan 2004 08:41:33 -0500


On Jan 23, 2004, at 7:36 AM, David Bateman wrote:
> ... In any case, the ziggurat code in
> this paper is for 2^32 bit seeds, which is a bit short for my liking.
> I'd prefer to target a 64 bit seed as does matlab, so a possible
> solution is just to rewrite the code from the paper for a 64 bit
> seed..

Do you want to generate 53 bits of precision in rand() or are you
happy with 32 bit?  It will of course slow down the code by almost
a factor of two.

Paul Kienzle
pkienzle at users dot sf dot net