From 8a9fa0086dfe1b2b9d828712ef41fc7dd4b8f94c Mon Sep 17 00:00:00 2001 From: Joseph Myers Date: Thu, 18 Feb 2016 21:40:39 +0000 Subject: Fix ldbl-128ibm ceill for non-default rounding modes (bug 19592). The ldbl-128ibm implementation of ceill is only correct in round-to-nearest mode (in other modes, there are incorrect results and overflow exceptions in some cases). It is also unnecessarily complicated, rounding both high and low parts to the nearest integer and then adjusting for the semantics of ceil, when it seems more natural to take the ceiling of the high part (__ceil optimized inline versions can be used), and that of the low part if the high part is an integer, as was done for floorl. This patch makes it use that simpler approach. Tested for powerpc. [BZ #19592] * sysdeps/ieee754/ldbl-128ibm/s_ceill.c (__ceill): Use __ceil on high and low parts then use ldbl_canonicalize_int if needed. --- ChangeLog | 4 ++++ 1 file changed, 4 insertions(+) (limited to 'ChangeLog') diff --git a/ChangeLog b/ChangeLog index 6f79779923..e622d9e14a 100644 --- a/ChangeLog +++ b/ChangeLog @@ -1,5 +1,9 @@ 2016-02-18 Joseph Myers + [BZ #19592] + * sysdeps/ieee754/ldbl-128ibm/s_ceill.c (__ceill): Use __ceil on + high and low parts then use ldbl_canonicalize_int if needed. + [BZ #17899] * sysdeps/ieee754/ldbl-128ibm/math_ldbl.h (ldbl_canonicalize_int): New function. -- cgit 1.4.1