Question

短篇小说：尝试修改在堆上运行时分配的大型3D数组。我相信修改下面显示的数组vcross的函数正在创建不会被破坏的内存。

更长的故事： 我有一个大的3D双阵列（~126000x3x3或大约8.6MB），我需要运行一些操作。我不知道这个数组的第一个维度在编译时有多大，所以我使用new和delete操作将内存分配给堆。

当我尝试将值存储到此数组时，我会遇到分段违规。这让我觉得在将数值存储到数组时，我会在某处浪费创建内存，最终填满堆。

代码编译正常，但在我运行时会遇到seg违规。

static void inpolyh(
    double (*f)[3],//pointer to array[3], treated as 2D array where I don't know the first dimension until run-time.
    double (*v)[3],
    double (*p)[3],
    size_t numF,
    size_t numP)
{

    /*Calculate the baseNormals*/
    //allocate memory on the heap
    double (*baseNormals)[3] = NULL;//pointer to an array[3]
    if ( !(baseNormals = new double[numF][3]) ) { out_of_memory(); }

    //store the vector cross products in each array[3] of baseNormals
    for (int i=0; i<numF; i++) {
         vcross(baseNormals[i],
             v[(int)f[i][0]],
             v[(int)f[i][1]],
             v[(int)f[i][2]]);
         //THIS WORKS
    }

    /*Calculate face normals of tetrahedron*/
    //allocate memory on the heap (THIS WORKS)
    double (*faceNormals)[3][3] = NULL; //pointer to an array[3] of arrays[3]
    if ( !(faceNormals = new double[numP][3][3]) ) { out_of_memory(); }

    //store vector cross products into each array[3] of faceNormals
    for (int i=0; i<numP; i++ ) {
        for (int j=0; j<3; j++ ) {
            vcross(faceNormals[i][j],
                p[i],
                v[ (int) f[i][j] ],
                v[ (int) f[i][ (j + 1) % 3 ] ] );
            //SEG VIOLATION at i=37560
        }
    }

    delete [] baseNormals;
    delete [] faceNormals;
}

这就是我相信罪魁祸首的地方。我认为这个函数会创建一个永远不会被破坏的内存。向量交叉乘积函数接受四个数组[3]参数，并为第一个输入参数赋值，并通过引用传递。

static void vcross(
    double (&n)[3],
    double a[3],
    double b[3],
    double c[3])
{
    n[0] = b[1] * c[2] - a[1] * c[2] + a[1] * b[2] - b[2]
                * c[1] + a[2] * c[1] - a[2] * b[1];
    n[1] = b[2] * c[0] - a[2] * c[0] + a[2] * b[0] - b[0]
                * c[2] + a[0] * c[2] - a[0] * b[2];
    n[2] = b[0] * c[1] - a[0] * c[1] + a[0] * b[1] - b[1]
                * c[0] + a[1] * c[0] - a[1] * b[0];
    return;
}

其他可能重要的细节：

这是一个在matlab中运行的mex函数。
默认编码：windows-1252
MATLAB Root：C：\ Program Files \ MATLAB \ R2012b
MATLAB版本：8.0.0.783（R2012b）
操作系统：Microsoft Windows 7
处理器ID：x86 Family 6 Model 58 Stepping 9，GenuineIntel
虚拟机：Java 1.6.0_17-b04与Sun Microsystems Inc. Java HotSpot（TM）64位服务器VM混合模式
Window System：版本6.1（Build 7601：Service Pack 1）

Answer 1

您确定(int)f[i][j]始终在[0，＆lt;维度v＆gt;]范围内吗？我要做的第一件事是在循环运行时打印出(int)f[i][j]的值。或者只需触发调试器即可查看崩溃时的值。

修改大型动态大小的3D数组mex / C ++

1 个答案: