OpenMP paralelo bucles for anidados vs paralelo interno para

Question

si uso paralelo bucles for anidados como esto:

#pragma omp parallel for schedule(dynamic,1) 
for (int x = 0; x < x_max; ++x) { 
    #pragma omp parallel for schedule(dynamic,1) 
    for (int y = 0; y < y_max; ++y) { 
    //parallelize this code here 
    } 
//IMPORTANT: no code in here 
} 

es esto equivale a:

for (int x = 0; x < x_max; ++x) { 
    #pragma omp parallel for schedule(dynamic,1) 
    for (int y = 0; y < y_max; ++y) { 
    //parallelize this code here 
    } 
//IMPORTANT: no code in here 
} 

es el paralelo externa para hacer otra cosa que creando una nueva tarea?

Answer 1

Si su compilador soporta OpenMP 3.0, puede utilizar la cláusula de collapse:

#pragma omp parallel for schedule(dynamic,1) collapse(2) 
for (int x = 0; x < x_max; ++x) { 
    for (int y = 0; y < y_max; ++y) { 
    //parallelize this code here 
    } 
//IMPORTANT: no code in here 
} 

Si no (por ejemplo, sólo OpenMP 2.5 es compatible), hay una solución sencilla:

#pragma omp parallel for schedule(dynamic,1) 
for (int xy = 0; xy < x_max*y_max; ++xy) { 
    int x = xy/y_max; 
    int y = xy % y_max; 
    //parallelize this code here 
} 

Puede habilitar el paralelismo anidado con omp_set_nested(1); y su código anidado omp parallel for funcionará, pero esa podría no ser la mejor idea.

Por cierto, ¿por qué la programación dinámica? ¿Se evalúa cada iteración de bucle en tiempo no constante?

Answer 2

NO.

El primer #pragma omp parallel creará un equipo de subprocesos paralelos y el segundo intentará crear para cada uno de los subprocesos originales otro equipo, es decir, un equipo de equipos. Sin embargo, en casi todas las implementaciones existentes, el segundo equipo tiene solo un hilo: la segunda región paralela esencialmente no se usa. Por lo tanto, su código es más como equivalente a

#pragma omp parallel for schedule(dynamic,1) 
for (int x = 0; x < x_max; ++x) { 
    // only one x per thread 
    for (int y = 0; y < y_max; ++y) { 
     // code here: each thread loops all y 
    } 
} 

Si no quiere eso, pero sólo paralelizar el bucle interno, se puede hacer esto:

#pragma omp parallel 
for (int x = 0; x < x_max; ++x) { 
    // each thread loops over all x 
#pragma omp for schedule(dynamic,1) 
    for (int y = 0; y < y_max; ++y) { 
     // code here, only one y per thread 
    } 
}