From 4cec31fc23b28271a2ed0b290e65699dffac3fee Mon Sep 17 00:00:00 2001 From: Szabolcs Nagy Date: Mon, 2 Sep 2013 23:33:20 +0000 Subject: math: ilogb cleanup * consistent code style * explicit union instead of typedef for double and float bit access * turn FENV_ACCESS ON to make 0/0.0f raise invalid flag * (untested) ld128 version of ilogbl (used by logbl which has ld128 support) --- src/math/ilogb.c | 14 ++++++++------ 1 file changed, 8 insertions(+), 6 deletions(-) (limited to 'src/math/ilogb.c') diff --git a/src/math/ilogb.c b/src/math/ilogb.c index 5a1819d8..64d40154 100644 --- a/src/math/ilogb.c +++ b/src/math/ilogb.c @@ -3,22 +3,24 @@ int ilogb(double x) { - union dshape u = {x}; - int e = u.bits>>52 & 0x7ff; + #pragma STDC FENV_ACCESS ON + union {double f; uint64_t i;} u = {x}; + uint64_t i = u.i; + int e = i>>52 & 0x7ff; if (!e) { - u.bits <<= 12; - if (u.bits == 0) { + i <<= 12; + if (i == 0) { FORCE_EVAL(0/0.0f); return FP_ILOGB0; } /* subnormal x */ - for (e = -0x3ff; u.bits < (uint64_t)1<<63; e--, u.bits<<=1); + for (e = -0x3ff; i>>63 == 0; e--, i<<=1); return e; } if (e == 0x7ff) { FORCE_EVAL(0/0.0f); - return u.bits<<12 ? FP_ILOGBNAN : INT_MAX; + return i<<12 ? FP_ILOGBNAN : INT_MAX; } return e - 0x3ff; } -- cgit v1.2.1