Annotation of /trunk/xvidcore/examples/xvid_bench.c

Revision 1804 - (view) (download)

1 :	edgomez	1382	/*****************************************************************************
2 :	Isibaar	225	*
3 :	edgomez	1382	* XVID MPEG-4 VIDEO CODEC
4 :			* - Unit tests and benches -
5 :	Isibaar	225	*
6 :	edgomez	1382	* Copyright(C) 2002 Pascal Massimino <skal@planet-d.net>
7 :	Isibaar	225	*
8 :	edgomez	1382	* This program is free software; you can redistribute it and/or modify
9 :			* it under the terms of the GNU General Public License as published by
10 :			* the Free Software Foundation; either version 2 of the License, or
11 :			* (at your option) any later version.
12 :	Isibaar	225	*
13 :	edgomez	1382	* This program is distributed in the hope that it will be useful,
14 :			* but WITHOUT ANY WARRANTY; without even the implied warranty of
15 :			* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
16 :			* GNU General Public License for more details.
17 :	Isibaar	225	*
18 :	edgomez	1382	* You should have received a copy of the GNU General Public License
19 :			* along with this program; if not, write to the Free Software
20 :			* Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
21 :			*
22 :	Isibaar	1804	* $Id: xvid_bench.c,v 1.38 2008-11-26 23:37:28 Isibaar Exp $
23 :	edgomez	1382	*
24 :			****************************************************************************/
25 :	Isibaar	225
26 :	edgomez	1382	/*****************************************************************************
27 :	Isibaar	225	*
28 :			* 'Reference' output is at the end of file.
29 :			*
30 :	edgomez	1382	* compiles with something like:
31 :			* gcc -o xvid_bench xvid_bench.c -I../src/ -lxvidcore -lm
32 :	Isibaar	225	*
33 :	edgomez	1382	****************************************************************************/
34 :	Isibaar	225
35 :			#include <stdio.h>
36 :			#include <stdlib.h>
37 :	edgomez	1382	#include <string.h> /* for memset */
38 :	Isibaar	225	#include <assert.h>
39 :
40 :	suxen_drol	860	#ifndef WIN32
41 :	edgomez	1382	#include <sys/time.h> /* for gettimeofday */
42 :	suxen_drol	860	#else
43 :			#include <time.h>
44 :			#endif
45 :
46 :
47 :	Isibaar	225	#include "xvid.h"
48 :
49 :	edgomez	851	// inner guts
50 :	Isibaar	225	#include "dct/idct.h"
51 :			#include "dct/fdct.h"
52 :			#include "image/colorspace.h"
53 :			#include "image/interpolate8x8.h"
54 :			#include "utils/mem_transfer.h"
55 :	edgomez	1382	#include "quant/quant.h"
56 :	Isibaar	225	#include "motion/sad.h"
57 :			#include "utils/emms.h"
58 :			#include "utils/timer.h"
59 :			#include "quant/quant_matrix.c"
60 :			#include "bitstream/cbp.h"
61 :	Skal	1746	#include "bitstream/bitstream.h"
62 :	Isibaar	225
63 :	Isibaar	262	#include <math.h>
64 :	suxen_drol	860
65 :			#ifndef M_PI
66 :			#define M_PI 3.14159265358979323846
67 :			#endif
68 :
69 :	Skal	1614	int speed_ref = 100; /* on slow machines, decrease this value */
70 :			int verbose = 0;
71 :			unsigned int cpu_mask;
72 :	Isibaar	225
73 :			/*********************************************************************
74 :			* misc
75 :			*********************************************************************/
76 :
77 :	edgomez	1382	/* returns time in micro-s*/
78 :	Isibaar	225	double gettime_usec()
79 :			{
80 :	suxen_drol	860	#ifndef WIN32
81 :	edgomez	1382	struct timeval tv;
82 :			gettimeofday(&tv, 0);
83 :			return tv.tv_sec*1.0e6 + tv.tv_usec;
84 :	suxen_drol	860	#else
85 :			clock_t clk;
86 :			clk = clock();
87 :	Skal	1613	return clk * 1000. / CLOCKS_PER_SEC; /* clock() returns time in Milliseconds */
88 :	suxen_drol	860	#endif
89 :	Isibaar	225	}
90 :
91 :	edgomez	1382	/* returns squared deviates (mean(vv)-mean(v)^2) of a 8x8 block /
92 :	Isibaar	225	double sqr_dev(uint8_t v[8*8])
93 :			{
94 :	edgomez	1382	double sum=0.;
95 :			double sum2=0.;
96 :			int n;
97 :			for (n=0;n<8*8;n++)
98 :			{
99 :			sum += v[n];
100 :			sum2 += v[n]*v[n];
101 :			}
102 :			sum2 /= n;
103 :			sum /= n;
104 :			return sum2-sum*sum;
105 :	Isibaar	225	}
106 :
107 :			/*********************************************************************
108 :			* cpu init
109 :			*********************************************************************/
110 :
111 :			typedef struct {
112 :	edgomez	1382	const char *name;
113 :			unsigned int cpu;
114 :	Isibaar	225	} CPU;
115 :
116 :	Skal	1621	CPU cpu_list[] = {
117 :			{ "PLAINC ", 0 },
118 :	Isibaar	1801	#if defined(ARCH_IS_IA32) \|\| defined(ARCH_IS_X86_64)
119 :	Skal	1621	{ "MMX ", XVID_CPU_MMX },
120 :			{ "MMXEXT ", XVID_CPU_MMXEXT \| XVID_CPU_MMX },
121 :			{ "SSE2 ", XVID_CPU_SSE2 \| XVID_CPU_MMX },
122 :	Isibaar	1794	{ "SSE3 ", XVID_CPU_SSE3 \| XVID_CPU_SSE2 \| XVID_CPU_MMX },
123 :			{ "SSE41 ", XVID_CPU_SSE41\| XVID_CPU_SSE3 \| XVID_CPU_SSE2 \| XVID_CPU_MMX },
124 :	Isibaar	1801	{ "3DNOW ", XVID_CPU_3DNOW },
125 :	Skal	1621	{ "3DNOWE ", XVID_CPU_3DNOW \| XVID_CPU_3DNOWEXT },
126 :	edgomez	1382	#endif
127 :	edgomez	1412	#ifdef ARCH_IS_PPC
128 :	Skal	1621	{ "ALTIVEC", XVID_CPU_ALTIVEC },
129 :	edgomez	1412	#endif
130 :	Skal	1621	#ifdef ARCH_IS_IA64
131 :			// { "IA64 ", XVID_CPU_IA64 },
132 :	edgomez	1382	#endif
133 :	Skal	1621	// { "TSC ", XVID_CPU_TSC },
134 :			{ 0, 0 }
135 :			};
136 :	Isibaar	225
137 :
138 :			int init_cpu(CPU *cpu)
139 :			{
140 :	edgomez	1382	xvid_gbl_info_t xinfo;
141 :	Isibaar	225
142 :	edgomez	1382	/* Get the available CPU flags */
143 :			memset(&xinfo, 0, sizeof(xinfo));
144 :			xinfo.version = XVID_VERSION;
145 :			xvid_global(NULL, XVID_GBL_INFO, &xinfo, NULL);
146 :
147 :			/* Are we trying to test a subset of the host CPU features */
148 :			if ((xinfo.cpu_flags & cpu->cpu) == cpu->cpu) {
149 :			int xerr;
150 :			xvid_gbl_init_t xinit;
151 :			memset(&xinit, 0, sizeof(xinit));
152 :			xinit.cpu_flags = cpu->cpu \| XVID_CPU_FORCE;
153 :			xinit.version = XVID_VERSION;
154 :			xerr = xvid_global(NULL, XVID_GBL_INIT, &xinit, NULL);
155 :			if (xerr==XVID_ERR_FAIL) {
156 :			/* libxvidcore failed to init */
157 :			return 0;
158 :			}
159 :			} else {
160 :			/* The host CPU doesn't support some required feature for this test */
161 :			return(0);
162 :			}
163 :			return 1;
164 :	Isibaar	225	}
165 :
166 :	edgomez	1382	#define CRC32_REMAINDER 0xCBF43926
167 :			#define CRC32_INITIAL 0xffffffff
168 :
169 :			#define DO1(c, crc) ((crc) = crc32tab[((unsigned int)((crc)>>24) ^ (*c++)) & 0xff] ^ ((crc) << 8))
170 :			#define DO2(c, crc) DO1(c, crc); DO1(c, crc);
171 :			#define DO4(c, crc) DO2(c, crc); DO2(c, crc);
172 :			#define DO8(c, crc) DO4(c, crc); DO4(c, crc);
173 :
174 :			/******************************************************************************
175 :			* Precomputed AAL5 CRC32 lookup table
176 :			******************************************************************************/
177 :
178 :			static unsigned long crc32tab[256] = {
179 :
180 :			0x00000000L, 0x04C11DB7L, 0x09823B6EL, 0x0D4326D9L,
181 :			0x130476DCL, 0x17C56B6BL, 0x1A864DB2L, 0x1E475005L,
182 :			0x2608EDB8L, 0x22C9F00FL, 0x2F8AD6D6L, 0x2B4BCB61L,
183 :			0x350C9B64L, 0x31CD86D3L, 0x3C8EA00AL, 0x384FBDBDL,
184 :			0x4C11DB70L, 0x48D0C6C7L, 0x4593E01EL, 0x4152FDA9L,
185 :			0x5F15ADACL, 0x5BD4B01BL, 0x569796C2L, 0x52568B75L,
186 :			0x6A1936C8L, 0x6ED82B7FL, 0x639B0DA6L, 0x675A1011L,
187 :			0x791D4014L, 0x7DDC5DA3L, 0x709F7B7AL, 0x745E66CDL,
188 :			0x9823B6E0L, 0x9CE2AB57L, 0x91A18D8EL, 0x95609039L,
189 :			0x8B27C03CL, 0x8FE6DD8BL, 0x82A5FB52L, 0x8664E6E5L,
190 :			0xBE2B5B58L, 0xBAEA46EFL, 0xB7A96036L, 0xB3687D81L,
191 :			0xAD2F2D84L, 0xA9EE3033L, 0xA4AD16EAL, 0xA06C0B5DL,
192 :			0xD4326D90L, 0xD0F37027L, 0xDDB056FEL, 0xD9714B49L,
193 :			0xC7361B4CL, 0xC3F706FBL, 0xCEB42022L, 0xCA753D95L,
194 :			0xF23A8028L, 0xF6FB9D9FL, 0xFBB8BB46L, 0xFF79A6F1L,
195 :			0xE13EF6F4L, 0xE5FFEB43L, 0xE8BCCD9AL, 0xEC7DD02DL,
196 :			0x34867077L, 0x30476DC0L, 0x3D044B19L, 0x39C556AEL,
197 :			0x278206ABL, 0x23431B1CL, 0x2E003DC5L, 0x2AC12072L,
198 :			0x128E9DCFL, 0x164F8078L, 0x1B0CA6A1L, 0x1FCDBB16L,
199 :			0x018AEB13L, 0x054BF6A4L, 0x0808D07DL, 0x0CC9CDCAL,
200 :			0x7897AB07L, 0x7C56B6B0L, 0x71159069L, 0x75D48DDEL,
201 :			0x6B93DDDBL, 0x6F52C06CL, 0x6211E6B5L, 0x66D0FB02L,
202 :			0x5E9F46BFL, 0x5A5E5B08L, 0x571D7DD1L, 0x53DC6066L,
203 :			0x4D9B3063L, 0x495A2DD4L, 0x44190B0DL, 0x40D816BAL,
204 :			0xACA5C697L, 0xA864DB20L, 0xA527FDF9L, 0xA1E6E04EL,
205 :			0xBFA1B04BL, 0xBB60ADFCL, 0xB6238B25L, 0xB2E29692L,
206 :			0x8AAD2B2FL, 0x8E6C3698L, 0x832F1041L, 0x87EE0DF6L,
207 :			0x99A95DF3L, 0x9D684044L, 0x902B669DL, 0x94EA7B2AL,
208 :			0xE0B41DE7L, 0xE4750050L, 0xE9362689L, 0xEDF73B3EL,
209 :			0xF3B06B3BL, 0xF771768CL, 0xFA325055L, 0xFEF34DE2L,
210 :			0xC6BCF05FL, 0xC27DEDE8L, 0xCF3ECB31L, 0xCBFFD686L,
211 :			0xD5B88683L, 0xD1799B34L, 0xDC3ABDEDL, 0xD8FBA05AL,
212 :			0x690CE0EEL, 0x6DCDFD59L, 0x608EDB80L, 0x644FC637L,
213 :			0x7A089632L, 0x7EC98B85L, 0x738AAD5CL, 0x774BB0EBL,
214 :			0x4F040D56L, 0x4BC510E1L, 0x46863638L, 0x42472B8FL,
215 :			0x5C007B8AL, 0x58C1663DL, 0x558240E4L, 0x51435D53L,
216 :			0x251D3B9EL, 0x21DC2629L, 0x2C9F00F0L, 0x285E1D47L,
217 :			0x36194D42L, 0x32D850F5L, 0x3F9B762CL, 0x3B5A6B9BL,
218 :			0x0315D626L, 0x07D4CB91L, 0x0A97ED48L, 0x0E56F0FFL,
219 :			0x1011A0FAL, 0x14D0BD4DL, 0x19939B94L, 0x1D528623L,
220 :			0xF12F560EL, 0xF5EE4BB9L, 0xF8AD6D60L, 0xFC6C70D7L,
221 :			0xE22B20D2L, 0xE6EA3D65L, 0xEBA91BBCL, 0xEF68060BL,
222 :			0xD727BBB6L, 0xD3E6A601L, 0xDEA580D8L, 0xDA649D6FL,
223 :			0xC423CD6AL, 0xC0E2D0DDL, 0xCDA1F604L, 0xC960EBB3L,
224 :			0xBD3E8D7EL, 0xB9FF90C9L, 0xB4BCB610L, 0xB07DABA7L,
225 :			0xAE3AFBA2L, 0xAAFBE615L, 0xA7B8C0CCL, 0xA379DD7BL,
226 :			0x9B3660C6L, 0x9FF77D71L, 0x92B45BA8L, 0x9675461FL,
227 :			0x8832161AL, 0x8CF30BADL, 0x81B02D74L, 0x857130C3L,
228 :			0x5D8A9099L, 0x594B8D2EL, 0x5408ABF7L, 0x50C9B640L,
229 :			0x4E8EE645L, 0x4A4FFBF2L, 0x470CDD2BL, 0x43CDC09CL,
230 :			0x7B827D21L, 0x7F436096L, 0x7200464FL, 0x76C15BF8L,
231 :			0x68860BFDL, 0x6C47164AL, 0x61043093L, 0x65C52D24L,
232 :			0x119B4BE9L, 0x155A565EL, 0x18197087L, 0x1CD86D30L,
233 :			0x029F3D35L, 0x065E2082L, 0x0B1D065BL, 0x0FDC1BECL,
234 :			0x3793A651L, 0x3352BBE6L, 0x3E119D3FL, 0x3AD08088L,
235 :			0x2497D08DL, 0x2056CD3AL, 0x2D15EBE3L, 0x29D4F654L,
236 :			0xC5A92679L, 0xC1683BCEL, 0xCC2B1D17L, 0xC8EA00A0L,
237 :			0xD6AD50A5L, 0xD26C4D12L, 0xDF2F6BCBL, 0xDBEE767CL,
238 :			0xE3A1CBC1L, 0xE760D676L, 0xEA23F0AFL, 0xEEE2ED18L,
239 :			0xF0A5BD1DL, 0xF464A0AAL, 0xF9278673L, 0xFDE69BC4L,
240 :			0x89B8FD09L, 0x8D79E0BEL, 0x803AC667L, 0x84FBDBD0L,
241 :			0x9ABC8BD5L, 0x9E7D9662L, 0x933EB0BBL, 0x97FFAD0CL,
242 :			0xAFB010B1L, 0xAB710D06L, 0xA6322BDFL, 0xA2F33668L,
243 :			0xBCB4666DL, 0xB8757BDAL, 0xB5365D03L, 0xB1F740B4L
244 :
245 :			};
246 :
247 :			uint32_t
248 :	Skal	1614	calc_crc(uint8_t *mem, int len, uint32_t crc)
249 :	edgomez	1382	{
250 :			while( len >= 8) {
251 :			DO8(mem, crc);
252 :			len -= 8;
253 :			}
254 :
255 :			while( len ) {
256 :			DO1(mem, crc);
257 :			len--;
258 :			}
259 :
260 :	Skal	1614	return crc;
261 :	edgomez	1382	}
262 :
263 :	Skal	1621	void byte_swap(uint8_t *mem, int len, int element_size) {
264 :			#ifdef ARCH_IS_BIG_ENDIAN
265 :			int i;
266 :
267 :			if(element_size == 1) {
268 :			/* No need to swap */
269 :			} else if(element_size == 2) {
270 :			uint8_t temp[2];
271 :
272 :			for(i=0; i < (len/2); i++ ) {
273 :			temp[0] = mem[0];
274 :			temp[1] = mem[1];
275 :			mem[0] = temp[1];
276 :			mem[1] = temp[0];
277 :
278 :			mem += 2;
279 :			}
280 :			} else if(element_size == 4) {
281 :			uint8_t temp[4];
282 :
283 :			for(i=0; i < (len/4); i++ ) {
284 :			temp[0] = mem[0];
285 :			temp[1] = mem[1];
286 :			temp[2] = mem[2];
287 :			temp[3] = mem[3];
288 :			mem[0] = temp[3];
289 :			mem[1] = temp[2];
290 :			mem[2] = temp[1];
291 :			mem[3] = temp[0];
292 :
293 :			mem += 4;
294 :			}
295 :			} else {
296 :			printf("ERROR: byte_swap unsupported element_size(%u)\n", element_size);
297 :			}
298 :			#endif
299 :			}
300 :
301 :	Isibaar	225	/*********************************************************************
302 :			* test DCT
303 :			*********************************************************************/
304 :
305 :			#define ABS(X) ((X)<0 ? -(X) : (X))
306 :
307 :			void test_dct()
308 :			{
309 :	edgomez	1382	const int nb_tests = 300*speed_ref;
310 :			int tst;
311 :			CPU *cpu;
312 :			int i;
313 :			DECLARE_ALIGNED_MATRIX(iDst0, 8, 8, short, 16);
314 :			DECLARE_ALIGNED_MATRIX(iDst, 8, 8, short, 16);
315 :			DECLARE_ALIGNED_MATRIX(fDst, 8, 8, short, 16);
316 :			double overhead;
317 :	Isibaar	225
318 :	edgomez	1382	printf( "\n ===== test fdct/idct =====\n" );
319 :	Isibaar	225
320 :	edgomez	1382	for(i=0; i<88; ++i) iDst0[i] = (i7-i*i) & 0x7f;
321 :			overhead = gettime_usec();
322 :			for(tst=0; tst<nb_tests; ++tst)
323 :			{
324 :			for(i=0; i<8*8; ++i) fDst[i] = iDst0[i];
325 :			for(i=0; i<8*8; ++i) iDst[i] = fDst[i];
326 :			}
327 :			overhead = gettime_usec() - overhead;
328 :	Isibaar	225
329 :	edgomez	1382	for(cpu = cpu_list; cpu->name!=0; ++cpu)
330 :			{
331 :			double t, PSNR, MSE;
332 :	Isibaar	225
333 :	edgomez	1382	if (!init_cpu(cpu))
334 :			continue;
335 :	Isibaar	225
336 :	edgomez	1382	t = gettime_usec();
337 :			emms();
338 :			for(tst=0; tst<nb_tests; ++tst)
339 :			{
340 :			for(i=0; i<8*8; ++i) fDst[i] = iDst0[i];
341 :			fdct(fDst);
342 :			for(i=0; i<8*8; ++i) iDst[i] = fDst[i];
343 :			idct(iDst);
344 :			}
345 :			emms();
346 :			t = (gettime_usec() - t - overhead) / nb_tests;
347 :			MSE = 0.;
348 :			for(i=0; i<8*8; ++i) {
349 :			double delta = 1.0*(iDst[i] - iDst0[i]);
350 :			MSE += delta*delta;
351 :			}
352 :			PSNR = (MSE==0.) ? 1.e6 : -4.3429448*log( MSE/64. );
353 :			printf( "%s - %.3f usec PSNR=%.3f MSE=%.3f %s\n",
354 :			cpu->name, t, PSNR, MSE,
355 :			(ABS(MSE)>=64)? "\| ERROR" :"");
356 :			}
357 :	Isibaar	225	}
358 :
359 :			/*********************************************************************
360 :			* test SAD
361 :			*********************************************************************/
362 :
363 :			void test_sad()
364 :			{
365 :	edgomez	1382	const int nb_tests = 2000*speed_ref;
366 :			int tst;
367 :			CPU *cpu;
368 :			int i;
369 :			DECLARE_ALIGNED_MATRIX(Cur, 16, 16, uint8_t, 16);
370 :			DECLARE_ALIGNED_MATRIX(Ref1, 16, 16, uint8_t, 16);
371 :			DECLARE_ALIGNED_MATRIX(Ref2, 16, 16, uint8_t, 16);
372 :	Isibaar	225
373 :	edgomez	1382	printf( "\n ====== test SAD ======\n" );
374 :			for(i=0; i<16*16;++i) {
375 :			Cur[i] = (i/5) ^ 0x05;
376 :			Ref1[i] = (i + 0x0b) & 0xff;
377 :			Ref2[i] = i ^ 0x76;
378 :			}
379 :	Isibaar	225
380 :	edgomez	1382	for(cpu = cpu_list; cpu->name!=0; ++cpu)
381 :			{
382 :			double t;
383 :			uint32_t s;
384 :			if (!init_cpu(cpu))
385 :			continue;
386 :	Isibaar	225
387 :	edgomez	1382	t = gettime_usec();
388 :			emms();
389 :			for(tst=0; tst<nb_tests; ++tst) s = sad8(Cur, Ref1, 16);
390 :			emms();
391 :			t = (gettime_usec() - t) / nb_tests;
392 :			printf("%s - sad8 %.3f usec sad=%d %s\n",
393 :			cpu->name, t, s,
394 :			(s!=3776)?"\| ERROR": "" );
395 :	Isibaar	225
396 :	edgomez	1382	t = gettime_usec();
397 :			emms();
398 :			for(tst=0; tst<nb_tests; ++tst) s = sad16(Cur, Ref1, 16, -1);
399 :			emms();
400 :			t = (gettime_usec() - t) / nb_tests;
401 :			printf("%s - sad16 %.3f usec sad=%d %s\n",
402 :			cpu->name, t, s,
403 :			(s!=27214)?"\| ERROR": "" );
404 :	Isibaar	225
405 :	edgomez	1382	t = gettime_usec();
406 :			emms();
407 :			for(tst=0; tst<nb_tests; ++tst) s = sad16bi(Cur, Ref1, Ref2, 16);
408 :			emms();
409 :			t = (gettime_usec() - t) / nb_tests;
410 :			printf( "%s - sad16bi %.3f usec sad=%d %s\n",
411 :			cpu->name, t, s,
412 :			(s!=26274)?"\| ERROR": "" );
413 :	chl	898
414 :	Isibaar	1804	t = gettime_usec();
415 :			emms();
416 :			for(tst=0; tst<nb_tests; ++tst) s = sad8bi(Cur, Ref1, Ref2, 8);
417 :			emms();
418 :			t = (gettime_usec() - t) / nb_tests;
419 :			printf( "%s - sad8bi %.3f usec sad=%d %s\n",
420 :			cpu->name, t, s,
421 :			(s!=4002)?"\| ERROR": "" );
422 :
423 :	edgomez	1382	t = gettime_usec();
424 :			emms();
425 :			for(tst=0; tst<nb_tests; ++tst) s = dev16(Cur, 16);
426 :			emms();
427 :			t = (gettime_usec() - t) / nb_tests;
428 :			printf( "%s - dev16 %.3f usec sad=%d %s\n",
429 :			cpu->name, t, s,
430 :			(s!=3344)?"\| ERROR": "" );
431 :	Isibaar	225
432 :	edgomez	1382	printf( " --- \n" );
433 :			}
434 :	Isibaar	225	}
435 :
436 :			/*********************************************************************
437 :			* test interpolation
438 :			*********************************************************************/
439 :
440 :			#define ENTER \
441 :	edgomez	1382	for(i=0; i<16*8; ++i) Dst[i] = 0; \
442 :			t = gettime_usec(); \
443 :			emms();
444 :	Isibaar	225
445 :			#define LEAVE \
446 :	edgomez	1382	emms(); \
447 :			t = (gettime_usec() - t) / nb_tests; \
448 :			iCrc = calc_crc((uint8_t*)Dst, sizeof(Dst), CRC32_INITIAL)
449 :	Isibaar	225
450 :			#define TEST_MB(FUNC, R) \
451 :	edgomez	1382	ENTER \
452 :			for(tst=0; tst<nb_tests; ++tst) (FUNC)(Dst, Src0, 16, (R)); \
453 :			LEAVE
454 :	Isibaar	225
455 :			#define TEST_MB2(FUNC) \
456 :	edgomez	1382	ENTER \
457 :			for(tst=0; tst<nb_tests; ++tst) (FUNC)(Dst, Src0, 16); \
458 :			LEAVE
459 :	Isibaar	225
460 :
461 :			void test_mb()
462 :			{
463 :	edgomez	1382	const int nb_tests = 2000*speed_ref;
464 :			CPU *cpu;
465 :			const uint8_t Src0[16*9] = {
466 :			/* try to have every possible combinaison of rounding... */
467 :			0, 0, 1, 0, 2, 0, 3, 0, 4 ,0,0,0, 0,0,0,0,
468 :			0, 1, 1, 1, 2, 1, 3, 1, 3 ,0,0,0, 0,0,0,0,
469 :			0, 2, 1, 2, 2, 2, 3, 2, 2 ,0,0,0, 0,0,0,0,
470 :			0, 3, 1, 3, 2, 3, 3, 3, 1 ,0,0,0, 0,0,0,0,
471 :			1, 3, 0, 2, 1, 0, 2, 3, 4 ,0,0,0, 0,0,0,0,
472 :			2, 2, 1, 2, 0, 1, 3, 5, 3 ,0,0,0, 0,0,0,0,
473 :			3, 1, 2, 3, 1, 2, 2, 6, 2 ,0,0,0, 0,0,0,0,
474 :			1, 0, 1, 3, 0, 3, 1, 6, 1 ,0,0,0, 0,0,0,0,
475 :			4, 3, 2, 1, 2, 3, 4, 0, 3 ,0,0,0, 0,0,0,0
476 :			};
477 :			uint8_t Dst[16*8] = {0};
478 :	Isibaar	225
479 :	edgomez	1382	printf( "\n === test block motion ===\n" );
480 :	Isibaar	225
481 :	edgomez	1382	for(cpu = cpu_list; cpu->name!=0; ++cpu)
482 :			{
483 :			double t;
484 :			int tst, i, iCrc;
485 :	Isibaar	225
486 :	edgomez	1382	if (!init_cpu(cpu))
487 :			continue;
488 :	Isibaar	225
489 :	edgomez	1382	TEST_MB(interpolate8x8_halfpel_h, 0);
490 :			printf("%s - interp- h-round0 %.3f usec crc32=0x%08x %s\n",
491 :			cpu->name, t, iCrc,
492 :			(iCrc!=0x115381ba)?"\| ERROR": "" );
493 :	Isibaar	225
494 :	edgomez	1382	TEST_MB(interpolate8x8_halfpel_h, 1);
495 :			printf("%s - round1 %.3f usec crc32=0x%08x %s\n",
496 :			cpu->name, t, iCrc,
497 :			(iCrc!=0x2b1f528f)?"\| ERROR": "" );
498 :	Isibaar	225
499 :
500 :	edgomez	1382	TEST_MB(interpolate8x8_halfpel_v, 0);
501 :			printf("%s - interp- v-round0 %.3f usec crc32=0x%08x %s\n",
502 :			cpu->name, t, iCrc,
503 :			(iCrc!=0x423cdcc7)?"\| ERROR": "" );
504 :	Isibaar	225
505 :	edgomez	1382	TEST_MB(interpolate8x8_halfpel_v, 1);
506 :			printf("%s - round1 %.3f usec crc32=0x%08x %s\n",
507 :			cpu->name, t, iCrc,
508 :			(iCrc!=0x42202efe)?"\| ERROR": "" );
509 :	Isibaar	225
510 :
511 :	edgomez	1382	TEST_MB(interpolate8x8_halfpel_hv, 0);
512 :			printf("%s - interp-hv-round0 %.3f usec crc32=0x%08x %s\n",
513 :			cpu->name, t, iCrc,
514 :			(iCrc!=0xd198d387)?"\| ERROR": "" );
515 :	Isibaar	225
516 :	edgomez	1382	TEST_MB(interpolate8x8_halfpel_hv, 1);
517 :			printf("%s - round1 %.3f usec crc32=0x%08x %s\n",
518 :			cpu->name, t, iCrc,
519 :			(iCrc!=0x9ecfd921)?"\| ERROR": "" );
520 :	Isibaar	225
521 :	Isibaar	262
522 :	edgomez	1382	/* this is a new function, as of 06.06.2002 */
523 :	Isibaar	262	#if 0
524 :	edgomez	1382	TEST_MB2(interpolate8x8_avrg);
525 :			printf("%s - interpolate8x8_c %.3f usec crc32=0x%08x %s\n",
526 :			cpu->name, t, iCrc,
527 :			(iCrc!=8107)?"\| ERROR": "" );
528 :	Isibaar	262	#endif
529 :
530 :	Skal	1651	/* New functions for field prediction by CK 1.10.2005 */
531 :			#pragma NEW8X4
532 :			TEST_MB(interpolate8x4_halfpel_h, 0);
533 :			printf("%s - interpfield-h -round0 %.3f usec crc32=0x%08x %s\n",
534 :			cpu->name, t, iCrc,
535 :			(iCrc!=0x9538d6df)?"\| ERROR": "" );
536 :
537 :			TEST_MB(interpolate8x4_halfpel_h, 1);
538 :			printf("%s - round1 %.3f usec crc32=0x%08x %s\n",
539 :			cpu->name, t, iCrc,
540 :			(iCrc!=0xde5f1db4)?"\| ERROR": "" );
541 :
542 :
543 :			TEST_MB(interpolate8x4_halfpel_v, 0);
544 :			printf("%s - interpfield- v-round0 %.3f usec crc32=0x%08x %s\n",
545 :			cpu->name, t, iCrc,
546 :			(iCrc!=0xea5a69ef)?"\| ERROR": "" );
547 :
548 :			TEST_MB(interpolate8x4_halfpel_v, 1);
549 :			printf("%s - round1 %.3f usec crc32=0x%08x %s\n",
550 :			cpu->name, t, iCrc,
551 :			(iCrc!=0x4f10ec0f)?"\| ERROR": "" );
552 :
553 :
554 :			TEST_MB(interpolate8x4_halfpel_hv, 0);
555 :			printf("%s - interpfield-hv-round0 %.3f usec crc32=0x%08x %s\n",
556 :			cpu->name, t, iCrc,
557 :			(iCrc!=0xf97ee367)?"\| ERROR": "" );
558 :
559 :			TEST_MB(interpolate8x4_halfpel_hv, 1);
560 :			printf("%s - round1 %.3f usec crc32=0x%08x %s\n",
561 :			cpu->name, t, iCrc,
562 :			(iCrc!=0xb6a9f581)?"\| ERROR": "" );
563 :			/* End of 8x4 functions */
564 :
565 :	edgomez	1382	printf( " --- \n" );
566 :			}
567 :	Isibaar	225	}
568 :
569 :	Skal	1707	#undef ENTER
570 :			#undef LEAVE
571 :			#undef TEST_MB
572 :			#undef TEST_MB2
573 :
574 :	Isibaar	225	/*********************************************************************
575 :			* test transfer
576 :			*********************************************************************/
577 :
578 :			#define INIT_TRANSFER \
579 :	edgomez	1382	for(i=0; i<8*32; ++i) { \
580 :			Src8[i] = i; Src16[i] = i; \
581 :			Dst8[i] = 0; Dst16[i] = 0; \
582 :			Ref1[i] = i^0x27; \
583 :			Ref2[i] = i^0x51; \
584 :			}
585 :	Isibaar	225
586 :			#define TEST_TRANSFER_BEGIN(DST) \
587 :	edgomez	1382	INIT_TRANSFER \
588 :			overhead = -gettime_usec(); \
589 :			for(tst=0; tst<nb_tests; ++tst) { \
590 :			for(i=0; i<8*32; ++i) (DST)[i] = i^0x6a;\
591 :			} \
592 :			overhead += gettime_usec(); \
593 :			t = gettime_usec(); \
594 :			emms(); \
595 :			for(tst=0; tst<nb_tests; ++tst) { \
596 :			for(i=0; i<8*32; ++i) (DST)[i] = i^0x6a;
597 :	Isibaar	225
598 :
599 :			#define TEST_TRANSFER_END(DST) \
600 :	edgomez	1382	} \
601 :			emms(); \
602 :			t = (gettime_usec()-t -overhead) / nb_tests;\
603 :	Skal	1621	byte_swap((uint8_t)(DST), 832*sizeof((DST)[0]), sizeof((DST)[0])); \
604 :	Skal	1620	s = calc_crc((uint8_t)(DST), 832*sizeof((DST)[0]), CRC32_INITIAL)
605 :	Isibaar	225
606 :			#define TEST_TRANSFER(FUNC, DST, SRC) \
607 :	edgomez	1382	TEST_TRANSFER_BEGIN(DST); \
608 :			(FUNC)((DST), (SRC), 32); \
609 :			TEST_TRANSFER_END(DST)
610 :	Isibaar	225
611 :
612 :			#define TEST_TRANSFER2_BEGIN(DST, SRC) \
613 :	edgomez	1382	INIT_TRANSFER \
614 :			overhead = -gettime_usec(); \
615 :			for(tst=0; tst<nb_tests; ++tst) { \
616 :			for(i=0; i<8*32; ++i) (DST)[i] = i^0x6a;\
617 :			for(i=0; i<8*32; ++i) (SRC)[i] = i^0x3e;\
618 :			} \
619 :			overhead += gettime_usec(); \
620 :			t = gettime_usec(); \
621 :			emms(); \
622 :			for(tst=0; tst<nb_tests; ++tst) { \
623 :			for(i=0; i<8*32; ++i) (DST)[i] = i^0x6a;\
624 :			for(i=0; i<8*32; ++i) (SRC)[i] = i^0x3e;
625 :	Isibaar	225
626 :			#define TEST_TRANSFER2_END(DST) \
627 :	edgomez	1382	} \
628 :			emms(); \
629 :			t = (gettime_usec()-t -overhead) / nb_tests;\
630 :	Skal	1621	byte_swap((uint8_t)(DST), 832*sizeof((DST)[0]), sizeof((DST)[0])); \
631 :	Skal	1620	s = calc_crc((uint8_t)(DST), 832*sizeof((DST)[0]), CRC32_INITIAL)
632 :	Isibaar	225
633 :			#define TEST_TRANSFER2(FUNC, DST, SRC, R1) \
634 :	edgomez	1382	TEST_TRANSFER2_BEGIN(DST,SRC); \
635 :			(FUNC)((DST), (SRC), (R1), 32); \
636 :			TEST_TRANSFER2_END(DST)
637 :	Isibaar	225
638 :			#define TEST_TRANSFER3(FUNC, DST, SRC, R1, R2)\
639 :	edgomez	1382	TEST_TRANSFER_BEGIN(DST); \
640 :			(FUNC)((DST), (SRC), (R1), (R2), 32); \
641 :			TEST_TRANSFER_END(DST)
642 :	Isibaar	225
643 :			void test_transfer()
644 :			{
645 :	edgomez	1382	const int nb_tests = 4000*speed_ref;
646 :			int i;
647 :			CPU *cpu;
648 :	Skal	1620	// uint8_t Src8[832], Dst8[832], Ref1[832], Ref2[832];
649 :			// int16_t Src16[832], Dst16[832];
650 :			DECLARE_ALIGNED_MATRIX(Src8, 8, 32, uint8_t, CACHE_LINE);
651 :			DECLARE_ALIGNED_MATRIX(Dst8, 8, 32, uint8_t, CACHE_LINE);
652 :			DECLARE_ALIGNED_MATRIX(Ref1, 8, 32, uint8_t, CACHE_LINE);
653 :			DECLARE_ALIGNED_MATRIX(Ref2, 8, 32, uint8_t, CACHE_LINE);
654 :			DECLARE_ALIGNED_MATRIX(Src16, 8, 32, uint16_t, CACHE_LINE);
655 :			DECLARE_ALIGNED_MATRIX(Dst16, 8, 32, uint16_t, CACHE_LINE);
656 :	Isibaar	225
657 :	edgomez	1382	printf( "\n === test transfer ===\n" );
658 :	Isibaar	225
659 :	edgomez	1382	for(cpu = cpu_list; cpu->name!=0; ++cpu)
660 :			{
661 :			double t, overhead;
662 :			int tst, s;
663 :	Isibaar	225
664 :	edgomez	1382	if (!init_cpu(cpu))
665 :			continue;
666 :	Isibaar	225
667 :	edgomez	1382	TEST_TRANSFER(transfer_8to16copy, Dst16, Src8);
668 :			printf("%s - 8to16 %.3f usec crc32=0x%08x %s\n",
669 :			cpu->name, t, s,
670 :			(s!=0x115814bb)?"\| ERROR": "");
671 :	Isibaar	225
672 :	edgomez	1382	TEST_TRANSFER(transfer_16to8copy, Dst8, Src16);
673 :			printf( "%s - 16to8 %.3f usec crc32=0x%08x %s\n",
674 :			cpu->name, t, s,
675 :			(s!=0xee7ccbb4)?"\| ERROR": "");
676 :	Isibaar	225
677 :	Skal	1651	/* New functions for field prediction by CK 1.10.2005 */
678 :			#pragma NEW8X4
679 :			TEST_TRANSFER(transfer8x4_copy, Dst8, Src8);
680 :			printf("%s - 8to4 %.3f usec crc32=0x%08x %s\n",
681 :			cpu->name, t, s,
682 :			(s!=0xbb9c3db5)?"\| ERROR": "");
683 :			/* End of new functions */
684 :
685 :	edgomez	1382	TEST_TRANSFER(transfer8x8_copy, Dst8, Src8);
686 :			printf("%s - 8to8 %.3f usec crc32=0x%08x %s\n",
687 :			cpu->name, t, s,
688 :			(s!=0xd37b3295)?"\| ERROR": "");
689 :	Isibaar	225
690 :	edgomez	1382	TEST_TRANSFER(transfer_16to8add, Dst8, Src16);
691 :			printf("%s - 16to8add %.3f usec crc32=0x%08x %s\n",
692 :			cpu->name, t, s,
693 :			(s!=0xdd817bf4)?"\| ERROR": "" );
694 :	Isibaar	225
695 :	edgomez	1382	TEST_TRANSFER2(transfer_8to16sub, Dst16, Src8, Ref1);
696 :			{
697 :			int s1, s2;
698 :	Skal	1620	s1 = calc_crc((uint8_t)Dst16, 832*sizeof(Dst16[0]), CRC32_INITIAL);
699 :			s2 = calc_crc((uint8_t)Src8, 832*sizeof(Src8[0]), CRC32_INITIAL);
700 :	edgomez	1382	printf("%s - 8to16sub %.3f usec crc32(1)=0x%08x crc32(2)=0x%08x %s %s\n",
701 :			cpu->name, t, s1, s2,
702 :			(s1!=0xa1e07163)?"\| ERROR1": "",
703 :			(s2!=0xd86c5d23)?"\| ERROR2": "" );
704 :			}
705 :
706 :			TEST_TRANSFER3(transfer_8to16sub2, Dst16, Src8, Ref1, Ref2);
707 :			printf("%s - 8to16sub2 %.3f usec crc32=0x%08x %s\n",
708 :			cpu->name, t, s,
709 :			(s!=0x99b6c4c7)?"\| ERROR": "" );
710 :
711 :			printf( " --- \n" );
712 :			}
713 :	Isibaar	225	}
714 :
715 :			/*********************************************************************
716 :			* test quantization
717 :			*********************************************************************/
718 :
719 :	Isibaar	262	#define TEST_QUANT(FUNC, DST, SRC) \
720 :	edgomez	1382	t = gettime_usec(); \
721 :			for(s=CRC32_INITIAL,qm=1; qm<=255; ++qm) { \
722 :			for(i=0; i<8*8; ++i) Quant[i] = qm; \
723 :			set_inter_matrix( mpeg_quant_matrices, Quant ); \
724 :			emms(); \
725 :			for(q=1; q<=max_Q; ++q) { \
726 :			for(tst=0; tst<nb_tests; ++tst) \
727 :			(FUNC)((DST), (SRC), q, mpeg_quant_matrices); \
728 :	Skal	1621	byte_swap((uint8_t)(DST), 64sizeof((DST)[0]), sizeof((DST)[0])); \
729 :			s = calc_crc((uint8_t)(DST), 64sizeof((DST)[0]), s); \
730 :	edgomez	1382	} \
731 :			emms(); \
732 :			} \
733 :			t = (gettime_usec()-t-overhead)/nb_tests/qm
734 :	Isibaar	225
735 :	Isibaar	262	#define TEST_QUANT2(FUNC, DST, SRC) \
736 :	edgomez	1382	t = gettime_usec(); \
737 :			for(s=CRC32_INITIAL,qm=1; qm<=255; ++qm) { \
738 :			for(i=0; i<8*8; ++i) Quant[i] = qm; \
739 :			set_intra_matrix( mpeg_quant_matrices, Quant ); \
740 :			emms(); \
741 :			for(q=1; q<=max_Q; ++q) { \
742 :	Isibaar	1801	init_intra_matrix( mpeg_quant_matrices, q ); \
743 :	edgomez	1382	for(tst=0; tst<nb_tests; ++tst) \
744 :			(FUNC)((DST), (SRC), q, q, mpeg_quant_matrices); \
745 :	Skal	1621	byte_swap((uint8_t)(DST), 64sizeof((DST)[0]), sizeof((DST)[0])); \
746 :			s = calc_crc((uint8_t)(DST), 64sizeof((DST)[0]), s); \
747 :	edgomez	1382	} \
748 :			emms(); \
749 :			} \
750 :			t = (gettime_usec()-t-overhead)/nb_tests/qm
751 :	Isibaar	225
752 :	edgomez	1398	#define TEST_INTRA(REFFUNC, NEWFUNC, RANGE) \
753 :	Isibaar	1801	{ int32_t i,q,s;\
754 :	edgomez	1398	DECLARE_ALIGNED_MATRIX(Src, 8, 8, int16_t, 16); \
755 :			DECLARE_ALIGNED_MATRIX(Dst, 8, 8, int16_t, 16); \
756 :			DECLARE_ALIGNED_MATRIX(Dst2,8, 8, int16_t, 16); \
757 :			for(q=1;q<=max_Q;q++) \
758 :			for(s=-RANGE;s<RANGE;s++) { \
759 :			for(i=0;i<64;i++) Src[i]=s; \
760 :			(REFFUNC)((Dst),(Src),q,q,mpeg_quant_matrices); \
761 :			(NEWFUNC)((Dst2),(Src),q,q,mpeg_quant_matrices); \
762 :			for(i=0;i<64;i++) \
763 :			if(Dst[i]!=Dst2[i]) printf("ERROR : " #NEWFUNC " i%d quant:%d input:%d C_result:%d ASM_result:%d\n",i,q,s,Dst[i],Dst2[i]); \
764 :			} \
765 :			}
766 :
767 :			#define TEST_INTER(REFFUNC, NEWFUNC, RANGE) \
768 :			{ int i,q,s; \
769 :			DECLARE_ALIGNED_MATRIX(Src, 8, 8, int16_t, 16); \
770 :			DECLARE_ALIGNED_MATRIX(Dst, 8, 8, int16_t, 16); \
771 :			DECLARE_ALIGNED_MATRIX(Dst2,8, 8, int16_t, 16); \
772 :			for(q=1;q<=max_Q;q++) \
773 :			for(s=-RANGE;s<RANGE;s++) { \
774 :			for(i=0;i<64;i++) Src[i]=s; \
775 :			(REFFUNC)((Dst),(Src),q,mpeg_quant_matrices); \
776 :			(NEWFUNC)((Dst2),(Src),q,mpeg_quant_matrices); \
777 :			emms(); \
778 :			for(i=0;i<64;i++) \
779 :			if(Dst[i]!=Dst2[i]) printf("ERROR : " #NEWFUNC " i%d quant:%d input:%d C_result:%d ASM_result:%d\n",i,q,s,Dst[i],Dst2[i]); \
780 :			} \
781 :			}
782 :
783 :	Isibaar	225	void test_quant()
784 :			{
785 :	Isibaar	1801	const int32_t nb_tests = 1*speed_ref;
786 :			const int32_t max_Q = 31;
787 :	edgomez	1382	DECLARE_ALIGNED_MATRIX(mpeg_quant_matrices, 8, 64, uint16_t, 16);
788 :	Isibaar	225
789 :	Isibaar	1801	int32_t i, qm;
790 :	edgomez	1382	CPU *cpu;
791 :			DECLARE_ALIGNED_MATRIX(Src, 8, 8, int16_t, 16);
792 :			DECLARE_ALIGNED_MATRIX(Dst, 8, 8, int16_t, 16);
793 :	edgomez	1398	DECLARE_ALIGNED_MATRIX(Dst2,8, 8, int16_t, 16);
794 :	edgomez	1382	uint8_t Quant[8*8];
795 :	Isibaar	225
796 :	edgomez	1382	printf( "\n ===== test quant =====\n" );
797 :	Isibaar	225
798 :	edgomez	1382	/* we deliberately enfringe the norm's specified range [-127,127], */
799 :			/* to test the robustness of the iquant module */
800 :			for(i=0; i<64; ++i) {
801 :			Src[i] = 1 + (i-32) * (i&6);
802 :			Dst[i] = 0;
803 :			}
804 :	Isibaar	225
805 :	edgomez	1382	for(cpu = cpu_list; cpu->name!=0; ++cpu)
806 :			{
807 :			double t, overhead;
808 :	Isibaar	1801	int32_t tst, q;
809 :	edgomez	1382	uint32_t s;
810 :	Isibaar	225
811 :	edgomez	1382	if (!init_cpu(cpu))
812 :	edgomez	1398	continue;
813 :	Isibaar	225
814 :	edgomez	1398	// exhaustive tests to compare against the (ref) C-version
815 :			TEST_INTRA(quant_h263_intra_c, quant_h263_intra, 2048);
816 :			TEST_INTRA(dequant_h263_intra_c, dequant_h263_intra , 512 );
817 :			TEST_INTER(quant_h263_inter_c, quant_h263_inter , 2048);
818 :			TEST_INTER(dequant_h263_inter_c, dequant_h263_inter , 512 );
819 :
820 :	edgomez	1382	overhead = -gettime_usec();
821 :			for(s=0,qm=1; qm<=255; ++qm) {
822 :			for(i=0; i<8*8; ++i) Quant[i] = qm;
823 :			set_inter_matrix(mpeg_quant_matrices, Quant );
824 :			for(q=1; q<=max_Q; ++q)
825 :			for(i=0; i<64; ++i) s+=Dst[i]^i^qm;
826 :			}
827 :			overhead += gettime_usec();
828 :	Isibaar	225
829 :	edgomez	1382	TEST_QUANT2(quant_mpeg_intra, Dst, Src);
830 :			printf("%s - quant_mpeg_intra %.3f usec crc32=0x%08x %s\n",
831 :			cpu->name, t, s,
832 :	Isibaar	1801	(s!=0x3b999af6)? "\| ERROR": "");
833 :	Isibaar	225
834 :	edgomez	1382	TEST_QUANT(quant_mpeg_inter, Dst, Src);
835 :			printf("%s - quant_mpeg_inter %.3f usec crc32=0x%08x %s\n",
836 :			cpu->name, t, s,
837 :			(s!=0xf6de7757)?"\| ERROR": "");
838 :	Isibaar	225
839 :	edgomez	1382	TEST_QUANT2(dequant_mpeg_intra, Dst, Src);
840 :			printf("%s - dequant_mpeg_intra %.3f usec crc32=0x%08x %s\n",
841 :			cpu->name, t, s,
842 :			(s!=0x2def7bc7)?"\| ERROR": "");
843 :	Isibaar	225
844 :	edgomez	1382	TEST_QUANT(dequant_mpeg_inter, Dst, Src);
845 :			printf("%s - dequant_mpeg_inter %.3f usec crc32=0x%08x %s\n",
846 :			cpu->name, t, s,
847 :			(s!=0xd878c722)?"\| ERROR": "");
848 :
849 :			TEST_QUANT2(quant_h263_intra, Dst, Src);
850 :			printf("%s - quant_h263_intra %.3f usec crc32=0x%08x %s\n",
851 :			cpu->name, t, s,
852 :			(s!=0x2eba9d43)?"\| ERROR": "");
853 :
854 :			TEST_QUANT(quant_h263_inter, Dst, Src);
855 :			printf("%s - quant_h263_inter %.3f usec crc32=0x%08x %s\n",
856 :			cpu->name, t, s,
857 :			(s!=0xbd315a7e)?"\| ERROR": "");
858 :
859 :			TEST_QUANT2(dequant_h263_intra, Dst, Src);
860 :			printf("%s - dequant_h263_intra %.3f usec crc32=0x%08x %s\n",
861 :			cpu->name, t, s,
862 :			(s!=0x9841212a)?"\| ERROR": "");
863 :
864 :			TEST_QUANT(dequant_h263_inter, Dst, Src);
865 :			printf("%s - dequant_h263_inter %.3f usec crc32=0x%08x %s\n",
866 :			cpu->name, t, s,
867 :			(s!=0xe7df8fba)?"\| ERROR": "");
868 :
869 :			printf( " --- \n" );
870 :			}
871 :	Isibaar	225	}
872 :
873 :			/*********************************************************************
874 :	edgomez	1424	* test distortion operators
875 :			*********************************************************************/
876 :
877 :			static void ieee_reseed(long s);
878 :			static long ieee_rand(int Min, int Max);
879 :
880 :			#define TEST_SSE(FUNCTION, SRC1, SRC2, STRIDE) \
881 :			do { \
882 :			t = gettime_usec(); \
883 :			tst = nb_tests; \
884 :			while((tst--)>0) sse = (FUNCTION)((SRC1), (SRC2), (STRIDE)); \
885 :			emms(); \
886 :			t = (gettime_usec() - t)/(double)nb_tests; \
887 :			} while(0)
888 :
889 :
890 :			void test_sse()
891 :			{
892 :			const int nb_tests = 100000*speed_ref;
893 :			int i;
894 :			CPU *cpu;
895 :			DECLARE_ALIGNED_MATRIX(Src1, 8, 8, int16_t, 16);
896 :			DECLARE_ALIGNED_MATRIX(Src2, 8, 8, int16_t, 16);
897 :			DECLARE_ALIGNED_MATRIX(Src3, 8, 8, int16_t, 16);
898 :			DECLARE_ALIGNED_MATRIX(Src4, 8, 8, int16_t, 16);
899 :
900 :			printf( "\n ===== test sse =====\n" );
901 :
902 :			ieee_reseed(1);
903 :			for(i=0; i<64; ++i) {
904 :			Src1[i] = ieee_rand(-2048, 2047);
905 :			Src2[i] = ieee_rand(-2048, 2047);
906 :			Src3[i] = ieee_rand(-2048, 2047);
907 :			Src4[i] = ieee_rand(-2048, 2047);
908 :			}
909 :
910 :			for(cpu = cpu_list; cpu->name!=0; ++cpu)
911 :			{
912 :			double t;
913 :			int tst, sse;
914 :
915 :			if (!init_cpu(cpu))
916 :			continue;
917 :
918 :			/* 16 bit element blocks */
919 :			TEST_SSE(sse8_16bit, Src1, Src2, 16);
920 :			printf("%s - sse8_16bit#1 %.3f usec sse=%d %s\n",
921 :			cpu->name, t, sse, (sse!=182013834)?"\| ERROR": "");
922 :			TEST_SSE(sse8_16bit, Src1, Src3, 16);
923 :			printf("%s - sse8_16bit#2 %.3f usec sse=%d %s\n",
924 :			cpu->name, t, sse, (sse!=142545203)?"\| ERROR": "");
925 :			TEST_SSE(sse8_16bit, Src1, Src4, 16);
926 :			printf("%s - sse8_16bit#3 %.3f usec sse=%d %s\n",
927 :			cpu->name, t, sse, (sse!=146340935)?"\| ERROR": "");
928 :			TEST_SSE(sse8_16bit, Src2, Src3, 16);
929 :			printf("%s - sse8_16bit#4 %.3f usec sse=%d %s\n",
930 :			cpu->name, t, sse, (sse!=130136661)?"\| ERROR": "");
931 :			TEST_SSE(sse8_16bit, Src2, Src4, 16);
932 :			printf("%s - sse8_16bit#5 %.3f usec sse=%d %s\n",
933 :			cpu->name, t, sse, (sse!=136870353)?"\| ERROR": "");
934 :			TEST_SSE(sse8_16bit, Src3, Src4, 16);
935 :			printf("%s - sse8_16bit#6 %.3f usec sse=%d %s\n",
936 :			cpu->name, t, sse, (sse!=164107772)?"\| ERROR": "");
937 :
938 :			/* 8 bit element blocks */
939 :			TEST_SSE(sse8_8bit, (int8_t)Src1, (int8_t)Src2, 8);
940 :			printf("%s - sse8_8bit#1 %.3f usec sse=%d %s\n",
941 :			cpu->name, t, sse, (sse!=1356423)?"\| ERROR": "");
942 :			TEST_SSE(sse8_8bit, (int8_t)Src1, (int8_t)Src3, 8);
943 :			printf("%s - sse8_8bit#2 %.3f usec sse=%d %s\n",
944 :			cpu->name, t, sse, (sse!=1173074)?"\| ERROR": "");
945 :			TEST_SSE(sse8_8bit, (int8_t)Src1, (int8_t)Src4, 8);
946 :			printf("%s - sse8_8bit#3 %.3f usec sse=%d %s\n",
947 :			cpu->name, t, sse, (sse!=1092357)?"\| ERROR": "");
948 :			TEST_SSE(sse8_8bit, (int8_t)Src2, (int8_t)Src3, 8);
949 :			printf("%s - sse8_8bit#4 %.3f usec sse=%d %s\n",
950 :			cpu->name, t, sse, (sse!=1360239)?"\| ERROR": "");
951 :			TEST_SSE(sse8_8bit, (int8_t)Src2, (int8_t)Src4, 8);
952 :			printf("%s - sse8_8bit#5 %.3f usec sse=%d %s\n",
953 :			cpu->name, t, sse, (sse!=1208414)?"\| ERROR": "");
954 :			TEST_SSE(sse8_8bit, (int8_t)Src3, (int8_t)Src4, 8);
955 :			printf("%s - sse8_8bit#6 %.3f usec sse=%d %s\n",
956 :			cpu->name, t, sse, (sse!=1099285)?"\| ERROR": "");
957 :
958 :			printf(" ---\n");
959 :			}
960 :			}
961 :
962 :			/*********************************************************************
963 :	Isibaar	225	* test non-zero AC counting
964 :			*********************************************************************/
965 :
966 :	Skal	1629	#define TEST_CBP(FUNC, SRC, NB) \
967 :	edgomez	1382	t = gettime_usec(); \
968 :			emms(); \
969 :	Skal	1629	for(tst=0; tst<NB; ++tst) { \
970 :	edgomez	1382	cbp = (FUNC)((SRC)); \
971 :			} \
972 :			emms(); \
973 :			t = (gettime_usec()-t ) / nb_tests;
974 :	Isibaar	225
975 :			void test_cbp()
976 :			{
977 :	edgomez	1382	const int nb_tests = 10000*speed_ref;
978 :	Skal	1629	int i, n, m;
979 :	edgomez	1382	CPU *cpu;
980 :			DECLARE_ALIGNED_MATRIX(Src1, 6, 64, int16_t, 16);
981 :			DECLARE_ALIGNED_MATRIX(Src2, 6, 64, int16_t, 16);
982 :			DECLARE_ALIGNED_MATRIX(Src3, 6, 64, int16_t, 16);
983 :			DECLARE_ALIGNED_MATRIX(Src4, 6, 64, int16_t, 16);
984 :	suxen_drol	1641	DECLARE_ALIGNED_MATRIX(Src5, 6, 64, int16_t, 16);
985 :	Isibaar	225
986 :	edgomez	1382	printf( "\n ===== test cbp =====\n" );
987 :	Isibaar	225
988 :	edgomez	1382	for(i=0; i<6*64; ++i) {
989 :			Src1[i] = (ii3/8192)&(i/64)&1; /* 'random' */
990 :			Src2[i] = (i<364); / half-full */
991 :			Src3[i] = ((i+32)>3*64);
992 :			Src4[i] = (i==(364+2) \|\| i==(564+9));
993 :	suxen_drol	1641	Src5[i] = ieee_rand(0,1) ? -1 : 1; /* +/- test */
994 :	edgomez	1382	}
995 :	Isibaar	225
996 :	edgomez	1382	for(cpu = cpu_list; cpu->name!=0; ++cpu)
997 :			{
998 :			double t;
999 :			int tst, cbp;
1000 :	Isibaar	225
1001 :	edgomez	1382	if (!init_cpu(cpu))
1002 :			continue;
1003 :
1004 :	Skal	1629	TEST_CBP(calc_cbp, Src1, nb_tests);
1005 :	edgomez	1424	printf("%s - calc_cbp#1 %.3f usec cbp=0x%02x %s\n",
1006 :	edgomez	1382	cpu->name, t, cbp, (cbp!=0x15)?"\| ERROR": "");
1007 :	Skal	1629	TEST_CBP(calc_cbp, Src2, nb_tests);
1008 :	edgomez	1424	printf("%s - calc_cbp#2 %.3f usec cbp=0x%02x %s\n",
1009 :	edgomez	1382	cpu->name, t, cbp, (cbp!=0x38)?"\| ERROR": "");
1010 :	Skal	1629	TEST_CBP(calc_cbp, Src3, nb_tests);
1011 :	edgomez	1424	printf("%s - calc_cbp#3 %.3f usec cbp=0x%02x %s\n",
1012 :	edgomez	1382	cpu->name, t, cbp, (cbp!=0x0f)?"\| ERROR": "" );
1013 :	Skal	1629	TEST_CBP(calc_cbp, Src4, nb_tests);
1014 :	edgomez	1424	printf("%s - calc_cbp#4 %.3f usec cbp=0x%02x %s\n",
1015 :	edgomez	1382	cpu->name, t, cbp, (cbp!=0x05)?"\| ERROR": "" );
1016 :	suxen_drol	1641	TEST_CBP(calc_cbp, Src5, nb_tests);
1017 :			printf("%s - calc_cbp#4 %.3f usec cbp=0x%02x %s\n",
1018 :			cpu->name, t, cbp, (cbp!=0x3f)?"\| ERROR": "" );
1019 :	edgomez	1382	printf( " --- \n" );
1020 :			}
1021 :	Skal	1629
1022 :			for(cpu = cpu_list; cpu->name!=0; ++cpu) /* bench suggested by Carlo (carlo dot bramix at libero dot it) */
1023 :			{
1024 :			double t;
1025 :			int tst, cbp, err;
1026 :
1027 :			if (!init_cpu(cpu))
1028 :			continue;
1029 :
1030 :			err = 0;
1031 :			for(n=0; n<6; ++n)
1032 :			{
1033 :			for(m=0; m<64; ++m)
1034 :			{
1035 :			for(i=0; i<6*64; ++i)
1036 :			Src1[i] = (i== (m + n*64));
1037 :
1038 :			TEST_CBP(calc_cbp, Src1, 1);
1039 :			if (cbp!= (((m!=0)<<(5-n))))
1040 :			{
1041 :			printf( "%s - calc_cbp#5: ERROR at pos %d / %d!\n", cpu->name, n, m);
1042 :			err = 1;
1043 :			break;
1044 :			}
1045 :			}
1046 :			}
1047 :			if (!err)
1048 :			printf( " %s - calc_cbp#5 : OK\n", cpu->name );
1049 :
1050 :			}
1051 :	Isibaar	225	}
1052 :
1053 :			/*********************************************************************
1054 :	Isibaar	262	* fdct/idct IEEE1180 compliance
1055 :			*********************************************************************/
1056 :
1057 :			typedef struct {
1058 :	edgomez	1382	long Errors[64];
1059 :			long Sqr_Errors[64];
1060 :			long Max_Errors[64];
1061 :			long Nb;
1062 :	Isibaar	262	} STATS_8x8;
1063 :
1064 :			void init_stats(STATS_8x8 *S)
1065 :			{
1066 :	edgomez	1382	int i;
1067 :			for(i=0; i<64; ++i) {
1068 :			S->Errors[i] = 0;
1069 :			S->Sqr_Errors[i] = 0;
1070 :			S->Max_Errors[i] = 0;
1071 :			}
1072 :			S->Nb = 0;
1073 :	Isibaar	262	}
1074 :
1075 :			void store_stats(STATS_8x8 *S, short Blk[64], short Ref[64])
1076 :			{
1077 :	edgomez	1382	int i;
1078 :			for(i=0; i<64; ++i)
1079 :			{
1080 :			short Err = Blk[i] - Ref[i];
1081 :			S->Errors[i] += Err;
1082 :			S->Sqr_Errors[i] += Err * Err;
1083 :			if (Err<0) Err = -Err;
1084 :			if (S->Max_Errors[i]<Err)
1085 :			S->Max_Errors[i] = Err;
1086 :			}
1087 :			S->Nb++;
1088 :	Isibaar	262	}
1089 :
1090 :			void print_stats(STATS_8x8 *S)
1091 :			{
1092 :	edgomez	1382	int i;
1093 :			double Norm;
1094 :	Isibaar	262
1095 :	edgomez	1382	assert(S->Nb>0);
1096 :			Norm = 1. / (double)S->Nb;
1097 :			printf("\n== Max absolute values of errors ==\n");
1098 :			for(i=0; i<64; i++) {
1099 :			printf(" %4ld", S->Max_Errors[i]);
1100 :			if ((i&7)==7) printf("\n");
1101 :			}
1102 :	Isibaar	262
1103 :	edgomez	1382	printf("\n== Mean square errors ==\n");
1104 :			for(i=0; i<64; i++)
1105 :			{
1106 :			double Err = Norm * (double)S->Sqr_Errors[i];
1107 :			printf(" %.3f", Err);
1108 :			if ((i&7)==7) printf("\n");
1109 :			}
1110 :	Isibaar	262
1111 :	edgomez	1382	printf("\n== Mean errors ==\n");
1112 :			for(i=0; i<64; i++)
1113 :			{
1114 :			double Err = Norm * (double)S->Errors[i];
1115 :			printf(" %.3f", Err);
1116 :			if ((i&7)==7) printf("\n");
1117 :			}
1118 :			printf("\n");
1119 :	Isibaar	262	}
1120 :
1121 :			static const char *CHECK(double v, double l) {
1122 :	edgomez	1382	if (fabs(v)<=l) return "ok";
1123 :			else return "FAIL!";
1124 :	Isibaar	262	}
1125 :
1126 :			void report_stats(STATS_8x8 S, const double Limits)
1127 :			{
1128 :	edgomez	1382	int i;
1129 :			double Norm, PE, PMSE, OMSE, PME, OME;
1130 :	Isibaar	262
1131 :	edgomez	1382	assert(S->Nb>0);
1132 :			Norm = 1. / (double)S->Nb;
1133 :			PE = 0.;
1134 :			for(i=0; i<64; i++) {
1135 :			if (PE<S->Max_Errors[i])
1136 :			PE = S->Max_Errors[i];
1137 :			}
1138 :	Isibaar	262
1139 :	edgomez	1382	PMSE = 0.;
1140 :			OMSE = 0.;
1141 :			for(i=0; i<64; i++)
1142 :			{
1143 :			double Err = Norm * (double)S->Sqr_Errors[i];
1144 :			OMSE += Err;
1145 :			if (PMSE < Err) PMSE = Err;
1146 :			}
1147 :			OMSE /= 64.;
1148 :	Isibaar	262
1149 :	edgomez	1382	PME = 0.;
1150 :			OME = 0.;
1151 :			for(i=0; i<64; i++)
1152 :			{
1153 :			double Err = Norm * (double)S->Errors[i];
1154 :			OME += Err;
1155 :			Err = fabs(Err);
1156 :			if (PME < Err) PME = Err;
1157 :			}
1158 :			OME /= 64.;
1159 :	Isibaar	262
1160 :	edgomez	1382	printf( "Peak error: %4.4f\n", PE );
1161 :			printf( "Peak MSE: %4.4f\n", PMSE );
1162 :			printf( "Overall MSE: %4.4f\n", OMSE );
1163 :			printf( "Peak ME: %4.4f\n", PME );
1164 :			printf( "Overall ME: %4.4f\n", OME );
1165 :
1166 :			if (Limits!=0)
1167 :			{
1168 :			printf( "[PE<=%.4f %s] ", Limits[0], CHECK(PE, Limits[0]) );
1169 :			printf( "\n" );
1170 :			printf( "[PMSE<=%.4f %s]", Limits[1], CHECK(PMSE, Limits[1]) );
1171 :			printf( "[OMSE<=%.4f %s]", Limits[2], CHECK(OMSE, Limits[2]) );
1172 :			printf( "\n" );
1173 :			printf( "[PME<=%.4f %s] ", Limits[3], CHECK(PME , Limits[3]) );
1174 :			printf( "[OME<=%.4f %s] ", Limits[4], CHECK(OME , Limits[4]) );
1175 :			printf( "\n" );
1176 :			}
1177 :	Isibaar	262	}
1178 :
1179 :	edgomez	1382	///* ////////////////////////////////////////////////////// */
1180 :	Isibaar	262	/* Pseudo-random generator specified by IEEE 1180 */
1181 :
1182 :			static long ieee_seed = 1;
1183 :			static void ieee_reseed(long s) {
1184 :	edgomez	1382	ieee_seed = s;
1185 :	Isibaar	262	}
1186 :			static long ieee_rand(int Min, int Max)
1187 :			{
1188 :	edgomez	1382	static double z = (double) 0x7fffffff;
1189 :	Isibaar	262
1190 :	edgomez	1382	long i,j;
1191 :			double x;
1192 :	Isibaar	262
1193 :	edgomez	1382	ieee_seed = (ieee_seed * 1103515245) + 12345;
1194 :			i = ieee_seed & 0x7ffffffe;
1195 :			x = ((double) i) / z;
1196 :			x *= (Max-Min+1);
1197 :			j = (long)x;
1198 :			j = j + Min;
1199 :			assert(j>=Min && j<=Max);
1200 :			return (short)j;
1201 :	Isibaar	262	}
1202 :
1203 :			#define CLAMP(x, M) (x) = ((x)<-(M)) ? (-(M)) : ((x)>=(M) ? ((M)-1) : (x))
1204 :
1205 :			static double Cos[8][8];
1206 :			static void init_ref_dct()
1207 :			{
1208 :	edgomez	1382	int i, j;
1209 :			for(i=0; i<8; i++)
1210 :			{
1211 :			double scale = (i == 0) ? sqrt(0.125) : 0.5;
1212 :			for (j=0; j<8; j++)
1213 :			Cos[i][j] = scalecos( (M_PI/8.0)i*(j + 0.5) );
1214 :			}
1215 :	Isibaar	262	}
1216 :
1217 :			void ref_idct(short *M)
1218 :			{
1219 :	edgomez	1382	int i, j, k;
1220 :			double Tmp[8][8];
1221 :	Isibaar	262
1222 :	edgomez	1382	for(i=0; i<8; i++) {
1223 :			for(j=0; j<8; j++)
1224 :			{
1225 :			double Sum = 0.0;
1226 :			for (k=0; k<8; k++) Sum += Cos[k][j]M[8i+k];
1227 :			Tmp[i][j] = Sum;
1228 :			}
1229 :			}
1230 :			for(i=0; i<8; i++) {
1231 :			for(j=0; j<8; j++) {
1232 :			double Sum = 0.0;
1233 :			for (k=0; k<8; k++) Sum += Cos[k][i]*Tmp[k][j];
1234 :			M[8*i+j] = (short)floor(Sum + .5);
1235 :			}
1236 :			}
1237 :	Isibaar	262	}
1238 :
1239 :			void ref_fdct(short *M)
1240 :			{
1241 :	edgomez	1382	int i, j, k;
1242 :			double Tmp[8][8];
1243 :	Isibaar	262
1244 :	edgomez	1382	for(i=0; i<8; i++) {
1245 :			for(j=0; j<8; j++)
1246 :			{
1247 :			double Sum = 0.0;
1248 :			for (k=0; k<8; k++) Sum += Cos[j][k]M[8i+k];
1249 :			Tmp[i][j] = Sum;
1250 :			}
1251 :			}
1252 :			for(i=0; i<8; i++) {
1253 :			for(j=0; j<8; j++) {
1254 :			double Sum = 0.0;
1255 :			for (k=0; k<8; k++) Sum += Cos[i][k]*Tmp[k][j];
1256 :			M[8*i+j] = (short)floor(Sum + 0.5);
1257 :			}
1258 :			}
1259 :	Isibaar	262	}
1260 :
1261 :			void test_IEEE1180_compliance(int Min, int Max, int Sign)
1262 :			{
1263 :	edgomez	1382	static const double ILimits[5] = { 1., 0.06, 0.02, 0.015, 0.0015 };
1264 :			int Loops = 10000;
1265 :			int i, m, n;
1266 :			DECLARE_ALIGNED_MATRIX(Blk0, 8, 8, short, 16); /* reference */
1267 :			DECLARE_ALIGNED_MATRIX(Blk, 8, 8, short, 16);
1268 :			DECLARE_ALIGNED_MATRIX(iBlk, 8, 8, short, 16);
1269 :			DECLARE_ALIGNED_MATRIX(Ref_FDCT, 8, 8, short, 16);
1270 :			DECLARE_ALIGNED_MATRIX(Ref_IDCT, 8, 8, short, 16);
1271 :	Isibaar	262
1272 :	edgomez	1382	STATS_8x8 FStats; /* forward dct stats */
1273 :			STATS_8x8 IStats; /* inverse dct stats */
1274 :	Isibaar	262
1275 :	edgomez	1382	CPU *cpu;
1276 :	Isibaar	262
1277 :	edgomez	1382	init_ref_dct();
1278 :	Isibaar	262
1279 :	edgomez	1382	for(cpu = cpu_list; cpu->name!=0; ++cpu)
1280 :			{
1281 :			if (!init_cpu(cpu))
1282 :			continue;
1283 :	Isibaar	262
1284 :	edgomez	1382	printf( "\n===== IEEE test for %s ==== (Min=%d Max=%d Sign=%d Loops=%d)\n",
1285 :			cpu->name, Min, Max, Sign, Loops);
1286 :	Isibaar	262
1287 :	edgomez	1382	init_stats(&IStats);
1288 :			init_stats(&FStats);
1289 :	Isibaar	262
1290 :	edgomez	1382	ieee_reseed(1);
1291 :			for(n=0; n<Loops; ++n)
1292 :			{
1293 :			for(i=0; i<64; ++i)
1294 :			Blk0[i] = (short)ieee_rand(Min,Max) * Sign;
1295 :	Isibaar	262
1296 :	edgomez	1382	/* hmm, I'm not quite sure this is exactly */
1297 :			/* the tests described in the norm. check... */
1298 :	Isibaar	262
1299 :	edgomez	1382	memcpy(Ref_FDCT, Blk0, 64*sizeof(short));
1300 :			ref_fdct(Ref_FDCT);
1301 :			for(i=0; i<64; i++) CLAMP( Ref_FDCT[i], 2048 );
1302 :	Isibaar	262
1303 :	edgomez	1382	memcpy(Blk, Blk0, 64*sizeof(short));
1304 :			emms(); fdct(Blk); emms();
1305 :			for(i=0; i<64; i++) CLAMP( Blk[i], 2048 );
1306 :	Isibaar	262
1307 :	edgomez	1382	store_stats(&FStats, Blk, Ref_FDCT);
1308 :	Isibaar	262
1309 :
1310 :	edgomez	1382	memcpy(Ref_IDCT, Ref_FDCT, 64*sizeof(short));
1311 :			ref_idct(Ref_IDCT);
1312 :			for (i=0; i<64; i++) CLAMP( Ref_IDCT[i], 256 );
1313 :	Isibaar	262
1314 :	edgomez	1382	memcpy(iBlk, Ref_FDCT, 64*sizeof(short));
1315 :			emms(); idct(iBlk); emms();
1316 :			for(i=0; i<64; i++) CLAMP( iBlk[i], 256 );
1317 :	Isibaar	262
1318 :	edgomez	1382	store_stats(&IStats, iBlk, Ref_IDCT);
1319 :			}
1320 :	Isibaar	262
1321 :
1322 :	edgomez	1382	printf( "\n -- FDCT report --\n" );
1323 :	edgomez	851	// print_stats(&FStats);
1324 :	edgomez	1382	report_stats(&FStats, 0); /* so far I know, IEEE1180 says nothing for fdct */
1325 :	Isibaar	262
1326 :	edgomez	1382	for(i=0; i<64; i++) Blk[i] = 0;
1327 :			emms(); fdct(Blk); emms();
1328 :			for(m=i=0; i<64; i++) if (Blk[i]!=0) m++;
1329 :			printf( "FDCT(0) == 0 ? %s\n", (m!=0) ? "NOPE!" : "yup." );
1330 :	Isibaar	262
1331 :	edgomez	1382	printf( "\n -- IDCT report --\n" );
1332 :	edgomez	851	// print_stats(&IStats);
1333 :	edgomez	1382	report_stats(&IStats, ILimits);
1334 :	Isibaar	262
1335 :
1336 :	edgomez	1382	for(i=0; i<64; i++) Blk[i] = 0;
1337 :			emms(); idct(Blk); emms();
1338 :			for(m=i=0; i<64; i++) if (Blk[i]!=0) m++;
1339 :			printf( "IDCT(0) == 0 ? %s\n", (m!=0) ? "NOPE!" : "yup." );
1340 :			}
1341 :	Isibaar	262	}
1342 :
1343 :
1344 :			void test_dct_saturation(int Min, int Max)
1345 :			{
1346 :	edgomez	1382	/* test behaviour on input range fringe */
1347 :	Isibaar	262
1348 :	edgomez	1382	int i, n, p;
1349 :			CPU *cpu;
1350 :			// const short IDCT_MAX = 2047; /* 12bits input */
1351 :	edgomez	851	// const short IDCT_MIN = -2048;
1352 :	edgomez	1382	// const short IDCT_OUT = 256; /* 9bits ouput */
1353 :			const int Partitions = 4;
1354 :			const int Loops = 10000 / Partitions;
1355 :	Isibaar	262
1356 :	edgomez	1382	init_ref_dct();
1357 :	Isibaar	262
1358 :	edgomez	1382	for(cpu = cpu_list; cpu->name!=0; ++cpu)
1359 :			{
1360 :			short Blk0[64], Blk[64];
1361 :			STATS_8x8 Stats;
1362 :	Isibaar	262
1363 :	edgomez	1382	if (!init_cpu(cpu))
1364 :			continue;
1365 :	Isibaar	262
1366 :	edgomez	1382	printf( "\n===== IEEE test for %s Min=%d Max=%d =====\n",
1367 :			cpu->name, Min, Max );
1368 :	Isibaar	262
1369 :	edgomez	1382	/* FDCT tests // */
1370 :	Isibaar	262
1371 :	edgomez	1382	init_stats(&Stats);
1372 :	Isibaar	262
1373 :	edgomez	1382	/* test each computation channels separately */
1374 :			for(i=0; i<64; i++) Blk[i] = Blk0[i] = ((i/8)==(i%8)) ? Max : 0;
1375 :			ref_fdct(Blk0);
1376 :			emms(); fdct(Blk); emms();
1377 :			store_stats(&Stats, Blk, Blk0);
1378 :	Isibaar	262
1379 :	edgomez	1382	for(i=0; i<64; i++) Blk[i] = Blk0[i] = ((i/8)==(i%8)) ? Min : 0;
1380 :			ref_fdct(Blk0);
1381 :			emms(); fdct(Blk); emms();
1382 :			store_stats(&Stats, Blk, Blk0);
1383 :	Isibaar	262
1384 :	edgomez	1382	/* randomly saturated inputs */
1385 :			for(p=0; p<Partitions; ++p)
1386 :			{
1387 :			for(n=0; n<Loops; ++n)
1388 :			{
1389 :			for(i=0; i<64; ++i)
1390 :			Blk0[i] = Blk[i] = (ieee_rand(0,Partitions)>=p)? Max : Min;
1391 :			ref_fdct(Blk0);
1392 :			emms(); fdct(Blk); emms();
1393 :			store_stats(&Stats, Blk, Blk0);
1394 :			}
1395 :			}
1396 :			printf( "\n -- FDCT saturation report --\n" );
1397 :			report_stats(&Stats, 0);
1398 :	Isibaar	262
1399 :
1400 :	edgomez	1382	/* IDCT tests // */
1401 :	Isibaar	262	#if 0
1402 :	edgomez	1382	/* no finished yet */
1403 :	Isibaar	262
1404 :	edgomez	1382	init_stats(&Stats);
1405 :	Isibaar	262
1406 :	edgomez	1382	/* test each computation channel separately */
1407 :			for(i=0; i<64; i++) Blk[i] = Blk0[i] = ((i/8)==(i%8)) ? IDCT_MAX : 0;
1408 :			ref_idct(Blk0);
1409 :			emms(); idct(Blk); emms();
1410 :			for(i=0; i<64; i++) { CLAMP(Blk0[i], IDCT_OUT); CLAMP(Blk[i], IDCT_OUT); }
1411 :			store_stats(&Stats, Blk, Blk0);
1412 :	Isibaar	262
1413 :	edgomez	1382	for(i=0; i<64; i++) Blk[i] = Blk0[i] = ((i/8)==(i%8)) ? IDCT_MIN : 0;
1414 :			ref_idct(Blk0);
1415 :			emms(); idct(Blk); emms();
1416 :			for(i=0; i<64; i++) { CLAMP(Blk0[i], IDCT_OUT); CLAMP(Blk[i], IDCT_OUT); }
1417 :			store_stats(&Stats, Blk, Blk0);
1418 :	Isibaar	262
1419 :	edgomez	1382	/* randomly saturated inputs */
1420 :			for(p=0; p<Partitions; ++p)
1421 :			{
1422 :			for(n=0; n<Loops; ++n)
1423 :			{
1424 :			for(i=0; i<64; ++i)
1425 :			Blk0[i] = Blk[i] = (ieee_rand(0,Partitions)>=p)? IDCT_MAX : IDCT_MIN;
1426 :			ref_idct(Blk0);
1427 :			emms(); idct(Blk); emms();
1428 :			for(i=0; i<64; i++) { CLAMP(Blk0[i],IDCT_OUT); CLAMP(Blk[i],IDCT_OUT); }
1429 :			store_stats(&Stats, Blk, Blk0);
1430 :			}
1431 :			}
1432 :	Isibaar	262
1433 :	edgomez	1382	printf( "\n -- IDCT saturation report --\n" );
1434 :			print_stats(&Stats);
1435 :			report_stats(&Stats, 0);
1436 :	Isibaar	262	#endif
1437 :	edgomez	1382	}
1438 :	Isibaar	262	}
1439 :
1440 :			/*********************************************************************
1441 :	Isibaar	225	* measure raw decoding speed
1442 :			*********************************************************************/
1443 :
1444 :	Skal	1614	void test_dec(const char *name, int width, int height, int ref_chksum)
1445 :	Isibaar	225	{
1446 :	edgomez	1382	FILE *f = 0;
1447 :			void *dechandle = 0;
1448 :			int xerr;
1449 :			xvid_gbl_init_t xinit;
1450 :			xvid_dec_create_t xparam;
1451 :			xvid_dec_frame_t xframe;
1452 :	Isibaar	225	double t = 0.;
1453 :			int nb = 0;
1454 :	edgomez	1382	uint8_t *buf = 0;
1455 :	Skal	1614	uint8_t *yuv_out = 0;
1456 :	edgomez	1382	int buf_size, pos;
1457 :			uint32_t chksum = 0;
1458 :	Skal	1614	int bps = (width+31) & ~31;
1459 :	Isibaar	225
1460 :	edgomez	1382	memset(&xinit, 0, sizeof(xinit));
1461 :	Skal	1614	xinit.cpu_flags = cpu_mask;
1462 :	edgomez	1382	xinit.version = XVID_VERSION;
1463 :			xvid_global(NULL, 0, &xinit, NULL);
1464 :	Isibaar	225
1465 :	edgomez	1382	memset(&xparam, 0, sizeof(xparam));
1466 :	Skal	1614	xparam.width = width;
1467 :	Isibaar	225	xparam.height = height;
1468 :	edgomez	1382	xparam.version = XVID_VERSION;
1469 :	Isibaar	225	xerr = xvid_decore(NULL, XVID_DEC_CREATE, &xparam, NULL);
1470 :	edgomez	1382	if (xerr==XVID_ERR_FAIL) {
1471 :	Skal	1614	printf("ERROR: can't init decoder (err=%d)\n", xerr);
1472 :	edgomez	1382	return;
1473 :	Isibaar	225	}
1474 :			dechandle = xparam.handle;
1475 :
1476 :
1477 :			f = fopen(name, "rb");
1478 :	edgomez	1382	if (f==0) {
1479 :	Skal	1614	printf( "ERROR: can't open file '%s'\n", name);
1480 :	edgomez	1382	return;
1481 :			}
1482 :			fseek(f, 0, SEEK_END);
1483 :			buf_size = ftell(f);
1484 :			fseek(f, 0, SEEK_SET);
1485 :			if (buf_size<=0) {
1486 :	Skal	1614	printf("ERROR: error while stating file\n");
1487 :	edgomez	1382	goto End;
1488 :			}
1489 :	Isibaar	225
1490 :	Skal	1614	buf = malloc(buf_size);
1491 :			yuv_out = calloc(1, bpsheight3/2 + 15);
1492 :			if (buf==0 \|\| yuv_out==0) {
1493 :			printf( "ERROR: malloc failed!\n" );
1494 :	edgomez	1382	goto End;
1495 :			}
1496 :	Isibaar	225
1497 :	edgomez	1382	if (fread(buf, buf_size, 1, f)!=1) {
1498 :	Skal	1614	printf( "ERROR: file-read failed\n" );
1499 :	edgomez	1382	goto End;
1500 :			}
1501 :	Isibaar	225
1502 :	edgomez	1382	nb = 0;
1503 :			pos = 0;
1504 :			t = -gettime_usec();
1505 :			while(1) {
1506 :	Skal	1614	int y;
1507 :
1508 :	edgomez	1382	memset(&xframe, 0, sizeof(xframe));
1509 :			xframe.version = XVID_VERSION;
1510 :			xframe.bitstream = buf + pos;
1511 :			xframe.length = buf_size - pos;
1512 :	Skal	1614	xframe.output.plane[0] = (uint8_t*)(((size_t)yuv_out + 15) & ~15);
1513 :	suxen_drol	1638	xframe.output.plane[1] = (uint8_t)xframe.output.plane[0] + bpsheight;
1514 :			xframe.output.plane[2] = (uint8_t*)xframe.output.plane[1] + bps/2;
1515 :	Skal	1614	xframe.output.stride[0] = bps;
1516 :			xframe.output.stride[1] = bps;
1517 :			xframe.output.stride[2] = bps;
1518 :			xframe.output.csp = XVID_CSP_I420;
1519 :	edgomez	1382	xerr = xvid_decore(dechandle, XVID_DEC_DECODE, &xframe, 0);
1520 :	Skal	1614	if (xerr<0) {
1521 :			printf("ERROR: decoding failed for frame #%d (err=%d)!\n", nb, xerr);
1522 :			break;
1523 :			}
1524 :			else if (xerr==0)
1525 :			break;
1526 :			else if (verbose>0) printf("#%d %d\n", nb, xerr );
1527 :
1528 :			pos += xerr;
1529 :	edgomez	1382	nb++;
1530 :	Skal	1614
1531 :			for(y=0; y<height/2; ++y) {
1532 :	suxen_drol	1638	chksum = calc_crc((uint8_t)xframe.output.plane[0] + (2y+0)*bps, width, chksum);
1533 :			chksum = calc_crc((uint8_t)xframe.output.plane[0] + (2y+1)*bps, width, chksum);
1534 :			chksum = calc_crc((uint8_t)xframe.output.plane[1] + ybps, width/2, chksum);
1535 :			chksum = calc_crc((uint8_t)xframe.output.plane[2] + ybps, width/2, chksum);
1536 :	edgomez	1382	}
1537 :			if (pos==buf_size)
1538 :			break;
1539 :			}
1540 :			t += gettime_usec();
1541 :	Skal	1614	if (ref_chksum==0) {
1542 :			if (t>0.)
1543 :			printf( "%d frames decoded in %.3f s -> %.1f FPS Checksum:0x%.8x\n", nb, t1.e-6f, (float)(nb1.e6f/t), chksum );
1544 :			}
1545 :			else {
1546 :			printf("FPS:%.1f Checksum: 0x%.8x Expected:0x%.8x \| %s\n",
1547 :			t>0. ? (float)(nb*1.e6f/t) : 0.f, chksum, ref_chksum, (chksum==ref_chksum) ? "OK" : "ERROR");
1548 :			}
1549 :	Isibaar	225
1550 :	edgomez	1382	End:
1551 :	Skal	1614	if (yuv_out!=0) free(yuv_out);
1552 :	edgomez	1382	if (buf!=0) free(buf);
1553 :			if (dechandle!=0) {
1554 :			xerr= xvid_decore(dechandle, XVID_DEC_DESTROY, NULL, NULL);
1555 :			if (xerr==XVID_ERR_FAIL)
1556 :	Skal	1614	printf("ERROR: destroy-decoder failed (err=%d)!\n", xerr);
1557 :	edgomez	1382	}
1558 :			if (f!=0) fclose(f);
1559 :	Isibaar	225	}
1560 :
1561 :			/*********************************************************************
1562 :			* non-regression tests
1563 :			*********************************************************************/
1564 :
1565 :			void test_bugs1()
1566 :			{
1567 :	edgomez	1382	CPU *cpu;
1568 :			uint16_t mpeg_quant_matrices[64*8];
1569 :	Isibaar	225
1570 :	edgomez	1382	printf( "\n ===== (de)quant4_intra saturation bug? =====\n" );
1571 :	Isibaar	225
1572 :	edgomez	1382	for(cpu = cpu_list; cpu->name!=0; ++cpu)
1573 :			{
1574 :			int i;
1575 :			int16_t Src[88], Dst[88];
1576 :	Isibaar	225
1577 :	edgomez	1382	if (!init_cpu(cpu))
1578 :			continue;
1579 :	Isibaar	225
1580 :	edgomez	1382	for(i=0; i<64; ++i) Src[i] = i-32;
1581 :			set_intra_matrix( mpeg_quant_matrices, get_default_intra_matrix() );
1582 :			dequant_mpeg_intra(Dst, Src, 31, 5, mpeg_quant_matrices);
1583 :			printf( "dequant_mpeg_intra with CPU=%s: ", cpu->name);
1584 :			printf( " Out[]= " );
1585 :			for(i=0; i<64; ++i) printf( "[%d]", Dst[i]);
1586 :			printf( "\n" );
1587 :			}
1588 :	Isibaar	225
1589 :	edgomez	1382	printf( "\n ===== (de)quant4_inter saturation bug? =====\n" );
1590 :	Isibaar	225
1591 :	edgomez	1382	for(cpu = cpu_list; cpu->name!=0; ++cpu)
1592 :			{
1593 :			int i;
1594 :			int16_t Src[88], Dst[88];
1595 :	Isibaar	225
1596 :	edgomez	1382	if (!init_cpu(cpu))
1597 :			continue;
1598 :	Isibaar	225
1599 :	edgomez	1382	for(i=0; i<64; ++i) Src[i] = i-32;
1600 :			set_inter_matrix( mpeg_quant_matrices, get_default_inter_matrix() );
1601 :			dequant_mpeg_inter(Dst, Src, 31, mpeg_quant_matrices);
1602 :			printf( "dequant_mpeg_inter with CPU=%s: ", cpu->name);
1603 :			printf( " Out[]= " );
1604 :			for(i=0; i<64; ++i) printf( "[%d]", Dst[i]);
1605 :			printf( "\n" );
1606 :			}
1607 :	Isibaar	225	}
1608 :
1609 :			void test_dct_precision_diffs()
1610 :			{
1611 :	edgomez	1382	CPU *cpu;
1612 :			DECLARE_ALIGNED_MATRIX(Blk, 8, 8, int16_t, 16);
1613 :			DECLARE_ALIGNED_MATRIX(Blk0, 8, 8, int16_t, 16);
1614 :	Isibaar	225
1615 :	edgomez	1382	printf( "\n ===== fdct/idct precision diffs =====\n" );
1616 :	Isibaar	225
1617 :	edgomez	1382	for(cpu = cpu_list; cpu->name!=0; ++cpu)
1618 :			{
1619 :			int i;
1620 :	Isibaar	225
1621 :	edgomez	1382	if (!init_cpu(cpu))
1622 :			continue;
1623 :	Isibaar	225
1624 :	edgomez	1382	for(i=0; i<8*8; ++i) {
1625 :			Blk0[i] = (i7-ii) & 0x7f;
1626 :			Blk[i] = Blk0[i];
1627 :			}
1628 :	Isibaar	225
1629 :	edgomez	1382	fdct(Blk);
1630 :			idct(Blk);
1631 :			printf( " fdct+idct diffs with CPU=%s: \n", cpu->name );
1632 :			for(i=0; i<8; ++i) {
1633 :			int j;
1634 :			for(j=0; j<8; ++j) printf( " %d ", Blk[i8+j]-Blk0[i8+j]);
1635 :			printf("\n");
1636 :			}
1637 :			printf("\n");
1638 :			}
1639 :	Isibaar	225	}
1640 :
1641 :	Isibaar	262	void test_quant_bug()
1642 :			{
1643 :	edgomez	1382	const int max_Q = 31;
1644 :			int i, n, qm, q;
1645 :			CPU *cpu;
1646 :			DECLARE_ALIGNED_MATRIX(Src, 8, 8, int16_t, 16);
1647 :			DECLARE_ALIGNED_MATRIX(Dst, 8, 8, int16_t, 16);
1648 :			uint8_t Quant[8*8];
1649 :			CPU cpu_bug_list[] = { { "PLAINC", 0 }, { "MMX ", XVID_CPU_MMX }, {0,0} };
1650 :			uint16_t Crcs_Inter[2][32];
1651 :			uint16_t Crcs_Intra[2][32];
1652 :			DECLARE_ALIGNED_MATRIX(mpeg_quant_matrices, 8, 64, uint16_t, 16);
1653 :	Isibaar	225
1654 :	edgomez	1382	printf( "\n ===== test MPEG4-quantize bug =====\n" );
1655 :	Isibaar	262
1656 :	edgomez	1382	for(i=0; i<64; ++i) Src[i] = 2048*(i-32)/32;
1657 :
1658 :	Isibaar	262	#if 1
1659 :	edgomez	1382	for(qm=1; qm<=255; ++qm)
1660 :			{
1661 :			for(i=0; i<8*8; ++i) Quant[i] = qm;
1662 :			set_inter_matrix( mpeg_quant_matrices, Quant );
1663 :	Isibaar	262
1664 :	edgomez	1382	for(n=0, cpu = cpu_bug_list; cpu->name!=0; ++cpu, ++n)
1665 :			{
1666 :			uint16_t s;
1667 :	Isibaar	262
1668 :	edgomez	1382	if (!init_cpu(cpu))
1669 :			continue;
1670 :	Isibaar	262
1671 :	edgomez	1382	for(q=1; q<=max_Q; ++q) {
1672 :			emms();
1673 :			quant_mpeg_inter( Dst, Src, q, mpeg_quant_matrices );
1674 :			emms();
1675 :			for(s=0, i=0; i<64; ++i) s+=((uint16_t)Dst[i])^i;
1676 :			Crcs_Inter[n][q] = s;
1677 :			}
1678 :			}
1679 :	Isibaar	262
1680 :	edgomez	1382	for(q=1; q<=max_Q; ++q)
1681 :			for(i=0; i<n-1; ++i)
1682 :			if (Crcs_Inter[i][q]!=Crcs_Inter[i+1][q])
1683 :			printf( "Discrepancy Inter: qm=%d, q=%d -> %d/%d !\n",
1684 :			qm, q, Crcs_Inter[i][q], Crcs_Inter[i+1][q]);
1685 :			}
1686 :	Isibaar	262	#endif
1687 :
1688 :			#if 1
1689 :	edgomez	1382	for(qm=1; qm<=255; ++qm)
1690 :			{
1691 :			for(i=0; i<8*8; ++i) Quant[i] = qm;
1692 :			set_intra_matrix( mpeg_quant_matrices, Quant );
1693 :	Isibaar	262
1694 :	edgomez	1382	for(n=0, cpu = cpu_bug_list; cpu->name!=0; ++cpu, ++n)
1695 :			{
1696 :			uint16_t s;
1697 :	Isibaar	262
1698 :	edgomez	1382	if (!init_cpu(cpu))
1699 :			continue;
1700 :	Isibaar	262
1701 :	edgomez	1382	for(q=1; q<=max_Q; ++q) {
1702 :			emms();
1703 :			quant_mpeg_intra( Dst, Src, q, q, mpeg_quant_matrices);
1704 :			emms();
1705 :			for(s=0, i=0; i<64; ++i) s+=((uint16_t)Dst[i])^i;
1706 :			Crcs_Intra[n][q] = s;
1707 :			}
1708 :			}
1709 :	Isibaar	262
1710 :	edgomez	1382	for(q=1; q<=max_Q; ++q)
1711 :			for(i=0; i<n-1; ++i)
1712 :			if (Crcs_Intra[i][q]!=Crcs_Intra[i+1][q])
1713 :			printf( "Discrepancy Intra: qm=%d, q=%d -> %d/%d!\n",
1714 :			qm, q, Crcs_Inter[i][q], Crcs_Inter[i+1][q]);
1715 :			}
1716 :	Isibaar	262	#endif
1717 :			}
1718 :	Skal	1707
1719 :			/*********************************************************************
1720 :			* test some YUV func
1721 :			*********************************************************************/
1722 :
1723 :			#define ENTER \
1724 :			for(i=0; i<(int)sizeof(Dst0); ++i) Dst0[0][i] = 0; \
1725 :			t = gettime_usec(); \
1726 :			emms();
1727 :
1728 :			#define LEAVE \
1729 :			emms(); \
1730 :			t = (gettime_usec() - t) / nb_tests; \
1731 :			iCrc = calc_crc((uint8_t*)Dst0, sizeof(Dst0), CRC32_INITIAL)
1732 :
1733 :	Isibaar	1801	#define TEST_YUYV(FUNC, S, FLIP) \
1734 :	Skal	1707	ENTER \
1735 :	Isibaar	1801	for(tst=0; tst<nb_tests; ++tst) (FUNC)(Dst0[0], S*WIDTH, Src0[0], Src0[1], Src0[2], WIDTH, WIDTH/2, WIDTH, HEIGHT, (FLIP)); \
1736 :	Skal	1707	LEAVE
1737 :
1738 :	Isibaar	1801	static const int yuv_CRCs[6][2] = {
1739 :			{0x0f4fb96b,0x780b6a68}
1740 :			, {0xa986b289,0x65e49b76}
1741 :			, {0x7f19c152,0xd539b86e}
1742 :			, {0x0f4fb96b,0x780b6a68}
1743 :			, {0xa986b289,0x65e49b76}
1744 :			, {0x36ab8b57,0x1cd92fee}
1745 :			};
1746 :
1747 :	Skal	1742	#define WIDTH 128
1748 :			#define HEIGHT 32
1749 :	Skal	1707	void test_yuv()
1750 :			{
1751 :			const int nb_tests = 200*speed_ref;
1752 :			CPU *cpu;
1753 :			uint8_t Src0[3][WIDTH*HEIGHT];
1754 :			uint8_t Dst0[4][WIDTH*HEIGHT];
1755 :	Isibaar	1801	int i, j, with_flip;
1756 :	Skal	1707	double t;
1757 :			int tst, iCrc;
1758 :
1759 :			colorspace_init();
1760 :			ieee_reseed(1);
1761 :			for(i=0; i<(int)sizeof(Src0); ++i) Src0[0][i] = ieee_rand(0,255);
1762 :	Skal	1742	for(i=0; i<(int)sizeof(Dst0); ++i) Dst0[0][i] = 0x5a;
1763 :	Isibaar	1801
1764 :			printf( "\n === test YUV ===\n" );
1765 :	Skal	1707
1766 :	Isibaar	1801	for(with_flip=0; with_flip<=1; ++with_flip) {
1767 :	Skal	1707
1768 :	Isibaar	1801	init_cpu(&cpu_list[0]);
1769 :			TEST_YUYV(yv12_to_yuyv_c, 4, with_flip);
1770 :			printf(" yv12_to_yuyv_c %.3f usec crc32=0x%08x %s\n",
1771 :			t, iCrc, (iCrc!=yuv_CRCs[0][with_flip])?"\| ERROR": "" );
1772 :			TEST_YUYV(yv12_to_uyvy_c, 4, with_flip);
1773 :			printf(" yv12_to_uyvy_c %.3f usec crc32=0x%08x %s\n",
1774 :			t, iCrc, (iCrc!=yuv_CRCs[1][with_flip])?"\| ERROR": "" );
1775 :
1776 :			TEST_YUYV(yv12_to_bgra_c, 4, with_flip);
1777 :			printf(" yv12_to_bgra_c %.3f usec crc32=0x%08x %s\n",
1778 :			t, iCrc, (iCrc!=yuv_CRCs[2][with_flip])?"\| ERROR": "" );
1779 :	Skal	1707
1780 :	Isibaar	1801	#if defined(ARCH_IS_IA32) \|\| defined(ARCH_IS_X86_64)
1781 :			init_cpu(&cpu_list[1]);
1782 :			TEST_YUYV(yv12_to_yuyv_mmx, 4, with_flip);
1783 :			printf(" yv12_to_yuyv_mmx %.3f usec crc32=0x%08x %s\n",
1784 :			t, iCrc, (iCrc!=yuv_CRCs[3][with_flip])?"\| ERROR": "" );
1785 :	Skal	1707
1786 :	Isibaar	1801	TEST_YUYV(yv12_to_uyvy_mmx, 4, with_flip);
1787 :			printf(" yv12_to_uyvy_mmx %.3f usec crc32=0x%08x %s\n",
1788 :			t, iCrc, (iCrc!=yuv_CRCs[4][with_flip])?"\| ERROR": "" );
1789 :
1790 :			TEST_YUYV(yv12_to_bgra_mmx, 4, with_flip);
1791 :			printf(" yv12_to_bgra_mmx %.3f usec crc32=0x%08x %s\n",
1792 :			t, iCrc, (iCrc!=yuv_CRCs[5][with_flip])?"\| ERROR": "" );
1793 :
1794 :	Skal	1707	#endif
1795 :
1796 :			#ifdef ARCH_IS_PPC
1797 :	Isibaar	1801	init_cpu(&cpu_list[1]);
1798 :			TEST_YUYV(yv12_to_yuyv_altivec_c, 4, with_flip);
1799 :			printf(" yv12_to_yuyv_altivec_c %.3f usec crc32=0x%08x %s\n",
1800 :			t, iCrc, (iCrc!=yuv_CRCs[3][with_flip])?"\| ERROR": "" );
1801 :	Skal	1707
1802 :	Isibaar	1801	TEST_YUYV(yv12_to_uyvy_altivec_c, 4, with_flip);
1803 :			printf(" yv12_to_uyvy_altivec_c %.3f usec crc32=0x%08x %s\n",
1804 :			t, iCrc, (iCrc!=yuv_CRCs[4][with_flip])?"\| ERROR": "" );
1805 :
1806 :			TEST_YUYV(yv12_to_bgra_altivec_c, 4, with_flip);
1807 :			printf(" yv12_to_bgra_altivec_c %.3f usec crc32=0x%08x %s\n",
1808 :			t, iCrc, (iCrc!=yuv_CRCs[5][with_flip])?"\| ERROR": "" );
1809 :	Skal	1707	#endif
1810 :	Isibaar	1801	}
1811 :	Skal	1707	printf( " --- \n" );
1812 :			}
1813 :	Skal	1742
1814 :			#define TEST_YV2(FUNC, WITH_UV, WITH_FLIP) \
1815 :			ENTER \
1816 :			for(tst=0; tst<nb_tests; ++tst) (FUNC)(Dst0[0], Dst0[1], Dst0[2], WIDTH, WIDTH, \
1817 :			Src0[0], (WITH_UV) ? Src0[1] : 0, (WITH_UV) ? Src0[2] : 0, WIDTH, WIDTH, \
1818 :			WIDTH-2, HEIGHT-2, WITH_FLIP); \
1819 :			LEAVE
1820 :
1821 :			#define PRINT_NxN(DATA,W,H,STR) { \
1822 :			int i,j; \
1823 :			for(j=0; j<(H); ++j) { \
1824 :			for(i=0; i<(W); ++i) printf( "0x%.2x ", (DATA)[i+j*(STR)] );\
1825 :			printf("\n"); \
1826 :			} \
1827 :			printf("---\n"); \
1828 :			}
1829 :
1830 :			static const int yv12_CRCs[2][2] = {
1831 :			{0x5cab7cf0,0xdab46541}
1832 :			, {0xe8bae865,0x1faf77b7}
1833 :			};
1834 :
1835 :			void test_yuv2()
1836 :			{
1837 :			const int nb_tests = 800*speed_ref;
1838 :			CPU *cpu;
1839 :			uint8_t Src0[3][WIDTH*HEIGHT];
1840 :			uint8_t Dst0[3][WIDTH*HEIGHT];
1841 :			int with_uv, with_flip;
1842 :			int i, j;
1843 :			double t;
1844 :			int tst, iCrc;
1845 :
1846 :			colorspace_init();
1847 :			ieee_reseed(1);
1848 :			for(i=0; i<(int)sizeof(Src0); ++i) Src0[0][i] = ieee_rand(0,255);
1849 :
1850 :			printf( "\n === test YV2 ===\n" );
1851 :			for(with_flip=0; with_flip<=1; ++with_flip) {
1852 :			for(with_uv=0; with_uv<=1; ++with_uv) {
1853 :			init_cpu(&cpu_list[0]);
1854 :			TEST_YV2(yv12_to_yv12_c, with_uv, with_flip);
1855 :			printf(" yv12_to_yv12_c %.3f usec \tcrc32=0x%08x %s\n",
1856 :			t, iCrc, (iCrc!=yv12_CRCs[with_flip][with_uv])?"\| ERROR": "" );
1857 :			/* if (!with_uv) PRINT_NxN(Dst0[1], WIDTH/2, HEIGHT/2, WIDTH ); */
1858 :
1859 :	Isibaar	1801	#if defined(ARCH_IS_IA32) \|\| defined(ARCH_IS_X86_64)
1860 :	Skal	1742	init_cpu(&cpu_list[1]);
1861 :			TEST_YV2(yv12_to_yv12_mmx, with_uv, with_flip);
1862 :			printf(" yv12_to_yv12_mmx %.3f usec \tcrc32=0x%08x %s\n",
1863 :			t, iCrc, (iCrc!=yv12_CRCs[with_flip][with_uv])?"\| ERROR": "" );
1864 :			/* if (!with_uv) PRINT_NxN(Dst0[1], WIDTH/2, HEIGHT/2, WIDTH ); */
1865 :
1866 :			TEST_YV2(yv12_to_yv12_xmm, with_uv, with_flip);
1867 :			printf(" yv12_to_yv12_xmm %.3f usec \tcrc32=0x%08x %s\n",
1868 :			t, iCrc, (iCrc!=yv12_CRCs[with_flip][with_uv])?"\| ERROR": "" );
1869 :			#endif
1870 :			}
1871 :
1872 :			printf( " --- \n" );
1873 :			}
1874 :			printf( " ===== \n" );
1875 :			}
1876 :
1877 :	Skal	1707	#undef WIDTH
1878 :			#undef HEIGHT
1879 :			#undef ENTER
1880 :			#undef LEAVE
1881 :
1882 :	Skal	1616	/*********************************************************************/
1883 :	Isibaar	262
1884 :	Skal	1616	static uint32_t __inline log2bin_v1(uint32_t value)
1885 :			{
1886 :			int n = 0;
1887 :			while (value) {
1888 :			value >>= 1;
1889 :			n++;
1890 :			}
1891 :			return n;
1892 :			}
1893 :
1894 :	Skal	1617	static const uint8_t log2_tab_16[16] = { 0, 1, 2, 2, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4, 4 };
1895 :	Skal	1616
1896 :			static uint32_t __inline log2bin_v2(uint32_t value)
1897 :			{
1898 :			int n = 0;
1899 :			if (value & 0xffff0000) {
1900 :			value >>= 16;
1901 :			n += 16;
1902 :			}
1903 :			if (value & 0xff00) {
1904 :			value >>= 8;
1905 :			n += 8;
1906 :			}
1907 :			if (value & 0xf0) {
1908 :			value >>= 4;
1909 :			n += 4;
1910 :			}
1911 :			return n + log2_tab_16[value];
1912 :			}
1913 :
1914 :			void test_log2bin()
1915 :			{
1916 :			const int nb_tests = 3000*speed_ref;
1917 :			int n, crc1=0, crc2=0;
1918 :			uint32_t s, s0;
1919 :			double t1, t2;
1920 :
1921 :			t1 = gettime_usec();
1922 :			s0 = (int)(t1*31.241);
1923 :			for(s=s0, n=0; n<nb_tests; ++n, s=(s*12363+31)&0x7fffffff)
1924 :			crc1 += log2bin_v1(s);
1925 :			t1 = (gettime_usec()-t1) / nb_tests;
1926 :
1927 :			t2 = gettime_usec();
1928 :			for(s=s0, n=0; n<nb_tests; ++n, s=(s*12363+31)&0x7fffffff)
1929 :			crc2 += log2bin_v2(s);
1930 :			t2 = (gettime_usec() - t2) / nb_tests;
1931 :	Skal	1617
1932 :	Skal	1616	printf( "log2bin_v1: %.3f sec crc=%d\n", t1, crc1 );
1933 :			printf( "log2bin_v2: %.3f sec crc=%d\n", t2, crc2 );
1934 :			if (crc1!=crc2) printf( " CRC ERROR !\n" );
1935 :			}
1936 :
1937 :	Skal	1617	/*********************************************************************/
1938 :
1939 :			static void __inline old_gcd(int num, int den)
1940 :			{
1941 :			int i = *num;
1942 :			while (i > 1) {
1943 :			if (num % i == 0 && den % i == 0) {
1944 :			*num /= i;
1945 :			*den /= i;
1946 :			i = *num;
1947 :			continue;
1948 :			}
1949 :			i--;
1950 :			}
1951 :			}
1952 :
1953 :			static uint32_t gcd(int num, int den)
1954 :			{
1955 :			int tmp;
1956 :			while( (tmp=num%den) ) { num = den; den = tmp; }
1957 :			return den;
1958 :			}
1959 :			static void __inline new_gcd(int num, int den)
1960 :			{
1961 :			const int div = gcd(num, den);
1962 :			if (num) {
1963 :			*num /= div;
1964 :			*den /= div;
1965 :			}
1966 :			}
1967 :
1968 :			void test_gcd()
1969 :			{
1970 :	Skal	1746	const int nb_tests = 10*speed_ref;
1971 :	Skal	1617	int i;
1972 :			uint32_t crc1=0, crc2=0;
1973 :			uint32_t n0, n, d0, d;
1974 :			double t1, t2;
1975 :
1976 :			t1 = gettime_usec();
1977 :			n0 = 0xfffff & (int)(t1*31.241);
1978 :			d0 = 0xfffff & (int)( ((n0*4123)%17) \| 1 );
1979 :			for(n=n0, d=d0, i=0; i<nb_tests; ++i) {
1980 :			old_gcd(&n, &d);
1981 :			crc1 = (((crc1>>4)^d) + ((crc1<<2)^n) ) & 0xffffff;
1982 :			n = d;
1983 :			d = (d*12363+31) & 0xffff;
1984 :			d \|= !d;
1985 :			}
1986 :			t1 = (gettime_usec()-t1) / nb_tests;
1987 :
1988 :			t2 = gettime_usec();
1989 :			for(n=n0, d=d0, i=0; i<nb_tests; ++i) {
1990 :			new_gcd(&n, &d);
1991 :			crc2 = (((crc2>>4)^d) + ((crc2<<2)^n) ) & 0xffffff;
1992 :			n = d;
1993 :			d = (d*12363+31) & 0xffff;
1994 :			d \|= !d;
1995 :			}
1996 :			t2 = (gettime_usec() - t2) / nb_tests;
1997 :
1998 :			printf( "old_gcd: %.3f sec crc=%d\n", t1, crc1 );
1999 :			printf( "new_gcd: %.3f sec crc=%d\n", t2, crc2 );
2000 :			if (crc1!=crc2) printf( " CRC ERROR !\n" );
2001 :			}
2002 :
2003 :	Isibaar	225	/*********************************************************************
2004 :	Skal	1651	* test compiler
2005 :			*********************************************************************/
2006 :
2007 :			void test_compiler() {
2008 :			int nb_err = 0;
2009 :			int32_t v;
2010 :			if (sizeof(uint16_t)<2) {
2011 :	Skal	1652	printf( "ERROR: sizeof(uint16_t)<2 !!\n" );
2012 :	Skal	1651	nb_err++;
2013 :			}
2014 :			if (sizeof(int16_t)<2) {
2015 :	Skal	1652	printf( "ERROR: sizeof(int16_t)<2 !!\n" );
2016 :	Skal	1651	nb_err++;
2017 :			}
2018 :			if (sizeof(uint8_t)!=1) {
2019 :	Skal	1652	printf( "ERROR: sizeof(uint8_t)!=1 !!\n" );
2020 :	Skal	1651	nb_err++;
2021 :			}
2022 :			if (sizeof(int8_t)!=1) {
2023 :	Skal	1652	printf( "ERROR: sizeof(int8_t)!=1 !!\n" );
2024 :	Skal	1651	nb_err++;
2025 :			}
2026 :			if (sizeof(uint32_t)<4) {
2027 :	Skal	1652	printf( "ERROR: sizeof(uint32_t)<4 !!\n" );
2028 :	Skal	1651	nb_err++;
2029 :			}
2030 :			if (sizeof(int32_t)<4) {
2031 :	Skal	1652	printf( "ERROR: sizeof(int32_t)<4 !!\n" );
2032 :	Skal	1651	nb_err++;
2033 :			}
2034 :			/* yes, i know, this test is silly. But better be safe than sorry. :) */
2035 :			for(v=1000; v>=0; v--) {
2036 :			if ( (v>>2) != v/4)
2037 :			nb_err++;
2038 :			}
2039 :			for(v=-1000; v!=-1; v++) {
2040 :			if ( (v>>2) != (v/4)-!!(v%4))
2041 :			nb_err++;
2042 :			}
2043 :			if (nb_err!=0) {
2044 :			printf( "ERROR! please post your platform/compiler specs to xvid-devel@xvid.org !\n" );
2045 :			}
2046 :			}
2047 :	Skal	1746
2048 :	Skal	1731	/*********************************************************************
2049 :			* test SSIM functions
2050 :			*********************************************************************/
2051 :	Skal	1651
2052 :	Skal	1731	typedef int (lumfunc)(uint8_t ptr, int stride);
2053 :			typedef void (csfunc)(uint8_t ptro, uint8_t* ptrc, int stride, int lumo, int lumc, int* pdevo, int* pdevc, int* pcorr);
2054 :
2055 :			extern int lum_8x8_c(uint8_t* ptr, int stride);
2056 :			extern int lum_8x8_mmx(uint8_t* ptr, int stride);
2057 :			extern int lum_2x8_c(uint8_t* ptr, int stride);
2058 :	Skal	1743	extern void consim_c(uint8_t* ptro, uint8_t* ptrc, int stride, int lumo, int lumc, int* pdevo, int* pdevc, int* pcorr);
2059 :	Skal	1731	extern void consim_mmx(uint8_t* ptro, uint8_t* ptrc, int stride, int lumo, int lumc, int* pdevo, int* pdevc, int* pcorr);
2060 :			extern void consim_sse2(uint8_t* ptro, uint8_t* ptrc, int stride, int lumo, int lumc, int* pdevo, int* pdevc, int* pcorr);
2061 :
2062 :			void test_SSIM()
2063 :			{
2064 :			const int nb_tests = 3000*speed_ref;
2065 :			int tst;
2066 :			CPU *cpu;
2067 :			int i;
2068 :			int devs[3];
2069 :			long lumo, lumc;
2070 :			DECLARE_ALIGNED_MATRIX(Ref1, 16, 16, uint8_t, 16);
2071 :			DECLARE_ALIGNED_MATRIX(Ref2, 16, 16, uint8_t, 16);
2072 :			lumfunc lum8x8;
2073 :			lumfunc lum2x8;
2074 :			csfunc csim;
2075 :
2076 :			ieee_reseed(1);
2077 :			printf( "\n ====== test SSIM ======\n" );
2078 :			for(i=0; i<16*16;++i) {
2079 :			long v1, v2;
2080 :			v1 = ieee_rand(-256, 511);
2081 :			v2 = ieee_rand(-256, 511);
2082 :			Ref1[i] = (v1<0) ? 0 : (v1>255) ? 255 : v1;
2083 :			Ref2[i] = (v2<0) ? 0 : (v2>255) ? 255 : v2;
2084 :			}
2085 :			lumc = ieee_rand(0, 255);
2086 :			lumo = ieee_rand(0, 255);
2087 :
2088 :			for(cpu = cpu_list; cpu->name!=0; ++cpu)
2089 :			{
2090 :			double t;
2091 :			int m;
2092 :			if (!init_cpu(cpu))
2093 :			continue;
2094 :			lum8x8 = lum_8x8_c;
2095 :			lum2x8 = lum_2x8_c;
2096 :	Skal	1743	csim = consim_c;
2097 :	Isibaar	1801	#if defined(ARCH_IS_IA32) \|\| defined(ARCH_IS_X86_64)
2098 :	Skal	1731	if (cpu->cpu & XVID_CPU_MMX){
2099 :			lum8x8 = lum_8x8_mmx;
2100 :			csim = consim_mmx;
2101 :			}
2102 :			if (cpu->cpu & XVID_CPU_MMX){
2103 :			csim = consim_sse2;
2104 :			}
2105 :	Isibaar	1767	#endif
2106 :	Skal	1731	t = gettime_usec();
2107 :			emms();
2108 :			for(tst=0; tst<nb_tests; ++tst) m = lum8x8(Ref1, 16);
2109 :			emms();
2110 :			t = (gettime_usec() - t) / nb_tests;
2111 :			printf("%s - ssim-lum8x8 %.3f usec m=%d %s\n",
2112 :			cpu->name, t, m,
2113 :			(m!=8230)?"\| ERROR": "" );
2114 :
2115 :			t = gettime_usec();
2116 :			emms();
2117 :	Skal	1738	for(tst=0; tst<nb_tests; ++tst) m = lum2x8(Ref1+8, 16);
2118 :	Skal	1731	emms();
2119 :			t = (gettime_usec() - t) / nb_tests;
2120 :			printf("%s - ssim-lum2x8 %.3f usec m=%d %s\n",
2121 :			cpu->name, t, m,
2122 :	Skal	1738	(m!=681)?"\| ERROR": "" );
2123 :	Skal	1731
2124 :			t = gettime_usec();
2125 :			emms();
2126 :			for(tst=0; tst<nb_tests; ++tst) csim(Ref1, Ref2, 16, lumo, lumc, devs+0, devs+1, devs+2);
2127 :			emms();
2128 :			t = (gettime_usec() - t) / nb_tests;
2129 :	Skal	1737	printf("%s - ssim-consim %.3f usec devs=[0x%x 0x%x 0x%x] %s\n",
2130 :	Skal	1731	cpu->name, t, devs[0], devs[1], devs[2],
2131 :	Skal	1743	(devs[0]!=0x1bdf0f \|\| devs[1]!=0x137258 \|\| devs[2]!=0xcdb13)?"\| ERROR": "" );
2132 :	Skal	1731	printf( " --- \n" );
2133 :			}
2134 :			}
2135 :
2136 :	Skal	1651	/*********************************************************************
2137 :	Skal	1746	* test bitstream functions
2138 :			*********************************************************************/
2139 :
2140 :			#define BIT_BUF_SIZE 2000
2141 :
2142 :			static void test_bits()
2143 :			{
2144 :			const int nb_tests = 50*speed_ref;
2145 :			int tst;
2146 :			uint32_t Crc;
2147 :			uint8_t Buf[BIT_BUF_SIZE];
2148 :			uint32_t Extracted[BIT_BUF_SIZE8]; / worst case: bits read 1 by 1 */
2149 :			int Lens[BIT_BUF_SIZE*8];
2150 :			double t1;
2151 :
2152 :
2153 :			printf( "\n === test bitstream ===\n" );
2154 :			ieee_reseed(1);
2155 :			Crc = 0;
2156 :
2157 :			t1 = gettime_usec();
2158 :			for(tst=0; tst<nb_tests; ++tst) {
2159 :			Bitstream bs;
2160 :			int m, m2, l, l2;
2161 :
2162 :			for(l=0; l<BIT_BUF_SIZE; ++l)
2163 :			Buf[l] = (uint8_t)ieee_rand(0,255);
2164 :
2165 :			l = BIT_BUF_SIZE - ieee_rand(1,BIT_BUF_SIZE/10);
2166 :			BitstreamInit(&bs, (void*)(Buf+BIT_BUF_SIZE-l), l);
2167 :
2168 :
2169 :			BitstreamReset(&bs);
2170 :			for(l2=l*8, m=0; l2>0; m++) {
2171 :			const int b = ieee_rand(1,32);
2172 :			Lens[m] = b;
2173 :			l2 -= b;
2174 :			if (l2<0) break;
2175 :			Extracted[m] = BitstreamShowBits(&bs, b);
2176 :			BitstreamSkip(&bs, b);
2177 :			// printf( "<= %d: %d 0x%x\n", m, b, Extracted[m]);
2178 :			}
2179 :
2180 :			BitstreamReset(&bs);
2181 :			for(m2=0; m2<m; ++m2) {
2182 :			const int b = Lens[m2];
2183 :			const uint32_t v = BitstreamGetBits(&bs, b);
2184 :			Crc \|= (v!=Extracted[m2]);
2185 :			// printf( "=> %d: %d 0x%x %c\n", m2, b, v, " *"[Crc]);
2186 :			}
2187 :			}
2188 :			t1 = (gettime_usec() - t1) / nb_tests;
2189 :			printf(" test_bits %.3f usec %s\n", t1, (Crc!=0)?"\| ERROR": "" );
2190 :			}
2191 :
2192 :			/*********************************************************************
2193 :	Isibaar	225	* main
2194 :			*********************************************************************/
2195 :
2196 :	Skal	1614	static void arg_missing(const char *opt)
2197 :	Isibaar	225	{
2198 :	Skal	1614	printf( "missing argument after option '%s'\n", opt);
2199 :			exit(-1);
2200 :			}
2201 :
2202 :			int main(int argc, const char *argv[])
2203 :			{
2204 :			int c, what = 0;
2205 :			int width, height;
2206 :			uint32_t chksum = 0;
2207 :	Skal	1746	const char * test_bitstream = 0;
2208 :	Skal	1614
2209 :			cpu_mask = 0; // default => will use autodectect
2210 :			for(c=1; c<argc; ++c)
2211 :			{
2212 :			if (!strcmp(argv[c], "-v")) verbose++;
2213 :			else if (!strcmp(argv[c], "-c")) cpu_mask = 0 /* PLAIN_C */ \| XVID_CPU_FORCE;
2214 :			else if (!strcmp(argv[c], "-mmx")) cpu_mask = XVID_CPU_MMX \| XVID_CPU_FORCE;
2215 :			else if (!strcmp(argv[c], "-mmxext")) cpu_mask = XVID_CPU_MMXEXT \| XVID_CPU_MMX \| XVID_CPU_FORCE;
2216 :	Skal	1618	else if (!strcmp(argv[c], "-sse2")) cpu_mask = XVID_CPU_SSE2 \| XVID_CPU_MMXEXT \| XVID_CPU_MMX \| XVID_CPU_FORCE;
2217 :	Isibaar	1794	else if (!strcmp(argv[c], "-sse3")) cpu_mask = XVID_CPU_SSE3 \| XVID_CPU_SSE2 \| XVID_CPU_MMXEXT \| XVID_CPU_MMX \| XVID_CPU_FORCE;
2218 :			else if (!strcmp(argv[c], "-sse4")) cpu_mask = XVID_CPU_SSE41 \| XVID_CPU_SSE3 \| XVID_CPU_SSE2 \| XVID_CPU_MMXEXT \| XVID_CPU_MMX \| XVID_CPU_FORCE;
2219 :			else if (!strcmp(argv[c], "-3dnow")) cpu_mask = XVID_CPU_3DNOW \| XVID_CPU_FORCE;
2220 :	Skal	1614	else if (!strcmp(argv[c], "-3dnowe")) cpu_mask = XVID_CPU_3DNOW \| XVID_CPU_3DNOWEXT \| XVID_CPU_FORCE;
2221 :			else if (!strcmp(argv[c], "-altivec")) cpu_mask = XVID_CPU_ALTIVEC \| XVID_CPU_FORCE;
2222 :			else if (!strcmp(argv[c], "-spd")) {
2223 :			if (++c==argc) arg_missing( argv[argc-1] );
2224 :			speed_ref = atoi(argv[c]);
2225 :			}
2226 :			else if (argv[c][0]!='-') {
2227 :			what = atoi(argv[c]);
2228 :			if (what==9) {
2229 :			if (c+4>argc) {
2230 :			printf("usage: %s %d bitstream width height (checksum)\n", argv[0], what);
2231 :			exit(-1);
2232 :			}
2233 :			test_bitstream = argv[++c];
2234 :			width = atoi(argv[++c]);
2235 :			height = atoi(argv[++c]);
2236 :			if (c+1<argc && argv[c+1][0]!='-') {
2237 :			if (sscanf(argv[c+1], "0x%x", &chksum)!=1) {
2238 :			printf( "can't read checksum value.\n" );
2239 :			exit(-1);
2240 :			}
2241 :			else c++;
2242 :			}
2243 :			// printf( "[%s] %dx%d (0x%.8x)\n", test_bitstream, width, height, chksum);
2244 :			}
2245 :			}
2246 :			else {
2247 :			printf( "unrecognized option '%s'\n", argv[c]);
2248 :			exit(-1);
2249 :			}
2250 :			}
2251 :
2252 :
2253 :	edgomez	1382	if (what==0 \|\| what==1) test_dct();
2254 :			if (what==0 \|\| what==2) test_mb();
2255 :			if (what==0 \|\| what==3) test_sad();
2256 :			if (what==0 \|\| what==4) test_transfer();
2257 :			if (what==0 \|\| what==5) test_quant();
2258 :			if (what==0 \|\| what==6) test_cbp();
2259 :	edgomez	1424	if (what==0 \|\| what==10) test_sse();
2260 :	Skal	1616	if (what==0 \|\| what==11) test_log2bin();
2261 :	Skal	1617	if (what==0 \|\| what==12) test_gcd();
2262 :	Skal	1651	if (what==0 \|\| what==13) test_compiler();
2263 :	Skal	1707	if (what==0 \|\| what==14) test_yuv();
2264 :	Skal	1731	if (what==0 \|\| what==15) test_SSIM();
2265 :	Skal	1742	if (what==0 \|\| what==16) test_yuv2();
2266 :	Skal	1746	if (what==0 \|\| what==17) test_bits();
2267 :	Isibaar	225
2268 :	edgomez	1382	if (what==7) {
2269 :			test_IEEE1180_compliance(-256, 255, 1);
2270 :			test_IEEE1180_compliance(-256, 255,-1);
2271 :			test_IEEE1180_compliance( -5, 5, 1);
2272 :			test_IEEE1180_compliance( -5, 5,-1);
2273 :			test_IEEE1180_compliance(-300, 300, 1);
2274 :			test_IEEE1180_compliance(-300, 300,-1);
2275 :			}
2276 :			if (what==8) test_dct_saturation(-256, 255);
2277 :	Isibaar	262
2278 :	Skal	1614	if (test_bitstream)
2279 :			test_dec(test_bitstream, width, height, chksum);
2280 :	edgomez	1382	if (what==-1) {
2281 :			test_dct_precision_diffs();
2282 :			test_bugs1();
2283 :			}
2284 :			if (what==-2)
2285 :			test_quant_bug();
2286 :	Isibaar	262
2287 :	edgomez	1424	if ((what >= 0 && what <= 6) \|\| what == 10) {
2288 :	edgomez	1382	printf("\n\n"
2289 :			"NB: If a function isn't optimised for a specific set of intructions,\n"
2290 :			" a C function is used instead. So don't panic if some functions\n"
2291 :			" may appear to be slow.\n");
2292 :			}
2293 :
2294 :	Isibaar	1801	#if defined(ARCH_IS_IA32) \|\| defined(ARCH_IS_X86_64)
2295 :	edgomez	1382	if (what == 0 \|\| what == 5) {
2296 :			printf("\n"
2297 :			"NB: MMX mpeg4 quantization is known to have very small errors (+/-1 magnitude)\n"
2298 :			" for 1 or 2 coefficients a block. This is mainly caused by the fact the unit\n"
2299 :			" test goes far behind the usual limits of real encoding. Please do not report\n"
2300 :			" this error to the developers.\n");
2301 :			}
2302 :			#endif
2303 :
2304 :			return 0;
2305 :	Isibaar	225	}
2306 :
2307 :	Skal	1617	/*********************************************************************/

No admin address has been configured	ViewVC Help
Powered by ViewVC 1.0.4