多维超大型数组
多维超大型数组
提问于 2018-11-27 17:49:23
嗨，我想使用多维非常大的数组。我试着跟踪代码。它会编译，但当我执行它时，它会给我分段错误。
'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.;
}
}
}'但是当我在下面的部分中工作时，我得到了分割错误，我的代码下一部分是包括二维数组的行。这个二维数组非常大，也许正因为如此，我得到了分割错误
'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;' 因此，工作和不工作之间的代码存在问题，导致分割错误。任何对此的解决方案。
回答 1
Stack Overflow用户
发布于 2018-11-29 12:29:53
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;
}
}
}
}页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接：
https://stackoverflow.com/questions/53496862
复制相关文章
相似问题
腾讯云开发者
