c - 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 文件中获取数字,它们的编码有点奇怪,所以。
否则我将提供所有其他信息,所以请不要关闭这个问题,我真的很想知道为什么会发生这种情况以及如何解决它。
这也是异常的另一个例子。 例子
我想知道的是,即使在这个例子中没有除以 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,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
所以问题不是由代码的结构引起的。它肯定来自一个特殊的算术问题。
推荐阅读
- spring-boot - OAuth2AuthenticationToken 与 OAuth2Authentication
- c# - 通过调用“IServiceCollection.AddHealthChecks”添加所有必需的服务
- java - 我想知道 system.out.println() 在这段代码中做了什么
- django - Django LiveServerTestCase 似乎在测试之间重新加载夹具
- javascript - Chrome 扩展顶部访问列表
- nginx - 通过查看 nginx 错误日志,我的服务器是否受到攻击?
- node.js - 如何使用 express 将请求响应从一条路由传递到另一条路由
- mongoid - Omniauth 会话在身份验证期间丢失
- azure-cognitive-search - 执行 indexClient.Documents.Index(batch) 时 index 方法失败
- python - Keras 支持对函数求导