#include <math.h>
#include <stdint.h>

float frexpf(float x, int* e) {
  union {
    float f;
    uint32_t i;
  } y = {x};
  int ee = y.i >> 23 & 0xff;

  if (!ee) {
    if (x) {
      x = frexpf(x * 0x1p64, e);
      *e -= 64;
    } else
      *e = 0;
    return x;
  } else if (ee == 0xff) {
    return x;
  }

  *e = ee - 0x7e;
  y.i &= 0x807ffffful;
  y.i |= 0x3f000000ul;
  return y.f;
}