From d79b56892651fcfc1ae0f44d3d29d0fa03538ec1 Mon Sep 17 00:00:00 2001 From: Joseph Myers Date: Tue, 6 Oct 2015 15:14:26 +0000 Subject: Fix ldbl-128ibm log1pl (-1) sign of infinity (bug 19076). The ldbl-128ibm implementation of log1pl produces an infinity with the wrong sign for log1pl (-1) in FE_DOWNWARD mode. This patch fixes this by changing a division (-1.0L / (x - x)) (incorrect in FE_DOWNWARD mode) to (-1.0L / 0.0L) (correct in all rounding modes). Tested for powerpc. [BZ #19076] * sysdeps/ieee754/ldbl-128ibm/s_log1pl.c (__log1pl): Divide by constant 0.0L when computing infinite result. --- ChangeLog | 6 ++++++ 1 file changed, 6 insertions(+) (limited to 'ChangeLog') diff --git a/ChangeLog b/ChangeLog index 6cb9b73dc6..17388a7927 100644 --- a/ChangeLog +++ b/ChangeLog @@ -1,3 +1,9 @@ +2015-10-06 Joseph Myers + + [BZ #19076] + * sysdeps/ieee754/ldbl-128ibm/s_log1pl.c (__log1pl): Divide by + constant 0.0L when computing infinite result. + 2015-10-06 Florian Weimer [BZ #10432] -- cgit 1.4.1