src/math/ilogbl.c

   1 #include <limits.h>
   2 #include "libm.h"
   3
   4 #if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024
   5 int ilogbl(long double x)
   6 {
   7         return ilogb(x);
   8 }
   9 #elif LDBL_MANT_DIG == 64 && LDBL_MAX_EXP == 16384
  10 int ilogbl(long double x)
  11 {
  12         union ldshape u = {x};
  13         uint64_t m = u.bits.m;
  14         int e = u.bits.exp;
  15
  16         if (!e) {
  17                 if (m == 0) {
  18                         FORCE_EVAL(0/0.0f);
  19                         return FP_ILOGB0;
  20                 }
  21                 /* subnormal x */
  22                 for (e = -0x3fff+1; m < (uint64_t)1<<63; e--, m<<=1);
  23                 return e;
  24         }
  25         if (e == 0x7fff) {
  26                 FORCE_EVAL(0/0.0f);
  27                 /* in ld80 msb is set in inf */
  28                 return m & (uint64_t)-1>>1 ? FP_ILOGBNAN : INT_MAX;
  29         }
  30         return e - 0x3fff;
  31 }
  32 #endif