<html><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Hi all,<div><br></div><div>I have a question concerning the Mersenne Twister random number generation in numpy: &nbsp;when I seed it with 0, I get a different sequence of numbers in numpy, compared to GSL.</div><div><br></div><div>In numpy:&nbsp;</div><div><br></div><div>r =&nbsp;numpy.Random.RandomState(seed=0)</div><div>r.uniform(size=5) &nbsp;----> &nbsp;array([ 0.5488135 , &nbsp;0.71518937, &nbsp;0.60276338, &nbsp;0.54488318, &nbsp;0.4236548 ])</div><div><br></div><div>whereas in GSL the first numbers are&nbsp;<span class="Apple-style-span" style="font-family: -webkit-monospace; font-size: 13px; white-space: pre; ">&nbsp;</span></div><div><span class="Apple-style-span" style="font-family: -webkit-monospace; font-size: 13px; white-space: pre; ">0.99974    0.16291   0.2826    0.94720&nbsp;   <span class="Apple-style-span" style="font-family: Helvetica; font-size: 12px; white-space: normal; "><span class="Apple-style-span" style="font-family: -webkit-monospace; font-size: 13px; white-space: pre; ">0.23166</span>&nbsp;</span></span></div><div><br></div><div>Matlab gives the same result as numpy...</div><div><br></div><div>I have translated some python code to c, and would like to debug it -- therefore, I would like to have exactly the same set of random numbers... &nbsp;How can I provoke this ?</div><div><br></div><div>Best.</div><div><br></div><div>Giovanni</div></body></html>