Как использовать 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; 
}