fortran读取文件中实数的精度问题
我写了一个利用fortran程序读取文件数据的程序，发现数据的实数小数点后的数据总是不对呢？请大家帮忙看看啊，先谢谢了！程序如下：
PROGRAM t1
IMPLICIT NONE
INTEGER :: number_nodes
REAL, DIMENSION(6) :: info_nodes
INTEGER :: ierror
INTEGER :: i,j
! open the file
OPEN (UNIT=8, FILE = 'nodes1.dat', STATUS='OLD', ACTION = 'READ', IOSTAT=ierror)
! readintegers
READ (8,*) number_nodes
! read nodes
DO i = 1,number_nodes
READ(8,*) info_nodes(i)
END DO
WRITE(*,*) 'number of nodes is', number_nodes
DO i = 1,number_nodes
WRITE(*,*) 'node ', i , 'is ', info_nodes(i)
END DO
CLOSE(UNIT=8)
END PROGRAM t1
其中文件 nodes1.dat的格式为：
6
0.0
0.2
0.4
0.6
0.8
1.0
利用ubuntu下的gfortran编译，得到的结果为：
number of nodes is 6
node 1 is 0.00000000
node 2 is 0.200000003
node 3 is 0.400000006
node 4 is 0.600000024
node 5 is 0.800000012
node 6 is 1.00000000
其中node2---nodes5后面为什么总有些莫名奇妙小数呢？
刚学fortran不就，请大家不吝赐教.
这是正常的。并不是读取问题，而是浮点数问题。
单精度浮点只有6-7位有效数字。
gfortran的表控格式，会用9位有效数字来输出单精度，这样就会出现显示误差。
你可以指定输出格式来避免显示误差。但你无法避免内存中的存储误差。 我明白了，谢谢！ 比如我读取的是0.2，真实数值也是0.2吧，是显示有问题！可以这样理解吧。 本帖最后由 楚香饭 于 2015-12-30 23:03 编辑
不，真实数值也不是 0.2，IEEE浮点数无法精确的表达 0.2
就好像十进制无法精确的表达三分之一一样。它只能用 0.3333333 近似 matliu 发表于 2015-12-30 22:53
比如我读取的是0.2，真实数值也是0.2吧，是显示有问题！可以这样理解吧。
一般而言，双精度的浮点数机器误差大约是2.2e-16，所以在计算机里面的双精度0.2实际范围在0.2-2.2e-16到0.2+2.2e-16之间 明白，谢谢，高手！:-)
页:
[1]