Subversion Repositories openrisc

! { 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