Fortran Coder

查看: 21797|回复: 6
打印 上一主题 下一主题

[文件读写] fortran读取文件中实数的精度问题

[复制链接]

6

帖子

2

主题

0

精华

入门

F 币
31 元
贡献
18 点
跳转到指定楼层
楼主
发表于 2015-12-30 15:07:32 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式
我写了一个利用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)
! read  integers
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不就,请大家不吝赐教.

分享到:  微信微信
收藏收藏 点赞点赞 点踩点踩

2033

帖子

12

主题

5

精华

论坛跑堂

臭石头雪球

F 币
1641 元
贡献
709 点

美女勋章热心勋章星光勋章新人勋章贡献勋章管理勋章帅哥勋章爱心勋章规矩勋章元老勋章水王勋章

沙发
发表于 2015-12-30 16:25:30 | 只看该作者
这是正常的。并不是读取问题,而是浮点数问题。
单精度浮点只有6-7位有效数字。

gfortran的表控格式,会用9位有效数字来输出单精度,这样就会出现显示误差。
你可以指定输出格式来避免显示误差。但你无法避免内存中的存储误差。

6

帖子

2

主题

0

精华

入门

F 币
31 元
贡献
18 点
板凳
 楼主| 发表于 2015-12-30 22:50:27 | 只看该作者
我明白了,谢谢!

6

帖子

2

主题

0

精华

入门

F 币
31 元
贡献
18 点
地板
 楼主| 发表于 2015-12-30 22:53:48 | 只看该作者
比如我读取的是0.2,真实数值也是0.2吧,是显示有问题!可以这样理解吧。

737

帖子

4

主题

0

精华

大师

农村外出务工人员

F 币
705 元
贡献
361 点

新人勋章爱心勋章水王勋章元老勋章热心勋章

5#
发表于 2015-12-30 23:02:55 | 只看该作者
本帖最后由 楚香饭 于 2015-12-30 23:03 编辑

不,真实数值也不是 0.2,IEEE浮点数无法精确的表达 0.2
就好像十进制无法精确的表达三分之一一样。它只能用 0.3333333 近似

490

帖子

4

主题

0

精华

大宗师

F 币
3298 元
贡献
1948 点

水王勋章元老勋章热心勋章

6#
发表于 2015-12-31 09:42:08 | 只看该作者
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之间

6

帖子

2

主题

0

精华

入门

F 币
31 元
贡献
18 点
7#
 楼主| 发表于 2015-12-31 15:25:02 | 只看该作者
明白,谢谢,高手!
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

捐赠本站|Archiver|关于我们 About Us|小黑屋|Fcode ( 京ICP备18005632-2号 )

GMT+8, 2024-12-25 09:46

Powered by Tencent X3.4

© 2013-2024 Tencent

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