RRO.EX2

unsigned int approx(float f) {
    float ax = fabs(f);
    unsigned int extra = 1;
    if (ax >= 128.f) {
        ax = 129.f;
        extra = 0;
    }
    double w;
    float ax2 = modf(ax, &w);
    unsigned int sign = (f < 0.0f) << 31;
    return sign | ((((int)ax) << 23) | (int)(ax2 * 8388607.f) + extra);
}