Cómo paralelizar correctamente un bucles for anidados

Question

que estoy trabajando con OpenMP para paralelizar un escalar anidado bucle for:

double P[N][N]; 
double x=0.0,y=0.0; 

for (int i=0; i<N; i++) 
{ 
    for (int j=0; j<N; j++) 
    { 
     P[i][j]=someLongFunction(x,y); 
     y+=1; 
    } 
    x+=1; 
} 

En este bucle, lo importante es que la matriz P debe ser la misma en ambos escalares y versiones paralelas:

todos mis posibles ensayos no tuvieron éxito ...

Answer 1

El problema aquí es que se ha añadido dependencias de iteración a iteración con:

x+=1; 
y+=1; 

Por lo tanto, como el código está parado ahora, no es paralelizable. Intentar hacerlo dará como resultado resultados incorrectos. (Como es probable que esté viendo)

Afortunadamente, en su caso, se puede calcular directamente sin la introducción de esta dependencia:

for (int i=0; i<N; i++) 
{ 
    for (int j=0; j<N; j++) 
    { 
     P[i][j]=someLongFunction((double)i, (double)N*i + j); 
    } 
} 

Ahora usted puede intentar tirar un pragma OpenMP sobre esto y ver si funciona:

#pragma omp parallel for 
for (int i=0; i<N; i++) 
{ 
    for (int j=0; j<N; j++) 
    { 
     P[i][j]=someLongFunction((double)i, (double)N*i + j); 
    } 
}