src/math/modff.c

   1 #include "libm.h"
   2
   3 float modff(float x, float *iptr)
   4 {
   5         union {float x; uint32_t n;} u = {x};
   6         uint32_t mask;
   7         int e;
   8
   9         e = (int)(u.n>>23 & 0xff) - 0x7f;
  10
  11         /* no fractional part */
  12         if (e >= 23) {
  13                 *iptr = x;
  14                 if (e == 0x80 && u.n<<9 != 0) { /* nan */
  15                         return x;
  16                 }
  17                 u.n &= 0x80000000;
  18                 return u.x;
  19         }
  20         /* no integral part */
  21         if (e < 0) {
  22                 u.n &= 0x80000000;
  23                 *iptr = u.x;
  24                 return x;
  25         }
  26
  27         mask = 0x007fffff>>e;
  28         if ((u.n & mask) == 0) {
  29                 *iptr = x;
  30                 u.n &= 0x80000000;
  31                 return u.x;
  32         }
  33         u.n &= ~mask;
  34         *iptr = u.x;
  35         STRICT_ASSIGN(float, x, x - *iptr);
  36         return x;
  37 }