nsz Git - musl/blob - src/math/ilogb.c

   1 #include <limits.h>
   2 #include "libm.h"
   3
   4 int ilogb(double x)
   5 {
   6         union dshape u = {x};
   7         int e = u.bits>>52 & 0x7ff;
   8
   9         if (!e) {
  10                 u.bits <<= 12;
  11                 if (u.bits == 0)
  12                         return FP_ILOGB0;
  13                 /* subnormal x */
  14                 for (e = -0x3ff; u.bits < (uint64_t)1<<63; e--, u.bits<<=1);
  15                 return e;
  16         }
  17         if (e == 0x7ff)
  18                 return u.bits<<12 ? FP_ILOGBNAN : INT_MAX;
  19         return e - 0x3ff;
  20 }