about summary refs log tree commit diff
path: root/sysdeps/i386/i686/fpu/multiarch
diff options
context:
space:
mode:
authorJoseph Myers <joseph@codesourcery.com>2016-03-24 01:32:52 +0000
committerJoseph Myers <joseph@codesourcery.com>2016-03-24 01:32:52 +0000
commitc898991d8bcfacc825097ba389ffccc5367c2b2d (patch)
treeef1f3fec885e5e5ed89048016ffe8b8b2c0e8393 /sysdeps/i386/i686/fpu/multiarch
parent7e1ff08c260ae105208edb4e778e50525dff05d3 (diff)
downloadglibc-c898991d8bcfacc825097ba389ffccc5367c2b2d.tar.gz
glibc-c898991d8bcfacc825097ba389ffccc5367c2b2d.tar.xz
glibc-c898991d8bcfacc825097ba389ffccc5367c2b2d.zip
Fix x86_64 / x86 powl inaccuracy for integer exponents (bug 19848).
Bug 19848 reports cases where powl on x86 / x86_64 has error
accumulation, for small integer exponents, larger than permitted by
glibc's accuracy goals, at least in some rounding modes.  This patch
further restricts the exponent range for which the
small-integer-exponent logic is used to limit the possible error
accumulation.

Tested for x86_64 and x86 and ulps updated accordingly.

	[BZ #19848]
	* sysdeps/i386/fpu/e_powl.S (p3): Rename to p2 and change value
	from 8 to 4.
	(__ieee754_powl): Compare integer exponent against 4 not 8.
	* sysdeps/x86_64/fpu/e_powl.S (p3): Rename to p2 and change value
	from 8 to 4.
	(__ieee754_powl): Compare integer exponent against 4 not 8.
	* math/auto-libm-test-in: Add more tests of pow.
	* math/auto-libm-test-out: Regenerated.
	* sysdeps/i386/i686/fpu/multiarch/libm-test-ulps: Update.
	* sysdeps/x86_64/fpu/libm-test-ulps: Likewise.
Diffstat (limited to 'sysdeps/i386/i686/fpu/multiarch')
-rw-r--r--sysdeps/i386/i686/fpu/multiarch/libm-test-ulps8
1 files changed, 4 insertions, 4 deletions
diff --git a/sysdeps/i386/i686/fpu/multiarch/libm-test-ulps b/sysdeps/i386/i686/fpu/multiarch/libm-test-ulps
index bbb644a591..585be7880d 100644
--- a/sysdeps/i386/i686/fpu/multiarch/libm-test-ulps
+++ b/sysdeps/i386/i686/fpu/multiarch/libm-test-ulps
@@ -1903,14 +1903,14 @@ ldouble: 4
 Function: "pow_towardzero":
 double: 1
 idouble: 1
-ildouble: 1
-ldouble: 1
+ildouble: 4
+ldouble: 4
 
 Function: "pow_upward":
 double: 1
 idouble: 1
-ildouble: 2
-ldouble: 2
+ildouble: 4
+ldouble: 4
 
 Function: "sin":
 ildouble: 1