3 ! This test checks whether the largest possible
4 ! floating-point number works. That's usually
5 ! REAL(16) -- either because the hardware supports it or
6 ! because of libquadmath. However, it can also be
10 use iso_fortran_env, only: real_kinds
12 integer, parameter :: QP = real_kinds(ubound(real_kinds,dim=1))
13 real(QP), parameter :: Z1 = 1,HALF_PI = asin(Z1),PI = HALF_PI+HALF_PI
14 real(QP) :: x = 0.124_QP
15 complex(QP) :: z = 0.124_QP
16 print *, 'kind = ', qp
19 print *, 16*atan(0.2_QP)-4*atan(Z1/239)
29 print *, PI*cmplx(0.0_qp, 1.0_qp)
30 print *, 16*atan(0.2_QP)-4*atan(Z1/239)
33 print *, sinh(z) ! asinh not implemented, cf. PR 46416