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 /lib/msun/ld128 | |
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
Diffstat (limited to 'lib/msun/ld128')
-rw-r--r-- | lib/msun/ld128/e_lgammal_r.c | 6 |
1 files changed, 5 insertions, 1 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) { |