The 80-bit extended precision version of tgammal(3) has the
tgamma(+-0) == +-Inf
corner case wrong.
ok?
Index: src/ld80/e_tgammal.c
===================================================================
RCS file: /home/cvs/src/lib/libm/src/ld80/e_tgammal.c,v
retrieving revision 1.2
diff -u -p -r1.2 e_tgammal.c
--- src/ld80/e_tgammal.c 20 Jul 2011 21:02:51 -0000 1.2
+++ src/ld80/e_tgammal.c 5 Aug 2013 10:14:28 -0000
@@ -229,6 +229,8 @@ if(x == INFINITY)
return(INFINITY);
if(x == -INFINITY)
return(x - x);
+if( x == 0.0L )
+ return( 1.0L / x );
q = fabsl(x);
if( q > 13.0L )
