From 25dbcb277a9dea5f241c0b888a30b9341310d941 Mon Sep 17 00:00:00 2001 From: Andreas Schwab Date: Fri, 25 May 2012 11:57:33 +0200 Subject: Optimize handling of denormals in logb/logbf/logbl --- sysdeps/ieee754/ldbl-128ibm/s_logbl.c | 7 ++----- 1 file changed, 2 insertions(+), 5 deletions(-) (limited to 'sysdeps/ieee754/ldbl-128ibm/s_logbl.c') diff --git a/sysdeps/ieee754/ldbl-128ibm/s_logbl.c b/sysdeps/ieee754/ldbl-128ibm/s_logbl.c index 678b6ca..92ce2c1 100644 --- a/sysdeps/ieee754/ldbl-128ibm/s_logbl.c +++ b/sysdeps/ieee754/ldbl-128ibm/s_logbl.c @@ -30,7 +30,7 @@ __logbl (long double x) GET_LDOUBLE_WORDS64 (hx, lx, x); hx &= 0x7fffffffffffffffLL; /* high |x| */ - if ((hx | (lx & 0x7fffffffffffffffLL)) == 0) + if (hx == 0) return -1.0 / fabs (x); if (hx >= 0x7ff0000000000000LL) return x * x; @@ -38,10 +38,7 @@ __logbl (long double x) { /* POSIX specifies that denormal number is treated as though it were normalized. */ - int m1 = (hx == 0) ? 0 : __builtin_clzll (hx); - int m2 = (lx == 0) ? 0 : __builtin_clzll (lx); - int ma = (m1 == 0) ? m2 + 64 : m1; - return -1022.0 + (long double)(11 - ma); + rhx -= __builtin_clzll (hx) - 12; } return (long double) (rhx - 1023); } -- cgit v1.1