Comment créer un `omp parallèle for` avec synchronisation (` barrier`) de tous les fils dans le milieu avec OpenMP

Question

I ont deux fonctions, do_step_one(i) et do_step_two(i), pour i de 0 à N-1.

Actuellement, j'ai ce code (séquentiel):

for(unsigned int i=0; i<N; i++) { 
    do_step_one(i); 
} 

for(unsigned int i=0; i<N; i++) { 
    do_step_two(i); 
} 

Chaque appel de do_step_one() et do_step2() peut être fait dans l'ordre et en parallèle, mais toute do_step_two() a besoin de la fin de tous les do_step_one() pour commencer (il utilise do_step_one() résultats).

J'ai essayé les éléments suivants:

#omp parallel for 
for(unsigned int i=0; i<N; i++) { 
    do_step_one(i); 

#omp barrier 

    do_step_two(i); 
} 

Mais gcc se plaint

convolve_slices.c: 21: avertissement: zone barrière ne peut pas être étroitement imbriquées à l'intérieur du partage du travail, critique, ordonné, master ou zone de tâche explicite.

Qu'est-ce que je ne comprends pas? Comment résoudre ce problème?

Answer 1

Un problème que je vois avec ce code, est que le code ne respecte pas la spécification :)

Si vous avez besoin tous do_step_one() est à la fin, vous aurez besoin de quelque chose comme ce qui suit:

#pragma omp parallel for 
for(unsigned int i=0; i<N; i++){ 
    do_step_one(i); 
} 

#pragma omp parallel for 
for(unsigned int i=0; i<N; i++){ 
    do_step_two(i); 
} 

Le résultat de ceci serait un parallélisme du premier pour, puis un parallélisme du second pour.

Answer 2

Juste une note de côté, si vous voulez vous assurer que les fils ne sont pas recréées, séparer la déclaration de parallèle et la déclaration de pour:

#pragma omp parallel 
{ 
    #pragma omp for 
    for(unsigned int i=0; i<N; i++){ 
    do_step_one(i); 
    } 
    //implicit barrier here 
    #pragma omp for 
    for(unsigned int i=0; i<N; i++){ 
    do_step_two(i); 
    } 
}