Я написал код Монте-Карло (MC) на Python с расширением Fortran (скомпилированным с помощью f2py). Поскольку это стохастическое интегрирование, алгоритм сильно зависит от случайных чисел, а именно я использую ~ 10^8 - 10^9 случайных чисел для типичного запуска. До сих пор я действительно не возражал против «качества» случайных чисел — однако это то, что я хочу проверить.
Мой вопрос: достаточно ли Mersenne-Twister, используемого numpy, или есть лучшие генераторы случайных чисел, которые следует (можно) использовать? (лучше в смысле времени выполнения, а также качества сгенерированной последовательности)
Любые предложения/опыты определенно приветствуются, спасибо!
