Thu, 29 Jul 2010 02:40:00 +0200, Sturla Molden wrote: >> Mon, 26 Jul 2010 23:58:11 +0800, Ralf Gommers wrote: >> Is the current algorithm in the trunk the ziggurat one, or the previous >> one? IIRC, the problem was that the ziggurat broke reproducibility of >> random numbers with a given seed. > > Ziggurat (in Enthought) did not break reproducibility but backwards > compatibility.
Exactly what I meant, it breaks reproducibility of old results with the same code. > It might be better to put a ziggurat module in scipy.random, just to > avoid confusion. Also it should be modified to use a Mersenne Twister > instead of the simple uniform generator in Marsaglia's code, as well as > getting rid of global variables. Want to submit a patch? ;) Cheers, Pauli _______________________________________________ NumPy-Discussion mailing list NumPy-Discussion@scipy.org http://mail.scipy.org/mailman/listinfo/numpy-discussion