Fortran Coder

查看: 18697|回复: 9

[数值问题] 请问迭代结果都是NaN是啥原因

[复制链接]

11

帖子

3

主题

0

精华

入门

F 币
44 元
贡献
41 点
发表于 2021-3-8 22:00:28 | 显示全部楼层 |阅读模式
有没有大佬帮看看,我写这代码最终输出结果如图所示,都是NaN,这是啥原因,新手上路,请多指教

输出结果

输出结果

陈亮3.8.1.f90

14.45 KB, 下载次数: 2

250

帖子

2

主题

0

精华

宗师

F 币
1730 元
贡献
872 点

规矩勋章

发表于 2021-3-8 23:03:58 | 显示全部楼层
pshr2=0,导致
flowep(1)=(sqrt(2.0)/3.0)*(1/sqrt(pshr2))*(2*pstrain(1)-pstrain(2)-pstrain(3))
这里有问题

11

帖子

3

主题

0

精华

入门

F 币
44 元
贡献
41 点
 楼主| 发表于 2021-3-9 10:50:40 | 显示全部楼层
necrohan 发表于 2021-3-8 23:03
pshr2=0,导致
flowep(1)=(sqrt(2.0)/3.0)*(1/sqrt(pshr2))*(2*pstrain(1)-pstrain(2)-pstrain(3))
这里有 ...

谢谢大佬,我先去看看

11

帖子

3

主题

0

精华

入门

F 币
44 元
贡献
41 点
 楼主| 发表于 2021-3-11 10:38:09 | 显示全部楼层
necrohan 发表于 2021-3-8 23:03
pshr2=0,导致
flowep(1)=(sqrt(2.0)/3.0)*(1/sqrt(pshr2))*(2*pstrain(1)-pstrain(2)-pstrain(3))
这里有 ...

大佬,我把这个公式换了之后,这个循环能进行一部分了,但是后面还是出现NaN,可能又有公式出现了问题,最近我看了下单步调试,但没咋看懂,请问这个断点的设置有什么规则么,怎么知道应该在哪设置断点?
(8]9~B}OIIPAX`82R4FDL87.png

陈亮程序3.10.1.f90

13.89 KB, 下载次数: 1

18

帖子

2

主题

0

精华

入门

F 币
80 元
贡献
50 点
发表于 2021-3-11 11:13:01 | 显示全部楼层
中间计算的数太大了,超过了显示范围,我之前计算的结果有出现过这样的情况

11

帖子

3

主题

0

精华

入门

F 币
44 元
贡献
41 点
 楼主| 发表于 2021-3-11 12:03:00 | 显示全部楼层
kellsi 发表于 2021-3-11 11:13
中间计算的数太大了,超过了显示范围,我之前计算的结果有出现过这样的情况 ...

可能逻辑出现了点问题 这个本来都是小数点后几位的数 突然变成了个位数,大哥你之前咋解决的

18

帖子

2

主题

0

精华

入门

F 币
80 元
贡献
50 点
发表于 2021-3-11 15:45:19 | 显示全部楼层
旅人 发表于 2021-3-11 12:03
可能逻辑出现了点问题 这个本来都是小数点后几位的数 突然变成了个位数,大哥你之前咋解决的 ...

如果是做案例的话,先把带入计算值设置成一个合理的小值去验算整个程序的代码对不对,之后再调整这个数值到正常值

11

帖子

3

主题

0

精华

入门

F 币
44 元
贡献
41 点
 楼主| 发表于 2021-3-11 18:45:13 | 显示全部楼层
kellsi 发表于 2021-3-11 15:45
如果是做案例的话,先把带入计算值设置成一个合理的小值去验算整个程序的代码对不对,之后再调整这个数值 ...

这是按一篇论文里面的计算流程来搞的,找不到论文原来的代码,只有论文里的图来参考一下最终的数据结果,我太难了

250

帖子

2

主题

0

精华

宗师

F 币
1730 元
贡献
872 点

规矩勋章

发表于 2021-3-11 20:47:29 | 显示全部楼层
旅人 发表于 2021-3-11 18:45
这是按一篇论文里面的计算流程来搞的,找不到论文原来的代码,只有论文里的图来参考一下最终的数据结果, ...

逻辑问题,不是专业的,爱莫能助

5

帖子

1

主题

0

精华

入门

F 币
42 元
贡献
21 点
发表于 2021-8-30 16:59:36 | 显示全部楼层
请问最后解决了吗?同样遇到了这种情况
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

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

GMT+8, 2024-3-29 17:04

Powered by Tencent X3.4

© 2013-2024 Tencent

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