Mercurial > libavcodec.hg

changeset 5733:c2f88af57c16 libavcodec

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
was computing one more autocorrelation coefficient that was actually used
author lorenm
date Sat, 29 Sep 2007 07:43:52 +0000
parents d6fc148d1a48
children 0590ad246eca
files flacenc.c
diffstat 1 files changed, 7 insertions(+), 4 deletions(-) [+]
line wrap: on
line diff
--- a/flacenc.c	Sat Sep 29 05:41:27 2007 +0000
+++ b/flacenc.c	Sat Sep 29 07:43:52 2007 +0000
@@ -608,13 +608,14 @@
                              double *autoc)
 {
     int i, j;
-    double tmp[len + lag];
+    double tmp[len + lag + 1];
     double *data1= tmp + lag;
 
     apply_welch_window(data, len, data1);
 
     for(j=0; j<lag; j++)
         data1[j-lag]= 0.0;
+    data1[len] = 0.0;
 
     for(j=0; j<lag; j+=2){
         double sum0 = 1.0, sum1 = 1.0;
@@ -628,8 +629,10 @@
 
     if(j==lag){
         double sum = 1.0;
-        for(i=0; i<len; i++)
-            sum += data1[i] * data1[i-j];
+        for(i=0; i<len; i+=2){
+            sum += data1[i  ] * data1[i-j  ]
+                 + data1[i+1] * data1[i-j+1];
+        }
         autoc[j] = sum;
     }
 }
@@ -757,7 +760,7 @@
     assert(max_order >= MIN_LPC_ORDER && max_order <= MAX_LPC_ORDER);
 
     if(use_lpc == 1){
-        compute_autocorr(samples, blocksize, max_order+1, autoc);
+        compute_autocorr(samples, blocksize, max_order, autoc);
 
         compute_lpc_coefs(autoc, max_order, lpc, ref);
     }else{