随机数还有蒙特C

Question

我用Python编写了一个montecarlo（MC）代码，扩展名为Fortran（用f2py编译）。由于这是一个随机积分，算法很大程度上依赖于随机数，即我使用~ 10^8 - 10^9随机数进行典型运行。到目前为止，我并不在意随机数的“质量”——不过，这是我想检查的东西。在

我的问题是：纽比使用的梅森捻线器是否足够，或者是否有更好的随机数生成器可以使用？（在运行时和生成序列的质量方面更好）

欢迎任何建议/经验，谢谢！在

Answer 1

我认为没有人可以告诉你这个算法是否足够，而不知道随机数是如何使用的。在

我要做的是用其他东西代替numpy随机数，当然还有其他模块提供了不同的算法。 如果您的模拟结果不受随机数生成器选择的影响，这已经是一个好迹象。在