OpenMP et std::vector ?

Bonjour.

Je viens juste de me mettre tout doucement à la programmation parallèle avec OpenMP et j'aurai quelques questions. En fait j'ai un programme pour l'instant monotâche que je cherche à paralléliser. Enfin plutôt je cherche à paralléliser 1 boucle pour le moment.

Schématiquement, j'ai un truc dans le genre :

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17

int main() { std::vector<unsigned int> resultat;   // // Bout de programme et données qui doivent rester "monotache" //   for(i=0;i<N;i++) { algorithme(i,N,resultat) }   // // Bout de programme et données qui doivent rester "monotache" // }

---

1) Déjà, en ce qui concernent les parties qui doivent rester "monotache" est-ce que j'ai besoin de changer quelque chose du point de vue d'OpenMP ou je n'y touche pas ?

2) Ensuite, au niveau de la boucle for, chaque itération est indépendante ... sauf que la fonction "algorithme(unsigned int i,unsigned int N,std::vector<unsigned int> & resultat)" modifie le vecteur "resultat" en faisant des push_back...

L'ordre dans lequel les éléments sont placés dans le vector importe peut puisqu'ensuite de toutes façon je fais un "sort". Mais je me dis que si deux thread veulent faire un push_back en même temps, cela risque de poser problème....

Bref, ça doit être un truc de base, mais comment s'en sortir ?

Merci :zoubi:


EDIT : je viens d'avoir une idée en lisant un tuto sur OpenMP et à mon avis c'est comme cela qu'il faut procéder (par contre je ne sais pas trop comment l'implémenter) :
- chaque thread travaille sur un std::vector privé
- une fois le travail effectué, on ajoute "séquentiellemment" (avec un "ordered" à mon avis) le contenu de chaque std::vector privé dans un grand std::vector global.

Salut,
Dans ton cas le plus simple et de decalre un tableau de vecteurs, en fait un par thread. Ainsi tu pourra paralléliser ta boucle sans souci. Ensuite à la fin tu concatène l'ensemble.
Il y a une fonction openMP qui permet de savoir combien de threads seront crée avant de rentrer dans ta boucle, cette info te permettra de créer ton tableau de vecteur. Je crois que c'est omp_get_num_thread(), enfin regarde la doc de openMP

Salut

Citation:

---

Envoyé par Kaluza

EDIT : je viens d'avoir une idée en lisant un tuto sur OpenMP et à mon avis c'est comme cela qu'il faut procéder (par contre je ne sais pas trop comment l'implémenter) :
- chaque thread travaille sur un std::vector privé
- une fois le travail effectué, on ajoute "séquentiellemment" (avec un "ordered" à mon avis) le contenu de chaque std::vector privé dans un grand std::vector global.

---

normalement tu peux faire.

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31

int main() { std::vector<unsigned int> resultat;   // // Bout de programme et données qui doivent rester "monotache" //     #pragma omp parrallel {   //variable locale à chaque thread std::vector<unsigned int> threadRes;   #pragma omp for for(i=0;i<N;i++) { algorithme(i,N,threadRes) }   #pragma omp critical { //merge les threadRes dans resultat. } }   // // Bout de programme et données qui doivent rester "monotache" // }

---

Une autre solution consisterait à faire un resize() de ton vector avant la section parallèle, et dans ta boucle d'y ajouter des éléments non pas avec des push_back() mais à la façon des tableaux (avec l'opérateur [] )