[BZ #4407]

2008-02-12  Pete Eberlein  <eberlein@us.ibm.com>
	[BZ #4407]
	* sysdeps/ieee754/ldbl-128/e_lgammal_r.c (__ieee754_lgammal_r):
	Preserve sign in signgamp when x is zero.

diff --git a/ChangeLog b/ChangeLog
index 247c79247fba28bfe76210f31bfe683dd5a50037..374553e507e4fba765b221497dac17682e377b51 100644 (file)
--- a/ChangeLog
+++ b/ChangeLog
@@ -1,3 +1,9 @@
+2008-02-12  Pete Eberlein  <eberlein@us.ibm.com>
+
+       [BZ #4407]
+       * sysdeps/ieee754/ldbl-128/e_lgammal_r.c (__ieee754_lgammal_r):
+       Preserve sign in signgamp when x is zero.
+
 2007-10-26  Steven Munroe  <sjmunroe@us.ibm.com>
 
        * sysdeps/powerpc/powerpc32/power4/fpu/w_sqrt.S: New File.

diff --git a/sysdeps/ieee754/ldbl-128/e_lgammal_r.c b/sysdeps/ieee754/ldbl-128/e_lgammal_r.c
index b9302974c00948394380738086ecae4c82b8be60..d08044847671e66d68feb2b74d15c6ddc93f5a41 100644 (file)
--- a/sysdeps/ieee754/ldbl-128/e_lgammal_r.c
+++ b/sysdeps/ieee754/ldbl-128/e_lgammal_r.c
@@ -772,6 +772,12 @@ __ieee754_lgammal_r (x, signgamp)
   if (! __finitel (x))
     return x * x;
 
+  if (x == 0.0L)
+    {
+      if (__signbitl (x))
+        *signgamp = -1;
+    }
+
   if (x < 0.0L)
     {
       q = -x;