在三个嵌套并行循环中使用OpenMP和O0在C代码中进行分段错误？

Question

我正在尝试使用以下代码并行化函数

        double*** calchr(int snr, int snf ,int snz, double*** shr,double*** shr_temp, double*** sef, double*** sez,double*** ser,double*** shf, double*** a1r_h,double*** a2r_h,double dr,double df,double dz,double dt, int singularity,int pmlthick)
    {
        int i,j,k,l;
        double temp,per_temp,sum1,sum2;
        int fi_pi_2,fi_3pi_2;

        for(i=0; i<=(snr-1); i++)
        {
            for(k=0; k<=(snz-1); k++)
                {
                    k+=pmlthick;
                    shr_temp[i][snf][k]=shr[i][snf][k];
                    #ifdef parfor
                    #pragma omp parallel for schedule(static) private(i, j, k, temp, fi_3pi_2, fi_pi_2, per_temp, sum1, sum2) shared(shr_temp, shr, a2r_h, sef, sez, a1r_h)  num_threads(THREADS) ordered
                    #endif                    
                    for(j=0; j<=(snf-1); j++)
                        {
                        temp=shr[i][j][k];
                        shr_temp[i][j][k]=temp;
                        if( i!=0)
                            shr[i][j][k]=(a2r_h[i][j][k]/dz)*(sef[i][j][k+1]-sef[i][j][k])-( (a2r_h[i][j][k]/(double(i)*dr*df))*(sez[i][j+1][k]-sez[i][j][k]) )+ a1r_h[i][j][k]*temp;
                        else
                        {  
                            {   
                            fi_pi_2=j+snf/4;
                            fi_3pi_2=j+(3*snf/4);
                            if(fi_pi_2>=snf)
                            fi_pi_2=fi_pi_2-snf;
                            if(fi_3pi_2>=snf)
                            fi_3pi_2=fi_3pi_2-snf;
                            shr[0][j][k] = temp + (0.5*dt/MI_0)*( ((sez[1][fi_3pi_2][k]-sez[1][fi_pi_2][k])/dr) + ((-1.0*ser[0][fi_3pi_2][k+1] + ser[0][fi_pi_2][k+1] - ser[0][fi_pi_2][k] + ser[0][fi_3pi_2][k])/dz));                         
                            }
                        }
                        shr[i][snf][k]=shr[i][0][k];
                        k-=pmlthick;
                        }
                }
        }
        return(shr,sef,sez);    
    }

num_threads（THREADS）作为带有适当标志-DTHREADS = 4的Makefile的参数（例如）

当我使用-O2标志编译代码时，不会出现分段错误。 当我使用-O0标志进行编译（无优化）时，出现分段错误。 该代码的串行版本在-O2和-O0优化级别上都可以正常运行。 所有运算都是加法和乘法，并且矩阵在这些函数之外动态分配。 是什么导致此问题？ 我已经读过了，在并行嵌套循环中，我必须将所有循环的索引声明为私有变量。 我做到了，但细分错误仍然存​​在。 在内部循环中，我放置了一条printf消息，并且看起来可以很好地运行到某个索引号。 我想知道是什么原因导致了这个问题。 任何想法将不胜感激