OpenMPの最小値の配列

Question

は、私は元のコードがあります。

min = INT_MAX; 
for (i=0;i<N;i++) 
    if (A[i]<min) 
    min = A[i]; 
for (i=0;i<N;i++) 
    A[i]=A[i]-min; 

を、私はこれの並列バージョンを取得したいと私はこれでした：

min = INT_MAX; 
#pragma omp parallel private(i){ 
minl = INT_MAX; 
#pragma omp for 
for (i=0;i<N;i++) 
    if (A[i]<minl) 
    minl=A[i]; 
#pragma omp critical{ 
if (minl<min) 
    min=minl; 
} 
#pragma omp for 
for (i=0;i<N;i++) 
    A[i]=A[i]-min; 
} 

は、並列コードの右か？ #pragma omp barrierを#pragma omp criticalの前に書く必要があるのか​​どうか疑問に思っていたので、グローバル最小値を計算する前にすべての最小値が計算されていることを確認してください。

Answer 1

コードは正しいです。 1つのスレッドがクリティカルセクションに入ったときにすべてmin_lを計算する必要がないため、#pragma omp barrierを追加する必要はありません。また、ループ領域の終わりに暗黙のバリアがあります。

さらに、ループ反復変数iを明示的に宣言する必要はありません。

あなたが代わりにminlのあなたの手動マージの削減を使用してコードを向上させることができます。

#pragma omp for reduction(min:min) 
for (i=0;i<N;i++) 
    if (A[i]<min) 
    min=A[i]; 

注：削減のためのminオペレータは、OpenMP 3.1以降で使用可能です。