Fortran Coder

查看: 72|回复: 2

[编译器] 求助fortran 高精度计算问题

[复制链接]

2

帖子

1

主题

0

精华

新人

F 币
19 元
贡献
8 点
发表于 2020-9-5 08:22:23 | 显示全部楼层 |阅读模式
大家好,向大家咨询一个关于fortran高精度的计算问题。在我的计算中,经常遇到双精度的复杂计算,随着计算时间变长,结果就不正确了,不知道为何,
比如:
real(8) :: a=0.12345678988
real(8) b
b=a+a
write(*,*) b

理论上计算b的值应该为:0.24691357976
计算机输出的结果是:     0.246913582086563

请问,如何修改能保证计算结果。
谢谢。


回复

使用道具 举报

57

帖子

2

主题

0

精华

专家

F 币
391 元
贡献
170 点

规矩勋章

发表于 2020-9-5 12:55:11 | 显示全部楼层
尽管变量a和b定义为64位(双精度),但是赋值的常数只有32位(单精度),改为如下:
real(8) :: a=0.12345678988_8

167

帖子

14

主题

0

精华

大宗师

F 币
4021 元
贡献
823 点
发表于 2020-9-8 10:44:00 | 显示全部楼层
风平老涡 发表于 2020-9-5 12:55
尽管变量a和b定义为64位(双精度),但是赋值的常数只有32位(单精度),改为如下:
real(8) :: a=0.123456 ...

感觉更常见易于理解的写法是 0.12345678988d0
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

QQ|捐赠本站|Archiver|关于我们 About Us|群聊|Fcode

GMT+8, 2020-9-27 08:08

Powered by Discuz! X3.2

© 2001-2017 Comsenz Inc.

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