OpenMP mit einer C++ Code

Question

I hav parallelisiert mein C++ Code openmp wie unten parallelisieren:

....//some code 
    double S_theta = 0.0, S_x = 0.0, S_y = 0.0; 
#pragma omp parallel for private(dx,dy,theta_new) reduction(+ : S_x,S_y,S_theta) 
    for(int j = 0; j < N; j++) 
     { 
     if (j==i) continue; 
      double delta_x = x[i]-x[j], 

     for(int ky = -1; ky<= 1; ky++) 
      { 
      for(int kx = -1; kx<= 1; kx++) 
       { 

       if (r_ij_square > l0_two) 
        { 
        double r_ij = sqrt(r_ij_square); 
        r_x_ij/= r_ij; 
        r_y_ij/= r_ij; 
        double rdote = r_x_ij * e_x[j] + r_y_ij * e_y[j]; 
        S_theta += e_dot_e * (e_cross_e - rdote * (e_x[i] * r_y_ij - e_y[i] * r_x_ij))/(r_ij_square * r_ij); 
        double S = rdote /r_ij_square; 
        S_x += r_x_ij * S; 
        S_y += r_y_ij * S; 
        } 
       } 
     } 

    double zeta = -eta/2.0 ; 
    theta_new[i] += zeta + dt * dlthet; 
    } 
    ....//some code 

Aber das Ergebnis ändert sich jedes Mal, wenn ich es laufen! Wenn ich diese Linie

#pragma omp parallel for reduction(+ : S_x,S_y,S_theta) 

entferne, ist das Ergebnis das gleiche in verschiedenen Läufen! Was ist falsch an meinem Code?

Answer 1

Ihre Variable S wird außerhalb der Schleife deklariert und von allen Threads geteilt, die Ihre Schleife ausführen. Wenn es verwendet wird, kann es innerhalb der Schleife, in der es verwendet wird, als double S = rdote/r_ij_square; deklariert werden.