1/**
 2 * This file is part of the mingw-w64 runtime package.
 3 * No warranty is given; refer to the file DISCLAIMER within this package.
 4 */
 5#define __CRT__NO_INLINE
 6#include <math.h>
 7
 8int __fpclassifyl (long double _x){
 9#if __SIZEOF_LONG_DOUBLE__ == __SIZEOF_DOUBLE__
10    return __fpclassify(_x);
11#elif defined(_AMD64_) || defined(__x86_64__)
12    __mingw_ldbl_type_t hlp;
13    unsigned int e;
14    hlp.x = _x;
15    e = hlp.lh.sign_exponent & 0x7fff;
16    if (!e)
17      {
18        unsigned int h = hlp.lh.high;
19        if (!(hlp.lh.low | h))
20          return FP_ZERO;
21        else if (!(h & 0x80000000))
22          return FP_SUBNORMAL;
23      }
24    else if (e == 0x7fff)
25      return (((hlp.lh.high & 0x7fffffff) | hlp.lh.low) == 0 ?
26	      FP_INFINITE : FP_NAN);
27    return FP_NORMAL;
28#elif defined(__i386__) || defined(_X86_)
29  unsigned short sw;
30  __asm__ __volatile__ (
31	"fxam; fstsw %%ax;"
32	: "=a" (sw)
33	: "t" (_x)
34	);
35  return sw & (FP_NAN | FP_NORMAL | FP_ZERO );
36#endif
37}