如何知道两种方法是否对相同数据的不同表示做同样的事情

Question

我有两种不同的矩阵表示。

矩阵保存为数组 float*:

| 1.0   2.0   3.0   4.0   5.0   6.0  7.0   8.0   9.0 |

作为数组的数组 float**

| 1.0  2.0  3.0 |
| 4.0  5.0  6.0 |
| 7.0  8.0  9.0 |

我有这种方法可以将矩阵用作数组：

float a(float* x, float* b, int m, int n) {
    int i, j;
    float e = 0, de, yp;

    for (i = 0; i < m; i++) {
        yp = 0;
        for (j = 0; j < n; j++)
            yp += x[i*(n+1)+j] * b[j];
        de = fabs(x[i*(n+1)+n]-yp);
        e += (de*de);
    }
    return e;
}

在我看来，这个应该做同样的事情：

float a(float **x, float **b, int m, int n) {
    int i, j;
    float e = 0, de, yp;

    for (i = 0; i < m; i++) {
        yp = 0;
        for (j = 0; j < n; j++)
            yp += x[i][j] * b[j][0];
        de = fabs(x[i][j]-yp);
        e += (de*de);
    }
    return e;
}

有没有一种聪明的方法来比较这两种方法，而无需尝试编写相同输入的两种不同表示并检查两者的输出是否相同？

注意： fis 方法中的de = fabs(x[i*(n+1)+n]-yp); 行似乎是一个错误，因为我看不到这条线应该做什么......我不是该方法的单维版本的作者，所以我不知道真实的作者的意图，我正在尝试调整代码以处理多维矩阵，因此我将此行解释为de = fabs(x[i*(n+1)+j]-yp);，但我不确定。

Answer 1

即使您已正确分配内存，您也没有按照您想要的方式索引第一个数组。看看这些行：

for (i = 0; i < m; i++) {
    yp = 0;
    for (j = 0; j < n; j++)
        yp += x[i*(n+1)+j] * b[j];

你希望这是一样的

  for (i = 0; i < m; i++) {
    yp = 0;
    for (j = 0; j < n-1; j++)
        yp += x[i][j] * b[j][0];

但是当你有一个一维数组时，如果每行有 n 个元素，则相当于索引一个二维数组 x[i][j] 是用

 yp += x[i*n + j];

没有+1...