У меня есть следующий код на C, использующий OpenMP:
#pragma omp parallel default(none) private(i, j, k) shared(lb0, ub0, lb1, ub1, lb2, ub2, g, current, update, diff, dg)
#pragma omp for collapse(3)
for( i = lb0; i <= ub0; i++ ) {
for( j = lb1; j <= ub1; j++ ) {
for( k = lb2; k <= ub2; k++ ) {
g->data[ update ][ i ][ j ][ k ] =
ONE_SIXTH * ( g->data[ current ][ i + 1 ][ j ][ k ] +
g->data[ current ][ i - 1 ][ j ][ k ] +
g->data[ current ][ i ][ j + 1 ][ k ] +
g->data[ current ][ i ][ j - 1 ][ k ] +
g->data[ current ][ i ][ j ][ k + 1 ] +
g->data[ current ][ i ][ j ][ k - 1 ] );
diff = fabs( g->data[ update ][ i ][ j ][ k ] - g->data[ current ][ i ][ j ][ k ] );
dg = dg > diff ? dg : diff;
}
}
}
Насколько я знаю, я следую правилам использования предложения collapse с директивой #pragma omp для директивы, но когда я компилируюэто я получаю следующую ошибку:
grid.c: In function ‘grid_update’:
grid.c:202: error: expected ‘#pragma omp’ clause before ‘collapse’
make: *** [grid.o] Error 1
Что я здесь не так делаю?Я пытался поместить #pragma omp collapse
в качестве отдельной директивы, но это тоже не работает.