Problems with GAMMA functions?

[Thomas Koenig]

On Tue, 2007-10-09 at 09:23 +0200, Angelo Graziosi wrote:

> > - relax the tolerances in gamma_5.f90 to something that doesn't fail
> >   on as many systems
> 
> This does not seem a good solution: it is too simple to relax the
> tollerance so that the test passes. Too simple!

Hello Angelo,

I have submitted this inaccuracy as a bug to the glibc people
as http://sources.redhat.com/bugzilla/show_bug.cgi?id=5159 .
Gcc, g++ and gfortran on glibc-based systems will all
benefit if a solution is found there.

If this is important for you, I would suggest that you contribute
on this bug report.  In the meantime, you are probably better
off using the external gamma function from CERNLIB.

Again, thanks a lot for your bug report!

	Thomas