diff options
author | Ben Cheng <bccheng@google.com> | 2013-03-28 11:14:20 -0700 |
---|---|---|
committer | Ben Cheng <bccheng@google.com> | 2013-03-28 12:40:33 -0700 |
commit | af0c51ac87ab2a87caa03fa108f0d164987a2764 (patch) | |
tree | 4b8b470f7c5b69642fdab8d0aa1fbc148d02196b /gcc-4.8/gcc/testsuite/gfortran.dg/gamma_5.f90 | |
parent | d87cae247d39ebf4f5a6bf25c932a14d2fdb9384 (diff) | |
download | toolchain_gcc-af0c51ac87ab2a87caa03fa108f0d164987a2764.tar.gz toolchain_gcc-af0c51ac87ab2a87caa03fa108f0d164987a2764.tar.bz2 toolchain_gcc-af0c51ac87ab2a87caa03fa108f0d164987a2764.zip |
[GCC 4.8] Initial check-in of GCC 4.8.0
Change-Id: I0719d8a6d0f69b367a6ab6f10eb75622dbf12771
Diffstat (limited to 'gcc-4.8/gcc/testsuite/gfortran.dg/gamma_5.f90')
-rw-r--r-- | gcc-4.8/gcc/testsuite/gfortran.dg/gamma_5.f90 | 42 |
1 files changed, 42 insertions, 0 deletions
diff --git a/gcc-4.8/gcc/testsuite/gfortran.dg/gamma_5.f90 b/gcc-4.8/gcc/testsuite/gfortran.dg/gamma_5.f90 new file mode 100644 index 000000000..467c57962 --- /dev/null +++ b/gcc-4.8/gcc/testsuite/gfortran.dg/gamma_5.f90 @@ -0,0 +1,42 @@ +! { dg-do run } +! PR 33683 - we used to pick up the wrong gamma function +! from the library on some systems. +program main + implicit none + integer, parameter :: n_max = 20 + double precision, dimension(0:n_max) :: c + double precision :: pi + integer :: n + double precision :: td, xd + real :: ts,xs + + pi = 4 * atan(1.d0) + c(0) = 1. + do n=1, n_max + c(n) = (2*n-1)*c(n-1)*0.5d0 + end do + + do n=1, n_max + xs = n + 0.5 + xd = n + 0.5d0 + td = c(n)*sqrt(pi) + ts = c(n)*sqrt(pi) + if (abs(gamma(xs)-ts)/ts > 9e-6) call abort + if (abs(gamma(xd)-td)/td > 5e-14) call abort + end do + call tst_s(2.3, gamma(2.3)) + call tst_s(3.7, gamma(3.7)) + call tst_s(5.5, gamma(5.5)) + call tst_d(4.2d0, gamma(4.2d0)) + call tst_d(8.1d0, gamma(8.1d0)) +contains + subroutine tst_s(a, b) + real :: a, b + if (abs(gamma(a) - b)/b > 1e-6) call abort + end subroutine tst_s + + subroutine tst_d(a, b) + double precision :: a,b + if (abs(gamma(a) - b)/b > 5e-14) call abort + end subroutine tst_d +end program main |