1

我目前正在寻找两个矩阵的平方距离之和，数据保存在 double* 数组中。它们中的第一个保持不变，而另一个通过使用返回两个索引之间的 32x32 数组的函数循环。

但是，当我尝试在“e”的第一次递增之后调用“getTile(d,e)”时，它会引发堆损坏异常：

double* Matrix::ssd(int i, int j, Matrix& rhs){
double sum = 0, val = 0; int g = 0, h=0;
double* bestMatch = new double[32*32]; double* sameTile = new double[32*32];    double* changeTile = new double[32*32]; 

for(int x = i-32; x <i; x++){
    for(int y = j-32; y <j; y++){
        sameTile[g*32+h] = data[x*N+y];
        h++;
    }g++; h = 0;
}

system("pause");

for(int d = 32; d<=512; d+=32){
    for(int e = 32; e<=512; e+=32){

        changeTile = rhs.getTile(d,e);

        for(int out = 0; out < 32; out++){
            for(int in = 0; in < 32; in++){
                val = sameTile[out*32+in] - changeTile[out*32+in];
                val = val*val;
                sum = sum + val;

            }       
        }
        cout << sum << endl;

        sum = 0; val = 0;

        system("pause");
    }   
}

getTile(int i, int j) 函数：

double* Matrix::getTile(int i, int j){
double* tile = new double[32*32]; int g = 0; int h = 0;
for(int x=i-32; x<i; x++){
    for(int y=j-32; y<j; y++){
        tile[g*32+h] = data[x*N+y];
        h++;
    }
    cout << endl;
    g++;
}
return tile;
}

我相信在 changeTile double* 中分配内存时会发生错误？

任何帮助将不胜感激。

4

1 回答 1

3

您的代码中有很多问题都与不正确地访问数组元素有关。

在第一个循环中：

 sameTile[g*32+h] = data[x*N+y];

至少使数据数组下溢。考虑如果 i=0、j=0 和 N=512，那么您尝试在循环的第一遍中访问 data[-16416]。

第二个问题是 getTile() 方法，您忘记在内部循环结束时重置h为 0（就像您在 ssd() 方法中所做的那样）。这会导致 tile[] 溢出

我还会仔细检查该行：

 changeTile = rhs.getTile(d, e);

和方法 getTile() 确保在 data[] 上不会发生数组溢出。

总的来说，我建议尽可能使用适当的 std:: 容器。正确使用它们应该可以完全消除此类错误。如果您确实需要使用原始指针/数组，那么除了在需要的地方检查边界之外，您还需要确保对它们的所有索引都尽可能清晰。

于 2012-04-25T02:14:19.840 回答