[Bug fortran/33698] FAIL: gfortran.dg/gamma_5.f90

[dave at hiauly1 dot hia dot nrc dot ca]

------- Comment #13 from dave at hiauly1 dot hia dot nrc dot ca  2007-11-03 19:37 -------
Subject: Re:  FAIL: gfortran.dg/gamma_5.f90

> This implements the fallback functions, but naturally
> doesn't do anything on my linux system (which has all tgamma*
> and lgamma* functions).

You might hack config.h to test.

> Can anybody test this?  John?

c99_functions.c didn't compile due to some typos, etc.  I've attached
ammended version of the patch which builds on hppa2.0w-hp-hpux11.11.

Now, the hard part.  The gfortran.dg/gamma_5.f90 test fails at n = 16.

Breakpoint 2, tgamma (z=16.5)
    at ../../../gcc/libgfortran/intrinsics/c99_functions.c:1730
1730      if (z < 0.0)
(gdb)
Continuing.

Breakpoint 1, _gfortran_abort ()
    at ../../../gcc/libgfortran/intrinsics/abort.c:38
38        close_units ();

It fails here:
     if (abs(gamma(xd)-td)/td > 5e-14) call abort

I hacked the test a bit to make it easier to see the return value
from gamma(xd)-td):

(gdb) p gd
$1 = 5189998453040.5889
(gdb) p td
$2 = 5189998453040.125
(gdb) p (gd-td)/td
$4 = 8.9377134058350696e-14

I guess the first thing to check is the coefficients.

Dave


------- Comment #14 from dave at hiauly1 dot hia dot nrc dot ca  2007-11-03 19:37 -------
Created an attachment (id=14479)
 --> (http://gcc.gnu.org/bugzilla/attachment.cgi?id=14479&action=view)


-- 


http://gcc.gnu.org/bugzilla/show_bug.cgi?id=33698