comparison ac3dec.c @ 5313:4bd60d3bab7e libavcodec

AC-3 decoder, soc revision 53, Aug 17 08:53:44 2006 UTC by cloud9 Converted the window from double to float. Now sound produced is as good as sound produced by liba52.

5312:2c875e6274d5

/**
 * Generate a Kaiser Window.
 */
static void
k_window_init(int alpha, double *window, int n, int iter)
{
    int j, k;
    double a, x;
    a = alpha * M_PI / n;
    a = a*a;
    for(k=0; k<n; k++) {
        x = k * (n - k) * a;
        window[k] = 1.0;

 * @param window array to fill with window values
 * @param n      length of the window
 * @param iter   number of iterations to use in BesselI0
 */
static void
kbd_window_init(int alpha, double *window, int n, int iter)
{
    int k, n2;
    double *kwindow;

    n2 = n >> 1;
    kwindow = &window[n2];
    k_window_init(alpha, kwindow, n2, iter);
    window[0] = kwindow[0];

}

static void ac3_tables_init(void)
{
    int i, j, k, l, v;
    float alpha;
    /* compute bndtab and masktab from bandsz */
    k = 0;
    l = 0;
    for(i=0;i<50;i++) {
        bndtab[i] = l;

    generate_quantizers_table_2(l11_quantizers_1, 11, 11, 11, 128);
    generate_quantizers_table_3(l11_quantizers_2, 11, 11, 11, 128);

    //for level-15 quantizers
    generate_quantizers_table(l15_quantizers, 15, 15);

    /* Twiddle Factors for IMDCT. */
    for(i = 0; i < N / 4; i++) {
        alpha = 2 * M_PI * (8 * i + 1) / (8 * N);
        x_cos1[i] = -cos(alpha);
        x_sin1[i] = -sin(alpha);
    }

    for (i = 0; i < N / 8; i++) {
        alpha = 2 * M_PI * (8 * i + 1) / (4 * N);
        x_cos2[i] = -cos(alpha);
        x_sin2[i] = -sin(alpha);
    }

    /* Kaiser-Bessel derived window. */
    kbd_window_init(5, window, 256, 100);
}

static int ac3_decode_frame(AVCodecContext * avctx, void *data, int *data_size, uint8_t *buf, int buf_size)
{
    AC3DecodeContext *ctx = (AC3DecodeContext *)avctx->priv_data;
    int frame_start;
    int16_t *out_samples = (int16_t *)data;
    int i, j, k, value, fs_58;
    uint16_t crc1;

    av_log(NULL, AV_LOG_INFO, "decoding frame %d buf_size = %d\n", frame_count++, buf_size);

    //Synchronize the frame.
    frame_start = ac3_synchronize(buf, buf_size);

5313:4bd60d3bab7e

/**
 * Generate a Kaiser Window.
 */
static void
k_window_init(int alpha, float *window, int n, int iter)
{
    int j, k;
    float a, x;
    a = alpha * M_PI / n;
    a = a*a;
    for(k=0; k<n; k++) {
        x = k * (n - k) * a;
        window[k] = 1.0;

 * @param window array to fill with window values
 * @param n      length of the window
 * @param iter   number of iterations to use in BesselI0
 */
static void
kbd_window_init(int alpha, float *window, int n, int iter)
{
    int k, n2;
    float *kwindow;

    n2 = n >> 1;
    kwindow = &window[n2];
    k_window_init(alpha, kwindow, n2, iter);
    window[0] = kwindow[0];

}

static void ac3_tables_init(void)
{
    int i, j, k, l, v;
    /* compute bndtab and masktab from bandsz */
    k = 0;
    l = 0;
    for(i=0;i<50;i++) {
        bndtab[i] = l;

    generate_quantizers_table_2(l11_quantizers_1, 11, 11, 11, 128);
    generate_quantizers_table_3(l11_quantizers_2, 11, 11, 11, 128);

    //for level-15 quantizers
    generate_quantizers_table(l15_quantizers, 15, 15);

    /* Kaiser-Bessel derived window. */
    kbd_window_init(5, window, 256, 100);
}

static int ac3_decode_frame(AVCodecContext * avctx, void *data, int *data_size, uint8_t *buf, int buf_size)
{
    AC3DecodeContext *ctx = (AC3DecodeContext *)avctx->priv_data;
    int frame_start;
    int16_t *out_samples = (int16_t *)data;
    int i, j, k, value;

    av_log(NULL, AV_LOG_INFO, "decoding frame %d buf_size = %d\n", frame_count++, buf_size);

    //Synchronize the frame.
    frame_start = ac3_synchronize(buf, buf_size);