Mercurial > libavcodec.hg

changeset 3908:045c650f0951 libavcodec

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Original Commit: r109 | ods15 | 2006-10-01 19:14:52 +0200 (Sun, 01 Oct 2006) | 2 lines 5% speed improovement by changing put_vector
author ods15
date Mon, 02 Oct 2006 06:09:41 +0000
parents 5f5525341673
children 3cac42cde8ea
files vorbis_enc.c
diffstat 1 files changed, 12 insertions(+), 10 deletions(-) [+]
line wrap: on
line diff
--- a/vorbis_enc.c	Mon Oct 02 06:09:39 2006 +0000
+++ b/vorbis_enc.c	Mon Oct 02 06:09:41 2006 +0000
@@ -22,6 +22,7 @@
  * @author Oded Shimon <ods15@ods15.dyndns.org>
  */
 
+#include <float.h>
 #include "avcodec.h"
 #include "dsputil.h"
 #include "vorbis.h"
@@ -40,6 +41,7 @@
     int lookup;
     int * quantlist;
     float * dimentions;
+    float * pow2;
 } codebook_t;
 
 typedef struct {
@@ -645,10 +647,11 @@
 
     ff_vorbis_len2vlc(cb->lens, cb->codewords, cb->nentries);
 
-    if (!cb->lookup) cb->dimentions = NULL;
+    if (!cb->lookup) cb->pow2 = cb->dimentions = NULL;
     else {
         int vals = cb_lookup_vals(cb->lookup, cb->ndimentions, cb->nentries);
         cb->dimentions = av_malloc(sizeof(float) * cb->nentries * cb->ndimentions);
+        cb->pow2 = av_mallocz(sizeof(float) * cb->nentries);
         for (i = 0; i < cb->nentries; i++) {
             float last = 0;
             int j;
@@ -660,8 +663,10 @@
 
                 cb->dimentions[i * cb->ndimentions + j] = last + cb->min + cb->quantlist[off] * cb->delta;
                 if (cb->seq_p) last = cb->dimentions[i * cb->ndimentions + j];
+                cb->pow2[i] += cb->dimentions[i * cb->ndimentions + j]*cb->dimentions[i * cb->ndimentions + j];
                 div *= vals;
             }
+            cb->pow2[i] /= 2.;
         }
     }
 
@@ -1218,19 +1223,15 @@
 }
 
 static float * put_vector(codebook_t * book, PutBitContext * pb, float * num) {
-    int i;
-    int entry = -1;
-    float distance = 0;
+    int i, entry = -1;
+    float distance = FLT_MAX;
     assert(book->dimentions);
     for (i = 0; i < book->nentries; i++) {
-        float d = 0.;
+        float * vec = book->dimentions + i * book->ndimentions, d = book->pow2[i];
         int j;
         if (!book->lens[i]) continue;
-        for (j = 0; j < book->ndimentions; j++) {
-            float a = (book->dimentions[i * book->ndimentions + j] - num[j]);
-            d += a*a;
-        }
-        if (entry == -1 || distance > d) {
+        for (j = 0; j < book->ndimentions; j++) d -= vec[j] * num[j];
+        if (distance > d) {
             entry = i;
             distance = d;
         }
@@ -1459,6 +1460,7 @@
         av_freep(&venc->codebooks[i].codewords);
         av_freep(&venc->codebooks[i].quantlist);
         av_freep(&venc->codebooks[i].dimentions);
+        av_freep(&venc->codebooks[i].pow2);
     }
     av_freep(&venc->codebooks);