# HG changeset patch
# User michael
# Date 1112650518 0
# Node ID b70b4b69960be4e38a75d416e553189e54bcfc5a
# Parent  65874d96c2dae06ae5bd5427cca1f0dcf915d4d6
merge vertical lifting steps and a little cleanup

diff -r 65874d96c2da -r b70b4b69960b snow.c
--- a/snow.c	Mon Apr 04 11:25:44 2005 +0000
+++ b/snow.c	Mon Apr 04 21:35:18 2005 +0000
@@ -1403,6 +1403,25 @@
     }
 }
 
+static void vertical_compose97i(DWTELEM *b0, DWTELEM *b1, DWTELEM *b2, DWTELEM *b3, DWTELEM *b4, DWTELEM *b5, int width){
+    int i;
+    
+    for(i=0; i<width; i++){
+        int r;
+        b4[i] -= (W_DM*(b3[i] + b5[i])+W_DO)>>W_DS;
+#ifdef lift5
+        b3[i] -= (W_CM*(b2[i] + b4[i])+W_CO)>>W_CS;
+#else
+        r= 3*(b2[i] + b4[i]);
+        r+= r>>4;
+        r+= r>>8;
+        b3[i] -= (r+W_CO)>>W_CS;
+#endif
+        b2[i] += (W_BM*(b1[i] + b3[i])+W_BO)>>W_BS;
+        b1[i] += (W_AM*(b0[i] + b2[i])+W_AO)>>W_AS;
+    }
+}
+
 static void spatial_compose97i_buffered_init(dwt_compose_t *cs, slice_buffer * sb, int height, int stride_line){
     cs->b0 = slice_buffer_get_line(sb, mirror(-3-1, height-1) * stride_line);
     cs->b1 = slice_buffer_get_line(sb, mirror(-3  , height-1) * stride_line);
@@ -1435,31 +1454,17 @@
     DWTELEM *b4= slice_buffer_get_line(sb, mirror4 * stride_line);
     DWTELEM *b5= slice_buffer_get_line(sb, mirror5 * stride_line);
         
-        if(stride_line == 1 && y+4 < height && 0){ 
-            int x;
-            for(x=0; x<width/2; x++)
-                b5[x] += 64*2;
-            for(; x<width; x++)
-                b5[x] += 169*2;
-        }
-        
-//        if(mirror3 <= mirror5 && mirror2 <= mirror4 && mirror1 <= mirror3 && mirror0 <= mirror2)
-//        {
-//{START_TIMER
-//            vertical_compose97_complete(b0, b1, b2, b3, b4, b5, width);
-//if(width>400){
-//STOP_TIMER("vertical_compose97i-NEW")}}
-//        }
-//        else
-//        {
 {START_TIMER
-            if(mirror3 <= mirror5) vertical_compose97iL1(b3, b4, b5, width);
-            if(mirror2 <= mirror4) vertical_compose97iH1(b2, b3, b4, width);
-            if(mirror1 <= mirror3) vertical_compose97iL0(b1, b2, b3, width);
-            if(mirror0 <= mirror2) vertical_compose97iH0(b0, b1, b2, width);
+    if(y>0 && y+4<height){
+        vertical_compose97i(b0, b1, b2, b3, b4, b5, width);
+    }else{
+        if(mirror3 <= mirror5) vertical_compose97iL1(b3, b4, b5, width);
+        if(mirror2 <= mirror4) vertical_compose97iH1(b2, b3, b4, width);
+        if(mirror1 <= mirror3) vertical_compose97iL0(b1, b2, b3, width);
+        if(mirror0 <= mirror2) vertical_compose97iH0(b0, b1, b2, width);
+    }
 if(width>400){
 STOP_TIMER("vertical_compose97i")}}
-//        }
 
 {START_TIMER
         if(y-1>=  0) horizontal_compose97i(b0, width);