Fix ldbl-128ibm sinhl inaccuracy near 0 (bug 18789).

ldbl-128ibm sinhl uses a too-big threshold to decide when to return the argument, resulting in large errors. This patch fixes it to use a more appropriate threshold. Tested for x86_64, x86 and powerpc. [BZ #18789] * sysdeps/ieee754/ldbl-128ibm/e_sinhl.c (__ieee754_sinhl): Use smaller threshold for returning the argument. * math/auto-libm-test-in: Add more tests of sinh. * math/auto-libm-test-out: Regenerated. * sysdeps/i386/fpu/libm-test-ulps: Update.
git-mirror · Aug 10, 2015 · d0649b2 · d0649b2
1 parent d2ff039
commit d0649b2
Show file tree

Hide file tree

Showing 6 changed files with 785 additions and 2 deletions.
diff --git a/ChangeLog b/ChangeLog
@@ -1,3 +1,12 @@
+2015-08-10  Joseph Myers  <joseph@codesourcery.com>
+
+	[BZ #18789]
+	* sysdeps/ieee754/ldbl-128ibm/e_sinhl.c (__ieee754_sinhl): Use
+	smaller threshold for returning the argument.
+	* math/auto-libm-test-in: Add more tests of sinh.
+	* math/auto-libm-test-out: Regenerated.
+	* sysdeps/i386/fpu/libm-test-ulps: Update.
+
 2015-08-09  H.J. Lu  <hongjiu.lu@intel.com>
 
 	[BZ #18674]

diff --git a/NEWS b/NEWS
@@ -10,7 +10,7 @@ Version 2.23
 * The following bugs are resolved with this release:
 
   16517, 16519, 17905, 18265, 18480, 18525, 18618, 18647, 18661, 18674,
-  18787.
+  18787, 18789.
 
 Version 2.22
 

diff --git a/math/auto-libm-test-in b/math/auto-libm-test-in
@@ -2508,6 +2508,28 @@ sinh 0
 sinh -0
 sinh 0.75
 sinh 0x8p-32
+sinh 0x1p-5
+sinh -0x1p-5
+sinh 0x1p-10
+sinh -0x1p-10
+sinh 0x1p-20
+sinh -0x1p-20
+sinh 0x1p-30
+sinh -0x1p-30
+sinh 0x1p-40
+sinh -0x1p-40
+sinh 0x1p-50
+sinh -0x1p-50
+sinh 0x1p-60
+sinh -0x1p-60
+sinh 0x1p-70
+sinh -0x1p-70
+sinh 0x1p-100
+sinh -0x1p-100
+sinh 0x1p-1000
+sinh -0x1p-1000
+sinh 0x1p-10000
+sinh -0x1p-10000
 sinh 22
 sinh 23
 sinh 24