Bonjour,

Le code suivant ne donne pas le m�me r�sultat si je d�sative les directives openMP. Alors, la diff�rence est mineures, 10-8, mais cela suffit � donner un comportement diff�rent dans le reste du code. Est-ce normal compte tenu de la parall�lisation ou cela indique-t-il un probl�me. Merci d'avance de vos lumi�res.

Code : S�lectionner tout - Visualiser dans une fen�tre � part
1
2
3
4
5
6
7
8
9
10
11
12
13
14
 
double X,Y,XP,YP,P,E;
X=Y=Xp=Yp=P=E=0;
int npoint0=1000000
int NN=(npoint0-1)*6;
#pragma omp parallel for  schedule(guided) reduction(+:X,Y,Xp,Yp,P,E)
for (int i=0;i<=NN;i=i+6) {
   X+=cord[i];
   Xp+=cord[i+1];
   Y+=cord[i+2];
   Yp+=cord[i+3];
   P+=cord[i+4];
   E+=cord[i+5];
}