random_number() 如何并行工作？

Question

如何random_number()与 OpenMP 并行工作？

如果我在没有并行化的情况下运行我的程序，我总是会得到相同的结果，但是通过并行化，我每次都会得到不同（但相似）的结果。

Answer 1

一般没有关于线程安全或线程性能的保证random_number。Fortran 标准根本不知道 OpenMP。

个别编译器可能会为您提供一些保证，但它们仅对特定编译器中存在的版本有效。例如，当前的 gfortran 版本提供了一个线程安全的随机数生成器和“请注意，在多线程程序中（例如使用 OpenMP 指令），每个线程都将有自己的随机数状态。” 其他编译器可能不同。值得注意的是，您的用户可能想要使用的编译器可能有所不同，而您可能对此一无所知。

有专用的并行随机数生成器可用。例如，我使用库的修改版本，该版本使用 Ziggurat 方法处理多个随机数分布，由 Gib Bogle 并行化，并添加了 xoroshiro128+ 的实现作为底层算法，类似于 Gfortran 使用的算法。还有其他类似算法的实现可用，标准 C++ 包含一些新的生成器，这些生成器实际上定义为使用特定算法，因此您可以调用它们。