+2007-08-28 Francois-Xavier Coudert <fxcoudert@gcc.gnu.org>
+
+ * gfortran.dg/gamma_1.f90: Remove large-real tests.
+ * gfortran.dg/gamma_2.f90: Remove large-real tests.
+ * gfortran.dg/gamma_4.f90: New: large-real tests.
+
2007-08-28 Jason Merrill <jason@redhat.com>
PR c++/28558
intrinsic :: gamma, lgamma
integer, parameter :: sp = kind(1.0)
integer, parameter :: dp = kind(1.0d0)
-integer, parameter :: qp = selected_real_kind(p=15,r=900)
real(sp) :: rsp
real(dp) :: rdp
-real(dp) :: rqp
if (abs(gamma(1.0_sp) - 1.0_sp) > tiny(1.0_sp)) call abort()
if (abs(gamma(1.0_dp) - 1.0_dp) > tiny(1.0_dp)) call abort()
-if (abs(gamma(1.0_qp) - 1.0_qp) > tiny(1.0_qp)) call abort()
if (abs(dgamma(1.0_dp) - 1.0_dp) > tiny(1.0_dp)) call abort()
if (abs(lgamma(1.0_sp)) > tiny(1.0_sp)) call abort()
if (abs(lgamma(1.0_dp)) > tiny(1.0_dp)) call abort()
-if (abs(lgamma(1.0_qp)) > tiny(1.0_qp)) call abort()
if (abs(algama(1.0_sp)) > tiny(1.0_sp)) call abort()
if (abs(dlgama(1.0_dp)) > tiny(1.0_dp)) call abort()
end program gamma_test
integer, parameter :: sp = kind(1.0)
integer, parameter :: dp = kind(1.0d0)
-integer, parameter :: qp = selected_real_kind(p=15,r=900)
real(sp) :: rsp = 1.0_sp
real(dp) :: rdp = 1.0_dp
-real(dp) :: rqp = 1.0_qp
rsp = gamma(rsp) ! FIXME: "is not included in the selected standard"
rdp = gamma(rdp) ! FIXME: "is not included in the selected standard"
-rqp = gamma(rqp) ! FIXME: "is not included in the selected standard"
rdp = dgamma(rdp) ! { dg-error "is not included in the selected standard" }
rsp = lgamma(rsp) ! FIXME: "is not included in the selected standard"
rdp = lgamma(rdp) ! FIXME: "is not included in the selected standard"
-rqp = lgamma(rqp) ! FIXME: "is not included in the selected standard"
rsp = algama(rsp) ! { dg-error "is not included in the selected standard" }
rdp = dlgama(rdp) ! { dg-error "is not included in the selected standard" }
end subroutine foo
--- /dev/null
+! { dg-do run }
+! { dg-require-effective-target fortran_large_real }
+!
+! Test the vendor intrinsic (d)gamma, lgamma and algama/dlgama
+! gamma is also part of the Fortran 2008 draft; lgamma is called
+! log_gamma in the Fortran 2008 draft.
+!
+! PR fortran/32980
+!
+program gamma_test
+implicit none
+intrinsic :: gamma, lgamma
+integer, parameter :: qp = selected_real_kind(p=15,r=900)
+
+real(qp) :: rqp
+
+if (abs(gamma(1.0_qp) - 1.0_qp) > tiny(1.0_qp)) call abort()
+if (abs(lgamma(1.0_qp)) > tiny(1.0_qp)) call abort()
+end program gamma_test
+