C++ openMP-对变量和原子的并发访问

我想与openMP并行化一个函数，该函数对一个盒子进行采样（在盒子中随机选择点，并在这些点上计算给定函数）。我编写了以下代码

//storing points
double** points_ = new double*[N-m];
for(int i=0;i<N-m;i++)
{
    points_[i]=new double[ndim];
}
double* evals_ = new double[N-m];

#pragma omp parallel for 
for(int i=0;i<N-m;i++)
{
    double* pt_ = randomPoint(lower,upper);
    for(int k=0;k<ndim;k++)
    {
        points_[i][k]=pt_[k];
    }
    evals_[i]=evalFunc(pt_);
    delete pt_;
}

//存储点
双**点=新双*[N-m]；
对于（inti=0；i，代码中不需要任何原子子句

原因是外循环在索引
i
上被分割，因此每个线程将从
点和
评估点获得一组元素，这些元素不会与另一个线程的工作集重叠

evals\是一个数组，因此每个线程将获得一个连续的子数组（由于隐式静态调度方案），例如

点
是一个二维矩阵，每个线程将获得一组连续的行：

   0 
   1
t1 2
   3
   .
   4
   5
t2 6
   7
   .
   8
   9
t3 10
   11
   ...

在第二种情况下，
k
的值似乎有重叠，因为每个线程的
k
范围相同，但更新的点落在不同的行上（索引
i
），如上图所示，不同线程之间不会重叠。
我同意！但是，现在我不明白为什么代码的另一部分（我用这个示例编辑了上面的文章）会崩溃，告诉我在
sum+=set[j][I]行中存在并发访问
。谢谢！@dlib:这不应该是
#pragma omp parallel for reduce…
？您似乎缺少
的
，这为计算提供了完全不同的语义。
void atomic_example(float *x, float *y, int *index, int n)
{
    int i;
    #pragma omp parallel for shared(x, y, index, n)
    for (i=0; i<n; i++) {
    #pragma omp atomic update
    x[index[i]] += work1(i);
    y[i] += work2(i);
    }
}

for(i=0;i<m_ndim;i++)
{
    double sum_=0;
    #pragma omp parallel reduction(+:sum_)
    for(j=0;j<m_npts;j++)
    {
        sum_ += set_[j][i];
    }
    Sum_[i] = sum_;
}

0 1 2 3 . 4 5 6 7 . 8 9 10 11...
   t1        t2        t3

   0 
   1
t1 2
   3
   .
   4
   5
t2 6
   7
   .
   8
   9
t3 10
   11
   ...