diff options
author | Joseph Myers <joseph@codesourcery.com> | 2015-02-26 21:06:34 +0000 |
---|---|---|
committer | Joseph Myers <joseph@codesourcery.com> | 2015-02-26 21:06:34 +0000 |
commit | 1d9ab20c14211abbf45514e63c190dfeec27fa5e (patch) | |
tree | cbe705e320799863f477d32c3062f6260c0c72a2 /NEWS | |
parent | ec0ce0d3be4ae21add09139c705852ed8bb114c7 (diff) | |
download | glibc-1d9ab20c14211abbf45514e63c190dfeec27fa5e.tar glibc-1d9ab20c14211abbf45514e63c190dfeec27fa5e.tar.gz glibc-1d9ab20c14211abbf45514e63c190dfeec27fa5e.tar.bz2 glibc-1d9ab20c14211abbf45514e63c190dfeec27fa5e.zip |
Fix ldbl-128/ldbl-128ibm acosl inaccuracy (bug 18038, bug 18039).
The ldbl-128 and ldbl-128ibm implementations of acosl have similar
bugs, using a threshold of 0x1p-57L to determine when they just return
pi/2. Since the result pi/2 - asinl (x) is roughly pi/2 - x for small
x, the relevant cut-off is actually x being < 0.5ulp of 1. This patch
fixes the implementations to use that cut-off and adds tests of small
acos arguments.
Tested for powerpc and mips64. Also tested for x86_64 and x86; no
ulps updates needed.
[BZ #18038]
[BZ #18039]
* sysdeps/ieee754/ldbl-128/e_acosl.c (__ieee754_acosl): Only
return pi/2 for arguments below 0x1p-113L.
* sysdeps/ieee754/ldbl-128ibm/e_acosl.c (__ieee754_acosl): Only
return pi/2 for arguments below 0x1p-106L.
* math/auto-libm-test-in: Add more tests of acos.
* math/auto-libm-test-out: Regenerated.
Diffstat (limited to 'NEWS')
-rw-r--r-- | NEWS | 2 |
1 files changed, 1 insertions, 1 deletions
@@ -12,7 +12,7 @@ Version 2.22 4719, 14841, 13064, 14094, 15319, 15467, 15790, 15969, 16351, 16560, 16783, 17269, 17523, 17569, 17588, 17792, 17836, 17912, 17916, 17932, 17944, 17949, 17964, 17965, 17967, 17969, 17978, 17987, 17991, 17996, - 17998, 17999, 18019, 18020, 18029, 18030, 18032. + 17998, 17999, 18019, 18020, 18029, 18030, 18032, 18038, 18039. * Character encoding and ctype tables were updated to Unicode 7.0.0, using new generator scripts contributed by Pravin Satpute and Mike FABIAN (Red |