我正在处理10个数组,其中一些是双1024x1392。
我试着用在堆上动态分配它们
double **x_array;
x_array = new double*[NUM_ROWS];
for(int i=0; i < NUM_ROWS; i++) {
x_array[i] = new double[NUM_COLS];
}
for(int ix=0; ix < NUM_COLS; ix++) {
for(int iy=0; iy < NUM_ROWS; iy++) {
x_array[ix][iy]=(x1y1*(ix+1) + x2y1*(iy+1) + x3y1);
//y_array[ix][iy]=(x1y2*(ix+1) + x2y2*(iy+1) + x3y2);
}
}
}
但我仍然会在说时出错
unhandled exception: System.Runtime.InteropServices.SEGException: External Component has thrown an exception. at line 106
106是我开始初始化上面代码中的数组的地方:
x_array = new double*[NUM_ROWS];
我真的没有空间了,还是我做错了什么?
您已将数组索引转置:
for(int ix=0; ix < NUM_COLS; ix++) {
for(int iy=0; iy < NUM_ROWS; iy++) {
x_array[ix][iy]=(x1y1*(ix+1) + x2y1*(iy+1) + x3y1);
应该是:
for(int iy=0; iy < NUM_ROWS; iy++) {
for(int ix=0; ix < NUM_COLS; ix++) {
x_array[iy][ix]=(x1y1*(ix+1) + x2y1*(iy+1) + x3y1);
或者如果你真的因为某种原因不得不保持缓存的敌对循环顺序:
for(int ix=0; ix < NUM_COLS; ix++) {
for(int iy=0; iy < NUM_ROWS; iy++) {
x_array[iy][ix]=(x1y1*(ix+1) + x2y1*(iy+1) + x3y1);