diff options
author | Joseph Myers <joseph@codesourcery.com> | 2015-06-29 17:51:32 +0000 |
---|---|---|
committer | Joseph Myers <joseph@codesourcery.com> | 2015-06-29 17:51:32 +0000 |
commit | 4aa10d01b654b7395c6e856b8e1e0959cfc80739 (patch) | |
tree | 01a52e164acecacf0e5b0e1851740d5defcd9d34 /sysdeps/ieee754/dbl-64/e_gamma_r.c | |
parent | 63dbe5f32238858c7b953b867ed0588c7808dd4f (diff) | |
download | glibc-4aa10d01b654b7395c6e856b8e1e0959cfc80739.tar.gz glibc-4aa10d01b654b7395c6e856b8e1e0959cfc80739.tar.xz glibc-4aa10d01b654b7395c6e856b8e1e0959cfc80739.zip |
Fix ldbl-128 j1l spurious underflows (bug 18612).
The ldbl-128 implementation of j1l produces spurious underflow exceptions for some small arguments, as a result of squaring the argument. This patch fixes it just to use a linear approximation for sufficiently small arguments, and then to force an underflow exception only in the cases where it is required. Tested for mips64. [BZ #18612] * sysdeps/ieee754/ldbl-128/e_j1l.c (__ieee754_j1l): For small arguments, just return 0.5 times the argument, with underflow forced as needed. * math/auto-libm-test-in: Add more tests of j1. * math/auto-libm-test-out: Regenerated.
Diffstat (limited to 'sysdeps/ieee754/dbl-64/e_gamma_r.c')
0 files changed, 0 insertions, 0 deletions