use alternate formula for acos asm to avoid loss of precision

[oweals/musl.git] / src / math / __fpclassifyf.c

diff --git a/src/math/__fpclassifyf.c b/src/math/__fpclassifyf.c
index bf59d0d4709f066ce612e6f5ac5448b1aeb2f157..8149087b0d3f4e80c843acbaffad2875c41d1947 100644 (file)
--- a/src/math/__fpclassifyf.c
+++ b/src/math/__fpclassifyf.c
@@ -1,14 +1,10 @@
-#include <stdint.h>
-#include <math.h>
+#include "libm.h"
 
-int __fpclassifyf(float __x)
+int __fpclassifyf(float x)
 {
-       union {
-               float __f;
-               __uint32_t __i;
-       } __y = { __x };
-       int __ee = __y.__i>>23 & 0xff;
-       if (!__ee) return __y.__i<<1 ? FP_SUBNORMAL : FP_ZERO;
-       if (__ee==0xff) return __y.__i<<9 ? FP_NAN : FP_INFINITE;
+       union fshape u = { x };
+       int e = u.bits>>23 & 0xff;
+       if (!e) return u.bits<<1 ? FP_SUBNORMAL : FP_ZERO;
+       if (e==0xff) return u.bits<<9 ? FP_NAN : FP_INFINITE;
        return FP_NORMAL;
 }