diff options
Diffstat (limited to 'src/math/__fpclassify.c')
-rw-r--r-- | src/math/__fpclassify.c | 14 |
1 files changed, 14 insertions, 0 deletions
diff --git a/src/math/__fpclassify.c b/src/math/__fpclassify.c new file mode 100644 index 0000000..1605110 --- /dev/null +++ b/src/math/__fpclassify.c @@ -0,0 +1,14 @@ +#include <stdint.h> +#include <math.h> + +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; + return FP_NORMAL; +} |