的Fortran 95：直列評價如果條件

Question

這裏的代碼，對於一個真實值返回的ε-（）一個小片段：

program epstest 
real :: eps=1.0, d 
do 
    d=1.0+eps 
    if (d==1.0) then 
    eps=eps*2 
    exit 
    else 
    eps=eps/2 
    end if 
end do 
write(*,*) eps, epsilon(d) 
pause 
end program 

現在，當我通過

if (1.0+eps==1.0) then 

更換如果條件

該程序應該有相同的回報，但不幸的是它不！我在Linux和Windows上使用g95的最新（快照）發行版進行了測試。

有人可以向我解釋這個問題嗎？

Answer 1

浮點算法有許多微妙的問題。一種形式的源代碼可以產生與另一個看起來幾乎相同的源代碼不同的機器指令。例如，「d」可能會被存儲到內存位置，而「1.0 + eps」可能僅通過寄存器進行評估......這會導致不同的精度。更一般地說，爲什麼不使用爲Fortran 95提供的內部函數來揭示reals特定精度的特性？