diff options
| author | Joseph Myers <joseph@codesourcery.com> | 2012-11-22 19:56:47 +0000 |
|---|---|---|
| committer | Joseph Myers <joseph@codesourcery.com> | 2012-11-22 19:56:47 +0000 |
| commit | 0a42601f272ea962b200d452a54e0ca374163f60 (patch) | |
| tree | 79f2d88d56d8c44a291ac0e235ace15b8a3cfe6d /sysdeps/ieee754/ldbl-128ibm/s_atanl.c | |
| parent | 79c9b9feb45d175d4a394edd9ecf89cde97890c4 (diff) | |
| download | glibc-0a42601f272ea962b200d452a54e0ca374163f60.zip glibc-0a42601f272ea962b200d452a54e0ca374163f60.tar.gz glibc-0a42601f272ea962b200d452a54e0ca374163f60.tar.bz2 | |
Fix ldbl-128ibm atanl spurious underflows (bug 14871).
Diffstat (limited to 'sysdeps/ieee754/ldbl-128ibm/s_atanl.c')
| -rw-r--r-- | sysdeps/ieee754/ldbl-128ibm/s_atanl.c | 16 |
1 files changed, 16 insertions, 0 deletions
diff --git a/sysdeps/ieee754/ldbl-128ibm/s_atanl.c b/sysdeps/ieee754/ldbl-128ibm/s_atanl.c index 779209d..2a36d16 100644 --- a/sysdeps/ieee754/ldbl-128ibm/s_atanl.c +++ b/sysdeps/ieee754/ldbl-128ibm/s_atanl.c @@ -199,6 +199,22 @@ __atanl (long double x) return atantbl[83]; } + if (k <= 0x3c800000) /* |x| <= 2**-55. */ + { + /* Raise inexact. */ + if (1e300L + x > 0.0) + return x; + } + + if (k >= 0x46c00000) /* |x| >= 2**109. */ + { + /* Saturate result to {-,+}pi/2. */ + if (sign) + return -atantbl[83]; + else + return atantbl[83]; + } + if (sign) x = -x; |