first commit of the new libm!

[musl] / src / math / __fpclassify.c

diff --git a/src/math/__fpclassify.c b/src/math/__fpclassify.c
index 1605110..c9dd027 100644 (file)
--- a/src/math/__fpclassify.c
+++ b/src/math/__fpclassify.c
@@ -1,14 +1,10 @@
-#include <stdint.h>
-#include <math.h>
+#include "libm.h"

 
 

-int __fpclassify(double __x)
+int __fpclassify(double x)

 {
 {

-       union {
-               double __d;
-               __uint64_t __i;
-       } __y = { __x };
-       int __ee = __y.__i>>52 & 0x7ff;
-       if (!__ee) return __y.__i<<1 ? FP_SUBNORMAL : FP_ZERO;
-       if (__ee==0x7ff) return __y.__i<<12 ? FP_NAN : FP_INFINITE;
+       union dshape u = { x };
+       int e = u.bits>>52 & 0x7ff;
+       if (!e) return u.bits<<1 ? FP_SUBNORMAL : FP_ZERO;
+       if (e==0x7ff) return u.bits<<12 ? FP_NAN : FP_INFINITE;

        return FP_NORMAL;
 }
        return FP_NORMAL;
 }