Fortran Coder

查看: 79|回复: 4

[求助] 学习过程中遇到的一个关于单精度变量的问题

[复制链接]

3

帖子

1

主题

0

精华

新人

F 币
28 元
贡献
13 点
发表于 2019-8-23 14:02:51 | 显示全部楼层 |阅读模式
[Fortran] 纯文本查看 复制代码
program ex0408
real a,b
a=1000
b=0.1
write(*,*)  a,"+",b,"=",a+b
stop
end



编译运行结果如下:
pcx@ubuntu:~/a/ab$ vim ex0408.f90
pcx@ubuntu:~/a/ab$ gfortran ex0408.f90
pcx@ubuntu:~/a/ab$ ./a.out
   1000.00000     +  0.100000001     =   1000.09998  


为什么我给b赋值0.1,但在运行的时候b变成了超过单精度浮点数有效位的0.100000001.
还望大佬们能给萌新解释一下。

回复

使用道具 举报

373

帖子

1

主题

0

精华

宗师

F 币
1967 元
贡献
1285 点
发表于 2019-8-23 15:02:56 | 显示全部楼层
本帖最后由 li913 于 2019-8-23 15:23 编辑

我用ivf,没这个情况,但gfortran的确是这样,搞不懂。
QQ截图20190823150210.png

3

帖子

1

主题

0

精华

新人

F 币
28 元
贡献
13 点
 楼主| 发表于 2019-8-27 21:09:39 | 显示全部楼层
li913 发表于 2019-8-23 15:02
我用ivf,没这个情况,但gfortran的确是这样,搞不懂。

感谢回复。可能是精度的问题,我查了一下,对gfortran编译器来说,0.1和0.1d0是两个不同的数字,我把这个程序 第4行改成:b=0.1d0以后,输出结果精度高了很多:   
1.0000000000000000      +  0.10000000000000001      =   1.1000000000000001      b=  0.10000000000000001  

因为导师让我在Ubuntu环境下编程,所以用不了ivf,感觉ivf的编译器更人性化一点。

708

帖子

0

主题

0

精华

大师

F 币
737 元
贡献
410 点

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

QQ
发表于 2019-8-28 08:04:21 | 显示全部楼层
本帖最后由 vvt 于 2019-8-28 08:14 编辑

精度肯定影响这个问题的。
不只是 0.1 和 0.1d0 不一样(常量的精度),变量也需要定义为相同的精度类型(kind=8之类的)

ivf和gfortran的不同在于:用表控格式输出write(*,*),对单精度(实际只有7位有效数字),gfortran会默认输出9个有效数字。但ivf只默认7个
[Fortran] 纯文本查看 复制代码
program ex0408
real a,b
real(8) a2,b2
a=1000
b=0.1
a2=1000d0
b2=0.1d0
write(*,*)  a,"+",b,"=",a+b
write(*,'(3(f10.5,a))')  a,"+",b,"=",a+b
write(*,'(3(es13.6,a))')  a,"+",b,"=",a+b
write(*,*)  a2,"+",b2,"=",a2+b2
end

所以你感觉ivf和gfortran精度不同而已。
实际上,你查看两者结果的二进制,是一样的。
或者用 es13.6 查看,也是一样的
(楼主可以用ivf和gfortran分别编译运行以上代码,对比查看)

3

帖子

1

主题

0

精华

新人

F 币
28 元
贡献
13 点
 楼主| 发表于 2019-9-2 08:42:09 | 显示全部楼层
vvt 发表于 2019-8-28 08:04
精度肯定影响这个问题的。
不只是 0.1 和 0.1d0 不一样(常量的精度),变量也需要定义为相同的精度类型(k ...

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

本版积分规则

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

GMT+8, 2019-9-21 17:09

Powered by Discuz! X3.2

© 2001-2017 Comsenz Inc.

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