summaryrefslogtreecommitdiff
path: root/lib/mlibc/options/ansi/musl-generic-math/scalbn.c
diff options
context:
space:
mode:
Diffstat (limited to 'lib/mlibc/options/ansi/musl-generic-math/scalbn.c')
-rw-r--r--lib/mlibc/options/ansi/musl-generic-math/scalbn.c33
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;
-}