openMP nidificato parallelo per cicli vs parallelo interno per

Question

se uso annidati parallelo per cicli come questo:

#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 
} 

è questo 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 
} 

È il parallelo esterno per fare qualcosa di diverso creando un nuovo compito?

Answer 1

Se il compilatore supporta OpenMP 3.0, è possibile utilizzare la clausola 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 
} 

Se non (ad esempio solo OpenMP 2.5 è supportato), c'è una soluzione semplice:

01.235.

È possibile abilitare il parallelismo nidificato con omp_set_nested(1); e il codice nidificato omp parallel for funzionerà ma potrebbe non essere l'idea migliore.

A proposito, perché la pianificazione dinamica? Ogni iterazione del ciclo viene valutata in un tempo non costante?

Answer 2

NO.

Il primo #pragma omp parallel creerà un team di thread paralleli e il secondo tenterà quindi di creare per ciascuno dei thread originali un altro team, ovvero una squadra di team. Tuttavia, su quasi tutte le implementazioni esistenti, la seconda squadra ha solo un thread: la seconda regione parallela non viene sostanzialmente utilizzata. Così, il codice è più simile equivale 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 
    } 
} 

Se non si vuole questo, ma parallelizzare solo il ciclo interno, si può fare questo:

#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 
    } 
}