diff options
author | Steve Kargl <kargl@FreeBSD.org> | 2014-09-17 19:01:22 +0000 |
---|---|---|
committer | Steve Kargl <kargl@FreeBSD.org> | 2014-09-17 19:01:22 +0000 |
commit | f382031d34cdcbaaa9a5f2fd17296574095b85ad (patch) | |
tree | cd812f7cced5f8cea889659130d00da95eac7fcd | |
parent | 21cffce593d4cc7eeade12627c06b65379144468 (diff) | |
download | src-f382031d34cdcbaaa9a5f2fd17296574095b85ad.tar.gz src-f382031d34cdcbaaa9a5f2fd17296574095b85ad.zip |
For targets that have a signed zero, lgamma_r(-0, &signgamp) should
set signgamp = -1.
Submitted by: enh at google dot com (e_lgamma[f]_r.c)
Notes
Notes:
svn path=/head/; revision=271719
-rw-r--r-- | lib/msun/ld128/e_lgammal_r.c | 6 | ||||
-rw-r--r-- | lib/msun/ld80/e_lgammal_r.c | 6 | ||||
-rw-r--r-- | lib/msun/src/e_lgamma_r.c | 6 | ||||
-rw-r--r-- | lib/msun/src/e_lgammaf_r.c | 6 |
4 files changed, 20 insertions, 4 deletions
diff --git a/lib/msun/ld128/e_lgammal_r.c b/lib/msun/ld128/e_lgammal_r.c index a612681d3d7d..13f8f75a46c6 100644 --- a/lib/msun/ld128/e_lgammal_r.c +++ b/lib/msun/ld128/e_lgammal_r.c @@ -241,7 +241,11 @@ lgammal_r(long double x, int *signgamp) *signgamp = 1; if((hx & 0x7fff) == 0x7fff) /* x is +-Inf or NaN */ return x*x; - if((hx==0||hx==0x8000)&&lx==0) return one/vzero; + if((hx==0||hx==0x8000)&&lx==0) { + if (hx&0x8000) + *signgamp = -1; + return one/vzero; + } /* purge off tiny and negative arguments */ if(fabsl(x)<0x1p-119L) { diff --git a/lib/msun/ld80/e_lgammal_r.c b/lib/msun/ld80/e_lgammal_r.c index bafa9a3b04ef..621b04e74022 100644 --- a/lib/msun/ld80/e_lgammal_r.c +++ b/lib/msun/ld80/e_lgammal_r.c @@ -261,7 +261,11 @@ lgammal_r(long double x, int *signgamp) *signgamp = 1; if((hx & 0x7fff) == 0x7fff) /* x is +-Inf or NaN */ return x*x; - if((hx==0||hx==0x8000)&&lx==0) return one/vzero; + if((hx==0||hx==0x8000)&&lx==0) { + if (hx&0x8000) + *signgamp = -1; + return one/vzero; + } ENTERI(); diff --git a/lib/msun/src/e_lgamma_r.c b/lib/msun/src/e_lgamma_r.c index ce2af94ca1c3..a6fd691ffd27 100644 --- a/lib/msun/src/e_lgamma_r.c +++ b/lib/msun/src/e_lgamma_r.c @@ -214,7 +214,11 @@ __ieee754_lgamma_r(double x, int *signgamp) *signgamp = 1; ix = hx&0x7fffffff; if(ix>=0x7ff00000) return x*x; - if((ix|lx)==0) return one/vzero; + if((ix|lx)==0) { + if(hx<0) + *signgamp = -1; + return one/vzero; + } if(ix<0x3b900000) { /* |x|<2**-70, return -log(|x|) */ if(hx<0) { *signgamp = -1; diff --git a/lib/msun/src/e_lgammaf_r.c b/lib/msun/src/e_lgammaf_r.c index f940845429d2..9d23053b1f3b 100644 --- a/lib/msun/src/e_lgammaf_r.c +++ b/lib/msun/src/e_lgammaf_r.c @@ -132,7 +132,11 @@ __ieee754_lgammaf_r(float x, int *signgamp) *signgamp = 1; ix = hx&0x7fffffff; if(ix>=0x7f800000) return x*x; - if(ix==0) return one/vzero; + if(ix==0) { + if(hx<0) + *signgamp = -1; + return one/vzero; + } if(ix<0x35000000) { /* |x|<2**-21, return -log(|x|) */ if(hx<0) { *signgamp = -1; |