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

jojo

[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.
还望大佬们能给萌新解释一下。

jojo

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的编译器更人性化一点。

jojo

vvt 发表于 2019-8-28 08:04
精度肯定影响这个问题的。
不只是 0.1 和 0.1d0 不一样（常量的精度），变量也需要定义为相同的精度类型（k ...

非常感谢