14 x8
= 1.9_8
; x4
= 1.9_4
16 call check(erfc_scaled(x8
), erfc_scaled(1.9_8
))
17 call check(erfc_scaled(x4
), erfc_scaled(1.9_4
))
20 subroutine check_r4 (a
, b
)
21 real(kind
=4), intent(in
) :: a
, b
22 if (abs(a
- b
) > 1.e
-5 * abs(b
)) call abort
24 subroutine check_r8 (a
, b
)
25 real(kind
=8), intent(in
) :: a
, b
26 if (abs(a
- b
) > 1.e
-7 * abs(b
)) call abort