OpenMPの崩壊が間違った結果をもたらす

Question

私は3D配列zを持っています。すべての要素の値は1です。

今、私が行います

#pragma omp parallel for collapse(3) shared(z) 
for (int i=0; i < SIZE; ++i) { 
    for (int j=0; j < SIZE; ++j) { 
     for (int k=0; k < SIZE; ++k) { 
      for (int n=0; n < ITERATIONS-1; ++n) { 
       z[i][j][k] += 1; 
      } 
     } 
    } 
} 

これは、各要素にITERATIONSを追加する必要がありますし、それはありません。 collapse(3)をcollapse(4)に変更すると（for-loopsが4つあるため）、私は正しい結果を得られません。

4つのループをすべて折りたたむことはできませんか？

Answer 1

第4ループは、第3ループと同じように並列化できないという問題があります。あなた自身を説得するために、最後のループだけを考えてみてください。それはなる：

int zz = z[i][j][k]; 
for (int n=0; n < ITERATIONS-1; ++n) { 
    zz += 1; 
} 
z[i][j][k] = zz; 

それを並列化するために、あなたは右、reduction(+:zz)ディレクティブを追加する必要があるでしょうか？
あなたのcollapse(4)についてはまあ同じ話です。私はわからないんだけど、そのすべての可能な場合でも、reduction(+:z)を追加すると、いくつかの問題を提起します：

CまたはC++での配列のための reduction句は、以降のOpenMP 4.5でサポートされている、と私は知らない

1. コンパイラはそれを現時点でサポートしています（私は確信していますが）。
2. 縮小アスペクトを管理する複雑な仕組みのため、おそらくコードは非常に遅くなるでしょう。

要するに、必要に応じてcollapse(3)以下に固執するか、ループを別々に並列化します。