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 перед критическим #pragma omp, чтобы я убедился, что все минимумы вычисляются перед вычислением глобального минимума.

Answer 1

Код верный. Нет необходимости добавлять #pragma omp barrier, потому что нет необходимости, чтобы все min_l вычислялись, когда один поток попадает в критический раздел. Существует также неявный барьер в конце области петли.

Кроме того, вам не обязательно явно указывать переменную итерации цикла 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.