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