comparison flacenc.c @ 3388:22c7b4c96c2d libavcodec

simplify compute_autocorr

3387:92fabe8c1422

 * A Welch window function is applied before calculation.
 */
static void compute_autocorr(const int32_t *data, int len, int lag,
                             double *autoc)
{
    int i;
    double *data1;
    int lag_ptr, ptr;

    data1 = av_malloc(len * sizeof(double));
    apply_welch_window(data, len, data1);

    for(i=0; i<lag; i++) autoc[i] = 1.0;

    ptr = 0;
    while(ptr <= lag) {
        lag_ptr = 0;
        while(lag_ptr <= ptr) {
            autoc[ptr-lag_ptr] += data1[ptr] * data1[lag_ptr];
            lag_ptr++;
        }
        ptr++;
    }
    while(ptr < len) {
        lag_ptr = ptr - lag;
        while(lag_ptr <= ptr) {
            autoc[ptr-lag_ptr] += data1[ptr] * data1[lag_ptr];
            lag_ptr++;
        }
        ptr++;
    }

    av_freep(&data1);
}

/**
 * Levinson-Durbin recursion.
 * Produces LPC coefficients from autocorrelation data.

3388:22c7b4c96c2d

 * A Welch window function is applied before calculation.
 */
static void compute_autocorr(const int32_t *data, int len, int lag,
                             double *autoc)
{
    int i, lag_ptr;
    double tmp[len + lag];
    double *data1= tmp + lag;

    apply_welch_window(data, len, data1);

    for(i=0; i<lag; i++){
        autoc[i] = 1.0;
        data1[i-lag]= 0.0;
    }

    for(i=0; i<len; i++){
        for(lag_ptr= i-lag; lag_ptr<=i; lag_ptr++){
            autoc[i-lag_ptr] += data1[i] * data1[lag_ptr];
        }
    }
}

/**
 * Levinson-Durbin recursion.
 * Produces LPC coefficients from autocorrelation data.