Hello, Le 30/03/2017 à 13:31, Pierre Haessig a écrit : > [....] > > But how come Julia is 4-5x faster since Numpy uses C implementation > for the entire process ? (Mersenne Twister -> uniform double -> > Box-Muller transform to get a Gaussian > https://github.com/numpy/numpy/blob/master/numpy/random/mtrand/randomkit.c). > Also I noticed that Julia uses a different algorithm (Ziggurat Method > from Marsaglia and Tsang , > https://github.com/JuliaLang/julia/blob/master/base/random.jl#L700) > but this doesn't explain the difference for uniform rng. > Any ideas?
Do you think Stackoverflow would be a better place for my question? best, Pierre
_______________________________________________ NumPy-Discussion mailing list NumPy-Discussion@python.org https://mail.python.org/mailman/listinfo/numpy-discussion