nsz Git - libm/blob - src/math/__fpclassifyl.c

   1 #include "libm.h"
   2
   3 #if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024
   4
   5 #elif LDBL_MANT_DIG == 64 && LDBL_MAX_EXP == 16384
   6 int __fpclassifyl(long double x)
   7 {
   8         union ldshape u = { x };
   9         int e = u.bits.exp;
  10         if (!e)
  11                 return u.bits.m ? FP_SUBNORMAL : FP_ZERO;
  12         if (e == 0x7fff)
  13                 return u.bits.m & ~((uint64_t)1<<63) ? FP_NAN : FP_INFINITE;
  14         return u.bits.m & (uint64_t)1<<63 ? FP_NORMAL : FP_NAN;
  15 }
  16 #elif LDBL_MANT_DIG == 113 && LDBL_MAX_EXP == 16384
  17 int __fpclassifyl(long double x)
  18 {
  19         union ldshape u = { x };
  20         int e = u.bits.exp;
  21         if (!e)
  22                 return u.bits.mlo | u.bits.mhi ? FP_SUBNORMAL : FP_ZERO;
  23         if (e == 0x7fff)
  24                 return u.bits.mlo | u.bits.mhi ? FP_NAN : FP_INFINITE;
  25         return FP_NORMAL;
  26 }
  27 #endif