增加双精度值

我现在正在为某个迭代运行一个程序。时间步长为0.01。我想在达到特定时间时写入一些信息。例如：增加双精度值

program abc 
    implicit none 
    double precision :: time,step,target 
    integer :: x 

    time = 0.d0 
    step = 0.01 
    target = 5.d0 

    do x = 1,6000 
    time = time + step 
    "some equations here to calculate the model parameters" 
    if(time.eq.target)then 
     write(*,*) "model parameters" 
     endif 
    enddo

然而，“时间”永不等于1.0或2.0或等它显示了如“0.999999866”而不是“1.0”和“1.99999845”，而不是“2.0”。

尽管我可以使用整数“x”来定义何时写信息，但我更喜欢使用时间步长。此外，我可能想要更改时间步长（0.01/0.02/0.05/etc）或目标（5.0/6.0/8.0/etc）。

有谁知道如何解决这个问题？谢谢你。

来源

2016-03-03 mee mee

'如果（x.eq.nint（目标/步骤））' – agentp

您现在已经发现了浮点算法！只要确保时间足够接近目标。

if(abs(time-target) < 0.5d0*step) then ...

应该做的伎俩。

来源

2016-03-03 06:46:24 RussF

谢谢RussF！现在它可以工作。你节省了我的时间。 –

浮点运算并不完美，根据变量的数字格式（32,64,128位），变量总是精确到某个机器错误。下面的例子说明本良好特性：

PROGRAM main 
    USE, INTRINSIC :: ISO_FORTRAN_ENV, qp => real128 
    IMPLICIT NONE 

    REAL(qp) :: a, b, c 

    a = 128._qp  
    b = a/120._qp + 1 
    c = 120._qp*(b-1) 

    PRINT*, "a = ", a 
    PRINT*, "c = ", c 
END PROGRAM main

这里是输出给此程序gfortran v.4.6.3：

a = 128.00000000000000000  
c = 127.99999999999999999

来源

2016-03-03 06:53:43 solalito

增加双精度值

回答

相关问题