EPSILON
copyright

green dotEPSILON

書式EPSILON(X)
機能 EPSILONは、引数Xと同じ型および同じ種別型パラメタの数体系における、 1に対してほとんど無視できる正の数(表2の数体系でb1-p)を 求める。
引数 X:実数型でスカラまたは配列
戻り値 型、種別型パラメタは引数と同じで、形状はスカラである。

green dotサンプルプログラム

HUGETINY、EPSILONの例
    INTEGER,PARAMETER :: k=SELECTED_INT_KIND(10)
    INTEGER,PARAMETER :: dp=SELECTED_REAL_KIND(10)
    INTEGER :: i ;  INTEGER(KIND=k) :: j
    REAL    :: x ;  REAL(KIND=dp)   :: y
    WRITE(*,*) HUGE(i)    ![例] 231-1  = 2147483647
    WRITE(*,*) HUGE(j)    ![例] 263-1  = 9223372036854775807
    WRITE(*,*) HUGE(x)    ![例] (1-2-24)×2128 = 0.340282…E+39
    WRITE(*,*) HUGE(y)    ![例] (1-2-53)×21024 = 0.179769…E+309
    WRITE(*,*) TINY(x)    ![例] 2-126  = 0.117549…E-37
    WRITE(*,*) TINY(y)    ![例] 2-1022 = 0.222507…E-307
    WRITE(*,*) EPSILON(x) ![例] 2-23   = 0.119209…E-06
    WRITE(*,*) EPSILON(y) ![例] 2-52   = 0.222044…E-15

(表2)実数の数体系
hyou2 b,p2以上の整数
fk 0以上b未満(f1≠0)
s+1または -1
e最小値の整数eminから最大値の
整数 emaxまでの間の整数

line-end