diff options
| author | Szabolcs Nagy <nsz@port70.net> | 2013-09-04 15:52:54 +0000 |
|---|---|---|
| committer | Szabolcs Nagy <nsz@port70.net> | 2013-09-05 11:30:08 +0000 |
| commit | 2eaed464e2080d8321d3903b71086a1ecfc4ee4a (patch) | |
| tree | b84e65faddf4f1042f50a7d50338ef0656c39996 /src/math/scalbnf.c | |
| parent | 34660d73bd0db29469d2758e1b48d2360edf3a2f (diff) | |
| download | musl-2eaed464e2080d8321d3903b71086a1ecfc4ee4a.tar.gz | |
math: use float_t and double_t in scalbnf and scalbn
remove STRICT_ASSIGN (c99 semantics is assumed) and use the conventional
union to prepare the scaling factor (so libm.h is no longer needed)
Diffstat (limited to 'src/math/scalbnf.c')
| -rw-r--r-- | src/math/scalbnf.c | 18 |
1 files changed, 10 insertions, 8 deletions
diff --git a/src/math/scalbnf.c b/src/math/scalbnf.c index c0eeaf84..0b62c3c7 100644 --- a/src/math/scalbnf.c +++ b/src/math/scalbnf.c @@ -1,29 +1,31 @@ -#include "libm.h" +#include <math.h> +#include <stdint.h> 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) 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) 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; } |