c - C中的算术异常

Question

算术异常

//m*n行列Aを用いてy = A*x +b を計算する
void fc(int m, int n, const float *x, const float *A, const float *b, float *y){
    int i, j;

    for (i = 0; i < m; i++){
        y[i] = b[i];
        for (j = 0; j < n; j++){
            y[i] += A[i * n + j] * x[j];
        }
    }
}

这是对矩阵进行 AX+b 计算的代码。但如照片所示，发生了算术异常。为什么会这样？即使它是乘法并且没有除以0。我该如何解决这个错误？

抱歉，我无法添加这些值，否则我将不得不在此处添加整个文件。这些是神经网络的参数，我必须在此处添加 .dat 文件，然后我还需要其他可以加载这些文件的代码。另外，我不知道如何只从 .dat 文件中获取数字，它们的编码有点奇怪，所以。

否则我将提供所有其他信息，所以请不要关闭这个问题，我真的很想知道为什么会发生这种情况以及如何解决它。

这也是异常的另一个例子。 例子

我想知道的是，即使在这个例子中没有除以 0，这怎么可能发生。我该如何解释这种情况。

Answer 1

#include <stdio.h>
void fc(int m, int n, const float *x, const float *A, const float *b, float *y){
    int i, j;

    for (i = 0; i < m; i++){
        y[i] = b[i];
        for (j = 0; j < n; j++){
            y[i] += A[i * n + j] * x[j];
        }
    }
}

int main()
{
    const float x[3]={1,1,1};  
    const float *xp=x;
    const float A[3][3]={{1,1,1},{1,1,1},{1,1,1}};
    const float b[3]={1,1,1};   
    const float *bp=b;
    float y[3]; float *yp=y;
    fc (3,3,xp,*A,bp,yp);
    printf("%f %f %f ",y[0],y[1],y[2]);

    return 0;
}

我已经用所有变量的虚值 1 和 3x3 和 3x1 的矩阵大小测试了程序。结果是正确的，没有错误。结果是

4.0000000 4.0000000 4.0000000 

所以问题不是由代码的结构引起的。它肯定来自一个特殊的算术问题。

Answer 2

根据您的图像，您的矩阵大小为 100x50 (m,n)，这意味着 5000 个项目。但是您输入A[j*m+i]了“j”等于 55 而“i”等于 0 的位置。这意味着访问不允许的数组的 5500 项。