Openmp for循环，种子随机数发生器在哪里？

Question

我读过，为了保证线程安全，它的方便种子并行区域内的RNG这样的：

int seedbase = 392872; 

#pragma omp parallel 
{ 
    srand(omp_get_thread_id * seedbase); 
    #pragma omp for 
    .... 
} 

但是，如果我的并行部分是另一种循环内？如果我有这样的情况：

int seedbase = 392872; 
for(int i=0; i<100; ++i) 
{ 
    #pragma omp parallel 
    { 
     srand(omp_get_thread_id * seedbase); 
     #pragma omp for 
     .... 
    } 
} 

我应该在哪里初始化我RNG？

Answer 1

我读过，为了保证线程安全，它的方便，播种RNG并行区域内是这样的：

便利尽管如此，你的技术是不是有效。不管您拨打srand()的哪个位置，标准rand()函数都不能成为线程安全的。它依赖于每次调用都会修改的内部静态数据，所以如果您从多个线程调用它（或srand()）而没有进行某种同步，则会因此创建数据竞争。

一次，POSIX定义了一个函数rand_r()。如果你有这个，那么它会符合你的目的，但它现在已被标记为过时。在OpenMP中正确使用rand_r()将涉及在并行区域内部建立私有（OpenMP感知）变量来保存种子。根据每个线程中是否需要相同的随机数序列，在每个线程中以不同或相同方式对其进行初始化。然后将一个指向该变量的指针作为参数传递给rand_r()。在每次调用之后，您可能想要将返回值或从中派生出来的东西作为新种子。