From facd1d8ef35edb863523a7165a64d55196cc2d36 Mon Sep 17 00:00:00 2001 From: Ulrich Drepper Date: Thu, 10 Apr 2008 04:58:03 +0000 Subject: [BZ #4407] 2008-02-12 Pete Eberlein [BZ #4407] * sysdeps/ieee754/ldbl-128/e_lgammal_r.c (__ieee754_lgammal_r): Preserve sign in signgamp when x is zero. --- sysdeps/ieee754/ldbl-128/e_lgammal_r.c | 6 ++++++ 1 file changed, 6 insertions(+) (limited to 'sysdeps/ieee754/ldbl-128/e_lgammal_r.c') diff --git a/sysdeps/ieee754/ldbl-128/e_lgammal_r.c b/sysdeps/ieee754/ldbl-128/e_lgammal_r.c index b9302974c0..d080448476 100644 --- 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; -- cgit v1.2.3-70-g09d2