程序运行中数值精度下降

jojo9

如题，求助论坛的各位大大。程序片段以及断点位置如图：


    iof%rie(1) = Velc + gam2m * iof%ac_c
    iof%rie(2) = iof%Vel8 - gam2m * 1.0

    iof%Ubar = ( iof%rie(1) + iof%rie(2) ) / 2.0


这里在调试下可知 iof%rie(1) = 5.00000，iof%rie(2) = -5.00000，都是单精度浮点数，但是二者求和除以2.0（或者2.d0）后等-2.3841858E-07？？


请问各位大大，造成这种现象的可能原因能有哪些呢？另外我还单独试了试，其他程序上单精度的（5.00000-5.00000）/2.0 就是等于0.0。

vvt

这是非常正常的事情。
单精度浮点数只有7位有效数字。

正如人类习惯用十进制。
计算 (1/3) = 0.3333333 ，然后 0.3333333*3=0.9999999
是吧？(1/3)*3竟然不等于1

计算机使用二进制，一样有误差的。使用双精度可以一定程度降低误差。