diff options
author | Ulrich Drepper <drepper@redhat.com> | 2001-02-19 23:16:05 +0000 |
---|---|---|
committer | Ulrich Drepper <drepper@redhat.com> | 2001-02-19 23:16:05 +0000 |
commit | f14bd805a73308dd187d9fc768358a07a8efc54a (patch) | |
tree | 3cd722357ad98e9af17ba75e95603bce20c067f6 /sysdeps/ieee754/dbl-64/e_pow.c | |
parent | d210ca02ec8a5ba034946c57b85c2f3e503a6017 (diff) | |
download | glibc-f14bd805a73308dd187d9fc768358a07a8efc54a.tar.gz |
Update.
* sysdeps/ieee754/flt-32/e_powf.c: Handle x == +-1 correctly.
* sysdeps/ieee754/dbl-64/e_pow.c: Likewise.
Diffstat (limited to 'sysdeps/ieee754/dbl-64/e_pow.c')
-rw-r--r-- | sysdeps/ieee754/dbl-64/e_pow.c | 4 |
1 files changed, 4 insertions, 0 deletions
diff --git a/sysdeps/ieee754/dbl-64/e_pow.c b/sysdeps/ieee754/dbl-64/e_pow.c index 1e1496f00d..73e8f471f6 100644 --- a/sysdeps/ieee754/dbl-64/e_pow.c +++ b/sysdeps/ieee754/dbl-64/e_pow.c @@ -149,6 +149,10 @@ C[] = { /* y==zero: x**0 = 1 */ if((iy|ly)==0) return C[1]; + /* x==+-1 */ + if(x == 1.0) return C[1]; + if(x == -1.0 && isinf(y)) return C[1]; + /* +-NaN return x+y */ if(ix > 0x7ff00000 || ((ix==0x7ff00000)&&(lx!=0)) || iy > 0x7ff00000 || ((iy==0x7ff00000)&&(ly!=0))) |