diff options
Diffstat (limited to 'lib/mlibc/options/ansi/musl-generic-math/scalbn.c')
-rw-r--r-- | lib/mlibc/options/ansi/musl-generic-math/scalbn.c | 33 |
1 files changed, 0 insertions, 33 deletions
diff --git a/lib/mlibc/options/ansi/musl-generic-math/scalbn.c b/lib/mlibc/options/ansi/musl-generic-math/scalbn.c deleted file mode 100644 index 182f561..0000000 --- a/lib/mlibc/options/ansi/musl-generic-math/scalbn.c +++ /dev/null @@ -1,33 +0,0 @@ -#include <math.h> -#include <stdint.h> - -double scalbn(double x, int n) -{ - union {double f; uint64_t i;} u; - double_t y = x; - - if (n > 1023) { - y *= 0x1p1023; - n -= 1023; - if (n > 1023) { - y *= 0x1p1023; - n -= 1023; - if (n > 1023) - n = 1023; - } - } else if (n < -1022) { - /* make sure final n < -53 to avoid double - rounding in the subnormal range */ - y *= 0x1p-1022 * 0x1p53; - n += 1022 - 53; - if (n < -1022) { - y *= 0x1p-1022 * 0x1p53; - n += 1022 - 53; - if (n < -1022) - n = -1022; - } - } - u.i = (uint64_t)(0x3ff+n)<<52; - x = y * u.f; - return x; -} |