[Numpy-discussion] Suggestion: Port Theano RNG implementation to NumPy

Sturla Molden sturla.molden at gmail.com
Tue Feb 18 11:47:36 EST 2014

Matthieu Brucher <matthieu.brucher at gmail.com> wrote:
> Hi,
> The main issue with PRNG and MT is that you don't know how to
> initialize all MT generators properly. A hash-based PRNG is much more
> efficient in that regard (see Random123 for a more detailed
> explanation).
>> From what I heard, if MT is indeed chosen for RNG in numerical world,
> in parallel world, it is not as obvious because of this pitfall.

It is possible to solve this by using a set of independent MT generators,
one per thread. Independence in this case means that the characteristic
polynomials are relatively prime to each other:


Undortunately the DCMT code was LGPL, not BSD, I don't know if this has


More information about the NumPy-Discussion mailing list