このループではどのようにopenMPを使用できますか？

Question

openMPを使ってこのループを並列にすることが可能かどうか疑問に思っています。

大まかでは、競争条件に問題があります。私は外側ループによって生成される内側ループのnと、D = A [n]との競合条件をどのように扱うのか不明です。これを並列化しようとするのが実用的だと思いますか？

for(n=0; n < 10000000; ++n) { 

    for (n2=0; n2< 100; ++n2) { 
     A[n]=A[n]+B[n2][n+C[n2]+200]; 

     } 

    D=D+A[n]; 

} 

Answer 1

はい、ポインタがエイリアス化されていないと仮定すると、これは実際には並列化可能です。

次のようにそれは実際にやや最適化することができ

int D = 0; // Or whatever the type is. 

#pragma omp parallel for reduction(+:D) private(n2) 
for (n=0; n < 10000000; ++n) { 

    for (n2 = 0; n2 < 100; ++n2) { 
     A[n] = A[n] + B[n2][n + C[n2] + 200]; 
    } 

    D += A[n]; 
} 

：

int D = 0; // Or whatever the type is. 

#pragma omp parallel for reduction(+:D) private(n2) 
for (n=0; n < 10000000; ++n) { 

    int tmp = A[n] 
    for (n2 = 0; n2 < 100; ++n2) { 
     tmp += B[n2][n + C[n2] + 200]; 
    } 

    A[n] = tmp; 
    D += tmp; 
}