comparison aac.c @ 9179:4172c9ef331a libavcodec

AAC: IEEE-754 type punning for 16-bit floating point rounding.

9178:c1e6aaea75f2

#include <assert.h>
#include <errno.h>
#include <math.h>
#include <string.h>

static VLC vlc_scalefactors;
static VLC vlc_spectral[11];

    }
    return 0;
}

static av_always_inline float flt16_round(float pf) {
    int exp;
    pf = frexpf(pf, &exp);
    pf = ldexpf(roundf(ldexpf(pf, 8)), exp-8);
    return pf;
}

static av_always_inline float flt16_even(float pf) {
    int exp;
    pf = frexpf(pf, &exp);
    pf = ldexpf(rintf(ldexpf(pf, 8)), exp-8);
    return pf;
}

static av_always_inline float flt16_trunc(float pf) {
    int exp;
    pf = frexpf(pf, &exp);
    pf = ldexpf(truncf(ldexpf(pf, 8)), exp-8);
    return pf;
}

static void predict(AACContext * ac, PredictorState * ps, float* coef, int output_enable) {
    const float a     = 0.953125; // 61.0/64
    const float alpha = 0.90625;  // 29.0/32

9179:4172c9ef331a

#include <assert.h>
#include <errno.h>
#include <math.h>
#include <string.h>

union float754 { float f; uint32_t i; };

static VLC vlc_scalefactors;
static VLC vlc_spectral[11];

    }
    return 0;
}

static av_always_inline float flt16_round(float pf) {
    union float754 tmp;
    tmp.f = pf;
    tmp.i = (tmp.i + 0x00008000U) & 0xFFFF0000U;
    return tmp.f;
}

static av_always_inline float flt16_even(float pf) {
    union float754 tmp;
    tmp.f = pf;
    tmp.i = (tmp.i + 0x00007FFFU + (tmp.i & 0x00010000U>>16)) & 0xFFFF0000U;
    return tmp.f;
}

static av_always_inline float flt16_trunc(float pf) {
    union float754 pun;
    pun.f = pf;
    pun.i &= 0xFFFF0000U;
    return pun.f;
}

static void predict(AACContext * ac, PredictorState * ps, float* coef, int output_enable) {
    const float a     = 0.953125; // 61.0/64
    const float alpha = 0.90625;  // 29.0/32