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原子具有不正確形式）althought我可以找到OpenMP的規格是examlpe，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); 
    } 
} 

和原子更新也是在這裏通過跟隨對數組的一種情感。

感謝和問候。

編輯--------

我同意鐸的回答。然而，似乎這個例子在這裏，在另一併行的代碼確實需要原子：行和_ + = ...，出現錯誤（併發訪問）

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_獲得一組元素，它們不會與另一個線程的工作集重疊。

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），如上所示，它們對於不同的線程不重疊。