X-Git-Url: http://nsz.repo.hu/git/?a=blobdiff_plain;f=src%2Fmath%2Fscalbnf.c;h=0b62c3c71f85d4b72a21a9b39752fb77ce0caf10;hb=04b8360adbb6487f61aa0c00e53ec3a90a5a0d29;hp=f94b5d5992bb0ad8ed67e19323e68853201df045;hpb=8bb181622222f2ee3462c8b021bcae4fcdbbd37a;p=musl diff --git a/src/math/scalbnf.c b/src/math/scalbnf.c index f94b5d59..0b62c3c7 100644 --- a/src/math/scalbnf.c +++ b/src/math/scalbnf.c @@ -1,33 +1,31 @@ -#include "libm.h" +#include +#include float scalbnf(float x, int n) { - float scale; + union {float f; uint32_t i;} u; + float_t y = x; if (n > 127) { - x *= 0x1p127f; + y *= 0x1p127f; n -= 127; if (n > 127) { - x *= 0x1p127f; + y *= 0x1p127f; n -= 127; - if (n > 127) { - STRICT_ASSIGN(float, x, x * 0x1p127f); - return x; - } + if (n > 127) + n = 127; } } else if (n < -126) { - x *= 0x1p-126f; + y *= 0x1p-126f; n += 126; if (n < -126) { - x *= 0x1p-126f; + y *= 0x1p-126f; n += 126; - if (n < -126) { - STRICT_ASSIGN(float, x, x * 0x1p-126f); - return x; - } + if (n < -126) + n = -126; } } - SET_FLOAT_WORD(scale, (uint32_t)(0x7f+n)<<23); - STRICT_ASSIGN(float, x, x * scale); + u.i = (uint32_t)(0x7f+n)<<23; + x = y * u.f; return x; }