blob: 2a961b1c1aaae0d74805c2638ac4ebbde7b4a5b8 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
|
#include <math.h> // for FP_INFINITE, FP_NAN, FP_NORMAL, FP_SUBNORMAL
#include <stdint.h> // for uint64_t
int __fpclassify(double x)
{
union {
double f;
uint64_t i;
} u = { x };
int e = u.i >> 52 & 0x7ff;
if (!e)
return u.i << 1 ? FP_SUBNORMAL : FP_ZERO;
if (e == 0x7ff)
return u.i << 12 ? FP_NAN : FP_INFINITE;
return FP_NORMAL;
}
|
