up to date regression tests
[ffmpeg.git] / tests / rotozoom.c
index d275d0e..bbf96e1 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)
@@ -139,6 +163,21 @@ int teta = 0;
 int h_cos [360];
 int h_sin [360];
 
+static int ipol(uint8_t *src, int x, int y){
+    int int_x= x>>16;
+    int int_y= y>>16;
+    int frac_x= x&0xFFFF;
+    int frac_y= y&0xFFFF;
+    int s00= src[ ( int_x   &255) + 256*( int_y   &255) ];
+    int s01= src[ ((int_x+1)&255) + 256*( int_y   &255) ];
+    int s10= src[ ( int_x   &255) + 256*((int_y+1)&255) ];
+    int s11= src[ ((int_x+1)&255) + 256*((int_y+1)&255) ];
+    int s0= (((1<<16) - frac_x)*s00 + frac_x*s01)>>8;
+    int s1= (((1<<16) - frac_x)*s10 + frac_x*s11)>>8;
+    
+    return (((1<<16) - frac_y)*s0 + frac_y*s1)>>24;
+}
+
 void gen_image(int num, int w, int h)
 {
   const int c = h_cos [teta];
@@ -149,8 +188,6 @@ void gen_image(int num, int w, int h)
   
   const int xj = -(h/2) * s;
   const int yj = -(h/2) * c;
-  
-  unsigned dep;
   int i,j;
   
   int x,y;
@@ -160,69 +197,63 @@ 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));
-      put_pixel(i, j, tab_r[dep], tab_g[dep], tab_b[dep]);
+#if 1
+      put_pixel(i, j, ipol(tab_r, x, y), ipol(tab_g, x, y), ipol(tab_b, x, y));
+#else
+      {
+          unsigned dep;
+          dep = ((x>>16)&255) + (((y>>16)&255)<<8);
+          put_pixel(i, j, tab_r[dep], tab_g[dep], tab_b[dep]);
+      }
+#endif
     }
   }
   teta = (teta+1) % 360;
 }
 
-void init_demo() {
+#define W 256
+#define H 256
+
+void init_demo(const char *filename) {
   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");
-  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];
-    }
-       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);
+  fichier = fopen(filename,"r");
+  if (!fichier) {
+      perror(filename);
+      exit(1);
   }
-  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];
+  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 + 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;
   }
 }
 
@@ -231,8 +262,8 @@ int main(int argc, char **argv)
     int w, h, i;
     char buf[1024];
 
-    if (argc != 2) {
-        printf("usage: %s directory/\n"
+    if (argc != 3) {
+        printf("usage: %s directory/ image.pnm\n"
                "generate a test video stream\n", argv[0]);
         exit(1);
     }
@@ -245,10 +276,10 @@ int main(int argc, char **argv)
     width = w;
     height = h;
 
-    init_demo();
+    init_demo(argv[2]);
 
     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);
     }