changeset 6158:2f43560f5dba libavcodec

simplify 4-point transform part a bit
author kostya
date Sat, 19 Jan 2008 11:04:02 +0000
parents 953c8efd5298
children eb2e3c3b7f78
files vc1dsp.c
diffstat 1 files changed, 25 insertions(+), 33 deletions(-) [+]
line wrap: on
line diff
--- a/vc1dsp.c	Sat Jan 19 06:52:29 2008 +0000
+++ b/vc1dsp.c	Sat Jan 19 11:04:02 2008 +0000
@@ -194,15 +194,13 @@
     for(i = 0; i < 8; i++){
         t1 = 17 * (src[ 0] + src[16]) + 64;
         t2 = 17 * (src[ 0] - src[16]) + 64;
-        t3 = 22 * src[ 8];
-        t4 = 22 * src[24];
-        t5 = 10 * src[ 8];
-        t6 = 10 * src[24];
+        t3 = 22 * src[ 8] + 10 * src[24];
+        t4 = 22 * src[24] - 10 * src[ 8];
 
-        dest[0*linesize] = cm[dest[0*linesize] + ((t1 + t3 + t6) >> 7)];
-        dest[1*linesize] = cm[dest[1*linesize] + ((t2 - t4 + t5) >> 7)];
-        dest[2*linesize] = cm[dest[2*linesize] + ((t2 + t4 - t5) >> 7)];
-        dest[3*linesize] = cm[dest[3*linesize] + ((t1 - t3 - t6) >> 7)];
+        dest[0*linesize] = cm[dest[0*linesize] + ((t1 + t3) >> 7)];
+        dest[1*linesize] = cm[dest[1*linesize] + ((t2 - t4) >> 7)];
+        dest[2*linesize] = cm[dest[2*linesize] + ((t2 + t4) >> 7)];
+        dest[3*linesize] = cm[dest[3*linesize] + ((t1 - t3) >> 7)];
 
         src ++;
         dest++;
@@ -223,15 +221,13 @@
     for(i = 0; i < 8; i++){
         t1 = 17 * (src[0] + src[2]) + 4;
         t2 = 17 * (src[0] - src[2]) + 4;
-        t3 = 22 * src[1];
-        t4 = 22 * src[3];
-        t5 = 10 * src[1];
-        t6 = 10 * src[3];
+        t3 = 22 * src[1] + 10 * src[3];
+        t4 = 22 * src[3] - 10 * src[1];
 
-        dst[0] = (t1 + t3 + t6) >> 3;
-        dst[1] = (t2 - t4 + t5) >> 3;
-        dst[2] = (t2 + t4 - t5) >> 3;
-        dst[3] = (t1 - t3 - t6) >> 3;
+        dst[0] = (t1 + t3) >> 3;
+        dst[1] = (t2 - t4) >> 3;
+        dst[2] = (t2 + t4) >> 3;
+        dst[3] = (t1 - t3) >> 3;
 
         src += 8;
         dst += 8;
@@ -273,7 +269,7 @@
 static void vc1_inv_trans_4x4_c(uint8_t *dest, int linesize, DCTELEM *block)
 {
     int i;
-    register int t1,t2,t3,t4,t5,t6;
+    register int t1,t2,t3,t4;
     DCTELEM *src, *dst;
     const uint8_t *cm = ff_cropTbl + MAX_NEG_CROP;
 
@@ -282,15 +278,13 @@
     for(i = 0; i < 4; i++){
         t1 = 17 * (src[0] + src[2]) + 4;
         t2 = 17 * (src[0] - src[2]) + 4;
-        t3 = 22 * src[1];
-        t4 = 22 * src[3];
-        t5 = 10 * src[1];
-        t6 = 10 * src[3];
+        t3 = 22 * src[1] + 10 * src[3];
+        t4 = 22 * src[3] - 10 * src[1];
 
-        dst[0] = (t1 + t3 + t6) >> 3;
-        dst[1] = (t2 - t4 + t5) >> 3;
-        dst[2] = (t2 + t4 - t5) >> 3;
-        dst[3] = (t1 - t3 - t6) >> 3;
+        dst[0] = (t1 + t3) >> 3;
+        dst[1] = (t2 - t4) >> 3;
+        dst[2] = (t2 + t4) >> 3;
+        dst[3] = (t1 - t3) >> 3;
 
         src += 8;
         dst += 8;
@@ -300,15 +294,13 @@
     for(i = 0; i < 4; i++){
         t1 = 17 * (src[ 0] + src[16]) + 64;
         t2 = 17 * (src[ 0] - src[16]) + 64;
-        t3 = 22 * src[ 8];
-        t4 = 22 * src[24];
-        t5 = 10 * src[ 8];
-        t6 = 10 * src[24];
+        t3 = 22 * src[ 8] + 10 * src[24];
+        t4 = 22 * src[24] - 10 * src[ 8];
 
-        dest[0*linesize] = cm[dest[0*linesize] + ((t1 + t3 + t6) >> 7)];
-        dest[1*linesize] = cm[dest[1*linesize] + ((t2 - t4 + t5) >> 7)];
-        dest[2*linesize] = cm[dest[2*linesize] + ((t2 + t4 - t5) >> 7)];
-        dest[3*linesize] = cm[dest[3*linesize] + ((t1 - t3 - t6) >> 7)];
+        dest[0*linesize] = cm[dest[0*linesize] + ((t1 + t3) >> 7)];
+        dest[1*linesize] = cm[dest[1*linesize] + ((t2 - t4) >> 7)];
+        dest[2*linesize] = cm[dest[2*linesize] + ((t2 + t4) >> 7)];
+        dest[3*linesize] = cm[dest[3*linesize] + ((t1 - t3) >> 7)];
 
         src ++;
         dest++;