並行delaunay三角測量樸素算法

Question

以下代碼（Pg.187，由Rourke編寫的C中的計算幾何）需要相同的時間以串行和並行方式運行（2 proc）。 請幫我找出問題所在。 這裏的平行部分

int chunk; 
chunk=10; 
#pragma omp parallel private(i,j,k,xn,yn,zn) 
{ 
#pragma omp for schedule(static,chunk) 
    for(i=0;i<n-2;i++) 
    { 
     for(j=i+1;j<n;j++) 
      for(k=i+1;k<n;k++) 
       if(j!=k) 
       { 
        xn=(y[j]-y[i])*(z[k]-z[i])-(y[k]-y[i])*(z[j]-z[i]); 
        yn=(x[k]-x[i])*(z[j]-z[i])-(x[j]-x[i])*(z[k]-z[i]); 
        zn=(x[j]-x[i])*(y[k]-y[i])-(x[k]-x[i])*(y[j]-y[i]); 
        if(flag=(zn<0)) 
         for(m=0;m<n;m++) 
          flag=flag && ((x[m]-x[i])*xn + (y[m]-y[i])*yn + (z[m]-z[i])*zn <=0); 
        if (flag) 
         printf("%d\t%d\t%d\n",i,j,k); 
       } 
    } 
} 

Answer 1

好的解決了它。 只需要讓剩下的變量-m和標誌private.Wow愚蠢的我!!!!!