c++ - 与相同引擎并行生成随机数失败

Question

我正在使用 C++11 提供的 RNG，我也在玩弄 OpenMP。我为每个线程分配了一个引擎，作为测试，我为每个引擎提供了相同的种子。这意味着我希望两个线程都能产生完全相同的随机生成数字序列。这是一个MWE：

#include <iostream>
#include <random>

using namespace std;


uniform_real_distribution<double> uni(0, 1);
normal_distribution<double> nor(0, 1);


int main()
{
    #pragma omp parallel
    {
        mt19937 eng(0); //GIVE EACH THREAD ITS OWN ENGINE
        vector<double> vec;

        #pragma omp for
        for(int i=0; i<5; i++)
        {
            nor(eng);
            vec.push_back(uni(eng));
        }
        #pragma omp critical
        cout << vec[0] << endl;
    }



    return 0;
}

大多数情况下，我得到了输出0.857946 0.857946，但有几次我得到了0.857946 0.592845。当两个线程具有相同的、不相关的引擎时，后一种结果怎么可能？！

Answer 1

你也必须把nor和uni放在omp parallel区域内。像这样：

#pragma omp parallel
{
    uniform_real_distribution<double> uni(0, 1);
    normal_distribution<double> nor(0, 1);
    mt19937 eng(0); //GIVE EACH THREAD ITS OWN ENGINE
    vector<double> vec;

否则每个线程只有一个副本，而实际上每个线程都需要自己的副本。

更新添加：我现在看到在这个 stackoverflow 线程中讨论了完全相同的问题 。