Change license headers to say 'FFmpeg' instead of 'this program/this library'
[ffmpeg.git] / libavcodec / motion_test.c
1 /*
2  * (c) 2001 Fabrice Bellard
3  *
4  * This file is part of FFmpeg.
5  *
6  * FFmpeg is free software; you can redistribute it and/or
7  * modify it under the terms of the GNU Lesser General Public
8  * License as published by the Free Software Foundation; either
9  * version 2.1 of the License, or (at your option) any later version.
10  *
11  * FFmpeg is distributed in the hope that it will be useful,
12  * but WITHOUT ANY WARRANTY; without even the implied warranty of
13  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
14  * Lesser General Public License for more details.
15  *
16  * You should have received a copy of the GNU Lesser General Public
17  * License along with FFmpeg; if not, write to the Free Software
18  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
19  */
20
21 /**
22  * @file motion_test.c
23  * motion test.
24  */
25
26 #include <stdlib.h>
27 #include <stdio.h>
28 #include <string.h>
29 #include <sys/time.h>
30 #include <unistd.h>
31
32 #include "dsputil.h"
33
34 #include "i386/mmx.h"
35
36 int pix_abs16x16_mmx(uint8_t *blk1, uint8_t *blk2, int lx);
37 int pix_abs16x16_mmx1(uint8_t *blk1, uint8_t *blk2, int lx);
38 int pix_abs16x16_x2_mmx(uint8_t *blk1, uint8_t *blk2, int lx);
39 int pix_abs16x16_x2_mmx1(uint8_t *blk1, uint8_t *blk2, int lx);
40 int pix_abs16x16_x2_c(uint8_t *blk1, uint8_t *blk2, int lx);
41 int pix_abs16x16_y2_mmx(uint8_t *blk1, uint8_t *blk2, int lx);
42 int pix_abs16x16_y2_mmx1(uint8_t *blk1, uint8_t *blk2, int lx);
43 int pix_abs16x16_y2_c(uint8_t *blk1, uint8_t *blk2, int lx);
44 int pix_abs16x16_xy2_mmx(uint8_t *blk1, uint8_t *blk2, int lx);
45 int pix_abs16x16_xy2_mmx1(uint8_t *blk1, uint8_t *blk2, int lx);
46 int pix_abs16x16_xy2_c(uint8_t *blk1, uint8_t *blk2, int lx);
47
48 typedef int motion_func(uint8_t *blk1, uint8_t *blk2, int lx);
49
50 #define WIDTH 64
51 #define HEIGHT 64
52
53 uint8_t img1[WIDTH * HEIGHT];
54 uint8_t img2[WIDTH * HEIGHT];
55
56 void fill_random(uint8_t *tab, int size)
57 {
58     int i;
59     for(i=0;i<size;i++) {
60 #if 1
61         tab[i] = random() % 256;
62 #else
63         tab[i] = i;
64 #endif
65     }
66 }
67
68 void help(void)
69 {
70     printf("motion-test [-h]\n"
71            "test motion implementations\n");
72     exit(1);
73 }
74
75 int64_t gettime(void)
76 {
77     struct timeval tv;
78     gettimeofday(&tv,NULL);
79     return (int64_t)tv.tv_sec * 1000000 + tv.tv_usec;
80 }
81
82 #define NB_ITS 500
83
84 int dummy;
85
86 void test_motion(const char *name,
87                  motion_func *test_func, motion_func *ref_func)
88 {
89     int x, y, d1, d2, it;
90     uint8_t *ptr;
91     int64_t ti;
92     printf("testing '%s'\n", name);
93
94     /* test correctness */
95     for(it=0;it<20;it++) {
96
97         fill_random(img1, WIDTH * HEIGHT);
98         fill_random(img2, WIDTH * HEIGHT);
99
100         for(y=0;y<HEIGHT-17;y++) {
101             for(x=0;x<WIDTH-17;x++) {
102                 ptr = img2 + y * WIDTH + x;
103                 d1 = test_func(img1, ptr, WIDTH);
104                 d2 = ref_func(img1, ptr, WIDTH);
105                 if (d1 != d2) {
106                     printf("error: mmx=%d c=%d\n", d1, d2);
107                 }
108             }
109         }
110     }
111     emms();
112
113     /* speed test */
114     ti = gettime();
115     d1 = 0;
116     for(it=0;it<NB_ITS;it++) {
117         for(y=0;y<HEIGHT-17;y++) {
118             for(x=0;x<WIDTH-17;x++) {
119                 ptr = img2 + y * WIDTH + x;
120                 d1 += test_func(img1, ptr, WIDTH);
121             }
122         }
123     }
124     emms();
125     dummy = d1; /* avoid optimisation */
126     ti = gettime() - ti;
127
128     printf("  %0.0f kop/s\n",
129            (double)NB_ITS * (WIDTH - 16) * (HEIGHT - 16) /
130            (double)(ti / 1000.0));
131 }
132
133
134 int main(int argc, char **argv)
135 {
136     int c;
137
138     for(;;) {
139         c = getopt(argc, argv, "h");
140         if (c == -1)
141             break;
142         switch(c) {
143         case 'h':
144             help();
145             break;
146         }
147     }
148
149     printf("ffmpeg motion test\n");
150
151     test_motion("mmx", pix_abs16x16_mmx, pix_abs16x16_c);
152     test_motion("mmx_x2", pix_abs16x16_x2_mmx, pix_abs16x16_x2_c);
153     test_motion("mmx_y2", pix_abs16x16_y2_mmx, pix_abs16x16_y2_c);
154     test_motion("mmx_xy2", pix_abs16x16_xy2_mmx, pix_abs16x16_xy2_c);
155     return 0;
156 }