diff options
author | Thomas Koenig <tkoenig@gcc.gnu.org> | 2007-10-07 21:36:09 +0000 |
---|---|---|
committer | Thomas Koenig <tkoenig@gcc.gnu.org> | 2007-10-07 21:36:09 +0000 |
commit | 4bf3289bd8cff84abfd38b4d4b78edf187f4a256 (patch) | |
tree | 26ffeef73c65d0bc34b9aa827a332404fe3efbc1 | |
parent | cbc9503d788ea852d7821d1baf6711ee619492ee (diff) | |
download | gcc-4bf3289bd8cff84abfd38b4d4b78edf187f4a256.zip gcc-4bf3289bd8cff84abfd38b4d4b78edf187f4a256.tar.gz gcc-4bf3289bd8cff84abfd38b4d4b78edf187f4a256.tar.bz2 |
re PR libfortran/33683 (calculating lgamma instead of gamma)
2007-10-07 Thomas Koenig <tkoenig@gcc.gnu.org>
PR libfortran/33683
* mathbuiltins.def (GAMMA): Change function name to
"tgamma" instad of "gamma".
2007-10-07 Thomas Koenig <tkoenig@gcc.gnu.org>
PR libfortran/33683
* gfortran.dg/gamma_5.f90: New test case
From-SVN: r129116
-rw-r--r-- | gcc/fortran/ChangeLog | 6 | ||||
-rw-r--r-- | gcc/fortran/mathbuiltins.def | 2 | ||||
-rw-r--r-- | gcc/testsuite/ChangeLog | 5 | ||||
-rw-r--r-- | gcc/testsuite/gfortran.dg/gamma_5.f90 | 42 |
4 files changed, 54 insertions, 1 deletions
diff --git a/gcc/fortran/ChangeLog b/gcc/fortran/ChangeLog index 28ed3f7..0f5758a 100644 --- a/gcc/fortran/ChangeLog +++ b/gcc/fortran/ChangeLog @@ -1,3 +1,9 @@ +2007-10-07 Thomas Koenig <tkoenig@gcc.gnu.org> + + PR libfortran/33683 + * mathbuiltins.def (GAMMA): Change function name to + "tgamma" instad of "gamma". + 2007-10-07 Tobias Schlüter <tobi@gcc.gnu.org> PR fortran/20851 diff --git a/gcc/fortran/mathbuiltins.def b/gcc/fortran/mathbuiltins.def index 33e87d1..37758e1a 100644 --- a/gcc/fortran/mathbuiltins.def +++ b/gcc/fortran/mathbuiltins.def @@ -30,5 +30,5 @@ DEFINE_MATH_BUILTIN (Y1, "y1", 0) DEFINE_MATH_BUILTIN (YN, "yn", 2) DEFINE_MATH_BUILTIN (ERF, "erf", 0) DEFINE_MATH_BUILTIN (ERFC, "erfc", 0) -DEFINE_MATH_BUILTIN (GAMMA, "gamma", 0) +DEFINE_MATH_BUILTIN (GAMMA, "tgamma", 0) DEFINE_MATH_BUILTIN (LGAMMA,"lgamma", 0) diff --git a/gcc/testsuite/ChangeLog b/gcc/testsuite/ChangeLog index 102c108..3997ac3 100644 --- a/gcc/testsuite/ChangeLog +++ b/gcc/testsuite/ChangeLog @@ -1,3 +1,8 @@ +2007-10-07 Thomas Koenig <tkoenig@gcc.gnu.org> + + PR libfortran/33683 + * gfortran.dg/gamma_5.f90: New test case + 2007-10-07 Tobias Schlüter <tobi@gcc.gnu.org> PR fortran/20851 diff --git a/gcc/testsuite/gfortran.dg/gamma_5.f90 b/gcc/testsuite/gfortran.dg/gamma_5.f90 new file mode 100644 index 0000000..6945cfd --- /dev/null +++ b/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 > 2e-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 |