diff options
Diffstat (limited to 'gcc/testsuite/gfortran.dg/gamma_1.f90')
| -rw-r--r-- | gcc/testsuite/gfortran.dg/gamma_1.f90 | 29 |
1 files changed, 29 insertions, 0 deletions
diff --git a/gcc/testsuite/gfortran.dg/gamma_1.f90 b/gcc/testsuite/gfortran.dg/gamma_1.f90 new file mode 100644 index 000000000..994616695 --- /dev/null +++ b/gcc/testsuite/gfortran.dg/gamma_1.f90 @@ -0,0 +1,29 @@ +! { dg-do run } +! +! 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, log_gamma +integer, parameter :: sp = kind(1.0) +integer, parameter :: dp = kind(1.0d0) + +real(sp) :: rsp +real(dp) :: rdp + +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(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(log_gamma(1.0_sp)) > tiny(1.0_sp)) call abort() +if (abs(log_gamma(1.0_dp)) > tiny(1.0_dp)) 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 + |