openMP - 對變量和原子的並發訪問

Question

我想與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_;
}

但是，我對此代碼沒有信心：evals_和points_在eauch線程中更新。 我想在那里添加一些原子語句：

#pragma omp parallel for 
for(int i=0;i<N-m;i++)
{
    double* pt_ = randomPoint(m_lower,m_upper);
    for(int k=0;k<m_ndim;k++)
    {
        #pragma omp atomic update
        points_[i][k]=pt_[k];
    }
    #pragma omp atomic update
    evals_[i]=evalFunc(pt_);
    delete pt_;
}

但我擔心這會非常低效：你有什么建議可以更准確地寫出來嗎？ 並且...這不是編譯...（錯誤：#pragma omp原子之后的表達式有不正確的形式）雖然我可以在openMP規范中找到該檢查，A22

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);
    }
}

並且原子更新也遵循對數組的影響。

感謝致敬。

編輯 - - - -

我同意都鐸的回答。 但是，似乎這個例子在另一個並行化的代碼片段中確實需要原子：on sum sum _ + = ...，發生錯誤（並發訪問）

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_;
}

為什么需要呢？ 還是別的錯了？

Answer 1

您的代碼中不需要任何原子句。

原因是外部循環在索引i被分割，因此每個線程將從points_和eval_獲得一組元素， points_元素不與另一個線程的工作集重疊。

evals_是一個數組，因此每個線程將獲得一個連續的子數組（由於隱式靜態調度方案），例如

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

points_是一個二維矩陣，每個線程將獲得一組連續的行：

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

在第二種情況下，看起來你的k值有重疊，因為每個線程的k范圍相同，但更新的點落在不同的行（索引i ）上，如上所示，不同線程的重疊。