[svn] / trunk / xvidcore / examples / xvid_bench.c Repository:
ViewVC logotype

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

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1398, Fri Apr 2 21:29:21 2004 UTC revision 1839, Tue Dec 2 13:44:55 2008 UTC
# Line 19  Line 19 
19   *  along with this program; if not, write to the Free Software   *  along with this program; if not, write to the Free Software
20   *  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA   *  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
21   *   *
22   * $Id: xvid_bench.c,v 1.12 2004-04-02 21:29:21 edgomez Exp $   * $Id: xvid_bench.c,v 1.39 2008-12-02 13:44:55 Isibaar Exp $
23   *   *
24   ****************************************************************************/   ****************************************************************************/
25    
26  /*****************************************************************************  /*****************************************************************************
27   *   *
28   *  'Reference' output is at the end of file.   *  'Reference' output is at the end of file.
  *  Don't take the checksums and crc too seriouly, they aren't  
  *  bullet-proof (should plug some .md5 here)...  
29   *   *
30   *   compiles with something like:   *   compiles with something like:
31   *   gcc -o xvid_bench xvid_bench.c  -I../src/ -lxvidcore -lm   *   gcc -o xvid_bench xvid_bench.c  -I../src/ -lxvidcore -lm
# Line 49  Line 47 
47  #include "xvid.h"  #include "xvid.h"
48    
49  // inner guts  // inner guts
50    #include "portab.h"
51  #include "dct/idct.h"  #include "dct/idct.h"
52  #include "dct/fdct.h"  #include "dct/fdct.h"
53  #include "image/colorspace.h"  #include "image/colorspace.h"
# Line 60  Line 59 
59  #include "utils/timer.h"  #include "utils/timer.h"
60  #include "quant/quant_matrix.c"  #include "quant/quant_matrix.c"
61  #include "bitstream/cbp.h"  #include "bitstream/cbp.h"
62    #include "bitstream/bitstream.h"
63    
64  #include <math.h>  #include <math.h>
65    
# Line 67  Line 67 
67  #define M_PI            3.14159265358979323846  #define M_PI            3.14159265358979323846
68  #endif  #endif
69    
70  const int speed_ref = 100;  /* on slow machines, decrease this value */  int speed_ref = 100;  /* on slow machines, decrease this value */
71    int verbose = 0;
72    unsigned int cpu_mask;
73    
74  /*********************************************************************  /*********************************************************************
75   * misc   * misc
# Line 83  Line 85 
85  #else  #else
86          clock_t clk;          clock_t clk;
87          clk = clock();          clk = clock();
88          return clk * 1000000 / CLOCKS_PER_SEC;          return clk * 1000. / CLOCKS_PER_SEC;  /* clock() returns time in Milliseconds */
89  #endif  #endif
90  }  }
91    
# Line 112  Line 114 
114          unsigned int cpu;          unsigned int cpu;
115  } CPU;  } CPU;
116    
117  CPU cpu_list[] =  CPU cpu_list[] = {
118  { { "PLAINC", 0 }          { "PLAINC ", 0 },
119  #ifdef ARCH_IS_IA32  #if defined(ARCH_IS_IA32) || defined(ARCH_IS_X86_64)
120    , { "MMX   ", XVID_CPU_MMX }          { "MMX    ", XVID_CPU_MMX },
121    , { "MMXEXT", XVID_CPU_MMXEXT | XVID_CPU_MMX }          { "MMXEXT ", XVID_CPU_MMXEXT | XVID_CPU_MMX },
122    , { "SSE2  ", XVID_CPU_SSE2 | XVID_CPU_MMX }          { "SSE2   ", XVID_CPU_SSE2 | XVID_CPU_MMX },
123    , { "3DNOW ", XVID_CPU_3DNOW }          { "SSE3   ", XVID_CPU_SSE3 | XVID_CPU_SSE2 | XVID_CPU_MMX },
124    , { "3DNOWE", XVID_CPU_3DNOW | XVID_CPU_3DNOWEXT }          { "SSE41  ", XVID_CPU_SSE41| XVID_CPU_SSE3 | XVID_CPU_SSE2 | XVID_CPU_MMX },
125            { "3DNOW  ", XVID_CPU_3DNOW },
126            { "3DNOWE ", XVID_CPU_3DNOW | XVID_CPU_3DNOWEXT },
127  #endif  #endif
128  //, { "IA64  ", XVID_CPU_IA64 }  #ifdef ARCH_IS_PPC
129  //, { "TSC   ", XVID_CPU_TSC }          { "ALTIVEC", XVID_CPU_ALTIVEC },
   , { 0, 0 } };  
   
 CPU  cpu_short_list[] =  
 { { "PLAINC", 0 }  
 #ifdef ARCH_IS_IA32  
   , { "MMX   ", XVID_CPU_MMX }  
 //, { "MMXEXT", XVID_CPU_MMXEXT | XVID_CPU_MMX }  
130  #endif  #endif
131  //, { "IA64  ", XVID_CPU_IA64 }  #ifdef ARCH_IS_IA64
132    , { 0, 0 } };  //      { "IA64   ", XVID_CPU_IA64 },
   
 CPU cpu_short_list2[] =  
 { { "PLAINC", 0 }  
 #ifdef ARCH_IS_IA32  
   , { "MMX   ", XVID_CPU_MMX }  
   , { "SSE2  ", XVID_CPU_SSE2 | XVID_CPU_MMX }  
133  #endif  #endif
134    , { 0, 0 } };  //      { "TSC    ", XVID_CPU_TSC },
135            { 0, 0 }
136    };
137    
138    
139  int init_cpu(CPU *cpu)  int init_cpu(CPU *cpu)
# Line 253  Line 246 
246  };  };
247    
248  uint32_t  uint32_t
249  calc_crc(uint8_t *mem, int len, uint32_t initial)  calc_crc(uint8_t *mem, int len, uint32_t crc)
250  {  {
   
         register unsigned int crc;  
   
         crc = initial;  
   
251          while( len >= 8) {          while( len >= 8) {
252                  DO8(mem, crc);                  DO8(mem, crc);
253                  len -= 8;                  len -= 8;
# Line 270  Line 258 
258                  len--;                  len--;
259          }          }
260    
261          return(crc);          return crc;
262    }
263    
264    void byte_swap(uint8_t *mem, int len, int element_size) {
265    #ifdef ARCH_IS_BIG_ENDIAN
266            int i;
267    
268            if(element_size == 1) {
269                    /* No need to swap */
270            } else if(element_size == 2) {
271                    uint8_t temp[2];
272    
273                    for(i=0; i < (len/2); i++ ) {
274                            temp[0] = mem[0];
275                            temp[1] = mem[1];
276                            mem[0] = temp[1];
277                            mem[1] = temp[0];
278    
279                            mem += 2;
280                    }
281            } else if(element_size == 4) {
282                    uint8_t temp[4];
283    
284                    for(i=0; i < (len/4); i++ ) {
285                            temp[0] = mem[0];
286                            temp[1] = mem[1];
287                            temp[2] = mem[2];
288                            temp[3] = mem[3];
289                            mem[0] = temp[3];
290                            mem[1] = temp[2];
291                            mem[2] = temp[1];
292                            mem[3] = temp[0];
293    
294                            mem += 4;
295                    }
296            } else {
297                    printf("ERROR: byte_swap unsupported element_size(%u)\n", element_size);
298            }
299    #endif
300  }  }
301    
302  /*********************************************************************  /*********************************************************************
# Line 389  Line 414 
414    
415                  t = gettime_usec();                  t = gettime_usec();
416                  emms();                  emms();
417                    for(tst=0; tst<nb_tests; ++tst) s = sad8bi(Cur, Ref1, Ref2, 8);
418                    emms();
419                    t = (gettime_usec() - t) / nb_tests;
420                    printf( "%s - sad8bi %.3f usec       sad=%d %s\n",
421                                    cpu->name, t, s,
422                                    (s!=4002)?"| ERROR": "" );
423    
424                    t = gettime_usec();
425                    emms();
426                  for(tst=0; tst<nb_tests; ++tst) s = dev16(Cur, 16);                  for(tst=0; tst<nb_tests; ++tst) s = dev16(Cur, 16);
427                  emms();                  emms();
428                  t = (gettime_usec() - t) / nb_tests;                  t = (gettime_usec() - t) / nb_tests;
# Line 494  Line 528 
528                             (iCrc!=8107)?"| ERROR": "" );                             (iCrc!=8107)?"| ERROR": "" );
529  #endif  #endif
530    
531        /* New functions for field prediction by CK 1.10.2005 */
532    #pragma NEW8X4
533                    TEST_MB(interpolate8x4_halfpel_h, 0);
534                    printf("%s - interpfield-h -round0 %.3f usec       crc32=0x%08x %s\n",
535                               cpu->name, t, iCrc,
536                               (iCrc!=0x9538d6df)?"| ERROR": "" );
537    
538                    TEST_MB(interpolate8x4_halfpel_h, 1);
539                    printf("%s -                round1 %.3f usec       crc32=0x%08x %s\n",
540                               cpu->name, t, iCrc,
541                               (iCrc!=0xde5f1db4)?"| ERROR": "" );
542    
543    
544                    TEST_MB(interpolate8x4_halfpel_v, 0);
545                    printf("%s - interpfield- v-round0 %.3f usec       crc32=0x%08x %s\n",
546                               cpu->name, t, iCrc,
547                               (iCrc!=0xea5a69ef)?"| ERROR": "" );
548    
549                    TEST_MB(interpolate8x4_halfpel_v, 1);
550                    printf("%s -                round1 %.3f usec       crc32=0x%08x %s\n",
551                               cpu->name, t, iCrc,
552                               (iCrc!=0x4f10ec0f)?"| ERROR": "" );
553    
554    
555                    TEST_MB(interpolate8x4_halfpel_hv, 0);
556                    printf("%s - interpfield-hv-round0 %.3f usec       crc32=0x%08x %s\n",
557                               cpu->name, t, iCrc,
558                               (iCrc!=0xf97ee367)?"| ERROR": "" );
559    
560                    TEST_MB(interpolate8x4_halfpel_hv, 1);
561                    printf("%s -                round1 %.3f usec       crc32=0x%08x %s\n",
562                               cpu->name, t, iCrc,
563                               (iCrc!=0xb6a9f581)?"| ERROR": "" );
564    /* End of 8x4 functions */
565    
566                  printf( " --- \n" );                  printf( " --- \n" );
567          }          }
568  }  }
569    
570    #undef ENTER
571    #undef LEAVE
572    #undef TEST_MB
573    #undef TEST_MB2
574    
575  /*********************************************************************  /*********************************************************************
576   * test transfer   * test transfer
577   *********************************************************************/   *********************************************************************/
# Line 527  Line 601 
601  }                                         \  }                                         \
602  emms();                                   \  emms();                                   \
603  t = (gettime_usec()-t -overhead) / nb_tests;\  t = (gettime_usec()-t -overhead) / nb_tests;\
604  s = calc_crc((uint8_t*)(DST), sizeof((DST)), CRC32_INITIAL)  byte_swap((uint8_t*)(DST), 8*32*sizeof((DST)[0]), sizeof((DST)[0]));  \
605    s = calc_crc((uint8_t*)(DST), 8*32*sizeof((DST)[0]), CRC32_INITIAL)
606    
607  #define TEST_TRANSFER(FUNC, DST, SRC)         \  #define TEST_TRANSFER(FUNC, DST, SRC)         \
608  TEST_TRANSFER_BEGIN(DST);                 \  TEST_TRANSFER_BEGIN(DST);                 \
# Line 553  Line 628 
628  }                                         \  }                                         \
629  emms();                                   \  emms();                                   \
630  t = (gettime_usec()-t -overhead) / nb_tests;\  t = (gettime_usec()-t -overhead) / nb_tests;\
631  s = calc_crc((uint8_t*)(DST), sizeof((DST)), CRC32_INITIAL)  byte_swap((uint8_t*)(DST), 8*32*sizeof((DST)[0]), sizeof((DST)[0]));  \
632    s = calc_crc((uint8_t*)(DST), 8*32*sizeof((DST)[0]), CRC32_INITIAL)
633    
634  #define TEST_TRANSFER2(FUNC, DST, SRC, R1)    \  #define TEST_TRANSFER2(FUNC, DST, SRC, R1)    \
635  TEST_TRANSFER2_BEGIN(DST,SRC);            \  TEST_TRANSFER2_BEGIN(DST,SRC);            \
# Line 570  Line 646 
646          const int nb_tests = 4000*speed_ref;          const int nb_tests = 4000*speed_ref;
647          int i;          int i;
648          CPU *cpu;          CPU *cpu;
649          uint8_t  Src8[8*32], Dst8[8*32], Ref1[8*32], Ref2[8*32];  //      uint8_t  Src8[8*32], Dst8[8*32], Ref1[8*32], Ref2[8*32];
650          int16_t Src16[8*32], Dst16[8*32];  //      int16_t Src16[8*32], Dst16[8*32];
651      DECLARE_ALIGNED_MATRIX(Src8, 8, 32, uint8_t, CACHE_LINE);
652      DECLARE_ALIGNED_MATRIX(Dst8, 8, 32, uint8_t, CACHE_LINE);
653      DECLARE_ALIGNED_MATRIX(Ref1, 8, 32, uint8_t, CACHE_LINE);
654      DECLARE_ALIGNED_MATRIX(Ref2, 8, 32, uint8_t, CACHE_LINE);
655      DECLARE_ALIGNED_MATRIX(Src16, 8, 32, uint16_t, CACHE_LINE);
656      DECLARE_ALIGNED_MATRIX(Dst16, 8, 32, uint16_t, CACHE_LINE);
657    
658          printf( "\n ===  test transfer ===\n" );          printf( "\n ===  test transfer ===\n" );
659    
# Line 593  Line 675 
675                                  cpu->name, t, s,                                  cpu->name, t, s,
676                                  (s!=0xee7ccbb4)?"| ERROR": "");                                  (s!=0xee7ccbb4)?"| ERROR": "");
677    
678        /* New functions for field prediction by CK 1.10.2005 */
679    #pragma NEW8X4
680                    TEST_TRANSFER(transfer8x4_copy, Dst8, Src8);
681                    printf("%s - 8to4      %.3f usec       crc32=0x%08x %s\n",
682                               cpu->name, t, s,
683                               (s!=0xbb9c3db5)?"| ERROR": "");
684    /* End of new functions */
685    
686                  TEST_TRANSFER(transfer8x8_copy, Dst8, Src8);                  TEST_TRANSFER(transfer8x8_copy, Dst8, Src8);
687                  printf("%s - 8to8      %.3f usec       crc32=0x%08x %s\n",                  printf("%s - 8to8      %.3f usec       crc32=0x%08x %s\n",
688                             cpu->name, t, s,                             cpu->name, t, s,
# Line 606  Line 696 
696                  TEST_TRANSFER2(transfer_8to16sub, Dst16, Src8, Ref1);                  TEST_TRANSFER2(transfer_8to16sub, Dst16, Src8, Ref1);
697                  {                  {
698                          int s1, s2;                          int s1, s2;
699                          s1 = calc_crc((uint8_t*)Dst16, sizeof(Dst16), CRC32_INITIAL);                          s1 = calc_crc((uint8_t*)Dst16, 8*32*sizeof(Dst16[0]), CRC32_INITIAL);
700                          s2 = calc_crc((uint8_t*)Src8, sizeof(Src8), CRC32_INITIAL);                          s2 = calc_crc((uint8_t*)Src8, 8*32*sizeof(Src8[0]), CRC32_INITIAL);
701                          printf("%s - 8to16sub  %.3f usec       crc32(1)=0x%08x crc32(2)=0x%08x %s %s\n",                          printf("%s - 8to16sub  %.3f usec       crc32(1)=0x%08x crc32(2)=0x%08x %s %s\n",
702                                     cpu->name, t, s1, s2,                                     cpu->name, t, s1, s2,
703                                     (s1!=0xa1e07163)?"| ERROR1": "",                                     (s1!=0xa1e07163)?"| ERROR1": "",
# Line 636  Line 726 
726    for(q=1; q<=max_Q; ++q) {                 \    for(q=1; q<=max_Q; ++q) {                 \
727          for(tst=0; tst<nb_tests; ++tst)         \          for(tst=0; tst<nb_tests; ++tst)         \
728            (FUNC)((DST), (SRC), q, mpeg_quant_matrices);              \            (FUNC)((DST), (SRC), q, mpeg_quant_matrices);              \
729          s = calc_crc((uint8_t*)(DST), 64*sizeof(int16_t), s); \          byte_swap((uint8_t*)(DST), 64*sizeof((DST)[0]), sizeof((DST)[0]));  \
730            s = calc_crc((uint8_t*)(DST), 64*sizeof((DST)[0]), s); \
731    }                                         \    }                                         \
732    emms();                                   \    emms();                                   \
733  }                                           \  }                                           \
# Line 649  Line 740 
740    set_intra_matrix( mpeg_quant_matrices, Quant );                \    set_intra_matrix( mpeg_quant_matrices, Quant );                \
741    emms();                                   \    emms();                                   \
742    for(q=1; q<=max_Q; ++q) {                 \    for(q=1; q<=max_Q; ++q) {                 \
743            init_intra_matrix( mpeg_quant_matrices, q ); \
744          for(tst=0; tst<nb_tests; ++tst)         \          for(tst=0; tst<nb_tests; ++tst)         \
745            (FUNC)((DST), (SRC), q, q, mpeg_quant_matrices);           \            (FUNC)((DST), (SRC), q, q, mpeg_quant_matrices);           \
746          s = calc_crc((uint8_t*)(DST), 64*sizeof(int16_t), s); \          byte_swap((uint8_t*)(DST), 64*sizeof((DST)[0]), sizeof((DST)[0]));  \
747            s = calc_crc((uint8_t*)(DST), 64*sizeof((DST)[0]), s); \
748    }                                         \    }                                         \
749    emms();                                   \    emms();                                   \
750  }                                           \  }                                           \
751  t = (gettime_usec()-t-overhead)/nb_tests/qm  t = (gettime_usec()-t-overhead)/nb_tests/qm
752    
753  #define TEST_INTRA(REFFUNC, NEWFUNC, RANGE)              \  #define TEST_INTRA(REFFUNC, NEWFUNC, RANGE)              \
754  { int i,q,s;\  { int32_t i,q,s;\
755          DECLARE_ALIGNED_MATRIX(Src, 8, 8, int16_t, 16); \          DECLARE_ALIGNED_MATRIX(Src, 8, 8, int16_t, 16); \
756    DECLARE_ALIGNED_MATRIX(Dst, 8, 8, int16_t, 16); \    DECLARE_ALIGNED_MATRIX(Dst, 8, 8, int16_t, 16); \
757    DECLARE_ALIGNED_MATRIX(Dst2,8, 8, int16_t, 16); \    DECLARE_ALIGNED_MATRIX(Dst2,8, 8, int16_t, 16); \
# Line 690  Line 783 
783    
784  void test_quant()  void test_quant()
785  {  {
786          const int nb_tests = 1*speed_ref;          const int32_t nb_tests = 1*speed_ref;
787          const int max_Q = 31;          const int32_t max_Q = 31;
788          DECLARE_ALIGNED_MATRIX(mpeg_quant_matrices, 8, 64, uint16_t, 16);          DECLARE_ALIGNED_MATRIX(mpeg_quant_matrices, 8, 64, uint16_t, 16);
789    
790          int i, qm;          int32_t i, qm;
791          CPU *cpu;          CPU *cpu;
792          DECLARE_ALIGNED_MATRIX(Src, 8, 8, int16_t, 16);          DECLARE_ALIGNED_MATRIX(Src, 8, 8, int16_t, 16);
793          DECLARE_ALIGNED_MATRIX(Dst, 8, 8, int16_t, 16);          DECLARE_ALIGNED_MATRIX(Dst, 8, 8, int16_t, 16);
# Line 713  Line 806 
806          for(cpu = cpu_list; cpu->name!=0; ++cpu)          for(cpu = cpu_list; cpu->name!=0; ++cpu)
807          {          {
808                  double t, overhead;                  double t, overhead;
809                  int tst, q;                  int32_t tst, q;
810                  uint32_t s;                  uint32_t s;
811    
812                  if (!init_cpu(cpu))                  if (!init_cpu(cpu))
# Line 737  Line 830 
830                  TEST_QUANT2(quant_mpeg_intra, Dst, Src);                  TEST_QUANT2(quant_mpeg_intra, Dst, Src);
831                  printf("%s -   quant_mpeg_intra %.3f usec       crc32=0x%08x %s\n",                  printf("%s -   quant_mpeg_intra %.3f usec       crc32=0x%08x %s\n",
832                             cpu->name, t, s,                             cpu->name, t, s,
833                             (s!=0xfd6a21a4)? "| ERROR": "");                             (s!=0x3b999af6)? "| ERROR": "");
834    
835                  TEST_QUANT(quant_mpeg_inter, Dst, Src);                  TEST_QUANT(quant_mpeg_inter, Dst, Src);
836                  printf("%s -   quant_mpeg_inter %.3f usec       crc32=0x%08x %s\n",                  printf("%s -   quant_mpeg_inter %.3f usec       crc32=0x%08x %s\n",
# Line 779  Line 872 
872  }  }
873    
874  /*********************************************************************  /*********************************************************************
875     * test distortion operators
876     *********************************************************************/
877    
878    static void ieee_reseed(long s);
879    static long ieee_rand(int Min, int Max);
880    
881    #define TEST_SSE(FUNCTION, SRC1, SRC2, STRIDE) \
882      do { \
883        t = gettime_usec(); \
884        tst = nb_tests; \
885        while((tst--)>0) sse = (FUNCTION)((SRC1), (SRC2), (STRIDE)); \
886        emms(); \
887        t = (gettime_usec() - t)/(double)nb_tests;  \
888      } while(0)
889    
890    
891    void test_sse()
892    {
893            const int nb_tests = 100000*speed_ref;
894            int i;
895            CPU *cpu;
896            DECLARE_ALIGNED_MATRIX(Src1, 8, 8, int16_t, 16);
897            DECLARE_ALIGNED_MATRIX(Src2, 8, 8, int16_t, 16);
898            DECLARE_ALIGNED_MATRIX(Src3, 8, 8, int16_t, 16);
899            DECLARE_ALIGNED_MATRIX(Src4, 8, 8, int16_t, 16);
900    
901            printf( "\n =====  test sse =====\n" );
902    
903            ieee_reseed(1);
904            for(i=0; i<64; ++i) {
905                    Src1[i] = ieee_rand(-2048, 2047);
906                    Src2[i] = ieee_rand(-2048, 2047);
907                    Src3[i] = ieee_rand(-2048, 2047);
908                    Src4[i] = ieee_rand(-2048, 2047);
909            }
910    
911            for(cpu = cpu_list; cpu->name!=0; ++cpu)
912            {
913                    double t;
914                    int tst, sse;
915    
916                    if (!init_cpu(cpu))
917                            continue;
918    
919                    /* 16 bit element blocks */
920                    TEST_SSE(sse8_16bit, Src1, Src2, 16);
921                    printf("%s -   sse8_16bit#1 %.3f usec       sse=%d %s\n",
922                               cpu->name, t, sse, (sse!=182013834)?"| ERROR": "");
923                    TEST_SSE(sse8_16bit, Src1, Src3, 16);
924                    printf("%s -   sse8_16bit#2 %.3f usec       sse=%d %s\n",
925                               cpu->name, t, sse, (sse!=142545203)?"| ERROR": "");
926                    TEST_SSE(sse8_16bit, Src1, Src4, 16);
927                    printf("%s -   sse8_16bit#3 %.3f usec       sse=%d %s\n",
928                               cpu->name, t, sse, (sse!=146340935)?"| ERROR": "");
929                    TEST_SSE(sse8_16bit, Src2, Src3, 16);
930                    printf("%s -   sse8_16bit#4 %.3f usec       sse=%d %s\n",
931                               cpu->name, t, sse, (sse!=130136661)?"| ERROR": "");
932                    TEST_SSE(sse8_16bit, Src2, Src4, 16);
933                    printf("%s -   sse8_16bit#5 %.3f usec       sse=%d %s\n",
934                               cpu->name, t, sse, (sse!=136870353)?"| ERROR": "");
935                    TEST_SSE(sse8_16bit, Src3, Src4, 16);
936                    printf("%s -   sse8_16bit#6 %.3f usec       sse=%d %s\n",
937                               cpu->name, t, sse, (sse!=164107772)?"| ERROR": "");
938    
939                    /* 8 bit element blocks */
940                    TEST_SSE(sse8_8bit, (int8_t*)Src1, (int8_t*)Src2, 8);
941                    printf("%s -    sse8_8bit#1 %.3f usec       sse=%d %s\n",
942                               cpu->name, t, sse, (sse!=1356423)?"| ERROR": "");
943                    TEST_SSE(sse8_8bit, (int8_t*)Src1, (int8_t*)Src3, 8);
944                    printf("%s -    sse8_8bit#2 %.3f usec       sse=%d %s\n",
945                               cpu->name, t, sse, (sse!=1173074)?"| ERROR": "");
946                    TEST_SSE(sse8_8bit, (int8_t*)Src1, (int8_t*)Src4, 8);
947                    printf("%s -    sse8_8bit#3 %.3f usec       sse=%d %s\n",
948                               cpu->name, t, sse, (sse!=1092357)?"| ERROR": "");
949                    TEST_SSE(sse8_8bit, (int8_t*)Src2, (int8_t*)Src3, 8);
950                    printf("%s -    sse8_8bit#4 %.3f usec       sse=%d %s\n",
951                               cpu->name, t, sse, (sse!=1360239)?"| ERROR": "");
952                    TEST_SSE(sse8_8bit, (int8_t*)Src2, (int8_t*)Src4, 8);
953                    printf("%s -    sse8_8bit#5 %.3f usec       sse=%d %s\n",
954                               cpu->name, t, sse, (sse!=1208414)?"| ERROR": "");
955                    TEST_SSE(sse8_8bit, (int8_t*)Src3, (int8_t*)Src4, 8);
956                    printf("%s -    sse8_8bit#6 %.3f usec       sse=%d %s\n",
957                               cpu->name, t, sse, (sse!=1099285)?"| ERROR": "");
958    
959                    printf(" ---\n");
960            }
961    }
962    
963    /*********************************************************************
964   * test non-zero AC counting   * test non-zero AC counting
965   *********************************************************************/   *********************************************************************/
966    
967  #define TEST_CBP(FUNC, SRC)               \  #define TEST_CBP(FUNC, SRC, NB)           \
968  t = gettime_usec();                       \  t = gettime_usec();                       \
969  emms();                                   \  emms();                                   \
970  for(tst=0; tst<nb_tests; ++tst) {         \  for(tst=0; tst<NB; ++tst) {         \
971    cbp = (FUNC)((SRC));                    \    cbp = (FUNC)((SRC));                    \
972  }                                         \  }                                         \
973  emms();                                   \  emms();                                   \
# Line 794  Line 976 
976  void test_cbp()  void test_cbp()
977  {  {
978          const int nb_tests = 10000*speed_ref;          const int nb_tests = 10000*speed_ref;
979          int i;          int i, n, m;
980          CPU *cpu;          CPU *cpu;
981          DECLARE_ALIGNED_MATRIX(Src1, 6, 64, int16_t, 16);          DECLARE_ALIGNED_MATRIX(Src1, 6, 64, int16_t, 16);
982          DECLARE_ALIGNED_MATRIX(Src2, 6, 64, int16_t, 16);          DECLARE_ALIGNED_MATRIX(Src2, 6, 64, int16_t, 16);
983          DECLARE_ALIGNED_MATRIX(Src3, 6, 64, int16_t, 16);          DECLARE_ALIGNED_MATRIX(Src3, 6, 64, int16_t, 16);
984          DECLARE_ALIGNED_MATRIX(Src4, 6, 64, int16_t, 16);          DECLARE_ALIGNED_MATRIX(Src4, 6, 64, int16_t, 16);
985      DECLARE_ALIGNED_MATRIX(Src5, 6, 64, int16_t, 16);
986    
987          printf( "\n =====  test cbp =====\n" );          printf( "\n =====  test cbp =====\n" );
988    
# Line 808  Line 991 
991                  Src2[i] = (i<3*64);               /* half-full */                  Src2[i] = (i<3*64);               /* half-full */
992                  Src3[i] = ((i+32)>3*64);                  Src3[i] = ((i+32)>3*64);
993                  Src4[i] = (i==(3*64+2) || i==(5*64+9));                  Src4[i] = (i==(3*64+2) || i==(5*64+9));
994        Src5[i] = ieee_rand(0,1) ? -1 : 1;  /* +/- test */
995          }          }
996    
997          for(cpu = cpu_list; cpu->name!=0; ++cpu)          for(cpu = cpu_list; cpu->name!=0; ++cpu)
# Line 818  Line 1002 
1002                  if (!init_cpu(cpu))                  if (!init_cpu(cpu))
1003                          continue;                          continue;
1004    
1005                  TEST_CBP(calc_cbp, Src1);                  TEST_CBP(calc_cbp, Src1, nb_tests);
1006                  printf("%s -   calc_cbp#1 %.3f usec       cbp=0x%02x\n",                  printf("%s -   calc_cbp#1 %.3f usec       cbp=0x%02x %s\n",
1007                             cpu->name, t, cbp, (cbp!=0x15)?"| ERROR": "");                             cpu->name, t, cbp, (cbp!=0x15)?"| ERROR": "");
1008                  TEST_CBP(calc_cbp, Src2);                  TEST_CBP(calc_cbp, Src2, nb_tests);
1009                  printf("%s -   calc_cbp#2 %.3f usec       cbp=0x%02x\n",                  printf("%s -   calc_cbp#2 %.3f usec       cbp=0x%02x %s\n",
1010                             cpu->name, t, cbp, (cbp!=0x38)?"| ERROR": "");                             cpu->name, t, cbp, (cbp!=0x38)?"| ERROR": "");
1011                  TEST_CBP(calc_cbp, Src3);                  TEST_CBP(calc_cbp, Src3, nb_tests);
1012                  printf("%s -   calc_cbp#3 %.3f usec       cbp=0x%02x\n",                  printf("%s -   calc_cbp#3 %.3f usec       cbp=0x%02x %s\n",
1013                             cpu->name, t, cbp, (cbp!=0x0f)?"| ERROR": "" );                             cpu->name, t, cbp, (cbp!=0x0f)?"| ERROR": "" );
1014                  TEST_CBP(calc_cbp, Src4);                  TEST_CBP(calc_cbp, Src4, nb_tests);
1015                  printf("%s -   calc_cbp#4 %.3f usec       cbp=0x%02x\n",                  printf("%s -   calc_cbp#4 %.3f usec       cbp=0x%02x %s\n",
1016                             cpu->name, t, cbp, (cbp!=0x05)?"| ERROR": "" );                             cpu->name, t, cbp, (cbp!=0x05)?"| ERROR": "" );
1017                    TEST_CBP(calc_cbp, Src5, nb_tests);
1018                    printf("%s -   calc_cbp#4 %.3f usec       cbp=0x%02x %s\n",
1019                               cpu->name, t, cbp, (cbp!=0x3f)?"| ERROR": "" );
1020                  printf( " --- \n" );                  printf( " --- \n" );
1021          }          }
1022    
1023            for(cpu = cpu_list; cpu->name!=0; ++cpu)  /* bench suggested by Carlo (carlo dot bramix at libero dot it) */
1024            {
1025                    double t;
1026                    int tst, cbp, err;
1027    
1028                    if (!init_cpu(cpu))
1029                            continue;
1030    
1031        err = 0;
1032        for(n=0; n<6; ++n)
1033        {
1034          for(m=0; m<64; ++m)
1035          {
1036            for(i=0; i<6*64; ++i)
1037              Src1[i] = (i== (m + n*64));
1038    
1039            TEST_CBP(calc_cbp, Src1, 1);
1040            if (cbp!= (((m!=0)<<(5-n))))
1041            {
1042              printf( "%s -   calc_cbp#5: ERROR at pos %d / %d!\n", cpu->name, n, m);
1043              err = 1;
1044              break;
1045            }
1046          }
1047        }
1048        if (!err)
1049          printf( " %s -    calc_cbp#5 : OK\n", cpu->name );
1050    
1051            }
1052  }  }
1053    
1054  /*********************************************************************  /*********************************************************************
# Line 1225  Line 1442 
1442   * measure raw decoding speed   * measure raw decoding speed
1443   *********************************************************************/   *********************************************************************/
1444    
1445  void test_dec(const char *name, int width, int height, int with_chksum)  void test_dec(const char *name, int width, int height, int ref_chksum)
1446  {  {
1447          FILE *f = 0;          FILE *f = 0;
1448          void *dechandle = 0;          void *dechandle = 0;
# Line 1236  Line 1453 
1453          double t = 0.;          double t = 0.;
1454          int nb = 0;          int nb = 0;
1455          uint8_t *buf = 0;          uint8_t *buf = 0;
1456          uint8_t *rgb_out = 0;          uint8_t *yuv_out = 0;
1457          int buf_size, pos;          int buf_size, pos;
1458          uint32_t chksum = 0;          uint32_t chksum = 0;
1459            int bps = (width+31) & ~31;
1460    
1461          memset(&xinit, 0, sizeof(xinit));          memset(&xinit, 0, sizeof(xinit));
1462          xinit.cpu_flags = XVID_CPU_MMX | XVID_CPU_FORCE;          xinit.cpu_flags = cpu_mask;
1463          xinit.version = XVID_VERSION;          xinit.version = XVID_VERSION;
1464          xvid_global(NULL, 0, &xinit, NULL);          xvid_global(NULL, 0, &xinit, NULL);
1465    
# Line 1251  Line 1469 
1469          xparam.version = XVID_VERSION;          xparam.version = XVID_VERSION;
1470          xerr = xvid_decore(NULL, XVID_DEC_CREATE, &xparam, NULL);          xerr = xvid_decore(NULL, XVID_DEC_CREATE, &xparam, NULL);
1471          if (xerr==XVID_ERR_FAIL) {          if (xerr==XVID_ERR_FAIL) {
1472                  printf("can't init decoder (err=%d)\n", xerr);                  printf("ERROR: can't init decoder (err=%d)\n", xerr);
1473                  return;                  return;
1474          }          }
1475          dechandle = xparam.handle;          dechandle = xparam.handle;
# Line 1259  Line 1477 
1477    
1478          f = fopen(name, "rb");          f = fopen(name, "rb");
1479          if (f==0) {          if (f==0) {
1480                  printf( "can't open file '%s'\n", name);                  printf( "ERROR: can't open file '%s'\n", name);
1481                  return;                  return;
1482          }          }
1483          fseek(f, 0, SEEK_END);          fseek(f, 0, SEEK_END);
1484          buf_size = ftell(f);          buf_size = ftell(f);
1485          fseek(f, 0, SEEK_SET);          fseek(f, 0, SEEK_SET);
1486          if (buf_size<=0) {          if (buf_size<=0) {
1487                  printf("error while stating file\n");                  printf("ERROR: error while stating file\n");
1488                  goto End;                  goto End;
1489          }          }
         else printf( "Input size: %d\n", buf_size);  
1490    
1491          buf = malloc(buf_size); /* should be enuf' */          buf = malloc(buf_size);
1492          rgb_out = calloc(4, width*height);  /* <-room for _RGB24 */          yuv_out = calloc(1, bps*height*3/2 + 15);
1493          if (buf==0 || rgb_out==0) {          if (buf==0 || yuv_out==0) {
1494                  printf( "malloc failed!\n" );                  printf( "ERROR: malloc failed!\n" );
1495                  goto End;                  goto End;
1496          }          }
1497    
1498          if (fread(buf, buf_size, 1, f)!=1) {          if (fread(buf, buf_size, 1, f)!=1) {
1499                  printf( "file-read failed\n" );                  printf( "ERROR: file-read failed\n" );
1500                  goto End;                  goto End;
1501          }          }
1502    
# Line 1287  Line 1504 
1504          pos = 0;          pos = 0;
1505          t = -gettime_usec();          t = -gettime_usec();
1506          while(1) {          while(1) {
1507              int y;
1508    
1509                  memset(&xframe, 0, sizeof(xframe));                  memset(&xframe, 0, sizeof(xframe));
1510                  xframe.version = XVID_VERSION;                  xframe.version = XVID_VERSION;
1511                  xframe.bitstream = buf + pos;                  xframe.bitstream = buf + pos;
1512                  xframe.length = buf_size - pos;                  xframe.length = buf_size - pos;
1513                  xframe.output.plane[0] = rgb_out;                  xframe.output.plane[0] = (uint8_t*)(((size_t)yuv_out + 15) & ~15);
1514                  xframe.output.stride[0] = width;                  xframe.output.plane[1] = (uint8_t*)xframe.output.plane[0] + bps*height;
1515                  xframe.output.csp = XVID_CSP_BGR;                  xframe.output.plane[2] = (uint8_t*)xframe.output.plane[1] + bps/2;
1516                    xframe.output.stride[0] = bps;
1517                    xframe.output.stride[1] = bps;
1518                    xframe.output.stride[2] = bps;
1519                    xframe.output.csp = XVID_CSP_I420;
1520                  xerr = xvid_decore(dechandle, XVID_DEC_DECODE, &xframe, 0);                  xerr = xvid_decore(dechandle, XVID_DEC_DECODE, &xframe, 0);
1521                    if (xerr<0) {
1522                            printf("ERROR: decoding failed for frame #%d (err=%d)!\n", nb, xerr);
1523                            break;
1524                    }
1525                    else if (xerr==0)
1526                      break;
1527        else if (verbose>0) printf("#%d %d\n", nb, xerr );
1528    
1529                    pos += xerr;
1530                  nb++;                  nb++;
1531                  pos += xframe.length;  
1532                  if (with_chksum) {      for(y=0; y<height/2; ++y) {
1533                          int k = width*height;                    chksum = calc_crc((uint8_t*)xframe.output.plane[0] + (2*y+0)*bps, width, chksum);
1534                          uint32_t *ptr = (uint32_t *)rgb_out;                          chksum = calc_crc((uint8_t*)xframe.output.plane[0] + (2*y+1)*bps, width, chksum);
1535                          while(k-->0) chksum += *ptr++;                          chksum = calc_crc((uint8_t*)xframe.output.plane[1] + y*bps, width/2, chksum);
1536                            chksum = calc_crc((uint8_t*)xframe.output.plane[2] + y*bps, width/2, chksum);
1537                  }                  }
1538                  if (pos==buf_size)                  if (pos==buf_size)
1539                          break;                          break;
                 if (xerr==XVID_ERR_FAIL) {  
                         printf("decoding failed for frame #%d (err=%d)!\n", nb, xerr);  
                         break;  
                 }  
1540          }          }
1541          t += gettime_usec();          t += gettime_usec();
1542            if (ref_chksum==0) {
1543          if (t>0.)          if (t>0.)
1544                  printf( "%d frames decoded in %.3f s -> %.1f FPS\n", nb, t*1.e-6f, (float)(nb*1.e6f/t) );                    printf( "%d frames decoded in %.3f s -> %.1f FPS   Checksum:0x%.8x\n", nb, t*1.e-6f, (float)(nb*1.e6f/t), chksum );
1545          if (with_chksum)    }
1546                  printf("checksum: 0x%.8x\n", chksum);    else {
1547                    printf("FPS:%.1f Checksum: 0x%.8x Expected:0x%.8x | %s\n",
1548                      t>0. ? (float)(nb*1.e6f/t) : 0.f, chksum, ref_chksum, (chksum==ref_chksum) ? "OK" : "ERROR");
1549      }
1550    
1551   End:   End:
1552          if (rgb_out!=0) free(rgb_out);          if (yuv_out!=0) free(yuv_out);
1553          if (buf!=0) free(buf);          if (buf!=0) free(buf);
1554          if (dechandle!=0) {          if (dechandle!=0) {
1555                  xerr= xvid_decore(dechandle, XVID_DEC_DESTROY, NULL, NULL);                  xerr= xvid_decore(dechandle, XVID_DEC_DESTROY, NULL, NULL);
1556                  if (xerr==XVID_ERR_FAIL)                  if (xerr==XVID_ERR_FAIL)
1557                          printf("destroy-decoder failed (err=%d)!\n", xerr);                          printf("ERROR: destroy-decoder failed (err=%d)!\n", xerr);
1558          }          }
1559          if (f!=0) fclose(f);          if (f!=0) fclose(f);
1560  }  }
# Line 1485  Line 1718 
1718  }  }
1719    
1720  /*********************************************************************  /*********************************************************************
1721     * test some YUV func
1722     *********************************************************************/
1723    
1724    #define ENTER \
1725    for(i=0; i<(int)sizeof(Dst0); ++i) Dst0[0][i] = 0;   \
1726    t = gettime_usec();                   \
1727    emms();
1728    
1729    #define LEAVE \
1730    emms();                             \
1731    t = (gettime_usec() - t) / nb_tests;  \
1732            iCrc = calc_crc((uint8_t*)Dst0, sizeof(Dst0), CRC32_INITIAL)
1733    
1734    #define TEST_YUYV(FUNC, S, FLIP)                \
1735    ENTER                               \
1736    for(tst=0; tst<nb_tests; ++tst) (FUNC)(Dst0[0], S*WIDTH, Src0[0], Src0[1], Src0[2], WIDTH, WIDTH/2, WIDTH, HEIGHT, (FLIP)); \
1737    LEAVE
1738    
1739    static const int yuv_CRCs[6][2] = {
1740            {0x0f4fb96b,0x780b6a68}
1741    ,       {0xa986b289,0x65e49b76}
1742    ,       {0x7f19c152,0xd539b86e}
1743    ,       {0x0f4fb96b,0x780b6a68}
1744    ,       {0xa986b289,0x65e49b76}
1745    ,       {0x36ab8b57,0x1cd92fee}
1746    };
1747    
1748    #define WIDTH 128
1749    #define HEIGHT 32
1750    void test_yuv()
1751    {
1752            const int nb_tests = 200*speed_ref;
1753            CPU *cpu;
1754            uint8_t Src0[3][WIDTH*HEIGHT];
1755            uint8_t Dst0[4][WIDTH*HEIGHT];
1756            int i, j, with_flip;
1757            double t;
1758            int tst, iCrc;
1759    
1760            colorspace_init();
1761            ieee_reseed(1);
1762            for(i=0; i<(int)sizeof(Src0); ++i) Src0[0][i] = ieee_rand(0,255);
1763            for(i=0; i<(int)sizeof(Dst0); ++i) Dst0[0][i] = 0x5a;
1764    
1765            printf( "\n ===  test YUV ===\n" );
1766    
1767            for(with_flip=0; with_flip<=1; ++with_flip) {
1768    
1769                    init_cpu(&cpu_list[0]);
1770                    TEST_YUYV(yv12_to_yuyv_c, 4, with_flip);
1771                    printf(" yv12_to_yuyv_c %.3f usec       crc32=0x%08x %s\n",
1772                               t, iCrc, (iCrc!=yuv_CRCs[0][with_flip])?"| ERROR": "" );
1773                    TEST_YUYV(yv12_to_uyvy_c, 4, with_flip);
1774                    printf(" yv12_to_uyvy_c %.3f usec       crc32=0x%08x %s\n",
1775                            t, iCrc, (iCrc!=yuv_CRCs[1][with_flip])?"| ERROR": "" );
1776    
1777                    TEST_YUYV(yv12_to_bgra_c, 4, with_flip);
1778                    printf(" yv12_to_bgra_c %.3f usec       crc32=0x%08x %s\n",
1779                            t, iCrc, (iCrc!=yuv_CRCs[2][with_flip])?"| ERROR": "" );
1780    
1781    #if defined(ARCH_IS_IA32) || defined(ARCH_IS_X86_64)
1782                    init_cpu(&cpu_list[1]);
1783                    TEST_YUYV(yv12_to_yuyv_mmx, 4, with_flip);
1784                    printf(" yv12_to_yuyv_mmx %.3f usec       crc32=0x%08x %s\n",
1785                            t, iCrc, (iCrc!=yuv_CRCs[3][with_flip])?"| ERROR": "" );
1786    
1787                    TEST_YUYV(yv12_to_uyvy_mmx, 4, with_flip);
1788                    printf(" yv12_to_uyvy_mmx %.3f usec       crc32=0x%08x %s\n",
1789                            t, iCrc, (iCrc!=yuv_CRCs[4][with_flip])?"| ERROR": "" );
1790    
1791                    TEST_YUYV(yv12_to_bgra_mmx, 4, with_flip);
1792                    printf(" yv12_to_bgra_mmx %.3f usec       crc32=0x%08x %s\n",
1793                            t, iCrc, (iCrc!=yuv_CRCs[5][with_flip])?"| ERROR": "" );
1794    
1795    #endif
1796    
1797    #ifdef ARCH_IS_PPC
1798                    init_cpu(&cpu_list[1]);
1799                    TEST_YUYV(yv12_to_yuyv_altivec_c, 4, with_flip);
1800                    printf(" yv12_to_yuyv_altivec_c %.3f usec       crc32=0x%08x %s\n",
1801                            t, iCrc, (iCrc!=yuv_CRCs[3][with_flip])?"| ERROR": "" );
1802    
1803                    TEST_YUYV(yv12_to_uyvy_altivec_c, 4, with_flip);
1804                    printf(" yv12_to_uyvy_altivec_c %.3f usec       crc32=0x%08x %s\n",
1805                            t, iCrc, (iCrc!=yuv_CRCs[4][with_flip])?"| ERROR": "" );
1806    
1807                    TEST_YUYV(yv12_to_bgra_altivec_c, 4, with_flip);
1808                    printf(" yv12_to_bgra_altivec_c %.3f usec       crc32=0x%08x %s\n",
1809                            t, iCrc, (iCrc!=yuv_CRCs[5][with_flip])?"| ERROR": "" );
1810    #endif
1811            }
1812            printf( " --- \n" );
1813    }
1814    
1815    #define TEST_YV2(FUNC, WITH_UV, WITH_FLIP)        \
1816    ENTER                               \
1817    for(tst=0; tst<nb_tests; ++tst) (FUNC)(Dst0[0], Dst0[1], Dst0[2], WIDTH, WIDTH, \
1818            Src0[0], (WITH_UV) ? Src0[1] : 0, (WITH_UV) ? Src0[2] : 0,  WIDTH, WIDTH, \
1819            WIDTH-2, HEIGHT-2, WITH_FLIP); \
1820    LEAVE
1821    
1822    #define PRINT_NxN(DATA,W,H,STR)   {   \
1823            int i,j; \
1824            for(j=0; j<(H); ++j) { \
1825                    for(i=0; i<(W); ++i) printf( "0x%.2x ", (DATA)[i+j*(STR)] );\
1826                    printf("\n"); \
1827            } \
1828            printf("---\n"); \
1829    }
1830    
1831    static const int yv12_CRCs[2][2] = {
1832            {0x5cab7cf0,0xdab46541}
1833    ,       {0xe8bae865,0x1faf77b7}
1834    };
1835    
1836    void test_yuv2()
1837    {
1838            const int nb_tests = 800*speed_ref;
1839            CPU *cpu;
1840            uint8_t Src0[3][WIDTH*HEIGHT];
1841            uint8_t Dst0[3][WIDTH*HEIGHT];
1842            int with_uv, with_flip;
1843            int i, j;
1844            double t;
1845            int tst, iCrc;
1846    
1847            colorspace_init();
1848            ieee_reseed(1);
1849            for(i=0; i<(int)sizeof(Src0); ++i) Src0[0][i] = ieee_rand(0,255);
1850    
1851            printf( "\n ===  test YV2 ===\n" );
1852            for(with_flip=0; with_flip<=1; ++with_flip) {
1853                    for(with_uv=0; with_uv<=1; ++with_uv) {
1854                            init_cpu(&cpu_list[0]);
1855                            TEST_YV2(yv12_to_yv12_c, with_uv, with_flip);
1856                            printf(" yv12_to_yv12_c   %.3f usec      \tcrc32=0x%08x %s\n",
1857                                    t, iCrc, (iCrc!=yv12_CRCs[with_flip][with_uv])?"| ERROR": "" );
1858                            /* if (!with_uv) PRINT_NxN(Dst0[1], WIDTH/2, HEIGHT/2, WIDTH ); */
1859    
1860    #if defined(ARCH_IS_IA32) || defined(ARCH_IS_X86_64)
1861                            init_cpu(&cpu_list[1]);
1862                            TEST_YV2(yv12_to_yv12_mmx, with_uv, with_flip);
1863                            printf(" yv12_to_yv12_mmx %.3f usec     \tcrc32=0x%08x %s\n",
1864                                    t, iCrc, (iCrc!=yv12_CRCs[with_flip][with_uv])?"| ERROR": "" );
1865                            /* if (!with_uv) PRINT_NxN(Dst0[1], WIDTH/2, HEIGHT/2, WIDTH ); */
1866    
1867                            TEST_YV2(yv12_to_yv12_xmm, with_uv, with_flip);
1868                            printf(" yv12_to_yv12_xmm %.3f usec     \tcrc32=0x%08x %s\n",
1869                                    t, iCrc, (iCrc!=yv12_CRCs[with_flip][with_uv])?"| ERROR": "" );
1870    #endif
1871                    }
1872    
1873                    printf( " --- \n" );
1874            }
1875            printf( " ===== \n" );
1876    }
1877    
1878    #undef WIDTH
1879    #undef HEIGHT
1880    #undef ENTER
1881    #undef LEAVE
1882    
1883    /*********************************************************************/
1884    
1885    static uint32_t __inline log2bin_v1(uint32_t value)
1886    {
1887      int n = 0;
1888      while (value) {
1889        value >>= 1;
1890        n++;
1891      }
1892      return n;
1893    }
1894    
1895    static const uint8_t log2_tab_16[16] =  { 0, 1, 2, 2, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4, 4 };
1896    
1897    static uint32_t __inline log2bin_v2(uint32_t value)
1898    {
1899      int n = 0;
1900      if (value & 0xffff0000) {
1901        value >>= 16;
1902        n += 16;
1903      }
1904      if (value & 0xff00) {
1905        value >>= 8;
1906        n += 8;
1907      }
1908      if (value & 0xf0) {
1909        value >>= 4;
1910        n += 4;
1911      }
1912     return n + log2_tab_16[value];
1913    }
1914    
1915    void test_log2bin()
1916    {
1917            const int nb_tests = 3000*speed_ref;
1918      int n, crc1=0, crc2=0;
1919      uint32_t s, s0;
1920      double t1, t2;
1921    
1922      t1 = gettime_usec();
1923      s0 = (int)(t1*31.241);
1924      for(s=s0, n=0; n<nb_tests; ++n, s=(s*12363+31)&0x7fffffff)
1925        crc1 += log2bin_v1(s);
1926      t1 = (gettime_usec()-t1) / nb_tests;
1927    
1928      t2 = gettime_usec();
1929      for(s=s0, n=0; n<nb_tests; ++n, s=(s*12363+31)&0x7fffffff)
1930        crc2 += log2bin_v2(s);
1931      t2 = (gettime_usec() - t2) / nb_tests;
1932    
1933      printf( "log2bin_v1: %.3f sec  crc=%d\n", t1, crc1 );
1934      printf( "log2bin_v2: %.3f sec  crc=%d\n", t2, crc2 );
1935      if (crc1!=crc2) printf( " CRC ERROR !\n" );
1936    }
1937    
1938    /*********************************************************************/
1939    
1940    static void __inline old_gcd(int *num, int *den)
1941    {
1942     int i = *num;
1943      while (i > 1) {
1944        if (*num % i == 0 && *den % i == 0) {
1945          *num /= i;
1946          *den /= i;
1947          i = *num;
1948          continue;
1949        }
1950        i--;
1951      }
1952    }
1953    
1954    static uint32_t gcd(int num, int den)
1955    {
1956      int tmp;
1957      while( (tmp=num%den) ) { num = den; den = tmp; }
1958      return den;
1959    }
1960    static void __inline new_gcd(int *num, int *den)
1961    {
1962      const int div = gcd(*num, *den);
1963      if (num) {
1964        *num /= div;
1965        *den /= div;
1966      }
1967    }
1968    
1969    void test_gcd()
1970    {
1971      const int nb_tests = 10*speed_ref;
1972      int i;
1973      uint32_t crc1=0, crc2=0;
1974      uint32_t n0, n, d0, d;
1975      double t1, t2;
1976    
1977      t1 = gettime_usec();
1978      n0 = 0xfffff & (int)(t1*31.241);
1979      d0 = 0xfffff & (int)( ((n0*4123)%17) | 1 );
1980      for(n=n0, d=d0, i=0; i<nb_tests; ++i) {
1981        old_gcd(&n, &d);
1982        crc1 = (((crc1>>4)^d) + ((crc1<<2)^n) ) & 0xffffff;
1983        n = d;
1984        d = (d*12363+31) & 0xffff;
1985        d |= !d;
1986      }
1987      t1 = (gettime_usec()-t1) / nb_tests;
1988    
1989      t2 = gettime_usec();
1990      for(n=n0, d=d0, i=0; i<nb_tests; ++i) {
1991        new_gcd(&n, &d);
1992        crc2 = (((crc2>>4)^d) + ((crc2<<2)^n) ) & 0xffffff;
1993        n = d;
1994        d = (d*12363+31) & 0xffff;
1995        d |= !d;
1996      }
1997      t2 = (gettime_usec() - t2) / nb_tests;
1998    
1999      printf( "old_gcd: %.3f sec  crc=%d\n", t1, crc1 );
2000      printf( "new_gcd: %.3f sec  crc=%d\n", t2, crc2 );
2001      if (crc1!=crc2) printf( " CRC ERROR !\n" );
2002    }
2003    
2004    /*********************************************************************
2005     * test compiler
2006     *********************************************************************/
2007    
2008    void test_compiler() {
2009      int nb_err = 0;
2010      int32_t v;
2011      if (sizeof(uint16_t)<2) {
2012        printf( "ERROR: sizeof(uint16_t)<2 !!\n" );
2013        nb_err++;
2014      }
2015      if (sizeof(int16_t)<2) {
2016        printf( "ERROR: sizeof(int16_t)<2 !!\n" );
2017        nb_err++;
2018      }
2019      if (sizeof(uint8_t)!=1) {
2020        printf( "ERROR: sizeof(uint8_t)!=1 !!\n" );
2021        nb_err++;
2022      }
2023      if (sizeof(int8_t)!=1) {
2024        printf( "ERROR: sizeof(int8_t)!=1 !!\n" );
2025        nb_err++;
2026      }
2027      if (sizeof(uint32_t)<4) {
2028        printf( "ERROR: sizeof(uint32_t)<4 !!\n" );
2029        nb_err++;
2030      }
2031      if (sizeof(int32_t)<4) {
2032        printf( "ERROR: sizeof(int32_t)<4 !!\n" );
2033        nb_err++;
2034      }
2035             /* yes, i know, this test is silly. But better be safe than sorry. :) */
2036      for(v=1000; v>=0; v--) {
2037        if ( (v>>2) != v/4)
2038          nb_err++;
2039      }
2040      for(v=-1000; v!=-1; v++) {
2041        if ( (v>>2) != (v/4)-!!(v%4))
2042          nb_err++;
2043      }
2044      if (nb_err!=0) {
2045        printf( "ERROR! please post your platform/compiler specs to xvid-devel@xvid.org !\n" );
2046      }
2047    }
2048    
2049    /*********************************************************************
2050     * test SSIM functions
2051     *********************************************************************/
2052    
2053    typedef int (*lumfunc)(uint8_t* ptr, int stride);
2054    typedef void (*csfunc)(uint8_t* ptro, uint8_t* ptrc, int stride, int lumo, int lumc, int* pdevo, int* pdevc, int* pcorr);
2055    
2056    extern int lum_8x8_c(uint8_t* ptr, int stride);
2057    extern int lum_8x8_mmx(uint8_t* ptr, int stride);
2058    extern int lum_2x8_c(uint8_t* ptr, int stride);
2059    extern void consim_c(uint8_t* ptro, uint8_t* ptrc, int stride, int lumo, int lumc, int* pdevo, int* pdevc, int* pcorr);
2060    extern void consim_mmx(uint8_t* ptro, uint8_t* ptrc, int stride, int lumo, int lumc, int* pdevo, int* pdevc, int* pcorr);
2061    extern void consim_sse2(uint8_t* ptro, uint8_t* ptrc, int stride, int lumo, int lumc, int* pdevo, int* pdevc, int* pcorr);
2062    
2063    void test_SSIM()
2064    {
2065            const int nb_tests = 3000*speed_ref;
2066            int tst;
2067            CPU *cpu;
2068            int i;
2069            int devs[3];
2070            long lumo, lumc;
2071            DECLARE_ALIGNED_MATRIX(Ref1, 16, 16, uint8_t, 16);
2072            DECLARE_ALIGNED_MATRIX(Ref2, 16, 16, uint8_t, 16);
2073            lumfunc lum8x8;
2074            lumfunc lum2x8;
2075            csfunc  csim;
2076    
2077            ieee_reseed(1);
2078            printf( "\n ======  test SSIM ======\n" );
2079            for(i=0; i<16*16;++i) {
2080                    long v1, v2;
2081                    v1 = ieee_rand(-256, 511);
2082                    v2 = ieee_rand(-256, 511);
2083                    Ref1[i] = (v1<0) ? 0 : (v1>255) ? 255 : v1;
2084                    Ref2[i] = (v2<0) ? 0 : (v2>255) ? 255 : v2;
2085            }
2086            lumc = ieee_rand(0, 255);
2087            lumo = ieee_rand(0, 255);
2088    
2089            for(cpu = cpu_list; cpu->name!=0; ++cpu)
2090            {
2091                    double t;
2092                    int m;
2093                    if (!init_cpu(cpu))
2094                            continue;
2095                    lum8x8 = lum_8x8_c;
2096                    lum2x8 = lum_2x8_c;
2097                    csim   = consim_c;
2098    #if defined(ARCH_IS_IA32) || defined(ARCH_IS_X86_64)
2099                    if (cpu->cpu & XVID_CPU_MMX){
2100                            lum8x8 = lum_8x8_mmx;
2101                            csim = consim_mmx;
2102                    }
2103                    if (cpu->cpu & XVID_CPU_MMX){
2104                            csim = consim_sse2;
2105                    }
2106    #endif
2107                    t = gettime_usec();
2108                    emms();
2109                    for(tst=0; tst<nb_tests; ++tst) m = lum8x8(Ref1, 16);
2110                    emms();
2111                    t = (gettime_usec() - t) / nb_tests;
2112                    printf("%s - ssim-lum8x8    %.3f usec       m=%d %s\n",
2113                               cpu->name, t, m,
2114                               (m!=8230)?"| ERROR": "" );
2115    
2116                    t = gettime_usec();
2117                    emms();
2118                    for(tst=0; tst<nb_tests; ++tst) m = lum2x8(Ref1+8, 16);
2119                    emms();
2120                    t = (gettime_usec() - t) / nb_tests;
2121                    printf("%s - ssim-lum2x8    %.3f usec       m=%d %s\n",
2122                               cpu->name, t, m,
2123                               (m!=681)?"| ERROR": "" );
2124    
2125                    t = gettime_usec();
2126                    emms();
2127                    for(tst=0; tst<nb_tests; ++tst) csim(Ref1, Ref2, 16, lumo, lumc, devs+0, devs+1, devs+2);
2128                    emms();
2129                    t = (gettime_usec() - t) / nb_tests;
2130                    printf("%s - ssim-consim    %.3f usec       devs=[0x%x 0x%x 0x%x] %s\n",
2131                               cpu->name, t, devs[0], devs[1], devs[2],
2132                               (devs[0]!=0x1bdf0f || devs[1]!=0x137258 ||  devs[2]!=0xcdb13)?"| ERROR": "" );
2133                    printf( " --- \n" );
2134            }
2135    }
2136    
2137    /*********************************************************************
2138     * test bitstream functions
2139     *********************************************************************/
2140    
2141    #define BIT_BUF_SIZE 2000
2142    
2143    static void test_bits()
2144    {
2145      const int nb_tests = 50*speed_ref;
2146      int tst;
2147      uint32_t Crc;
2148      uint8_t Buf[BIT_BUF_SIZE];
2149      uint32_t Extracted[BIT_BUF_SIZE*8]; /* worst case: bits read 1 by 1 */
2150      int Lens[BIT_BUF_SIZE*8];
2151      double t1;
2152    
2153    
2154      printf( "\n ===  test bitstream ===\n" );
2155      ieee_reseed(1);
2156      Crc = 0;
2157    
2158      t1 = gettime_usec();
2159      for(tst=0; tst<nb_tests; ++tst) {
2160            Bitstream bs;
2161            int m, m2, l, l2;
2162    
2163            for(l=0; l<BIT_BUF_SIZE; ++l)
2164                    Buf[l] = (uint8_t)ieee_rand(0,255);
2165    
2166            l = BIT_BUF_SIZE - ieee_rand(1,BIT_BUF_SIZE/10);
2167            BitstreamInit(&bs, (void*)(Buf+BIT_BUF_SIZE-l), l);
2168    
2169    
2170            BitstreamReset(&bs);
2171            for(l2=l*8, m=0; l2>0; m++) {
2172                    const int b = ieee_rand(1,32);
2173                    Lens[m] = b;
2174                    l2 -= b;
2175                    if (l2<0) break;
2176                    Extracted[m] = BitstreamShowBits(&bs, b);
2177                    BitstreamSkip(&bs, b);
2178    //              printf( "<= %d: %d 0x%x\n", m, b, Extracted[m]);
2179            }
2180    
2181            BitstreamReset(&bs);
2182            for(m2=0; m2<m; ++m2) {
2183                    const int b = Lens[m2];
2184                    const uint32_t v = BitstreamGetBits(&bs, b);
2185                    Crc |= (v!=Extracted[m2]);
2186    //              printf( "=> %d: %d 0x%x %c\n", m2, b, v, " *"[Crc]);
2187            }
2188      }
2189      t1 = (gettime_usec() - t1) / nb_tests;
2190      printf(" test_bits   %.3f usec   %s\n", t1, (Crc!=0)?"| ERROR": "" );
2191    }
2192    
2193    /*********************************************************************
2194   * main   * main
2195   *********************************************************************/   *********************************************************************/
2196    
2197  int main(int argc, char *argv[])  static void arg_missing(const char *opt)
2198    {
2199      printf( "missing argument after option '%s'\n", opt);
2200      exit(-1);
2201    }
2202    
2203    int main(int argc, const char *argv[])
2204    {
2205            int c, what = 0;
2206            int width, height;
2207            uint32_t chksum = 0;
2208            const char * test_bitstream = 0;
2209    #if defined(WIN32) && defined(ARCH_IS_X86_64)
2210            DECLARE_ALIGNED_MATRIX(xmm_save, 2, 4, uint64_t, 16);
2211            // assumes xmm6 and xmm7 won't be falsely preserved by C code
2212            for(c=0;c<4;c++)
2213                    xmm_save[c] = read_counter();
2214            prime_xmm(xmm_save);
2215    #endif
2216    
2217            cpu_mask = 0;  // default => will use autodectect
2218            for(c=1; c<argc; ++c)
2219  {  {
2220          int what = 0;            if (!strcmp(argv[c], "-v")) verbose++;
2221          if (argc>1) what = atoi(argv[1]);            else if (!strcmp(argv[c], "-c"))      cpu_mask = 0 /* PLAIN_C */ | XVID_CPU_FORCE;
2222              else if (!strcmp(argv[c], "-mmx"))    cpu_mask = XVID_CPU_MMX    | XVID_CPU_FORCE;
2223              else if (!strcmp(argv[c], "-mmxext")) cpu_mask = XVID_CPU_MMXEXT | XVID_CPU_MMX | XVID_CPU_FORCE;
2224              else if (!strcmp(argv[c], "-sse2"))   cpu_mask = XVID_CPU_SSE2   | XVID_CPU_MMXEXT | XVID_CPU_MMX | XVID_CPU_FORCE;
2225              else if (!strcmp(argv[c], "-sse3"))   cpu_mask = XVID_CPU_SSE3   | XVID_CPU_SSE2 | XVID_CPU_MMXEXT | XVID_CPU_MMX | XVID_CPU_FORCE;
2226              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;
2227          else if (!strcmp(argv[c], "-3dnow"))  cpu_mask = XVID_CPU_3DNOW  | XVID_CPU_FORCE;
2228              else if (!strcmp(argv[c], "-3dnowe")) cpu_mask = XVID_CPU_3DNOW  | XVID_CPU_3DNOWEXT | XVID_CPU_FORCE;
2229              else if (!strcmp(argv[c], "-altivec")) cpu_mask = XVID_CPU_ALTIVEC | XVID_CPU_FORCE;
2230              else if (!strcmp(argv[c], "-spd")) {
2231          if (++c==argc) arg_missing( argv[argc-1] );
2232          speed_ref = atoi(argv[c]);
2233        }
2234              else if (argv[c][0]!='-') {
2235                what = atoi(argv[c]);
2236                if (what==9) {
2237                  if (c+4>argc) {
2238                    printf("usage: %s %d bitstream width height (checksum)\n", argv[0], what);
2239                    exit(-1);
2240            }
2241            test_bitstream = argv[++c];
2242                  width  = atoi(argv[++c]);
2243                  height = atoi(argv[++c]);
2244                  if (c+1<argc && argv[c+1][0]!='-') {
2245                    if (sscanf(argv[c+1], "0x%x", &chksum)!=1) {
2246                      printf( "can't read checksum value.\n" );
2247                      exit(-1);
2248              }
2249              else c++;
2250            }
2251    //        printf( "[%s] %dx%d (0x%.8x)\n", test_bitstream, width, height, chksum);
2252          }
2253        }
2254        else {
2255          printf( "unrecognized option '%s'\n", argv[c]);
2256          exit(-1);
2257        }
2258      }
2259    
2260    
2261          if (what==0 || what==1) test_dct();          if (what==0 || what==1) test_dct();
2262          if (what==0 || what==2) test_mb();          if (what==0 || what==2) test_mb();
2263          if (what==0 || what==3) test_sad();          if (what==0 || what==3) test_sad();
2264          if (what==0 || what==4) test_transfer();          if (what==0 || what==4) test_transfer();
2265          if (what==0 || what==5) test_quant();          if (what==0 || what==5) test_quant();
2266          if (what==0 || what==6) test_cbp();          if (what==0 || what==6) test_cbp();
2267            if (what==0 || what==10) test_sse();
2268            if (what==0 || what==11) test_log2bin();
2269            if (what==0 || what==12) test_gcd();
2270            if (what==0 || what==13) test_compiler();
2271            if (what==0 || what==14) test_yuv();
2272            if (what==0 || what==15) test_SSIM();
2273            if (what==0 || what==16) test_yuv2();
2274            if (what==0 || what==17) test_bits();
2275    
2276          if (what==7) {          if (what==7) {
2277                  test_IEEE1180_compliance(-256, 255, 1);                  test_IEEE1180_compliance(-256, 255, 1);
# Line 1509  Line 2283 
2283          }          }
2284          if (what==8) test_dct_saturation(-256, 255);          if (what==8) test_dct_saturation(-256, 255);
2285    
2286          if (what==9) {          if (test_bitstream)
2287                  int width, height;            test_dec(test_bitstream, width, height, chksum);
                 if (argc<5) {  
                         printf("usage: %s %d [bitstream] [width] [height]\n", argv[0], what);  
                         return 1;  
                 }  
                 width = atoi(argv[3]);  
                 height = atoi(argv[4]);  
                 test_dec(argv[2], width, height, (argc>5));  
         }  
   
2288          if (what==-1) {          if (what==-1) {
2289                  test_dct_precision_diffs();                  test_dct_precision_diffs();
2290                  test_bugs1();                  test_bugs1();
# Line 1527  Line 2292 
2292          if (what==-2)          if (what==-2)
2293                  test_quant_bug();                  test_quant_bug();
2294    
2295          if (what >= 0 && what <= 6) {  #if defined(WIN32) && defined(ARCH_IS_X86_64)
2296            get_xmm(xmm_save+4);
2297            if (memcmp(xmm_save, xmm_save+4, 4*sizeof(int64_t))) {
2298                    printf("\nWIN64 ERROR: XMM6 and XMM7 contents not preserved!\n"
2299                           "        XMM6                             XMM7\n"
2300                           "Before: %.16I64X%.16I64X %.16I64X%.16I64X\n"
2301                           "After:  %.16I64X%.16I64X %.16I64X%.16I64X",
2302                            xmm_save[0],xmm_save[1],xmm_save[2],xmm_save[3],
2303                            xmm_save[4],xmm_save[5],xmm_save[6],xmm_save[7]);
2304            } else {
2305                    printf("\nWIN64: XMM6 and XMM7 contents preserved correctly.\n");
2306            }
2307    #endif
2308    
2309            if ((what >= 0 && what <= 6) || what == 10) {
2310                  printf("\n\n"                  printf("\n\n"
2311                             "NB: If a function isn't optimised for a specific set of intructions,\n"                             "NB: If a function isn't optimised for a specific set of intructions,\n"
2312                             "    a C function is used instead. So don't panic if some functions\n"                             "    a C function is used instead. So don't panic if some functions\n"
2313                             "    may appear to be slow.\n");                             "    may appear to be slow.\n");
2314          }          }
2315    
2316  #ifdef ARCH_IS_IA32  #if defined(ARCH_IS_IA32) || defined(ARCH_IS_X86_64)
2317          if (what == 0 || what == 5) {          if (what == 0 || what == 5) {
2318                  printf("\n"                  printf("\n"
2319                             "NB: MMX mpeg4 quantization is known to have very small errors (+/-1 magnitude)\n"                             "NB: MMX mpeg4 quantization is known to have very small errors (+/-1 magnitude)\n"
# Line 1547  Line 2326 
2326          return 0;          return 0;
2327  }  }
2328    
2329  /*********************************************************************  /*********************************************************************/
  * 'Reference' output (except for timing) on an Athlon XP 2200+  
  *********************************************************************/  
   
 /* as of 2002-01-07, there's a problem with MMX mpeg4-quantization */  
 /* as of 2003-11-30, the problem is still here */  
   
 /*********************************************************************  
   
   
  ===== test fdct/idct =====  
 PLAINC -  2.867 usec       PSNR=13.291  MSE=3.000  
 MMX    -  -0.211 usec       PSNR=9.611  MSE=7.000  
 MMXEXT -  -0.256 usec       PSNR=9.611  MSE=7.000  
 3DNOW  -  2.855 usec       PSNR=13.291  MSE=3.000  
 3DNOWE -  1.429 usec       PSNR=13.291  MSE=3.000  
   
  ===  test block motion ===  
 PLAINC - interp- h-round0 0.538 usec       crc32=0x115381ba  
 PLAINC -           round1 0.527 usec       crc32=0x2b1f528f  
 PLAINC - interp- v-round0 0.554 usec       crc32=0x423cdcc7  
 PLAINC -           round1 0.551 usec       crc32=0x42202efe  
 PLAINC - interp-hv-round0 1.041 usec       crc32=0xd198d387  
 PLAINC -           round1 1.038 usec       crc32=0x9ecfd921  
  ---  
 MMX    - interp- h-round0 0.051 usec       crc32=0x115381ba  
 MMX    -           round1 0.053 usec       crc32=0x2b1f528f  
 MMX    - interp- v-round0 0.048 usec       crc32=0x423cdcc7  
 MMX    -           round1 0.048 usec       crc32=0x42202efe  
 MMX    - interp-hv-round0 0.074 usec       crc32=0xd198d387  
 MMX    -           round1 0.073 usec       crc32=0x9ecfd921  
  ---  
 MMXEXT - interp- h-round0 0.020 usec       crc32=0x115381ba  
 MMXEXT -           round1 0.025 usec       crc32=0x2b1f528f  
 MMXEXT - interp- v-round0 0.016 usec       crc32=0x423cdcc7  
 MMXEXT -           round1 0.024 usec       crc32=0x42202efe  
 MMXEXT - interp-hv-round0 0.037 usec       crc32=0xd198d387  
 MMXEXT -           round1 0.037 usec       crc32=0x9ecfd921  
  ---  
 3DNOW  - interp- h-round0 0.020 usec       crc32=0x115381ba  
 3DNOW  -           round1 0.029 usec       crc32=0x2b1f528f  
 3DNOW  - interp- v-round0 0.016 usec       crc32=0x423cdcc7  
 3DNOW  -           round1 0.024 usec       crc32=0x42202efe  
 3DNOW  - interp-hv-round0 0.038 usec       crc32=0xd198d387  
 3DNOW  -           round1 0.039 usec       crc32=0x9ecfd921  
  ---  
 3DNOWE - interp- h-round0 0.020 usec       crc32=0x115381ba  
 3DNOWE -           round1 0.024 usec       crc32=0x2b1f528f  
 3DNOWE - interp- v-round0 0.016 usec       crc32=0x423cdcc7  
 3DNOWE -           round1 0.021 usec       crc32=0x42202efe  
 3DNOWE - interp-hv-round0 0.037 usec       crc32=0xd198d387  
 3DNOWE -           round1 0.036 usec       crc32=0x9ecfd921  
  ---  
   
  ======  test SAD ======  
 PLAINC - sad8    0.505 usec       sad=3776  
 PLAINC - sad16   1.941 usec       sad=27214  
 PLAINC - sad16bi 4.925 usec       sad=26274  
 PLAINC - dev16   4.254 usec       sad=3344  
  ---  
 MMX    - sad8    0.036 usec       sad=3776  
 MMX    - sad16   0.107 usec       sad=27214  
 MMX    - sad16bi 0.259 usec       sad=26274  
 MMX    - dev16   0.187 usec       sad=3344  
  ---  
 MMXEXT - sad8    0.016 usec       sad=3776  
 MMXEXT - sad16   0.050 usec       sad=27214  
 MMXEXT - sad16bi 0.060 usec       sad=26274  
 MMXEXT - dev16   0.086 usec       sad=3344  
  ---  
 3DNOW  - sad8    0.506 usec       sad=3776  
 3DNOW  - sad16   1.954 usec       sad=27214  
 3DNOW  - sad16bi 0.119 usec       sad=26274  
 3DNOW  - dev16   4.252 usec       sad=3344  
  ---  
 3DNOWE - sad8    0.017 usec       sad=3776  
 3DNOWE - sad16   0.038 usec       sad=27214  
 3DNOWE - sad16bi 0.052 usec       sad=26274  
 3DNOWE - dev16   0.067 usec       sad=3344  
  ---  
   
  ===  test transfer ===  
 PLAINC - 8to16     0.603 usec       crc32=0x115814bb  
 PLAINC - 16to8     1.077 usec       crc32=0xee7ccbb4  
 PLAINC - 8to8      0.679 usec       crc32=0xd37b3295  
 PLAINC - 16to8add  1.341 usec       crc32=0xdd817bf4  
 PLAINC - 8to16sub  1.566 usec       crc32(1)=0xa1e07163 crc32(2)=0xd86c5d23  
 PLAINC - 8to16sub2 2.206 usec       crc32=0x99b6c4c7  
  ---  
 MMX    - 8to16     -0.025 usec       crc32=0x115814bb  
 MMX    - 16to8     -0.049 usec       crc32=0xee7ccbb4  
 MMX    - 8to8      0.014 usec       crc32=0xd37b3295  
 MMX    - 16to8add  0.011 usec       crc32=0xdd817bf4  
 MMX    - 8to16sub  0.108 usec       crc32(1)=0xa1e07163 crc32(2)=0xd86c5d23  
 MMX    - 8to16sub2 0.164 usec       crc32=0x99b6c4c7  
  ---  
 MMXEXT - 8to16     -0.054 usec       crc32=0x115814bb  
 MMXEXT - 16to8     0.010 usec       crc32=0xee7ccbb4  
 MMXEXT - 8to8      0.015 usec       crc32=0xd37b3295  
 MMXEXT - 16to8add  0.008 usec       crc32=0xdd817bf4  
 MMXEXT - 8to16sub  0.263 usec       crc32(1)=0xa1e07163 crc32(2)=0xd86c5d23  
 MMXEXT - 8to16sub2 0.178 usec       crc32=0x99b6c4c7  
  ---  
 3DNOW  - 8to16     0.666 usec       crc32=0x115814bb  
 3DNOW  - 16to8     1.078 usec       crc32=0xee7ccbb4  
 3DNOW  - 8to8      0.665 usec       crc32=0xd37b3295  
 3DNOW  - 16to8add  1.365 usec       crc32=0xdd817bf4  
 3DNOW  - 8to16sub  1.356 usec       crc32(1)=0xa1e07163 crc32(2)=0xd86c5d23  
 3DNOW  - 8to16sub2 2.098 usec       crc32=0x99b6c4c7  
  ---  
 3DNOWE - 8to16     -0.024 usec       crc32=0x115814bb  
 3DNOWE - 16to8     0.010 usec       crc32=0xee7ccbb4  
 3DNOWE - 8to8      0.014 usec       crc32=0xd37b3295  
 3DNOWE - 16to8add  0.016 usec       crc32=0xdd817bf4  
 3DNOWE - 8to16sub  -0.000 usec       crc32(1)=0xa1e07163 crc32(2)=0xd86c5d23  
 3DNOWE - 8to16sub2 -0.031 usec       crc32=0x99b6c4c7  
  ---  
   
  =====  test quant =====  
 PLAINC -   quant_mpeg_intra 98.631 usec       crc32=0xfd6a21a4  
 PLAINC -   quant_mpeg_inter 104.876 usec       crc32=0xf6de7757  
 PLAINC - dequant_mpeg_intra 50.285 usec       crc32=0x2def7bc7  
 PLAINC - dequant_mpeg_inter 58.316 usec       crc32=0xd878c722  
 PLAINC -   quant_h263_intra 33.803 usec       crc32=0x2eba9d43  
 PLAINC -   quant_h263_inter 45.411 usec       crc32=0xbd315a7e  
 PLAINC - dequant_h263_intra 39.302 usec       crc32=0x9841212a  
 PLAINC - dequant_h263_inter 44.124 usec       crc32=0xe7df8fba  
  ---  
 MMX    -   quant_mpeg_intra 4.273 usec       crc32=0xdacabdb6 | ERROR  
 MMX    -   quant_mpeg_inter 3.576 usec       crc32=0x72883ab6 | ERROR  
 MMX    - dequant_mpeg_intra 3.793 usec       crc32=0x2def7bc7  
 MMX    - dequant_mpeg_inter 4.808 usec       crc32=0xd878c722  
 MMX    -   quant_h263_intra 2.881 usec       crc32=0x2eba9d43  
 MMX    -   quant_h263_inter 2.550 usec       crc32=0xbd315a7e  
 MMX    - dequant_h263_intra 2.974 usec       crc32=0x9841212a  
 MMX    - dequant_h263_inter 2.906 usec       crc32=0xe7df8fba  
  ---  
 MMXEXT -   quant_mpeg_intra 4.221 usec       crc32=0xfd6a21a4  
 MMXEXT -   quant_mpeg_inter 4.339 usec       crc32=0xf6de7757  
 MMXEXT - dequant_mpeg_intra 3.802 usec       crc32=0x2def7bc7  
 MMXEXT - dequant_mpeg_inter 4.821 usec       crc32=0xd878c722  
 MMXEXT -   quant_h263_intra 2.884 usec       crc32=0x2eba9d43  
 MMXEXT -   quant_h263_inter 2.554 usec       crc32=0xbd315a7e  
 MMXEXT - dequant_h263_intra 2.728 usec       crc32=0x9841212a  
 MMXEXT - dequant_h263_inter 2.611 usec       crc32=0xe7df8fba  
  ---  
 3DNOW  -   quant_mpeg_intra 98.512 usec       crc32=0xfd6a21a4  
 3DNOW  -   quant_mpeg_inter 104.873 usec       crc32=0xf6de7757  
 3DNOW  - dequant_mpeg_intra 50.219 usec       crc32=0x2def7bc7  
 3DNOW  - dequant_mpeg_inter 58.254 usec       crc32=0xd878c722  
 3DNOW  -   quant_h263_intra 33.778 usec       crc32=0x2eba9d43  
 3DNOW  -   quant_h263_inter 41.998 usec       crc32=0xbd315a7e  
 3DNOW  - dequant_h263_intra 39.344 usec       crc32=0x9841212a  
 3DNOW  - dequant_h263_inter 43.607 usec       crc32=0xe7df8fba  
  ---  
 3DNOWE -   quant_mpeg_intra 98.490 usec       crc32=0xfd6a21a4  
 3DNOWE -   quant_mpeg_inter 104.889 usec       crc32=0xf6de7757  
 3DNOWE - dequant_mpeg_intra 3.277 usec       crc32=0x2def7bc7  
 3DNOWE - dequant_mpeg_inter 4.485 usec       crc32=0xd878c722  
 3DNOWE -   quant_h263_intra 1.882 usec       crc32=0x2eba9d43  
 3DNOWE -   quant_h263_inter 2.246 usec       crc32=0xbd315a7e  
 3DNOWE - dequant_h263_intra 3.457 usec       crc32=0x9841212a  
 3DNOWE - dequant_h263_inter 3.275 usec       crc32=0xe7df8fba  
  ---  
   
  =====  test cbp =====  
 PLAINC -   calc_cbp#1 0.168 usec       cbp=0x15  
 PLAINC -   calc_cbp#2 0.168 usec       cbp=0x38  
 PLAINC -   calc_cbp#3 0.157 usec       cbp=0x0f  
 PLAINC -   calc_cbp#4 0.235 usec       cbp=0x05  
  ---  
 MMX    -   calc_cbp#1 0.070 usec       cbp=0x15  
 MMX    -   calc_cbp#2 0.062 usec       cbp=0x38  
 MMX    -   calc_cbp#3 0.062 usec       cbp=0x0f  
 MMX    -   calc_cbp#4 0.061 usec       cbp=0x05  
  ---  
 MMXEXT -   calc_cbp#1 0.062 usec       cbp=0x15  
 MMXEXT -   calc_cbp#2 0.061 usec       cbp=0x38  
 MMXEXT -   calc_cbp#3 0.061 usec       cbp=0x0f  
 MMXEXT -   calc_cbp#4 0.061 usec       cbp=0x05  
  ---  
 3DNOW  -   calc_cbp#1 0.168 usec       cbp=0x15  
 3DNOW  -   calc_cbp#2 0.168 usec       cbp=0x38  
 3DNOW  -   calc_cbp#3 0.157 usec       cbp=0x0f  
 3DNOW  -   calc_cbp#4 0.238 usec       cbp=0x05  
  ---  
 3DNOWE -   calc_cbp#1 0.049 usec       cbp=0x15  
 3DNOWE -   calc_cbp#2 0.049 usec       cbp=0x38  
 3DNOWE -   calc_cbp#3 0.049 usec       cbp=0x0f  
 3DNOWE -   calc_cbp#4 0.049 usec       cbp=0x05  
  ---  
   
   
 NB: If a function isn't optimised for a specific set of intructions,  
     a C function is used instead. So don't panic if some functions  
     may appear to be slow.  
   
 NB: MMX mpeg4 quantization is known to have very small errors (+/-1 magnitude)  
     for 1 or 2 coefficients a block. This is mainly caused by the fact the unit  
     test goes far behind the usual limits of real encoding. Please do not report  
     this error to the developers  
   
 *********************************************************************/  
Legend:
Removed from v.1398  
changed lines
  Added in v.1839
No admin address has been configured
 ViewVC Help
Powered by ViewVC 1.0.4