Fortran Coder

查看: 7818|回复: 2
打印 上一主题 下一主题

[数值问题] 使用矩阵求逆程序出现NaN

[复制链接]

954

帖子

0

主题

0

精华

大师

F 币
184 元
贡献
75 点

规矩勋章元老勋章新人勋章水王勋章热心勋章

QQ
楼主
发表于 2021-1-6 14:28:23 | 显示全部楼层
本帖最后由 vvt 于 2021-1-6 14:30 编辑

mateq函数中,第147,148 行。
记住,列主元是找绝对值大的,而不是数值大的。

本例中,你遇到了 -234,-5 和 0。然后你就找到 0 的那列了。然后算比例因子 temp=Ab(i,k)/Ab(k,k),得,除数为 0 ,结果为 NaN

147      if (dabs(Ab(i,k))>elmax) then
148         elmax=Ab(i,k)
改为
147      if (dabs(Ab(i,k))>elmax) then
148         elmax=abs(Ab(i,k))


评分

参与人数 1F 币 +3 收起 理由
fcode + 3

查看全部评分

您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

捐赠本站|Archiver|关于我们 About Us|小黑屋|Fcode ( 京ICP备18005632-2号 )

GMT+8, 2024-5-15 16:47

Powered by Tencent X3.4

© 2013-2024 Tencent

快速回复 返回顶部 返回列表