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