first commit of the new libm!

[musl] / src / math / frexpl.c

diff --git a/src/math/frexpl.c b/src/math/frexpl.c
new file mode 100644 (file)
index 0000000..f9d90a6
--- /dev/null
+++ b/src/math/frexpl.c
@@ -0,0 +1,37 @@
+#include <math.h>
+#include <stdint.h>
+#include <float.h>
+
+#if LDBL_MANT_DIG == 64 && LDBL_MAX_EXP == 16384
+
+/* This version is for 80-bit little endian long double */
+
+long double frexpl(long double x, int *e)
+{
+       union { long double ld; uint16_t hw[5]; } y = { x };
+       int ee = y.hw[4]&0x7fff;
+
+       if (!ee) {
+               if (x) {
+                       x = frexpl(x*0x1p64, e);
+                       *e -= 64;
+               } else *e = 0;
+               return x;
+       } else if (ee == 0x7fff) {
+               return x;
+       }
+
+       *e = ee - 0x3ffe;
+       y.hw[4] &= 0x8000;
+       y.hw[4] |= 0x3ffe;
+       return y.ld;
+}
+
+#else
+
+long double frexpl(long double x, int *e)
+{
+       return frexp(x, e);
+}
+
+#endif