comparison sipr.c @ 12465:1921e24d5886 libavcodec

Move lsp2lpc_sipr() function to common code so it can be reused in a AMRWB decoder. Patch by Marcelo Galv���o P���voa.

12464:2dd67ed2f947

        parms->gc_index[i]          = get_bits(pgb, p->gc_index_bits);
    }
}

static void lsp2lpc_sipr(const double *lsp, float *Az)
{
    int lp_half_order = LP_FILTER_ORDER >> 1;
    double buf[(LP_FILTER_ORDER >> 1) + 1];
    double pa[(LP_FILTER_ORDER >> 1) + 1];
    double *qa = buf + 1;
    int i,j;

    qa[-1] = 0.0;

    ff_lsp2polyf(lsp    , pa, lp_half_order    );
    ff_lsp2polyf(lsp + 1, qa, lp_half_order - 1);

    for (i = 1, j = LP_FILTER_ORDER - 1; i < lp_half_order; i++, j--) {
        double paf =  pa[i]            * (1 + lsp[LP_FILTER_ORDER - 1]);
        double qaf = (qa[i] - qa[i-2]) * (1 - lsp[LP_FILTER_ORDER - 1]);
        Az[i-1]  = (paf + qaf) * 0.5;
        Az[j-1]  = (paf - qaf) * 0.5;
    }

    Az[lp_half_order - 1] = (1.0 + lsp[LP_FILTER_ORDER - 1]) *
        pa[lp_half_order] * 0.5;

    Az[LP_FILTER_ORDER - 1] = lsp[LP_FILTER_ORDER - 1];
}

static void sipr_decode_lp(float *lsfnew, const float *lsfold, float *Az,
                           int num_subfr)
{
    double lsfint[LP_FILTER_ORDER];
    int i,j;

    t = t0 * 0.5;
    for (i = 0; i < num_subfr; i++) {
        for (j = 0; j < LP_FILTER_ORDER; j++)
            lsfint[j] = lsfold[j] * (1 - t) + t * lsfnew[j];

        lsp2lpc_sipr(lsfint, Az);
        Az += LP_FILTER_ORDER;
        t += t0;
    }
}

12465:1921e24d5886

        parms->gc_index[i]          = get_bits(pgb, p->gc_index_bits);
    }
}

static void sipr_decode_lp(float *lsfnew, const float *lsfold, float *Az,
                           int num_subfr)
{
    double lsfint[LP_FILTER_ORDER];
    int i,j;

    t = t0 * 0.5;
    for (i = 0; i < num_subfr; i++) {
        for (j = 0; j < LP_FILTER_ORDER; j++)
            lsfint[j] = lsfold[j] * (1 - t) + t * lsfnew[j];

        ff_amrwb_lsp2lpc(lsfint, Az, LP_FILTER_ORDER);
        Az += LP_FILTER_ORDER;
        t += t0;
    }
}