作者sjgau (sjgau)
看板Fortran
标题Re: [问题] 关於误差估计的问题
时间Thu May 20 11:25:58 2010
换成 real*8 以後，
(1.0D0 + 1.0D-8) 的结果，可以被计算出来，
有效位数，可以达到
2^52 --> 大概是 15位数
大家可以试试看
program VF1035
implicit none
real*8 x, y, z, ct1
x= 4.0D0*datan(1.0D0)
ct1= 0.0D0
print *, 'ct1, x= ', ct1, x
! pause
y= x + 1.0D0
! y > x
do while (y .GT. x)
x= x*(1.0D0 + 1.0D-8)
y= x + 1.0D0
ct1= ct1 + 1.0D0
end do
! y <= x
print *, 'x, y, ct1= ', x, y, ct1
end program VF1035
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 61.20.174.109