src/math/__fpclassifyl.c


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16

#include <stdint.h>
#include <math.h>

/* FIXME: move this to arch-specific file */
int __fpclassifyl(long double __x)
{
	union {
		long double __ld;
		__uint16_t __hw[5];
		int64_t __m;
	} __y = { __x };
	int __ee = __y.__hw[4]&0x7fff;
	if (!__ee) return __y.__m ? FP_SUBNORMAL : FP_ZERO;
	if (__ee==0x7fff) return __y.__m ? FP_NAN : FP_INFINITE;
	return __y.__m < 0 ? FP_NORMAL : FP_NAN;
}