C 中的算术异常

问题描述

//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 文件中提取数字，它们的编码有点奇怪，所以。

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

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

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

解决方法

#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}};
    const float b[3]={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

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

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

arithmeticexception c debugging debugging debugging