1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
|
! { 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
|
