1
2
3
4
5
6
7
8
9
10
11
12
13
|
program intrinsics_34
integer, parameter :: dp=kind(0d0)
real :: x = 3.143
real(dp) :: y = 2.33
print *, epsilon(x)
print *, epsilon(y)
print *, epsilon(1._dp) ** 0.5 !Part of Minpack
! Below numbers are corresponding output of gfortran and not magic numbers.
if (abs(epsilon(x) - 1.19209290E-07) > 1e-7) error stop
if (abs(epsilon(y) - 2.2204460492503131E-016_dp) > 1e-15_dp) error stop
if (abs((epsilon(1._dp) ** 0.5_dp) - 1.4901161193847656E-008_dp) > 1e-15_dp) error stop
end program
|
