[svn] / branches / dev-api-4 / xvidcore / examples / xvid_bench.c Repository:
ViewVC logotype

Annotation of /branches/dev-api-4/xvidcore/examples/xvid_bench.c

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1197 - (view) (download)

1 : edgomez 1054 /*****************************************************************************
2 : Isibaar 225 *
3 : edgomez 1054 * XVID MPEG-4 VIDEO CODEC
4 :     * - Unit tests and benches -
5 : Isibaar 225 *
6 : edgomez 1054 * Copyright(C) 2002 Pascal Massimino <skal@planet-d.net>
7 : Isibaar 225 *
8 : edgomez 1054 * 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 1054 * 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 1054 * 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 : edgomez 1197 * $Id: xvid_bench.c,v 1.9.2.6 2003-11-02 23:02:52 edgomez Exp $
23 : edgomez 1054 *
24 :     ****************************************************************************/
25 : Isibaar 225
26 : edgomez 1054 /*****************************************************************************
27 : Isibaar 225 *
28 :     * 'Reference' output is at the end of file.
29 :     * Don't take the checksums and crc too seriouly, they aren't
30 : Isibaar 262 * bullet-proof (should plug some .md5 here)...
31 : Isibaar 225 *
32 :     * compiles with something like:
33 :     * gcc -o xvid_bench xvid_bench.c -I../src/ -lxvidcore -lm
34 :     *
35 : edgomez 1054 ****************************************************************************/
36 : Isibaar 225
37 :     #include <stdio.h>
38 :     #include <stdlib.h>
39 : edgomez 1053 #include <string.h> /* for memset */
40 : Isibaar 225 #include <assert.h>
41 :    
42 : suxen_drol 860 #ifndef WIN32
43 : edgomez 1053 #include <sys/time.h> /* for gettimeofday */
44 : suxen_drol 860 #else
45 :     #include <time.h>
46 :     #endif
47 :    
48 :    
49 : Isibaar 225 #include "xvid.h"
50 :    
51 : edgomez 851 // inner guts
52 : Isibaar 225 #include "dct/idct.h"
53 :     #include "dct/fdct.h"
54 :     #include "image/colorspace.h"
55 :     #include "image/interpolate8x8.h"
56 :     #include "utils/mem_transfer.h"
57 : edgomez 1176 #include "quant/quant.h"
58 : Isibaar 225 #include "motion/sad.h"
59 :     #include "utils/emms.h"
60 :     #include "utils/timer.h"
61 :     #include "quant/quant_matrix.c"
62 :     #include "bitstream/cbp.h"
63 :    
64 : Isibaar 262 #include <math.h>
65 : suxen_drol 860
66 :     #ifndef M_PI
67 :     #define M_PI 3.14159265358979323846
68 :     #endif
69 :    
70 : edgomez 1053 const int speed_ref = 100; /* on slow machines, decrease this value */
71 : Isibaar 225
72 :     /*********************************************************************
73 :     * misc
74 :     *********************************************************************/
75 :    
76 : edgomez 1060 /* returns time in micro-s*/
77 : Isibaar 225 double gettime_usec()
78 :     {
79 : suxen_drol 860 #ifndef WIN32
80 : edgomez 1060 struct timeval tv;
81 :     gettimeofday(&tv, 0);
82 :     return tv.tv_sec*1.0e6 + tv.tv_usec;
83 : suxen_drol 860 #else
84 :     clock_t clk;
85 :     clk = clock();
86 :     return clk * 1000000 / CLOCKS_PER_SEC;
87 :     #endif
88 : Isibaar 225 }
89 :    
90 : edgomez 1060 /* returns squared deviates (mean(v*v)-mean(v)^2) of a 8x8 block */
91 : Isibaar 225 double sqr_dev(uint8_t v[8*8])
92 :     {
93 : edgomez 1060 double sum=0.;
94 :     double sum2=0.;
95 :     int n;
96 :     for (n=0;n<8*8;n++)
97 :     {
98 :     sum += v[n];
99 :     sum2 += v[n]*v[n];
100 :     }
101 :     sum2 /= n;
102 :     sum /= n;
103 :     return sum2-sum*sum;
104 : Isibaar 225 }
105 :    
106 :     /*********************************************************************
107 :     * cpu init
108 :     *********************************************************************/
109 :    
110 :     typedef struct {
111 : edgomez 1060 const char *name;
112 :     unsigned int cpu;
113 : Isibaar 225 } CPU;
114 :    
115 :     CPU cpu_list[] =
116 :     { { "PLAINC", 0 }
117 : edgomez 1122 #ifdef ARCH_IS_IA32
118 : edgomez 1060 , { "MMX ", XVID_CPU_MMX }
119 :     , { "MMXEXT", XVID_CPU_MMXEXT | XVID_CPU_MMX }
120 :     , { "SSE2 ", XVID_CPU_SSE2 | XVID_CPU_MMX }
121 :     , { "3DNOW ", XVID_CPU_3DNOW }
122 :     , { "3DNOWE", XVID_CPU_3DNOWEXT }
123 : edgomez 1122 #endif
124 : edgomez 1060 //, { "IA64 ", XVID_CPU_IA64 }
125 : edgomez 851 //, { "TSC ", XVID_CPU_TSC }
126 : edgomez 1060 , { 0, 0 } };
127 : Isibaar 225
128 : edgomez 1060 CPU cpu_short_list[] =
129 : Isibaar 225 { { "PLAINC", 0 }
130 : edgomez 1122 #ifdef ARCH_IS_IA32
131 : edgomez 1060 , { "MMX ", XVID_CPU_MMX }
132 : edgomez 851 //, { "MMXEXT", XVID_CPU_MMXEXT | XVID_CPU_MMX }
133 : edgomez 1122 #endif
134 : edgomez 1060 //, { "IA64 ", XVID_CPU_IA64 }
135 :     , { 0, 0 } };
136 : Isibaar 225
137 : edgomez 1060 CPU cpu_short_list2[] =
138 : Isibaar 225 { { "PLAINC", 0 }
139 : edgomez 1122 #ifdef ARCH_IS_IA32
140 : edgomez 1060 , { "MMX ", XVID_CPU_MMX }
141 :     , { "SSE2 ", XVID_CPU_SSE2 | XVID_CPU_MMX }
142 : edgomez 1122 #endif
143 : edgomez 1060 , { 0, 0 } };
144 : Isibaar 225
145 :    
146 :     int init_cpu(CPU *cpu)
147 :     {
148 : edgomez 1197 xvid_gbl_info_t xinfo;
149 : Isibaar 225
150 : edgomez 1197 /* Get the available CPU flags */
151 :     memset(&xinfo, 0, sizeof(xinfo));
152 :     xinfo.version = XVID_VERSION;
153 :     xvid_global(NULL, XVID_GBL_INFO, &xinfo, NULL);
154 :    
155 :     /* Are we trying to test a subset of the host CPU features */
156 :     if ((xinfo.cpu_flags & cpu->cpu) == cpu->cpu) {
157 :     int xerr;
158 :     xvid_gbl_init_t xinit;
159 :     memset(&xinit, 0, sizeof(xinit));
160 :     xinit.cpu_flags = cpu->cpu | XVID_CPU_FORCE;
161 :     xinit.version = XVID_VERSION;
162 :     xerr = xvid_global(NULL, XVID_GBL_INIT, &xinit, NULL);
163 :     if (xerr==XVID_ERR_FAIL) {
164 :     /* libxvidcore failed to init */
165 :     return 0;
166 :     }
167 :     } else {
168 :     /* The host CPU doesn't support some required feature for this test */
169 :     return(0);
170 : edgomez 1060 }
171 :     return 1;
172 : Isibaar 225 }
173 :    
174 :     /*********************************************************************
175 :     * test DCT
176 :     *********************************************************************/
177 :    
178 :     #define ABS(X) ((X)<0 ? -(X) : (X))
179 :    
180 :     void test_dct()
181 :     {
182 : edgomez 1060 const int nb_tests = 300*speed_ref;
183 :     int tst;
184 :     CPU *cpu;
185 :     int i;
186 :     short iDst0[8*8], iDst[8*8], fDst[8*8];
187 :     double overhead;
188 : Isibaar 225
189 : edgomez 1060 printf( "\n ===== test fdct/idct =====\n" );
190 : Isibaar 225
191 : edgomez 1060 for(i=0; i<8*8; ++i) iDst0[i] = (i*7-i*i) & 0x7f;
192 :     overhead = gettime_usec();
193 :     for(tst=0; tst<nb_tests; ++tst)
194 :     {
195 :     for(i=0; i<8*8; ++i) fDst[i] = iDst0[i];
196 :     for(i=0; i<8*8; ++i) iDst[i] = fDst[i];
197 :     }
198 :     overhead = gettime_usec() - overhead;
199 : Isibaar 225
200 : edgomez 1060 for(cpu = cpu_list; cpu->name!=0; ++cpu)
201 :     {
202 :     double t, PSNR, MSE;
203 : Isibaar 225
204 : edgomez 1060 if (!init_cpu(cpu))
205 :     continue;
206 : Isibaar 225
207 : edgomez 1060 t = gettime_usec();
208 :     emms();
209 :     for(tst=0; tst<nb_tests; ++tst)
210 :     {
211 :     for(i=0; i<8*8; ++i) fDst[i] = iDst0[i];
212 :     fdct(fDst);
213 :     for(i=0; i<8*8; ++i) iDst[i] = fDst[i];
214 :     idct(iDst);
215 :     }
216 :     emms();
217 :     t = (gettime_usec() - t - overhead) / nb_tests;
218 :     MSE = 0.;
219 :     for(i=0; i<8*8; ++i) {
220 :     double delta = 1.0*(iDst[i] - iDst0[i]);
221 :     MSE += delta*delta;
222 :     }
223 :     PSNR = (MSE==0.) ? 1.e6 : -4.3429448*log( MSE/64. );
224 :     printf( "%s - %.3f usec PSNR=%.3f MSE=%.3f\n",
225 :     cpu->name, t, PSNR, MSE );
226 :     if (ABS(MSE)>=64) printf( "*** CRC ERROR! ***\n" );
227 :     }
228 : Isibaar 225 }
229 :    
230 :     /*********************************************************************
231 :     * test SAD
232 :     *********************************************************************/
233 :    
234 :     void test_sad()
235 :     {
236 : edgomez 1060 const int nb_tests = 2000*speed_ref;
237 :     int tst;
238 :     CPU *cpu;
239 :     int i;
240 :     uint8_t Cur[16*16], Ref1[16*16], Ref2[16*16];
241 : Isibaar 225
242 : edgomez 1060 printf( "\n ====== test SAD ======\n" );
243 :     for(i=0; i<16*16;++i) {
244 :     Cur[i] = (i/5) ^ 0x05;
245 :     Ref1[i] = (i + 0x0b) & 0xff;
246 :     Ref2[i] = i ^ 0x76;
247 :     }
248 : Isibaar 225
249 : edgomez 1060 for(cpu = cpu_list; cpu->name!=0; ++cpu)
250 :     {
251 :     double t;
252 :     uint32_t s;
253 :     if (!init_cpu(cpu))
254 :     continue;
255 : Isibaar 225
256 : edgomez 1060 t = gettime_usec();
257 :     emms();
258 :     for(tst=0; tst<nb_tests; ++tst) s = sad8(Cur, Ref1, 16);
259 :     emms();
260 :     t = (gettime_usec() - t) / nb_tests;
261 :     printf( "%s - sad8 %.3f usec sad=%d\n", cpu->name, t, s );
262 :     if (s!=3776) printf( "*** CRC ERROR! ***\n" );
263 : Isibaar 225
264 : edgomez 1060 t = gettime_usec();
265 :     emms();
266 :     for(tst=0; tst<nb_tests; ++tst) s = sad16(Cur, Ref1, 16, -1);
267 :     emms();
268 :     t = (gettime_usec() - t) / nb_tests;
269 :     printf( "%s - sad16 %.3f usec sad=%d\n", cpu->name, t, s );
270 :     if (s!=27214) printf( "*** CRC ERROR! ***\n" );
271 : Isibaar 225
272 : edgomez 1060 t = gettime_usec();
273 :     emms();
274 :     for(tst=0; tst<nb_tests; ++tst) s = sad16bi(Cur, Ref1, Ref2, 16);
275 :     emms();
276 :     t = (gettime_usec() - t) / nb_tests;
277 :     printf( "%s - sad16bi %.3f usec sad=%d\n", cpu->name, t, s );
278 :     if (s!=26274) printf( "*** CRC ERROR! ***\n" );
279 : Isibaar 225
280 : edgomez 1060 t = gettime_usec();
281 :     emms();
282 :     for(tst=0; tst<nb_tests; ++tst) s = dev16(Cur, 16);
283 :     emms();
284 :     t = (gettime_usec() - t) / nb_tests;
285 :     printf( "%s - dev16 %.3f usec sad=%d\n", cpu->name, t, s );
286 :     if (s!=3344) printf( "*** CRC ERROR! ***\n" );
287 : Isibaar 225
288 : edgomez 1060 printf( " --- \n" );
289 :     }
290 : Isibaar 225 }
291 :    
292 :     /*********************************************************************
293 :     * test interpolation
294 :     *********************************************************************/
295 :    
296 :     #define ENTER \
297 : edgomez 1060 for(i=0; i<16*8; ++i) Dst[i] = 0; \
298 :     t = gettime_usec(); \
299 :     emms();
300 : Isibaar 225
301 :     #define LEAVE \
302 : edgomez 1060 emms(); \
303 :     t = (gettime_usec() - t) / nb_tests; \
304 :     iCrc = 0; \
305 :     for(i=0; i<16*8; ++i) { iCrc += Dst[i]^i; }
306 : Isibaar 225
307 :     #define TEST_MB(FUNC, R) \
308 : edgomez 1060 ENTER \
309 :     for(tst=0; tst<nb_tests; ++tst) (FUNC)(Dst, Src0, 16, (R)); \
310 :     LEAVE
311 : Isibaar 225
312 :     #define TEST_MB2(FUNC) \
313 : edgomez 1060 ENTER \
314 :     for(tst=0; tst<nb_tests; ++tst) (FUNC)(Dst, Src0, 16); \
315 :     LEAVE
316 : Isibaar 225
317 :    
318 :     void test_mb()
319 :     {
320 : edgomez 1060 const int nb_tests = 2000*speed_ref;
321 :     CPU *cpu;
322 :     const uint8_t Src0[16*9] = {
323 :     /* try to have every possible combinaison of rounding... */
324 :     0, 0, 1, 0, 2, 0, 3, 0, 4 ,0,0,0, 0,0,0,0
325 :     , 0, 1, 1, 1, 2, 1, 3, 1, 3 ,0,0,0, 0,0,0,0
326 :     , 0, 2, 1, 2, 2, 2, 3, 2, 2 ,0,0,0, 0,0,0,0
327 :     , 0, 3, 1, 3, 2, 3, 3, 3, 1 ,0,0,0, 0,0,0,0
328 :     , 1, 3, 0, 2, 1, 0, 2, 3, 4 ,0,0,0, 0,0,0,0
329 :     , 2, 2, 1, 2, 0, 1, 3, 5, 3 ,0,0,0, 0,0,0,0
330 :     , 3, 1, 2, 3, 1, 2, 2, 6, 2 ,0,0,0, 0,0,0,0
331 :     , 1, 0, 1, 3, 0, 3, 1, 6, 1 ,0,0,0, 0,0,0,0
332 :     , 4, 3, 2, 1, 2, 3, 4, 0, 3 ,0,0,0, 0,0,0,0
333 :     };
334 :     uint8_t Dst[16*8] = {0};
335 : Isibaar 225
336 : edgomez 1060 printf( "\n === test block motion ===\n" );
337 : Isibaar 225
338 : edgomez 1060 for(cpu = cpu_list; cpu->name!=0; ++cpu)
339 :     {
340 :     double t;
341 :     int tst, i, iCrc;
342 : Isibaar 225
343 : edgomez 1060 if (!init_cpu(cpu))
344 :     continue;
345 : Isibaar 225
346 : edgomez 1060 TEST_MB(interpolate8x8_halfpel_h, 0);
347 :     printf( "%s - interp- h-round0 %.3f usec iCrc=%d\n", cpu->name, t, iCrc );
348 :     if (iCrc!=8107) printf( "*** CRC ERROR! ***\n" );
349 : Isibaar 225
350 : edgomez 1060 TEST_MB(interpolate8x8_halfpel_h, 1);
351 :     printf( "%s - round1 %.3f usec iCrc=%d\n", cpu->name, t, iCrc );
352 :     if (iCrc!=8100) printf( "*** CRC ERROR! ***\n" );
353 : Isibaar 225
354 :    
355 : edgomez 1060 TEST_MB(interpolate8x8_halfpel_v, 0);
356 :     printf( "%s - interp- v-round0 %.3f usec iCrc=%d\n", cpu->name, t, iCrc );
357 :     if (iCrc!=8108) printf( "*** CRC ERROR! ***\n" );
358 : Isibaar 225
359 : edgomez 1060 TEST_MB(interpolate8x8_halfpel_v, 1);
360 :     printf( "%s - round1 %.3f usec iCrc=%d\n", cpu->name, t, iCrc );
361 :     if (iCrc!=8105) printf( "*** CRC ERROR! ***\n" );
362 : Isibaar 225
363 :    
364 : edgomez 1060 TEST_MB(interpolate8x8_halfpel_hv, 0);
365 :     printf( "%s - interp-hv-round0 %.3f usec iCrc=%d\n", cpu->name, t, iCrc );
366 :     if (iCrc!=8112) printf( "*** CRC ERROR! ***\n" );
367 : Isibaar 225
368 : edgomez 1060 TEST_MB(interpolate8x8_halfpel_hv, 1);
369 :     printf( "%s - round1 %.3f usec iCrc=%d\n", cpu->name, t, iCrc );
370 :     if (iCrc!=8103) printf( "*** CRC ERROR! ***\n" );
371 : Isibaar 225
372 : Isibaar 262
373 : edgomez 1060 /* this is a new function, as of 06.06.2002 */
374 : Isibaar 262 #if 0
375 : edgomez 1060 TEST_MB2(interpolate8x8_avrg);
376 :     printf( "%s - interpolate8x8_c %.3f usec iCrc=%d\n", cpu->name, t, iCrc );
377 :     if (iCrc!=8107) printf( "*** CRC ERROR! ***\n" );
378 : Isibaar 262 #endif
379 :    
380 : edgomez 1060 printf( " --- \n" );
381 :     }
382 : Isibaar 225 }
383 :    
384 :     /*********************************************************************
385 :     * test transfer
386 :     *********************************************************************/
387 :    
388 :     #define INIT_TRANSFER \
389 : edgomez 1060 for(i=0; i<8*32; ++i) { \
390 :     Src8[i] = i; Src16[i] = i; \
391 :     Dst8[i] = 0; Dst16[i] = 0; \
392 :     Ref1[i] = i^0x27; \
393 :     Ref2[i] = i^0x51; \
394 :     }
395 : Isibaar 225
396 :     #define TEST_TRANSFER_BEGIN(DST) \
397 : edgomez 1060 INIT_TRANSFER \
398 :     overhead = -gettime_usec(); \
399 :     for(tst=0; tst<nb_tests; ++tst) { \
400 :     for(i=0; i<8*32; ++i) (DST)[i] = i^0x6a;\
401 :     } \
402 :     overhead += gettime_usec(); \
403 :     t = gettime_usec(); \
404 :     emms(); \
405 :     for(tst=0; tst<nb_tests; ++tst) { \
406 :     for(i=0; i<8*32; ++i) (DST)[i] = i^0x6a;
407 : Isibaar 225
408 :    
409 :     #define TEST_TRANSFER_END(DST) \
410 : edgomez 1060 } \
411 :     emms(); \
412 :     t = (gettime_usec()-t -overhead) / nb_tests;\
413 :     s = 0; for(i=0; i<8*32; ++i) { s += (DST)[i]^i; }
414 : Isibaar 225
415 :     #define TEST_TRANSFER(FUNC, DST, SRC) \
416 : edgomez 1060 TEST_TRANSFER_BEGIN(DST); \
417 :     (FUNC)((DST), (SRC), 32); \
418 :     TEST_TRANSFER_END(DST)
419 : Isibaar 225
420 :    
421 :     #define TEST_TRANSFER2_BEGIN(DST, SRC) \
422 : edgomez 1060 INIT_TRANSFER \
423 :     overhead = -gettime_usec(); \
424 :     for(tst=0; tst<nb_tests; ++tst) { \
425 :     for(i=0; i<8*32; ++i) (DST)[i] = i^0x6a;\
426 :     for(i=0; i<8*32; ++i) (SRC)[i] = i^0x3e;\
427 :     } \
428 :     overhead += gettime_usec(); \
429 :     t = gettime_usec(); \
430 :     emms(); \
431 :     for(tst=0; tst<nb_tests; ++tst) { \
432 :     for(i=0; i<8*32; ++i) (DST)[i] = i^0x6a;\
433 :     for(i=0; i<8*32; ++i) (SRC)[i] = i^0x3e;
434 : Isibaar 225
435 :     #define TEST_TRANSFER2_END(DST) \
436 : edgomez 1060 } \
437 :     emms(); \
438 :     t = (gettime_usec()-t -overhead) / nb_tests;\
439 :     s = 0; for(i=0; i<8*32; ++i) { s += (DST)[i]; }
440 : Isibaar 225
441 :     #define TEST_TRANSFER2(FUNC, DST, SRC, R1) \
442 : edgomez 1060 TEST_TRANSFER2_BEGIN(DST,SRC); \
443 :     (FUNC)((DST), (SRC), (R1), 32); \
444 :     TEST_TRANSFER2_END(DST)
445 : Isibaar 225
446 :     #define TEST_TRANSFER3(FUNC, DST, SRC, R1, R2)\
447 : edgomez 1060 TEST_TRANSFER_BEGIN(DST); \
448 :     (FUNC)((DST), (SRC), (R1), (R2), 32); \
449 :     TEST_TRANSFER_END(DST)
450 : Isibaar 225
451 :     void test_transfer()
452 :     {
453 : edgomez 1060 const int nb_tests = 4000*speed_ref;
454 :     int i;
455 :     CPU *cpu;
456 :     uint8_t Src8[8*32], Dst8[8*32], Ref1[8*32], Ref2[8*32];
457 :     int16_t Src16[8*32], Dst16[8*32];
458 : Isibaar 225
459 : edgomez 1060 printf( "\n === test transfer ===\n" );
460 : Isibaar 225
461 : edgomez 1176 for(cpu = cpu_list; cpu->name!=0; ++cpu)
462 : edgomez 1060 {
463 :     double t, overhead;
464 :     int tst, s;
465 : Isibaar 225
466 : edgomez 1060 if (!init_cpu(cpu))
467 :     continue;
468 : Isibaar 225
469 : edgomez 1060 TEST_TRANSFER(transfer_8to16copy, Dst16, Src8);
470 :     printf( "%s - 8to16 %.3f usec crc=%d\n", cpu->name, t, s );
471 :     if (s!=28288) printf( "*** CRC ERROR! ***\n" );
472 : Isibaar 225
473 : edgomez 1060 TEST_TRANSFER(transfer_16to8copy, Dst8, Src16);
474 :     printf( "%s - 16to8 %.3f usec crc=%d\n", cpu->name, t, s );
475 :     if (s!=28288) printf( "*** CRC ERROR! ***\n" );
476 : Isibaar 225
477 : edgomez 1060 TEST_TRANSFER(transfer8x8_copy, Dst8, Src8);
478 :     printf( "%s - 8to8 %.3f usec crc=%d\n", cpu->name, t, s );
479 :     if (s!=20352) printf( "*** CRC ERROR! ***\n" );
480 : Isibaar 225
481 : edgomez 1060 TEST_TRANSFER(transfer_16to8add, Dst8, Src16);
482 :     printf( "%s - 16to8add %.3f usec crc=%d\n", cpu->name, t, s );
483 :     if (s!=25536) printf( "*** CRC ERROR! ***\n" );
484 : Isibaar 225
485 : edgomez 1060 TEST_TRANSFER2(transfer_8to16sub, Dst16, Src8, Ref1);
486 :     printf( "%s - 8to16sub %.3f usec crc1=%d ", cpu->name, t, s );
487 :     if (s!=28064) printf( "*** CRC ERROR! ***\n" );
488 :     s = 0; for(i=0; i<8*32; ++i) { s += (Src8[i]-Ref1[i])&i; }
489 :     printf( "crc2=%d\n", s);
490 :     if (s!=16256) printf( "*** CRC ERROR! ***\n" );
491 : Isibaar 262 #if 1
492 : edgomez 1060 TEST_TRANSFER3(transfer_8to16sub2, Dst16, Src8, Ref1, Ref2);
493 :     printf( "%s - 8to16sub2 %.3f usec crc=%d\n", cpu->name, t, s );
494 :     if (s!=20384) printf( "*** CRC ERROR! ***\n" );
495 : edgomez 851 // for(i=0; i<64; ++i) printf( "[%d]", Dst16[i]);
496 :     // printf("\n");
497 : Isibaar 262 #endif
498 : edgomez 1060 printf( " --- \n" );
499 :     }
500 : Isibaar 225 }
501 :    
502 :     /*********************************************************************
503 :     * test quantization
504 :     *********************************************************************/
505 :    
506 : Isibaar 262 #define TEST_QUANT(FUNC, DST, SRC) \
507 : edgomez 1060 t = gettime_usec(); \
508 :     for(s=0,qm=1; qm<=255; ++qm) { \
509 :     for(i=0; i<8*8; ++i) Quant[i] = qm; \
510 :     set_inter_matrix( Quant ); \
511 :     emms(); \
512 :     for(q=1; q<=max_Q; ++q) { \
513 :     for(tst=0; tst<nb_tests; ++tst) \
514 :     (FUNC)((DST), (SRC), q); \
515 :     for(i=0; i<64; ++i) s+=(DST)[i]^i^qm; \
516 :     } \
517 :     emms(); \
518 :     } \
519 :     t = (gettime_usec()-t-overhead)/nb_tests/qm;\
520 :     s = (s&0xffff)^(s>>16)
521 : Isibaar 225
522 : Isibaar 262 #define TEST_QUANT2(FUNC, DST, SRC) \
523 : edgomez 1060 t = gettime_usec(); \
524 :     for(s=0,qm=1; qm<=255; ++qm) { \
525 :     for(i=0; i<8*8; ++i) Quant[i] = qm; \
526 :     set_intra_matrix( Quant ); \
527 :     emms(); \
528 :     for(q=1; q<=max_Q; ++q) { \
529 :     for(tst=0; tst<nb_tests; ++tst) \
530 :     (FUNC)((DST), (SRC), q, q); \
531 :     for(i=0; i<64; ++i) s+=(DST)[i]^i^qm; \
532 :     } \
533 :     emms(); \
534 :     } \
535 :     t = (gettime_usec()-t-overhead)/nb_tests/qm;\
536 :     s = (s&0xffff)^(s>>16)
537 : Isibaar 225
538 :     void test_quant()
539 :     {
540 : edgomez 1060 const int nb_tests = 1*speed_ref;
541 :     const int max_Q = 31;
542 :     int i, qm;
543 :     CPU *cpu;
544 :     int16_t Src[8*8], Dst[8*8];
545 :     uint8_t Quant[8*8];
546 : Isibaar 225
547 : edgomez 1060 printf( "\n ===== test quant =====\n" );
548 : Isibaar 225
549 : edgomez 1060 /* we deliberately enfringe the norm's specified range [-127,127], */
550 :     /* to test the robustness of the iquant module */
551 :     for(i=0; i<64; ++i) {
552 :     Src[i] = 1 + (i-32) * (i&6);
553 :     Dst[i] = 0;
554 :     }
555 : Isibaar 225
556 : edgomez 1176 for(cpu = cpu_list; cpu->name!=0; ++cpu)
557 : edgomez 1060 {
558 :     double t, overhead;
559 :     int tst, q;
560 :     uint32_t s;
561 : Isibaar 225
562 : edgomez 1060 if (!init_cpu(cpu))
563 :     continue;
564 : Isibaar 225
565 : edgomez 1060 overhead = -gettime_usec();
566 :     for(s=0,qm=1; qm<=255; ++qm) {
567 :     for(i=0; i<8*8; ++i) Quant[i] = qm;
568 :     set_inter_matrix( Quant );
569 :     for(q=1; q<=max_Q; ++q)
570 :     for(i=0; i<64; ++i) s+=Dst[i]^i^qm;
571 :     }
572 :     overhead += gettime_usec();
573 : Isibaar 225
574 : Isibaar 262 #if 1
575 : edgomez 1176 TEST_QUANT2(quant_mpeg_intra, Dst, Src);
576 :     printf( "%s - quant_mpeg_intra %.3f usec crc=%d\n", cpu->name, t, s );
577 : edgomez 1060 if (s!=29809) printf( "*** CRC ERROR! ***\n" );
578 : Isibaar 225
579 : edgomez 1176 TEST_QUANT(quant_mpeg_inter, Dst, Src);
580 :     printf( "%s - quant_mpeg_inter %.3f usec crc=%d\n", cpu->name, t, s );
581 : edgomez 1060 if (s!=12574) printf( "*** CRC ERROR! ***\n" );
582 : Isibaar 262 #endif
583 :     #if 1
584 : edgomez 1176 TEST_QUANT2(dequant_mpeg_intra, Dst, Src);
585 :     printf( "%s - dequant_mpeg_intra %.3f usec crc=%d\n", cpu->name, t, s );
586 : edgomez 1060 if (s!=24052) printf( "*** CRC ERROR! ***\n" );
587 : Isibaar 225
588 : edgomez 1176 TEST_QUANT(dequant_mpeg_inter, Dst, Src);
589 :     printf( "%s - dequant_mpeg_inter %.3f usec crc=%d\n", cpu->name, t, s );
590 : edgomez 1060 if (s!=63847) printf( "*** CRC ERROR! ***\n" );
591 : Isibaar 262 #endif
592 :     #if 1
593 : edgomez 1176 TEST_QUANT2(quant_h263_intra, Dst, Src);
594 :     printf( "%s - quant_h263_intra %.3f usec crc=%d\n", cpu->name, t, s );
595 : edgomez 1060 if (s!=25662) printf( "*** CRC ERROR! ***\n" );
596 : Isibaar 225
597 : edgomez 1176 TEST_QUANT(quant_h263_inter, Dst, Src);
598 :     printf( "%s - quant_h263_inter %.3f usec crc=%d\n", cpu->name, t, s );
599 : edgomez 1060 if (s!=23972) printf( "*** CRC ERROR! ***\n" );
600 : Isibaar 262 #endif
601 :     #if 1
602 : edgomez 1176 TEST_QUANT2(dequant_h263_intra, Dst, Src);
603 :     printf( "%s - dequant_h263_intra %.3f usec crc=%d\n", cpu->name, t, s );
604 : edgomez 1060 if (s!=49900) printf( "*** CRC ERROR! ***\n" );
605 : Isibaar 225
606 : edgomez 1176 TEST_QUANT(dequant_h263_inter, Dst, Src);
607 :     printf( "%s - dequant_h263_inter %.3f usec crc=%d\n", cpu->name, t, s );
608 : edgomez 1060 if (s!=48899) printf( "*** CRC ERROR! ***\n" );
609 : Isibaar 262 #endif
610 : edgomez 1060 printf( " --- \n" );
611 :     }
612 : Isibaar 225 }
613 :    
614 :     /*********************************************************************
615 :     * test non-zero AC counting
616 :     *********************************************************************/
617 :    
618 :     #define TEST_CBP(FUNC, SRC) \
619 : edgomez 1060 t = gettime_usec(); \
620 :     emms(); \
621 :     for(tst=0; tst<nb_tests; ++tst) { \
622 :     cbp = (FUNC)((SRC)); \
623 :     } \
624 :     emms(); \
625 :     t = (gettime_usec()-t ) / nb_tests;
626 : Isibaar 225
627 :     void test_cbp()
628 :     {
629 : edgomez 1060 const int nb_tests = 10000*speed_ref;
630 :     int i;
631 :     CPU *cpu;
632 :     int16_t Src1[6*64], Src2[6*64], Src3[6*64], Src4[6*64];
633 : Isibaar 225
634 : edgomez 1060 printf( "\n ===== test cbp =====\n" );
635 : Isibaar 225
636 : edgomez 1060 for(i=0; i<6*64; ++i) {
637 :     Src1[i] = (i*i*3/8192)&(i/64)&1; /* 'random' */
638 :     Src2[i] = (i<3*64); /* half-full */
639 :     Src3[i] = ((i+32)>3*64);
640 :     Src4[i] = (i==(3*64+2) || i==(5*64+9));
641 :     }
642 : Isibaar 225
643 : edgomez 1176 for(cpu = cpu_list; cpu->name!=0; ++cpu)
644 : edgomez 1060 {
645 :     double t;
646 :     int tst, cbp;
647 : Isibaar 225
648 : edgomez 1060 if (!init_cpu(cpu))
649 :     continue;
650 :    
651 :     TEST_CBP(calc_cbp, Src1);
652 :     printf( "%s - calc_cbp#1 %.3f usec cbp=0x%x\n", cpu->name, t, cbp );
653 :     if (cbp!=0x15) printf( "*** CRC ERROR! ***\n" );
654 :     TEST_CBP(calc_cbp, Src2);
655 :     printf( "%s - calc_cbp#2 %.3f usec cbp=0x%x\n", cpu->name, t, cbp );
656 :     if (cbp!=0x38) printf( "*** CRC ERROR! ***\n" );
657 :     TEST_CBP(calc_cbp, Src3);
658 :     printf( "%s - calc_cbp#3 %.3f usec cbp=0x%x\n", cpu->name, t, cbp );
659 :     if (cbp!=0x0f) printf( "*** CRC ERROR! ***\n" );
660 :     TEST_CBP(calc_cbp, Src4);
661 :     printf( "%s - calc_cbp#4 %.3f usec cbp=0x%x\n", cpu->name, t, cbp );
662 :     if (cbp!=0x05) printf( "*** CRC ERROR! ***\n" );
663 :     printf( " --- \n" );
664 :     }
665 : Isibaar 225 }
666 :    
667 :     /*********************************************************************
668 : Isibaar 262 * fdct/idct IEEE1180 compliance
669 :     *********************************************************************/
670 :    
671 :     typedef struct {
672 : edgomez 1060 long Errors[64];
673 :     long Sqr_Errors[64];
674 :     long Max_Errors[64];
675 :     long Nb;
676 : Isibaar 262 } STATS_8x8;
677 :    
678 :     void init_stats(STATS_8x8 *S)
679 :     {
680 : edgomez 1060 int i;
681 :     for(i=0; i<64; ++i) {
682 :     S->Errors[i] = 0;
683 :     S->Sqr_Errors[i] = 0;
684 :     S->Max_Errors[i] = 0;
685 :     }
686 :     S->Nb = 0;
687 : Isibaar 262 }
688 :    
689 :     void store_stats(STATS_8x8 *S, short Blk[64], short Ref[64])
690 :     {
691 : edgomez 1060 int i;
692 :     for(i=0; i<64; ++i)
693 :     {
694 :     short Err = Blk[i] - Ref[i];
695 :     S->Errors[i] += Err;
696 :     S->Sqr_Errors[i] += Err * Err;
697 :     if (Err<0) Err = -Err;
698 :     if (S->Max_Errors[i]<Err)
699 :     S->Max_Errors[i] = Err;
700 :     }
701 :     S->Nb++;
702 : Isibaar 262 }
703 :    
704 :     void print_stats(STATS_8x8 *S)
705 :     {
706 : edgomez 1060 int i;
707 :     double Norm;
708 : Isibaar 262
709 : edgomez 1060 assert(S->Nb>0);
710 :     Norm = 1. / (double)S->Nb;
711 :     printf("\n== Max absolute values of errors ==\n");
712 :     for(i=0; i<64; i++) {
713 :     printf(" %4ld", S->Max_Errors[i]);
714 :     if ((i&7)==7) printf("\n");
715 :     }
716 : Isibaar 262
717 : edgomez 1060 printf("\n== Mean square errors ==\n");
718 :     for(i=0; i<64; i++)
719 :     {
720 :     double Err = Norm * (double)S->Sqr_Errors[i];
721 :     printf(" %.3f", Err);
722 :     if ((i&7)==7) printf("\n");
723 :     }
724 : Isibaar 262
725 : edgomez 1060 printf("\n== Mean errors ==\n");
726 :     for(i=0; i<64; i++)
727 :     {
728 :     double Err = Norm * (double)S->Errors[i];
729 :     printf(" %.3f", Err);
730 :     if ((i&7)==7) printf("\n");
731 :     }
732 :     printf("\n");
733 : Isibaar 262 }
734 :    
735 :     static const char *CHECK(double v, double l) {
736 : edgomez 1060 if (fabs(v)<=l) return "ok";
737 :     else return "FAIL!";
738 : Isibaar 262 }
739 :    
740 :     void report_stats(STATS_8x8 *S, const double *Limits)
741 :     {
742 : edgomez 1060 int i;
743 :     double Norm, PE, PMSE, OMSE, PME, OME;
744 : Isibaar 262
745 : edgomez 1060 assert(S->Nb>0);
746 :     Norm = 1. / (double)S->Nb;
747 :     PE = 0.;
748 :     for(i=0; i<64; i++) {
749 :     if (PE<S->Max_Errors[i])
750 :     PE = S->Max_Errors[i];
751 :     }
752 : Isibaar 262
753 : edgomez 1060 PMSE = 0.;
754 :     OMSE = 0.;
755 :     for(i=0; i<64; i++)
756 :     {
757 :     double Err = Norm * (double)S->Sqr_Errors[i];
758 :     OMSE += Err;
759 :     if (PMSE < Err) PMSE = Err;
760 :     }
761 :     OMSE /= 64.;
762 : Isibaar 262
763 : edgomez 1060 PME = 0.;
764 :     OME = 0.;
765 :     for(i=0; i<64; i++)
766 :     {
767 :     double Err = Norm * (double)S->Errors[i];
768 :     OME += Err;
769 :     Err = fabs(Err);
770 :     if (PME < Err) PME = Err;
771 :     }
772 :     OME /= 64.;
773 : Isibaar 262
774 : edgomez 1060 printf( "Peak error: %4.4f\n", PE );
775 :     printf( "Peak MSE: %4.4f\n", PMSE );
776 :     printf( "Overall MSE: %4.4f\n", OMSE );
777 :     printf( "Peak ME: %4.4f\n", PME );
778 :     printf( "Overall ME: %4.4f\n", OME );
779 :    
780 :     if (Limits!=0)
781 :     {
782 :     printf( "[PE<=%.4f %s] ", Limits[0], CHECK(PE, Limits[0]) );
783 :     printf( "\n" );
784 :     printf( "[PMSE<=%.4f %s]", Limits[1], CHECK(PMSE, Limits[1]) );
785 :     printf( "[OMSE<=%.4f %s]", Limits[2], CHECK(OMSE, Limits[2]) );
786 :     printf( "\n" );
787 :     printf( "[PME<=%.4f %s] ", Limits[3], CHECK(PME , Limits[3]) );
788 :     printf( "[OME<=%.4f %s] ", Limits[4], CHECK(OME , Limits[4]) );
789 :     printf( "\n" );
790 :     }
791 : Isibaar 262 }
792 :    
793 : edgomez 1053 ///* ////////////////////////////////////////////////////// */
794 : Isibaar 262 /* Pseudo-random generator specified by IEEE 1180 */
795 :    
796 :     static long ieee_seed = 1;
797 :     static void ieee_reseed(long s) {
798 : edgomez 1060 ieee_seed = s;
799 : Isibaar 262 }
800 :     static long ieee_rand(int Min, int Max)
801 :     {
802 : edgomez 1060 static double z = (double) 0x7fffffff;
803 : Isibaar 262
804 : edgomez 1060 long i,j;
805 :     double x;
806 : Isibaar 262
807 : edgomez 1060 ieee_seed = (ieee_seed * 1103515245) + 12345;
808 :     i = ieee_seed & 0x7ffffffe;
809 :     x = ((double) i) / z;
810 :     x *= (Max-Min+1);
811 :     j = (long)x;
812 :     j = j + Min;
813 :     assert(j>=Min && j<=Max);
814 :     return (short)j;
815 : Isibaar 262 }
816 :    
817 :     #define CLAMP(x, M) (x) = ((x)<-(M)) ? (-(M)) : ((x)>=(M) ? ((M)-1) : (x))
818 :    
819 :     static double Cos[8][8];
820 :     static void init_ref_dct()
821 :     {
822 : edgomez 1060 int i, j;
823 :     for(i=0; i<8; i++)
824 :     {
825 :     double scale = (i == 0) ? sqrt(0.125) : 0.5;
826 :     for (j=0; j<8; j++)
827 :     Cos[i][j] = scale*cos( (M_PI/8.0)*i*(j + 0.5) );
828 :     }
829 : Isibaar 262 }
830 :    
831 :     void ref_idct(short *M)
832 :     {
833 : edgomez 1060 int i, j, k;
834 :     double Tmp[8][8];
835 : Isibaar 262
836 : edgomez 1060 for(i=0; i<8; i++) {
837 :     for(j=0; j<8; j++)
838 :     {
839 :     double Sum = 0.0;
840 :     for (k=0; k<8; k++) Sum += Cos[k][j]*M[8*i+k];
841 :     Tmp[i][j] = Sum;
842 :     }
843 :     }
844 :     for(i=0; i<8; i++) {
845 :     for(j=0; j<8; j++) {
846 :     double Sum = 0.0;
847 :     for (k=0; k<8; k++) Sum += Cos[k][i]*Tmp[k][j];
848 :     M[8*i+j] = (short)floor(Sum + .5);
849 :     }
850 :     }
851 : Isibaar 262 }
852 :    
853 :     void ref_fdct(short *M)
854 :     {
855 : edgomez 1060 int i, j, k;
856 :     double Tmp[8][8];
857 : Isibaar 262
858 : edgomez 1060 for(i=0; i<8; i++) {
859 :     for(j=0; j<8; j++)
860 :     {
861 :     double Sum = 0.0;
862 :     for (k=0; k<8; k++) Sum += Cos[j][k]*M[8*i+k];
863 :     Tmp[i][j] = Sum;
864 :     }
865 :     }
866 :     for(i=0; i<8; i++) {
867 :     for(j=0; j<8; j++) {
868 :     double Sum = 0.0;
869 :     for (k=0; k<8; k++) Sum += Cos[i][k]*Tmp[k][j];
870 :     M[8*i+j] = (short)floor(Sum + 0.5);
871 :     }
872 :     }
873 : Isibaar 262 }
874 :    
875 :     void test_IEEE1180_compliance(int Min, int Max, int Sign)
876 :     {
877 : edgomez 1060 static const double ILimits[5] = { 1., 0.06, 0.02, 0.015, 0.0015 };
878 :     int Loops = 10000;
879 :     int i, m, n;
880 :     short Blk0[64]; /* reference */
881 :     short Blk[64], iBlk[64];
882 :     short Ref_FDCT[64];
883 :     short Ref_IDCT[64];
884 : Isibaar 262
885 : edgomez 1060 STATS_8x8 FStats; /* forward dct stats */
886 :     STATS_8x8 IStats; /* inverse dct stats */
887 : Isibaar 262
888 : edgomez 1060 CPU *cpu;
889 : Isibaar 262
890 : edgomez 1060 init_ref_dct();
891 : Isibaar 262
892 : edgomez 1060 for(cpu = cpu_list; cpu->name!=0; ++cpu)
893 :     {
894 :     if (!init_cpu(cpu))
895 :     continue;
896 : Isibaar 262
897 : edgomez 1060 printf( "\n===== IEEE test for %s ==== (Min=%d Max=%d Sign=%d Loops=%d)\n",
898 :     cpu->name, Min, Max, Sign, Loops);
899 : Isibaar 262
900 : edgomez 1060 init_stats(&IStats);
901 :     init_stats(&FStats);
902 : Isibaar 262
903 : edgomez 1060 ieee_reseed(1);
904 :     for(n=0; n<Loops; ++n)
905 :     {
906 :     for(i=0; i<64; ++i)
907 :     Blk0[i] = (short)ieee_rand(Min,Max) * Sign;
908 : Isibaar 262
909 : edgomez 1060 /* hmm, I'm not quite sure this is exactly */
910 :     /* the tests described in the norm. check... */
911 : Isibaar 262
912 : edgomez 1060 memcpy(Ref_FDCT, Blk0, 64*sizeof(short));
913 :     ref_fdct(Ref_FDCT);
914 :     for(i=0; i<64; i++) CLAMP( Ref_FDCT[i], 2048 );
915 : Isibaar 262
916 : edgomez 1060 memcpy(Blk, Blk0, 64*sizeof(short));
917 :     emms(); fdct(Blk); emms();
918 :     for(i=0; i<64; i++) CLAMP( Blk[i], 2048 );
919 : Isibaar 262
920 : edgomez 1060 store_stats(&FStats, Blk, Ref_FDCT);
921 : Isibaar 262
922 :    
923 : edgomez 1060 memcpy(Ref_IDCT, Ref_FDCT, 64*sizeof(short));
924 :     ref_idct(Ref_IDCT);
925 :     for (i=0; i<64; i++) CLAMP( Ref_IDCT[i], 256 );
926 : Isibaar 262
927 : edgomez 1060 memcpy(iBlk, Ref_FDCT, 64*sizeof(short));
928 :     emms(); idct(iBlk); emms();
929 :     for(i=0; i<64; i++) CLAMP( iBlk[i], 256 );
930 : Isibaar 262
931 : edgomez 1060 store_stats(&IStats, iBlk, Ref_IDCT);
932 :     }
933 : Isibaar 262
934 :    
935 : edgomez 1060 printf( "\n -- FDCT report --\n" );
936 : edgomez 851 // print_stats(&FStats);
937 : edgomez 1060 report_stats(&FStats, 0); /* so far I know, IEEE1180 says nothing for fdct */
938 : Isibaar 262
939 : edgomez 1060 for(i=0; i<64; i++) Blk[i] = 0;
940 :     emms(); fdct(Blk); emms();
941 :     for(m=i=0; i<64; i++) if (Blk[i]!=0) m++;
942 :     printf( "FDCT(0) == 0 ? %s\n", (m!=0) ? "NOPE!" : "yup." );
943 : Isibaar 262
944 : edgomez 1060 printf( "\n -- IDCT report --\n" );
945 : edgomez 851 // print_stats(&IStats);
946 : edgomez 1060 report_stats(&IStats, ILimits);
947 : Isibaar 262
948 :    
949 : edgomez 1060 for(i=0; i<64; i++) Blk[i] = 0;
950 :     emms(); idct(Blk); emms();
951 :     for(m=i=0; i<64; i++) if (Blk[i]!=0) m++;
952 :     printf( "IDCT(0) == 0 ? %s\n", (m!=0) ? "NOPE!" : "yup." );
953 :     }
954 : Isibaar 262 }
955 :    
956 :    
957 :     void test_dct_saturation(int Min, int Max)
958 :     {
959 : edgomez 1060 /* test behaviour on input range fringe */
960 : Isibaar 262
961 : edgomez 1060 int i, n, p;
962 :     CPU *cpu;
963 : edgomez 1053 // const short IDCT_MAX = 2047; /* 12bits input */
964 : edgomez 851 // const short IDCT_MIN = -2048;
965 : edgomez 1053 // const short IDCT_OUT = 256; /* 9bits ouput */
966 : edgomez 1060 const int Partitions = 4;
967 :     const int Loops = 10000 / Partitions;
968 : Isibaar 262
969 : edgomez 1060 init_ref_dct();
970 : Isibaar 262
971 : edgomez 1060 for(cpu = cpu_list; cpu->name!=0; ++cpu)
972 :     {
973 :     short Blk0[64], Blk[64];
974 :     STATS_8x8 Stats;
975 : Isibaar 262
976 : edgomez 1060 if (!init_cpu(cpu))
977 :     continue;
978 : Isibaar 262
979 : edgomez 1060 printf( "\n===== IEEE test for %s Min=%d Max=%d =====\n",
980 :     cpu->name, Min, Max );
981 : Isibaar 262
982 : edgomez 1060 /* FDCT tests // */
983 : Isibaar 262
984 : edgomez 1060 init_stats(&Stats);
985 : Isibaar 262
986 : edgomez 1060 /* test each computation channels separately */
987 :     for(i=0; i<64; i++) Blk[i] = Blk0[i] = ((i/8)==(i%8)) ? Max : 0;
988 :     ref_fdct(Blk0);
989 :     emms(); fdct(Blk); emms();
990 :     store_stats(&Stats, Blk, Blk0);
991 : Isibaar 262
992 : edgomez 1060 for(i=0; i<64; i++) Blk[i] = Blk0[i] = ((i/8)==(i%8)) ? Min : 0;
993 :     ref_fdct(Blk0);
994 :     emms(); fdct(Blk); emms();
995 :     store_stats(&Stats, Blk, Blk0);
996 : Isibaar 262
997 : edgomez 1060 /* randomly saturated inputs */
998 :     for(p=0; p<Partitions; ++p)
999 :     {
1000 :     for(n=0; n<Loops; ++n)
1001 :     {
1002 :     for(i=0; i<64; ++i)
1003 :     Blk0[i] = Blk[i] = (ieee_rand(0,Partitions)>=p)? Max : Min;
1004 :     ref_fdct(Blk0);
1005 :     emms(); fdct(Blk); emms();
1006 :     store_stats(&Stats, Blk, Blk0);
1007 :     }
1008 :     }
1009 :     printf( "\n -- FDCT saturation report --\n" );
1010 :     report_stats(&Stats, 0);
1011 : Isibaar 262
1012 :    
1013 : edgomez 1060 /* IDCT tests // */
1014 : Isibaar 262 #if 0
1015 : edgomez 1060 /* no finished yet */
1016 : Isibaar 262
1017 : edgomez 1060 init_stats(&Stats);
1018 : Isibaar 262
1019 : edgomez 1060 /* test each computation channel separately */
1020 :     for(i=0; i<64; i++) Blk[i] = Blk0[i] = ((i/8)==(i%8)) ? IDCT_MAX : 0;
1021 :     ref_idct(Blk0);
1022 :     emms(); idct(Blk); emms();
1023 :     for(i=0; i<64; i++) { CLAMP(Blk0[i], IDCT_OUT); CLAMP(Blk[i], IDCT_OUT); }
1024 :     store_stats(&Stats, Blk, Blk0);
1025 : Isibaar 262
1026 : edgomez 1060 for(i=0; i<64; i++) Blk[i] = Blk0[i] = ((i/8)==(i%8)) ? IDCT_MIN : 0;
1027 :     ref_idct(Blk0);
1028 :     emms(); idct(Blk); emms();
1029 :     for(i=0; i<64; i++) { CLAMP(Blk0[i], IDCT_OUT); CLAMP(Blk[i], IDCT_OUT); }
1030 :     store_stats(&Stats, Blk, Blk0);
1031 : Isibaar 262
1032 : edgomez 1060 /* randomly saturated inputs */
1033 :     for(p=0; p<Partitions; ++p)
1034 :     {
1035 :     for(n=0; n<Loops; ++n)
1036 :     {
1037 :     for(i=0; i<64; ++i)
1038 :     Blk0[i] = Blk[i] = (ieee_rand(0,Partitions)>=p)? IDCT_MAX : IDCT_MIN;
1039 :     ref_idct(Blk0);
1040 :     emms(); idct(Blk); emms();
1041 :     for(i=0; i<64; i++) { CLAMP(Blk0[i],IDCT_OUT); CLAMP(Blk[i],IDCT_OUT); }
1042 :     store_stats(&Stats, Blk, Blk0);
1043 :     }
1044 :     }
1045 : Isibaar 262
1046 : edgomez 1060 printf( "\n -- IDCT saturation report --\n" );
1047 :     print_stats(&Stats);
1048 :     report_stats(&Stats, 0);
1049 : Isibaar 262 #endif
1050 : edgomez 1060 }
1051 : Isibaar 262 }
1052 :    
1053 :     /*********************************************************************
1054 : Isibaar 225 * measure raw decoding speed
1055 :     *********************************************************************/
1056 :    
1057 :     void test_dec(const char *name, int width, int height, int with_chksum)
1058 :     {
1059 : edgomez 1060 FILE *f = 0;
1060 :     void *dechandle = 0;
1061 :     int xerr;
1062 :     xvid_gbl_init_t xinit;
1063 :     xvid_dec_create_t xparam;
1064 :     xvid_dec_frame_t xframe;
1065 : Isibaar 225 double t = 0.;
1066 :     int nb = 0;
1067 : edgomez 1060 uint8_t *buf = 0;
1068 :     uint8_t *rgb_out = 0;
1069 :     int buf_size, pos;
1070 :     uint32_t chksum = 0;
1071 : Isibaar 225
1072 : edgomez 1060 memset(&xinit, 0, sizeof(xinit));
1073 : Isibaar 262 xinit.cpu_flags = XVID_CPU_MMX | XVID_CPU_FORCE;
1074 : edgomez 1060 xinit.version = XVID_VERSION;
1075 :     xvid_global(NULL, 0, &xinit, NULL);
1076 : Isibaar 225
1077 : edgomez 1060 memset(&xparam, 0, sizeof(xparam));
1078 : Isibaar 225 xparam.width = width;
1079 :     xparam.height = height;
1080 : edgomez 1060 xparam.version = XVID_VERSION;
1081 : Isibaar 225 xerr = xvid_decore(NULL, XVID_DEC_CREATE, &xparam, NULL);
1082 : edgomez 1060 if (xerr==XVID_ERR_FAIL) {
1083 :     printf("can't init decoder (err=%d)\n", xerr);
1084 :     return;
1085 : Isibaar 225 }
1086 :     dechandle = xparam.handle;
1087 :    
1088 :    
1089 :     f = fopen(name, "rb");
1090 : edgomez 1060 if (f==0) {
1091 :     printf( "can't open file '%s'\n", name);
1092 :     return;
1093 :     }
1094 :     fseek(f, 0, SEEK_END);
1095 :     buf_size = ftell(f);
1096 :     fseek(f, 0, SEEK_SET);
1097 :     if (buf_size<=0) {
1098 :     printf("error while stating file\n");
1099 :     goto End;
1100 :     }
1101 :     else printf( "Input size: %d\n", buf_size);
1102 : Isibaar 225
1103 : edgomez 1060 buf = malloc(buf_size); /* should be enuf' */
1104 :     rgb_out = calloc(4, width*height); /* <-room for _RGB24 */
1105 :     if (buf==0 || rgb_out==0) {
1106 :     printf( "malloc failed!\n" );
1107 :     goto End;
1108 :     }
1109 : Isibaar 225
1110 : edgomez 1060 if (fread(buf, buf_size, 1, f)!=1) {
1111 :     printf( "file-read failed\n" );
1112 :     goto End;
1113 :     }
1114 : Isibaar 225
1115 : edgomez 1060 nb = 0;
1116 :     pos = 0;
1117 :     t = -gettime_usec();
1118 :     while(1) {
1119 :     memset(&xframe, 0, sizeof(xframe));
1120 :     xframe.version = XVID_VERSION;
1121 :     xframe.bitstream = buf + pos;
1122 :     xframe.length = buf_size - pos;
1123 :     xframe.output.plane[0] = rgb_out;
1124 :     xframe.output.stride[0] = width;
1125 :     xframe.output.csp = XVID_CSP_BGR;
1126 :     xerr = xvid_decore(dechandle, XVID_DEC_DECODE, &xframe, 0);
1127 :     nb++;
1128 :     pos += xframe.length;
1129 :     if (with_chksum) {
1130 :     int k = width*height;
1131 :     uint32_t *ptr = (uint32_t *)rgb_out;
1132 :     while(k-->0) chksum += *ptr++;
1133 :     }
1134 :     if (pos==buf_size)
1135 :     break;
1136 :     if (xerr==XVID_ERR_FAIL) {
1137 :     printf("decoding failed for frame #%d (err=%d)!\n", nb, xerr);
1138 :     break;
1139 :     }
1140 :     }
1141 :     t += gettime_usec();
1142 :     if (t>0.)
1143 :     printf( "%d frames decoded in %.3f s -> %.1f FPS\n", nb, t*1.e-6f, (float)(nb*1.e6f/t) );
1144 :     if (with_chksum)
1145 :     printf("checksum: 0x%.8x\n", chksum);
1146 : Isibaar 225
1147 : edgomez 1060 End:
1148 :     if (rgb_out!=0) free(rgb_out);
1149 :     if (buf!=0) free(buf);
1150 :     if (dechandle!=0) {
1151 :     xerr= xvid_decore(dechandle, XVID_DEC_DESTROY, NULL, NULL);
1152 :     if (xerr==XVID_ERR_FAIL)
1153 :     printf("destroy-decoder failed (err=%d)!\n", xerr);
1154 :     }
1155 :     if (f!=0) fclose(f);
1156 : Isibaar 225 }
1157 :    
1158 :     /*********************************************************************
1159 :     * non-regression tests
1160 :     *********************************************************************/
1161 :    
1162 :     void test_bugs1()
1163 :     {
1164 : edgomez 1060 CPU *cpu;
1165 : Isibaar 225
1166 : edgomez 1060 printf( "\n ===== (de)quant4_intra saturation bug? =====\n" );
1167 : Isibaar 225
1168 : edgomez 1176 for(cpu = cpu_list; cpu->name!=0; ++cpu)
1169 : edgomez 1060 {
1170 :     int i;
1171 :     int16_t Src[8*8], Dst[8*8];
1172 : Isibaar 225
1173 : edgomez 1060 if (!init_cpu(cpu))
1174 :     continue;
1175 : Isibaar 225
1176 : edgomez 1060 for(i=0; i<64; ++i) Src[i] = i-32;
1177 :     set_intra_matrix( get_default_intra_matrix() );
1178 : edgomez 1176 dequant_mpeg_intra(Dst, Src, 31, 5);
1179 :     printf( "dequant_mpeg_intra with CPU=%s: ", cpu->name);
1180 : edgomez 1060 printf( " Out[]= " );
1181 :     for(i=0; i<64; ++i) printf( "[%d]", Dst[i]);
1182 :     printf( "\n" );
1183 :     }
1184 : Isibaar 225
1185 : edgomez 1060 printf( "\n ===== (de)quant4_inter saturation bug? =====\n" );
1186 : Isibaar 225
1187 : edgomez 1176 for(cpu = cpu_list; cpu->name!=0; ++cpu)
1188 : edgomez 1060 {
1189 :     int i;
1190 :     int16_t Src[8*8], Dst[8*8];
1191 : Isibaar 225
1192 : edgomez 1060 if (!init_cpu(cpu))
1193 :     continue;
1194 : Isibaar 225
1195 : edgomez 1060 for(i=0; i<64; ++i) Src[i] = i-32;
1196 :     set_inter_matrix( get_default_inter_matrix() );
1197 : edgomez 1176 dequant_mpeg_inter(Dst, Src, 31);
1198 :     printf( "dequant_mpeg_inter with CPU=%s: ", cpu->name);
1199 : edgomez 1060 printf( " Out[]= " );
1200 :     for(i=0; i<64; ++i) printf( "[%d]", Dst[i]);
1201 :     printf( "\n" );
1202 :     }
1203 : Isibaar 225 }
1204 :    
1205 :     void test_dct_precision_diffs()
1206 :     {
1207 : edgomez 1060 CPU *cpu;
1208 :     short Blk[8*8], Blk0[8*8];
1209 : Isibaar 225
1210 : edgomez 1060 printf( "\n ===== fdct/idct precision diffs =====\n" );
1211 : Isibaar 225
1212 : edgomez 1176 for(cpu = cpu_list; cpu->name!=0; ++cpu)
1213 : edgomez 1060 {
1214 :     int i;
1215 : Isibaar 225
1216 : edgomez 1060 if (!init_cpu(cpu))
1217 :     continue;
1218 : Isibaar 225
1219 : edgomez 1060 for(i=0; i<8*8; ++i) {
1220 :     Blk0[i] = (i*7-i*i) & 0x7f;
1221 :     Blk[i] = Blk0[i];
1222 :     }
1223 : Isibaar 225
1224 : edgomez 1060 fdct(Blk);
1225 :     idct(Blk);
1226 :     printf( " fdct+idct diffs with CPU=%s: \n", cpu->name );
1227 :     for(i=0; i<8; ++i) {
1228 :     int j;
1229 :     for(j=0; j<8; ++j) printf( " %d ", Blk[i*8+j]-Blk0[i*8+j]);
1230 :     printf("\n");
1231 :     }
1232 :     printf("\n");
1233 :     }
1234 : Isibaar 225 }
1235 :    
1236 : Isibaar 262 void test_quant_bug()
1237 :     {
1238 : edgomez 1060 const int max_Q = 31;
1239 :     int i, n, qm, q;
1240 :     CPU *cpu;
1241 :     int16_t Src[8*8], Dst[8*8];
1242 :     uint8_t Quant[8*8];
1243 :     CPU cpu_bug_list[] = { { "PLAINC", 0 }, { "MMX ", XVID_CPU_MMX }, {0,0} };
1244 :     uint16_t Crcs_Inter[2][32];
1245 :     uint16_t Crcs_Intra[2][32];
1246 :     printf( "\n ===== test MPEG4-quantize bug =====\n" );
1247 : Isibaar 225
1248 : edgomez 1060 for(i=0; i<64; ++i) Src[i] = 2048*(i-32)/32;
1249 : Isibaar 262
1250 :     #if 1
1251 : edgomez 1060 for(qm=1; qm<=255; ++qm)
1252 :     {
1253 :     for(i=0; i<8*8; ++i) Quant[i] = qm;
1254 :     set_inter_matrix( Quant );
1255 : Isibaar 262
1256 : edgomez 1060 for(n=0, cpu = cpu_bug_list; cpu->name!=0; ++cpu, ++n)
1257 :     {
1258 :     uint16_t s;
1259 : Isibaar 262
1260 : edgomez 1060 if (!init_cpu(cpu))
1261 :     continue;
1262 : Isibaar 262
1263 : edgomez 1060 for(q=1; q<=max_Q; ++q) {
1264 :     emms();
1265 : edgomez 1176 quant_mpeg_inter( Dst, Src, q );
1266 : edgomez 1060 emms();
1267 :     for(s=0, i=0; i<64; ++i) s+=((uint16_t)Dst[i])^i;
1268 :     Crcs_Inter[n][q] = s;
1269 :     }
1270 :     }
1271 : Isibaar 262
1272 : edgomez 1060 for(q=1; q<=max_Q; ++q)
1273 :     for(i=0; i<n-1; ++i)
1274 :     if (Crcs_Inter[i][q]!=Crcs_Inter[i+1][q])
1275 :     printf( "Discrepancy Inter: qm=%d, q=%d -> %d/%d !\n",
1276 :     qm, q, Crcs_Inter[i][q], Crcs_Inter[i+1][q]);
1277 :     }
1278 : Isibaar 262 #endif
1279 :    
1280 :     #if 1
1281 : edgomez 1060 for(qm=1; qm<=255; ++qm)
1282 :     {
1283 :     for(i=0; i<8*8; ++i) Quant[i] = qm;
1284 :     set_intra_matrix( Quant );
1285 : Isibaar 262
1286 : edgomez 1060 for(n=0, cpu = cpu_bug_list; cpu->name!=0; ++cpu, ++n)
1287 :     {
1288 :     uint16_t s;
1289 : Isibaar 262
1290 : edgomez 1060 if (!init_cpu(cpu))
1291 :     continue;
1292 : Isibaar 262
1293 : edgomez 1060 for(q=1; q<=max_Q; ++q) {
1294 :     emms();
1295 : edgomez 1176 quant_mpeg_intra( Dst, Src, q, q);
1296 : edgomez 1060 emms();
1297 :     for(s=0, i=0; i<64; ++i) s+=((uint16_t)Dst[i])^i;
1298 :     Crcs_Intra[n][q] = s;
1299 :     }
1300 :     }
1301 : Isibaar 262
1302 : edgomez 1060 for(q=1; q<=max_Q; ++q)
1303 :     for(i=0; i<n-1; ++i)
1304 :     if (Crcs_Intra[i][q]!=Crcs_Intra[i+1][q])
1305 :     printf( "Discrepancy Intra: qm=%d, q=%d -> %d/%d!\n",
1306 :     qm, q, Crcs_Inter[i][q], Crcs_Inter[i+1][q]);
1307 :     }
1308 : Isibaar 262 #endif
1309 :     }
1310 :    
1311 : Isibaar 225 /*********************************************************************
1312 :     * main
1313 :     *********************************************************************/
1314 :    
1315 :     int main(int argc, char *argv[])
1316 :     {
1317 : edgomez 1060 int what = 0;
1318 :     if (argc>1) what = atoi(argv[1]);
1319 :     if (what==0 || what==1) test_dct();
1320 :     if (what==0 || what==2) test_mb();
1321 :     if (what==0 || what==3) test_sad();
1322 :     if (what==0 || what==4) test_transfer();
1323 :     if (what==0 || what==5) test_quant();
1324 :     if (what==0 || what==6) test_cbp();
1325 : Isibaar 225
1326 : edgomez 1060 if (what==7) {
1327 :     test_IEEE1180_compliance(-256, 255, 1);
1328 :     test_IEEE1180_compliance(-256, 255,-1);
1329 :     test_IEEE1180_compliance( -5, 5, 1);
1330 :     test_IEEE1180_compliance( -5, 5,-1);
1331 :     test_IEEE1180_compliance(-300, 300, 1);
1332 :     test_IEEE1180_compliance(-300, 300,-1);
1333 :     }
1334 :     if (what==8) test_dct_saturation(-256, 255);
1335 : Isibaar 262
1336 : edgomez 1060 if (what==9) {
1337 :     int width, height;
1338 :     if (argc<5) {
1339 :     printf("usage: %s %d [bitstream] [width] [height]\n", argv[0], what);
1340 :     return 1;
1341 :     }
1342 :     width = atoi(argv[3]);
1343 :     height = atoi(argv[4]);
1344 :     test_dec(argv[2], width, height, (argc>5));
1345 :     }
1346 : Isibaar 225
1347 : edgomez 1060 if (what==-1) {
1348 :     test_dct_precision_diffs();
1349 :     test_bugs1();
1350 :     }
1351 :     if (what==-2)
1352 :     test_quant_bug();
1353 : Isibaar 262
1354 : edgomez 1060 return 0;
1355 : Isibaar 225 }
1356 :    
1357 :     /*********************************************************************
1358 :     * 'Reference' output (except for timing) on a PIII 1.13Ghz/linux
1359 :     *********************************************************************/
1360 : Isibaar 262
1361 : edgomez 1060 /* as of 07/01/2002, there's a problem with mpeg4-quantization */
1362 : Isibaar 225 /*
1363 :    
1364 : edgomez 1060 ===== test fdct/idct =====
1365 : Isibaar 262 PLAINC - 3.312 usec PSNR=13.291 MSE=3.000
1366 :     MMX - 0.591 usec PSNR=13.291 MSE=3.000
1367 :     MMXEXT - 0.577 usec PSNR=13.291 MSE=3.000
1368 :     SSE2 - 0.588 usec PSNR=13.291 MSE=3.000
1369 : Isibaar 225 3DNOW - skipped...
1370 :     3DNOWE - skipped...
1371 :    
1372 : edgomez 1060 === test block motion ===
1373 : Isibaar 262 PLAINC - interp- h-round0 0.911 usec iCrc=8107
1374 :     PLAINC - round1 0.863 usec iCrc=8100
1375 :     PLAINC - interp- v-round0 0.860 usec iCrc=8108
1376 :     PLAINC - round1 0.857 usec iCrc=8105
1377 :     PLAINC - interp-hv-round0 2.103 usec iCrc=8112
1378 :     PLAINC - round1 2.050 usec iCrc=8103
1379 : edgomez 1060 ---
1380 : Isibaar 225 MMX - interp- h-round0 0.105 usec iCrc=8107
1381 : Isibaar 262 MMX - round1 0.106 usec iCrc=8100
1382 : Isibaar 225 MMX - interp- v-round0 0.106 usec iCrc=8108
1383 : Isibaar 262 MMX - round1 0.106 usec iCrc=8105
1384 : Isibaar 225 MMX - interp-hv-round0 0.145 usec iCrc=8112
1385 :     MMX - round1 0.145 usec iCrc=8103
1386 : edgomez 1060 ---
1387 : Isibaar 262 MMXEXT - interp- h-round0 0.028 usec iCrc=8107
1388 : Isibaar 225 MMXEXT - round1 0.041 usec iCrc=8100
1389 :     MMXEXT - interp- v-round0 0.027 usec iCrc=8108
1390 : Isibaar 262 MMXEXT - round1 0.041 usec iCrc=8105
1391 :     MMXEXT - interp-hv-round0 0.066 usec iCrc=8112
1392 :     MMXEXT - round1 0.065 usec iCrc=8103
1393 : edgomez 1060 ---
1394 : Isibaar 262 SSE2 - interp- h-round0 0.109 usec iCrc=8107
1395 : Isibaar 225 SSE2 - round1 0.105 usec iCrc=8100
1396 :     SSE2 - interp- v-round0 0.106 usec iCrc=8108
1397 : Isibaar 262 SSE2 - round1 0.109 usec iCrc=8105
1398 : Isibaar 225 SSE2 - interp-hv-round0 0.145 usec iCrc=8112
1399 :     SSE2 - round1 0.145 usec iCrc=8103
1400 : edgomez 1060 ---
1401 : Isibaar 225 3DNOW - skipped...
1402 :     3DNOWE - skipped...
1403 :    
1404 : edgomez 1060 ====== test SAD ======
1405 : Isibaar 262 PLAINC - sad8 0.251 usec sad=3776
1406 :     PLAINC - sad16 1.601 usec sad=27214
1407 :     PLAINC - sad16bi 2.371 usec sad=26274
1408 :     PLAINC - dev16 1.564 usec sad=3344
1409 : edgomez 1060 ---
1410 : Isibaar 225 MMX - sad8 0.057 usec sad=3776
1411 : Isibaar 262 MMX - sad16 0.182 usec sad=27214
1412 :     MMX - sad16bi 2.462 usec sad=26274
1413 :     MMX - dev16 0.311 usec sad=3344
1414 : edgomez 1060 ---
1415 : Isibaar 225 MMXEXT - sad8 0.036 usec sad=3776
1416 : Isibaar 262 MMXEXT - sad16 0.109 usec sad=27214
1417 :     MMXEXT - sad16bi 0.143 usec sad=26274
1418 :     MMXEXT - dev16 0.192 usec sad=3344
1419 : edgomez 1060 ---
1420 : Isibaar 225 SSE2 - sad8 0.057 usec sad=3776
1421 : Isibaar 262 SSE2 - sad16 0.179 usec sad=27214
1422 :     SSE2 - sad16bi 2.456 usec sad=26274
1423 :     SSE2 - dev16 0.321 usec sad=3344
1424 : edgomez 1060 ---
1425 : Isibaar 225 3DNOW - skipped...
1426 :     3DNOWE - skipped...
1427 :    
1428 : edgomez 1060 === test transfer ===
1429 : Isibaar 262 PLAINC - 8to16 0.151 usec crc=28288
1430 :     PLAINC - 16to8 1.113 usec crc=28288
1431 :     PLAINC - 8to8 0.043 usec crc=20352
1432 :     PLAINC - 16to8add 1.069 usec crc=25536
1433 :     PLAINC - 8to16sub 0.631 usec crc1=28064 crc2=16256
1434 :     PLAINC - 8to16sub2 0.597 usec crc=20384
1435 : edgomez 1060 ---
1436 : Isibaar 262 MMX - 8to16 0.032 usec crc=28288
1437 :     MMX - 16to8 0.024 usec crc=28288
1438 :     MMX - 8to8 0.020 usec crc=20352
1439 :     MMX - 16to8add 0.043 usec crc=25536
1440 :     MMX - 8to16sub 0.066 usec crc1=28064 crc2=16256
1441 :     MMX - 8to16sub2 0.111 usec crc=20384
1442 : edgomez 1060 ---
1443 : Isibaar 225
1444 : edgomez 1060 ===== test quant =====
1445 : Isibaar 262 PLAINC - quant4_intra 74.248 usec crc=29809
1446 :     PLAINC - quant4_inter 70.850 usec crc=12574
1447 :     PLAINC - dequant4_intra 40.628 usec crc=24052
1448 :     PLAINC - dequant4_inter 45.691 usec crc=63847
1449 :     PLAINC - quant_intra 43.357 usec crc=25662
1450 :     PLAINC - quant_inter 33.410 usec crc=23972
1451 :     PLAINC - dequant_intra 36.384 usec crc=49900
1452 :     PLAINC - dequant_inter 48.930 usec crc=48899
1453 : edgomez 1060 ---
1454 : Isibaar 262 MMX - quant4_intra 7.445 usec crc=3459
1455 :     *** CRC ERROR! ***
1456 :     MMX - quant4_inter 5.384 usec crc=51072
1457 :     *** CRC ERROR! ***
1458 :     MMX - dequant4_intra 5.515 usec crc=24052
1459 :     MMX - dequant4_inter 7.745 usec crc=63847
1460 :     MMX - quant_intra 4.661 usec crc=25662
1461 :     MMX - quant_inter 4.406 usec crc=23972
1462 :     MMX - dequant_intra 4.928 usec crc=49900
1463 :     MMX - dequant_inter 4.532 usec crc=48899
1464 : edgomez 1060 ---
1465 : Isibaar 225
1466 : edgomez 1060 ===== test cbp =====
1467 : Isibaar 262 PLAINC - calc_cbp#1 0.371 usec cbp=0x15
1468 :     PLAINC - calc_cbp#2 0.432 usec cbp=0x38
1469 :     PLAINC - calc_cbp#3 0.339 usec cbp=0xf
1470 :     PLAINC - calc_cbp#4 0.506 usec cbp=0x5
1471 : edgomez 1060 ---
1472 : Isibaar 225 MMX - calc_cbp#1 0.136 usec cbp=0x15
1473 : Isibaar 262 MMX - calc_cbp#2 0.134 usec cbp=0x38
1474 :     MMX - calc_cbp#3 0.138 usec cbp=0xf
1475 : Isibaar 225 MMX - calc_cbp#4 0.135 usec cbp=0x5
1476 : edgomez 1060 ---
1477 : Isibaar 262 SSE2 - calc_cbp#1 0.136 usec cbp=0x15
1478 :     SSE2 - calc_cbp#2 0.133 usec cbp=0x38
1479 :     SSE2 - calc_cbp#3 0.133 usec cbp=0xf
1480 :     SSE2 - calc_cbp#4 0.141 usec cbp=0x5
1481 : edgomez 1060 ---
1482 : Isibaar 262
1483 : Isibaar 225 */

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