From 613c92b3b59df6a06784cde1d4f410cef0b6da96 Mon Sep 17 00:00:00 2001 From: Joseph Myers Date: Wed, 9 Mar 2016 00:30:59 +0000 Subject: Fix ldbl-128ibm nearbyintl in non-default rounding modes (bug 19790). The ldbl-128ibm implementation of nearbyintl uses logic that only works in round-to-nearest mode. This contrasts with rintl, which works in all rounding modes. Now, arguably nearbyintl could simply be aliased to rintl, given that spurious "inexact" is generally allowed for ldbl-128ibm, even for the underlying arithmetic operations. But given that the only point of nearbyintl is to avoid "inexact", this patch follows the more conservative approach of adding conditionals to the rintl implementation to make it suitable for use to implement nearbyintl, then builds it for nearbyintl with USE_AS_NEARBYINTL defined. The test test-nearbyint-except-2 shows up issues when traps on "inexact" are enabled, which turn out to be problems with the powerpc fenv_private.h implementation (two functions that should disable exception traps potentially failing to do so in some cases); this patch duly fixes that as well (I don't see any other existing cases where this would be user-visible; there isn't much use of *_NOEX, *hold* etc. in libm that requires exceptions to be discarded and not trapped on). Tested for powerpc. [BZ #19790] * sysdeps/ieee754/ldbl-128ibm/s_rintl.c [USE_AS_NEARBYINTL] (rintl): Define as macro. [USE_AS_NEARBYINTL] (__rintl): Likewise. (__rintl) [USE_AS_NEARBYINTL]: Use SET_RESTORE_ROUND_NOEX instead of fesetround. Ensure results are evaluated before end of scope. * sysdeps/ieee754/ldbl-128ibm/s_nearbyintl.c: Define USE_AS_NEARBYINTL and include s_rintl.c. * sysdeps/powerpc/fpu/fenv_private.h (libc_feholdsetround_ppc): Disable exception traps in new environment. (libc_feholdsetround_ppc_ctx): Likewise. --- sysdeps/ieee754/ldbl-128ibm/s_nearbyintl.c | 109 +---------------------------- 1 file changed, 2 insertions(+), 107 deletions(-) (limited to 'sysdeps/ieee754/ldbl-128ibm/s_nearbyintl.c') diff --git a/sysdeps/ieee754/ldbl-128ibm/s_nearbyintl.c b/sysdeps/ieee754/ldbl-128ibm/s_nearbyintl.c index 08134ed..dfdefe3 100644 --- a/sysdeps/ieee754/ldbl-128ibm/s_nearbyintl.c +++ b/sysdeps/ieee754/ldbl-128ibm/s_nearbyintl.c @@ -17,110 +17,5 @@ License along with the GNU C Library; if not, see . */ -/* This has been coded in assembler because GCC makes such a mess of it - when it's coded in C. */ - -#include -#include -#include -#include -#include -#include - - -long double -__nearbyintl (long double x) -{ - fenv_t env; - static const long double TWO52 = 4503599627370496.0L; - union ibm_extended_long_double u; - u.ld = x; - - if (!isfinite (u.d[0].d)) - return x; - else if (fabs (u.d[0].d) < TWO52) - { - double xh = u.d[0].d; - double high = u.d[0].d; - feholdexcept (&env); - if (high > 0.0) - { - high += TWO52; - high -= TWO52; - if (high == -0.0) high = 0.0; - } - else if (high < 0.0) - { - high -= TWO52; - high += TWO52; - if (high == 0.0) high = -0.0; - } - if (u.d[1].d > 0.0 && (xh - high == 0.5)) - high += 1.0; - else if (u.d[1].d < 0.0 && (-(xh - high) == 0.5)) - high -= 1.0; - u.d[0].d = high; - u.d[1].d = 0.0; - math_force_eval (u.d[0]); - math_force_eval (u.d[1]); - fesetenv (&env); - } - else if (fabs (u.d[1].d) < TWO52 && u.d[1].d != 0.0) - { - double high = u.d[0].d, low = u.d[1].d, tau; - /* In this case we have to round the low double and handle any - adjustment to the high double that may be caused by rounding - (up). This is complicated by the fact that the high double - may already be rounded and the low double may have the - opposite sign to compensate. */ - feholdexcept (&env); - if (u.d[0].d > 0.0) - { - if (u.d[1].d > 0.0) - { - /* If the high/low doubles are the same sign then simply - round the low double. */ - } - else if (u.d[1].d < 0.0) - { - /* Else the high double is pre rounded and we need to - adjust for that. */ - - tau = __nextafter (u.d[0].d, 0.0); - tau = (u.d[0].d - tau) * 2.0; - high -= tau; - low += tau; - } - low += TWO52; - low -= TWO52; - } - else if (u.d[0].d < 0.0) - { - if (u.d[1].d < 0.0) - { - /* If the high/low doubles are the same sign then simply - round the low double. */ - } - else if (u.d[1].d > 0.0) - { - /* Else the high double is pre rounded and we need to - adjust for that. */ - tau = __nextafter (u.d[0].d, 0.0); - tau = (u.d[0].d - tau) * 2.0; - high -= tau; - low += tau; - } - low = TWO52 - low; - low = -(low - TWO52); - } - u.d[0].d = high + low; - u.d[1].d = high - u.d[0].d + low; - math_force_eval (u.d[0]); - math_force_eval (u.d[1]); - fesetenv (&env); - } - - return u.ld; -} - -long_double_symbol (libm, __nearbyintl, nearbyintl); +#define USE_AS_NEARBYINTL +#include "s_rintl.c" -- cgit v1.1