C++コード

Question

を並列化するOPENMPを使用して、私は以下のようにOpenMPを使用して、私のC++コードを並列化肝炎：

....//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 

しかし、結果は、私はそれを実行するたびに変わります！この行を削除したとき

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

の結果は、異なる実行で同じです！ 私のコードで何が問題になっていますか？

Answer 1

変数Sはループ外で宣言され、ループを実行するすべてのスレッド間で共有されます。どのように使用されているかを考えれば、ループの内側ではdouble S = rdote/r_ij_square;と宣言できます。