#include "libm.h"

float floorf(float x) {
  union {
    float f;
    uint32_t i;
  } u = {x};
  int e = (int)(u.i >> 23 & 0xff) - 0x7f;
  uint32_t m;

  if (e >= 23)
    return x;
  if (e >= 0) {
    m = 0x007fffff >> e;
    if ((u.i & m) == 0)
      return x;
    FORCE_EVAL(x + 0x1p120f);
    if (u.i >> 31)
      u.i += m;
    u.i &= ~m;
  } else {
    FORCE_EVAL(x + 0x1p120f);
    if (u.i >> 31 == 0)
      u.i = 0;
    else if (u.i << 1)
      u.f = -1.0;
  }
  return u.f;
}