! { dg-do link }
!
! This test checks whether the largest possible
! floating-point number works. That's usually
! REAL(16) -- either because the hardware supports it or
! because of libquadmath. However, it can also be
! REAL(10) or REAL(8)
!
program test_qp
   use iso_fortran_env, only: real_kinds
   implicit none
   integer, parameter :: QP = real_kinds(ubound(real_kinds,dim=1))
   real(QP), parameter :: Z1 = 1,HALF_PI = asin(Z1),PI = HALF_PI+HALF_PI
   real(QP) :: x = 0.124_QP
   complex(QP) :: z = 0.124_QP
   print *, 'kind = ', qp
   print *, x
   print *, PI
   print *, 16*atan(0.2_QP)-4*atan(Z1/239)
   print *, sin(PI)
   print *, cos(HALF_PI)
   print *, asinh(PI)
   print *, erfc(Z1)
   print *, epsilon(x)
   print *, precision(x)
   print *, digits(x)

   print *, z
   print *, PI*cmplx(0.0_qp, 1.0_qp)
! Disable the complex functions as not all "long-double" systems have
! a libm with those C99 functions. (libquadmath had), cf. PR 46584
!   print *, 16*atan(0.2_QP)-4*atan(Z1/239)
!   print *, sin(z)
!   print *, cos(z)
!   print *, sinh(z) ! asinh not implemented in libquadmath, cf. PR 46416
   print *, precision(z)
end program test_qp