Question

嗨，我想使用多维非常大的数组。我尝试了以下代码。它可以编译，但是当我执行它时，会给我分段错误。

'int NT = 35; int NX = 25; int NY = 25; int NZ = 25;
double dt = 0.1; double dx = 0.5; double dy = 0.5; double dz = 0.5;

double PosT[NT]; 
double PosX[NX]; double PosY[NY]; double PosZ[NZ];

for(int i=0;i<NT;i++)
    PosT[i] = i*dt+dt;
for(int i=0; i<NX;i++)
    PosX[i] = dx*i;
for(int i=0; i<NY;i++)
    PosY[i] = dy*i;
for(int i=0; i<NZ;i++)
    PosZ[i] = dz*i;

double* b_x=(double*)malloc(NX*NY*NZ*sizeof(double));
double* b_y=(double*)malloc(NX*NY*NZ*sizeof(double));
double** B=(double**)malloc(NX*NY*NZ*NT*sizeof(double*));

if(b_x==NULL||b_y==NULL){
    cout<<"Malloc space error!"<<endl;
    return 0;
}
    for(int ix=0;ix<NX;ix++){
        for(int iy=0;iy<NY;iy++){
            for(int iz=0;iz<NZ;iz++){
                int position=ix*NY*NZ+iy*NZ+iz; 
                b_x[position] =0.;
                b_y[position] =0.;
            }
        }
    }'

但是当我在下面的部分工作时，我遇到了分割错误，我的代码下一部分是包含2d数组的以下行。这个二维数组非常大，也许由于这个原因，我遇到了细分错误

    'if(B==NULL){
    cout<<"Malloc space error!"<<endl;
    return 0;
    }
    cout<<"work"<<endl;
        for(int ix=0;ix<NX;ix++){
        for(int iy=0;iy<NY;iy++){
            for(int iz=0;iz<NZ;iz++){
                int position=ix*NY*NZ+iy*NZ+iz;
                for(int it=0;it<NT;it++){               
                B[position][it]=0.;
            }
        }
    }
}
 cout<<"not working"<<endl;'

因此，工作与不工作之间的代码存在问题，导致分段错误。任何解决方案。

Answer 1

int NT = 35; int NX = 25; int NY = 25; int NZ = 25;

为简单起见，让我们将所有这些都更改为NT=NX=NY=NZ=2。这行：

double** B=(double**)malloc(NX*NY*NZ*NT*sizeof(double*));

然后将为16个指针分配空间。在循环的第一次迭代中，此行：

B[position][it]=0.;

将等同于：

double *tmp = B[0];  // Load uninitialized pointer from B[0]
tmp[0] = 0.0;        // Dereference uninitialized pointer to store something.

此代码生成SIGSEGV并不奇怪。

您可能的意思是：

double *B = malloc(NX*NY*NZ*NT*sizeof(double));
for(int ix = 0; ix < NX; ix++) {
  for(int iy = 0; iy < NY; iy++) {
    for(int iz = 0; iz < NZ; iz++) {
       for(int it = 0; it < NT; it++) {
         int position = NT * (NZ * (NY * ix + iy) + iz) + it;
         B[position] = 0.0;
       }
    }
  }
}

多维非常大的数组

1 个答案: