Почему этот код openmp дает ошибку сегментации?

Question

Я пытаюсь выполнить код OpenMP и был успешным в этом. Однако у меня есть сомнения относительно заявления #pragma omp parallel.

Рассмотрим эти два фрагмента кода:

#pragma omp parallel firstprivate(my_amax) 
{ 
    for (i=0; i<MATRIX_DIM; i++) { 
     #pragma omp for 
     for (j=0; j<MATRIX_DIM; j++) { 
      my_amax = abs_max(my_amax, A[i][j], B[i][j]); 
      #pragma omp critical 
      { 
       if(fabs(amax)<fabs(my_amax)) 
        amax=my_amax; 
      } 
     } 
    } 
} 

И

for (i=0; i<MATRIX_DIM; i++) { 
    #pragma omp parallel firstprivate(my_amax) 
    { 
     #pragma omp for 
     for (j=0; j<MATRIX_DIM; j++) { 
      my_amax = abs_max(my_amax, A[i][j], B[i][j]); 
      #pragma omp critical 
      { 
      if (fabs(amax)<fabs(my_amax)) 
       amax=my_amax; 
      } 
     } 
    } 
} 

Единственное различие в коде положение параллельной части. Первый код всегда дает мне ошибку сегментации, а второй - отлично. Почему это так?

Я знаю, что #pragam omp parallel порождает необходимые темы, но так как следующий i цикл не объявлен как параллельно, это не должно быть проблемой, т.е. i часть должна получить выполняются последовательно в то время как j итераций, которые на самом деле распараллеленные будет выполнять в параллели. Что именно происходит в первом случае с итерациями i?

Answer 1

Для того, что я вижу, вы забыли объявить i частным в первом случае. Поэтому i обновляется довольно случайным образом различными потоками, выполняющими соответствующий цикл, что приводит к отсутствию связанного доступа к массивам A и B.

Попробуйте добавить private(i) в свою директиву #pragma omp parallel firstprivate(my_amax) и посмотреть, что произойдет.