finished integrating rotozoom test in ffmpeg
[libav.git] / tests / rotozoom.c
index d275d0e..71d369f 100644 (file)
@@ -5,7 +5,31 @@
  */
 #include <stdlib.h>
 #include <stdio.h>
-#include <math.h>
+#include <inttypes.h>
+
+#define FIXP (1<<16)
+#define MY_PI 205887 //(M_PI*FIX)
+
+static int64_t int_pow(int64_t a, int p){
+    int64_t v= FIXP;
+
+    for(; p; p--){
+        v*= a;
+        v/= FIXP;
+    }
+
+    return v;
+}
+
+static int64_t int_sin(int64_t a){
+    if(a<0) a= MY_PI-a; // 0..inf
+    a %= 2*MY_PI;       // 0..2PI
+
+    if(a>=MY_PI*3/2) a -= 2*MY_PI;  // -PI/2 .. 3PI/2
+    if(a>=MY_PI/2  ) a = MY_PI - a; // -PI/2 ..  PI/2
+   
+    return a - int_pow(a, 3)/6 + int_pow(a, 5)/120 - int_pow(a, 7)/5040;
+}
 
 #define SCALEBITS 8
 #define ONE_HALF  (1 << (SCALEBITS - 1))
@@ -78,7 +102,7 @@ static void rgb24_to_yuv420p(UINT8 *lum, UINT8 *cb, UINT8 *cr,
 /* cif format */
 #define DEFAULT_WIDTH   352
 #define DEFAULT_HEIGHT  288
-#define DEFAULT_NB_PICT 360
+#define DEFAULT_NB_PICT 50
 
 void pgmyuv_save(const char *filename, int w, int h,
                  unsigned char *rgb_tab)
@@ -160,69 +184,51 @@ void gen_image(int num, int w, int h)
 
   for (j=0;j<h;j++) {
 
-    x = xprime + xi;
+    x = xprime + xi + FIXP*w/2;
     xprime += s;
 
-    y = yprime + yi;
+    y = yprime + yi + FIXP*h/2;
     yprime += c;
       
     for ( i=0 ; i<w ; i++ ) {
       x += c;
       y -= s;
-      dep = ((x>>8)&255) + (y&(255<<8));
+      dep = ((x>>16)&255) + (((y>>16)&255)<<8);
       put_pixel(i, j, tab_r[dep], tab_g[dep], tab_b[dep]);
     }
   }
   teta = (teta+1) % 360;
 }
 
+#define W 256
+#define H 256
+
 void init_demo() {
   int i,j;
-  double h;
-  double radian;
-  char line[3 * 128];
+  int h;
+  int radian;
+  char line[3 * W];
 
   FILE *fichier;
 
-
-  fichier = fopen("rotozoom-ffmpeg.pnm","r");
+  fichier = fopen("lena.pnm","r");
   fread(line, 1, 15, fichier);
-  for (i=0;i<128;i++) {
-    fread(line,1,3*128,fichier);
-    for (j=0;j<128;j++) {
-         tab_r[256*i+j] = line[3*j    ];
-         tab_g[256*i+j] = line[3*j + 1];
-         tab_b[256*i+j] = line[3*j + 2];
+  for (i=0;i<H;i++) {
+    fread(line,1,3*W,fichier);
+    for (j=0;j<W;j++) {
+         tab_r[W*i+j] = line[3*j    ];
+         tab_g[W*i+j] = line[3*j + 1];
+         tab_b[W*i+j] = line[3*j + 2];
     }
-       memcpy(tab_r + 257*128 + i*256, tab_r + i*256, 128);
-       memcpy(tab_g + 257*128 + i*256, tab_g + i*256, 128);
-       memcpy(tab_b + 257*128 + i*256, tab_b + i*256, 128);
   }
   fclose(fichier);
 
-  fichier = fopen("rotozoom-tux.pnm","r");
-  fread(line, 1, 15, fichier);
-  for (i=0;i<128;i++) {
-    fread(line,1,3*128,fichier);
-    for (j=0;j<128;j++) {
-         tab_r[128 + 256*i+j] = line[3*j    ];
-         tab_g[128 + 256*i+j] = line[3*j + 1];
-         tab_b[128 + 256*i+j] = line[3*j + 2];
-    }
-       memcpy(tab_r + 256*128 + i*256, tab_r + 128 + i*256, 128);
-       memcpy(tab_g + 256*128 + i*256, tab_g + 128 + i*256, 128);
-       memcpy(tab_b + 256*128 + i*256, tab_b + 128 + i*256, 128);
-  }
-  fclose(fichier);
-
-
-
   /* tables sin/cos */
   for (i=0;i<360;i++) {
-    radian = 2*i*M_PI/360;
-    h = 2 + cos (radian);      
-    h_cos[i] = 256 * ( h * cos (radian) );
-    h_sin[i] = 256 * ( h * sin (radian) );
+    radian = 2*i*MY_PI/360;
+    h = 2*FIXP + int_sin (radian);
+    h_cos[i] = ( h * int_sin (radian + MY_PI/2) )/2/FIXP;
+    h_sin[i] = ( h * int_sin (radian          ) )/2/FIXP;
   }
 }
 
@@ -248,7 +254,7 @@ int main(int argc, char **argv)
     init_demo();
 
     for(i=0;i<DEFAULT_NB_PICT;i++) {
-        snprintf(buf, sizeof(buf), "%s%03d.pgm", argv[1], i);
+        snprintf(buf, sizeof(buf), "%s%d.pgm", argv[1], i);
         gen_image(i, w, h);
         pgmyuv_save(buf, w, h, rgb_tab);
     }