Boucle Do parallele OpenMP

Bonjour,

J'ai un programme qui fonctionne bien avec un seul processeur.
Cependant, son temps d'execution est asse long alors j'aimerais utiliser les 2 processeurs de mon IBM. J'ai mis des instructions !$OMP DO / !$OMP END DO, cela fonctionne, les 2 cores tournent a 100% mais j'ai un probleme.
Voici un code pour illustrer mes propos :

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23

Program test integer i, compte real f,a(1024)   compte=0 !$OMP PARALLELE !$OMP DO do i=1, 1024 a(compte+1)=f(i) compte=compte+1 enddo !$OMP END DO !$OMP END PARALLELE   end program   Function f(nb) integer nb real f   f=(nb)**(1/2) return end function

---

Mon probleme est que pour certaines valeurs, dans le tableau a() ne sont pas remplis (a la limite entre les taches gerées par chacun des procs?) et pourtant la valeur final de compt est bien 1024.

Pouriez vous m'aider a resoudre ce probleme, avez vous des idées?

Merci

Bonjour,

une question : que ce passe-t-il si au lieu d'utiliser compte pour indicer ton tableau tu utilises i ?

Bonjour,

Il ne se passe rien, le problème reste le meme.
Dans mon tableau final, là où j'ai un problème, une valeur (environ) de -6E28 apparait, la meme valeur qui est dans chacune des lignes du tableau avant que je ne l'utilise.
Cette valeur de provient pas de la fonction qui retournerait cette valeur car j'ai mis un "if" pour afficher un "Erreur" à l'ecran si la fonction dans la boucle do renvoit une valeur négative.
Le problème vient d'un mauvais partage du travail entre les 2 cores.

Merci de t'être intéressé à mon problème.

Quel pourrait etre la cause de ces erreurs?

Merci

Bonsoir,

j'ai repris ton ecm (très bonne idée ;-)),

en modifiant, un peu ton code, j'obtiens :

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

Program test integer i real f,a(1024) !$OMP PARALLEL !$OMP DO do i=1, 1024 a(i)=f(i) enddo !$OMP END DO !$OMP END PARALLEL do i=1, 1024 write(*,*)i," ",a(i) enddo end program Function f(nb) integer nb real f f=(nb)**(1.0/2.0) return end function

---

Que je compile comme ça :

Code:

---

1 2	ifort -static -openmp Test.f90 -o test

---

et que je test comme ça :

Code:

---

1 2	./test | gawk '{print $1 "==" $2*$2}'

---

Bon du coup, il semble que pour l'ECM ça devrait marché (a priori si mon OS arrive bien à utiliser les deux cores)
=> Si cet ECM ne marche pas chez toi regarde du coté de la compilation
=> Sinon , tes problèmes venait peut-être d'autre choses (parties en rouge dans mon ecm)

Salut
Merci pour ta reponse. Et desole pour le retard.

En fait mon bout de code etait la pour illustrer mais ce n'est pas le code que je cherche a paralleliser.
Dans mon 'vrai' code (que je ne met pas parce qu'il est indigeste) je cherchais a paralleliser plusieur DO imbriques les uns dans les autres...
C'est ce que OMP a du mal a gerer (dans mon cas, peut etre existe-il une solution).
Donc je me suis debrouille pour reduire les 4 DO en un seul et ca semble fonctionner.
Merci pour ton aide.

(Tu compiles sous Linux? quels sont les avantages?)