[cvs] / xvidcore / src / xvid.c Repository:
ViewVC logotype

Diff of /xvidcore/src/xvid.c

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

revision 1.33.2.19, Mon Dec 30 10:49:17 2002 UTC revision 1.45.2.4, Sat May 17 13:26:42 2003 UTC
# Line 3  Line 3 
3   *  XVID MPEG-4 VIDEO CODEC   *  XVID MPEG-4 VIDEO CODEC
4   *  - Native API implementation  -   *  - Native API implementation  -
5   *   *
  *  This program is an implementation of a part of one or more MPEG-4  
  *  Video tools as specified in ISO/IEC 14496-2 standard.  Those intending  
  *  to use this software module in hardware or software products are  
  *  advised that its use may infringe existing patents or copyrights, and  
  *  any such use would be at such party's own risk.  The original  
  *  developer of this software module and his/her company, and subsequent  
  *  editors and their companies, will have no liability for use of this  
  *  software or modifications or derivatives thereof.  
  *  
6   *  This program is free software ; you can redistribute it and/or modify   *  This program is free software ; you can redistribute it and/or modify
7   *  it under the terms of the GNU General Public License as published by   *  it under the terms of the GNU General Public License as published by
8   *  the Free Software Foundation ; either version 2 of the License, or   *  the Free Software Foundation ; either version 2 of the License, or
# Line 26  Line 17 
17   *  along with this program ; if not, write to the Free Software   *  along with this program ; if not, write to the Free Software
18   *  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA   *  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
19   *   *
  ****************************************************************************/  
   
 /*****************************************************************************  
  *  
  *  History  
  *  
  *      - 23.06.2002    added XVID_CPU_CHKONLY  
  *  - 17.03.2002        Added interpolate8x8_halfpel_hv_xmm  
  *  - 22.12.2001  API change: added xvid_init() - Isibaar  
  *  - 16.12.2001        inital version; (c)2001 peter ross <pross@cs.rmit.edu.au>  
  *  
20   *  $Id$   *  $Id$
21   *   *
22   ****************************************************************************/   ****************************************************************************/
23    
24    #include <stdio.h>
25    #include <stdlib.h>
26    #include <string.h>
27    #include <time.h>
28    
29  #include "xvid.h"  #include "xvid.h"
30  #include "decoder.h"  #include "decoder.h"
31  #include "encoder.h"  #include "encoder.h"
# Line 60  Line 45 
45  #include "utils/timer.h"  #include "utils/timer.h"
46  #include "bitstream/mbcoding.h"  #include "bitstream/mbcoding.h"
47    
48  #if defined(ARCH_X86) && defined(EXPERIMENTAL_SSE2_CODE)  #if defined(_DEBUG)
49    unsigned int xvid_debug = 0; /* xvid debug mask */
50    #endif
51    
52    #if defined(ARCH_IS_IA32)
53    
54  #ifdef WIN32  #if defined(_MSC_VER)
55  #include <windows.h>  #include <windows.h>
56  #else  #else
57  #include <signal.h>  #include <signal.h>
58  #include <setjmp.h>  #include <setjmp.h>
 #endif  
   
   
 #ifndef WIN32  
59    
60  static jmp_buf mark;  static jmp_buf mark;
61    
# Line 83  Line 68 
68    
69    
70  /*  /*
71  calls the funcptr, and returns whether SIGILL (illegal instruction) was signalled   * Calls the funcptr, and returns whether SIGILL (illegal instruction) was
72  return values:   * signalled
73  -1 : could not determine   *
74  0  : SIGILL was *not* signalled   * Return values:
75  1  : SIGILL was signalled   *  -1 : could not determine
76     *   0 : SIGILL was *not* signalled
77     *   1 : SIGILL was signalled
78  */  */
79    
80  int  int
81  sigill_check(void (*func)())  sigill_check(void (*func)())
82  {  {
83  #ifdef WIN32  #if defined(_MSC_VER)
84          _try {          _try {
85                  func();                  func();
86          }          }
# Line 127  Line 114 
114  }  }
115  #endif  #endif
116    
117    
118    /* detect cpu flags  */
119    static unsigned int
120    detect_cpu_flags()
121    {
122            /* enable native assembly optimizations by default */
123            unsigned int cpu_flags = XVID_CPU_ASM;
124    
125    #if defined(ARCH_IS_IA32)
126            cpu_flags |= check_cpu_features();
127            if ((cpu_flags & XVID_CPU_SSE) && sigill_check(sse_os_trigger))
128                    cpu_flags &= ~XVID_CPU_SSE;
129    
130            if ((cpu_flags & XVID_CPU_SSE2) && sigill_check(sse2_os_trigger))
131                    cpu_flags &= ~XVID_CPU_SSE2;
132    #endif
133    
134    #if defined(ARCH_IS_PPC)
135    #if defined(ARCH_IS_PPC_ALTIVEC)
136            cpu_flags |= XVID_CPU_ALTIVEC;
137    #endif
138    #endif
139    
140            return cpu_flags;
141    }
142    
143    
144  /*****************************************************************************  /*****************************************************************************
145   * XviD Init Entry point   * XviD Init Entry point
146   *   *
# Line 143  Line 157 
157    
158    
159  static  static
160  int xvid_init_init(XVID_INIT_PARAM * init_param)  int xvid_gbl_init(xvid_gbl_init_t * init)
 {  
         int cpu_flags;  
   
         /* Inform the client the API version */  
         init_param->api_version = API_VERSION;  
   
         /* Inform the client the core build - unused because we're still alpha */  
         init_param->core_build = 1000;  
   
         /* Do we have to force CPU features  ? */  
         if ((init_param->cpu_flags & XVID_CPU_FORCE)) {  
   
                 cpu_flags = init_param->cpu_flags;  
   
         } else {  
   
                 cpu_flags = check_cpu_features();  
   
 #if defined(ARCH_X86) && defined(EXPERIMENTAL_SSE2_CODE)  
                 if ((cpu_flags & XVID_CPU_SSE) && sigill_check(sse_os_trigger))  
                         cpu_flags &= ~XVID_CPU_SSE;  
   
                 if ((cpu_flags & XVID_CPU_SSE2) && sigill_check(sse2_os_trigger))  
                         cpu_flags &= ~XVID_CPU_SSE2;  
 #endif  
         }  
   
         if ((init_param->cpu_flags & XVID_CPU_CHKONLY))  
161          {          {
162                  init_param->cpu_flags = cpu_flags;          unsigned int cpu_flags;
                 return XVID_ERR_OK;  
         }  
163    
164          init_param->cpu_flags = cpu_flags;          if (XVID_MAJOR(init->version) != 1) /* v1.x.x */
165                    return XVID_ERR_VERSION;
166    
167            cpu_flags = (init->cpu_flags & XVID_CPU_FORCE) ? init->cpu_flags : detect_cpu_flags();
168    
169          /* Initialize the function pointers */          /* Initialize the function pointers */
170          idct_int32_init();          idct_int32_init();
# Line 209  Line 195 
195          transfer_8to16copy = transfer_8to16copy_c;          transfer_8to16copy = transfer_8to16copy_c;
196          transfer_16to8copy = transfer_16to8copy_c;          transfer_16to8copy = transfer_16to8copy_c;
197          transfer_8to16sub  = transfer_8to16sub_c;          transfer_8to16sub  = transfer_8to16sub_c;
198            transfer_8to16subro  = transfer_8to16subro_c;
199          transfer_8to16sub2 = transfer_8to16sub2_c;          transfer_8to16sub2 = transfer_8to16sub2_c;
200          transfer_16to8add  = transfer_16to8add_c;          transfer_16to8add  = transfer_16to8add_c;
201          transfer8x8_copy   = transfer8x8_copy_c;          transfer8x8_copy   = transfer8x8_copy_c;
# Line 236  Line 223 
223          interpolate8x8_avg4 = interpolate8x8_avg4_c;          interpolate8x8_avg4 = interpolate8x8_avg4_c;
224    
225          /* reduced resoltuion */          /* reduced resoltuion */
   
226          copy_upsampled_8x8_16to8 = xvid_Copy_Upsampled_8x8_16To8_C;          copy_upsampled_8x8_16to8 = xvid_Copy_Upsampled_8x8_16To8_C;
227          add_upsampled_8x8_16to8 = xvid_Add_Upsampled_8x8_16To8_C;          add_upsampled_8x8_16to8 = xvid_Add_Upsampled_8x8_16To8_C;
 #ifdef ARCH_X86  
         vfilter_31 = xvid_VFilter_31_x86;  
         hfilter_31 = xvid_HFilter_31_x86;  
 #else  
228          vfilter_31 = xvid_VFilter_31_C;          vfilter_31 = xvid_VFilter_31_C;
229          hfilter_31 = xvid_HFilter_31_C;          hfilter_31 = xvid_HFilter_31_C;
 #endif  
230          filter_18x18_to_8x8 = xvid_Filter_18x18_To_8x8_C;          filter_18x18_to_8x8 = xvid_Filter_18x18_To_8x8_C;
231          filter_diff_18x18_to_8x8 = xvid_Filter_Diff_18x18_To_8x8_C;          filter_diff_18x18_to_8x8 = xvid_Filter_Diff_18x18_To_8x8_C;
232    
# Line 301  Line 282 
282          dev16    = dev16_c;          dev16    = dev16_c;
283          sad16v   = sad16v_c;          sad16v   = sad16v_c;
284    
285  //      Halfpel8_Refine = Halfpel8_Refine_c;  /*      Halfpel8_Refine = Halfpel8_Refine_c; */
286    
287    #if defined(ARCH_IS_IA32)
288    
289  #ifdef ARCH_X86          if ((cpu_flags & XVID_CPU_ASM))
290            {
291                    vfilter_31 = xvid_VFilter_31_x86;
292                    hfilter_31 = xvid_HFilter_31_x86;
293            }
294    
295          if ((cpu_flags & XVID_CPU_MMX) || (cpu_flags & XVID_CPU_MMXEXT) ||          if ((cpu_flags & XVID_CPU_MMX) || (cpu_flags & XVID_CPU_MMXEXT) ||
296                  (cpu_flags & XVID_CPU_3DNOW) || (cpu_flags & XVID_CPU_3DNOWEXT) ||                  (cpu_flags & XVID_CPU_3DNOW) || (cpu_flags & XVID_CPU_3DNOWEXT) ||
# Line 313  Line 300 
300                  emms = emms_mmx;                  emms = emms_mmx;
301          }          }
302    
303          if ((cpu_flags & XVID_CPU_MMX) > 0) {          if ((cpu_flags & XVID_CPU_MMX)) {
304    
305                  /* Forward and Inverse Discrete Cosine Transformation functions */                  /* Forward and Inverse Discrete Cosine Transformation functions */
306                  fdct = fdct_mmx;                  fdct = fdct_mmx;
# Line 334  Line 321 
321                  transfer_8to16copy = transfer_8to16copy_mmx;                  transfer_8to16copy = transfer_8to16copy_mmx;
322                  transfer_16to8copy = transfer_16to8copy_mmx;                  transfer_16to8copy = transfer_16to8copy_mmx;
323                  transfer_8to16sub  = transfer_8to16sub_mmx;                  transfer_8to16sub  = transfer_8to16sub_mmx;
324                    transfer_8to16subro  = transfer_8to16subro_mmx;
325                  transfer_8to16sub2 = transfer_8to16sub2_mmx;                  transfer_8to16sub2 = transfer_8to16sub2_mmx;
326                  transfer_16to8add  = transfer_16to8add_mmx;                  transfer_16to8add  = transfer_16to8add_mmx;
327                  transfer8x8_copy   = transfer8x8_copy_mmx;                  transfer8x8_copy   = transfer8x8_copy_mmx;
# Line 383  Line 371 
371                  sad8bi  = sad8bi_mmx;                  sad8bi  = sad8bi_mmx;
372                  dev16    = dev16_mmx;                  dev16    = dev16_mmx;
373                  sad16v   = sad16v_mmx;                  sad16v   = sad16v_mmx;
   
374          }          }
375    
376          /* these 3dnow functions are faster than mmx, but slower than xmm. */          /* these 3dnow functions are faster than mmx, but slower than xmm. */
377          if ((cpu_flags & XVID_CPU_3DNOW) > 0) {          if ((cpu_flags & XVID_CPU_3DNOW)) {
378    
379                    emms = emms_3dn;
380    
381                  /* ME functions */                  /* ME functions */
382                  sad16bi = sad16bi_3dn;                  sad16bi = sad16bi_3dn;
# Line 398  Line 387 
387          }          }
388    
389    
390          if ((cpu_flags & XVID_CPU_MMXEXT) > 0) {          if ((cpu_flags & XVID_CPU_MMXEXT)) {
391    
392                  /* Inverse DCT */                  /* Inverse DCT */
393                  idct = idct_xmm;                  idct = idct_xmm;
# Line 436  Line 425 
425                  sad16v   = sad16v_xmm;                  sad16v   = sad16v_xmm;
426          }          }
427    
428          if ((cpu_flags & XVID_CPU_3DNOW) > 0) {          if ((cpu_flags & XVID_CPU_3DNOW)) {
429    
430                  /* Interpolation */                  /* Interpolation */
431                  interpolate8x8_halfpel_h = interpolate8x8_halfpel_h_3dn;                  interpolate8x8_halfpel_h = interpolate8x8_halfpel_h_3dn;
# Line 444  Line 433 
433                  interpolate8x8_halfpel_hv = interpolate8x8_halfpel_hv_3dn;                  interpolate8x8_halfpel_hv = interpolate8x8_halfpel_hv_3dn;
434          }          }
435    
436          if ((cpu_flags & XVID_CPU_3DNOWEXT) > 0) {          if ((cpu_flags & XVID_CPU_3DNOWEXT)) {
437    
438                  /* Inverse DCT */                  /* Inverse DCT */
439                  idct =  idct_3dne;                  idct =  idct_3dne;
# Line 453  Line 442 
442                  transfer_8to16copy =  transfer_8to16copy_3dne;                  transfer_8to16copy =  transfer_8to16copy_3dne;
443                  transfer_16to8copy = transfer_16to8copy_3dne;                  transfer_16to8copy = transfer_16to8copy_3dne;
444                  transfer_8to16sub =  transfer_8to16sub_3dne;                  transfer_8to16sub =  transfer_8to16sub_3dne;
445                    transfer_8to16subro =  transfer_8to16subro_3dne;
446                  transfer_8to16sub2 =  transfer_8to16sub2_3dne;                  transfer_8to16sub2 =  transfer_8to16sub2_3dne;
447                  transfer_16to8add = transfer_16to8add_3dne;                  transfer_16to8add = transfer_16to8add_3dne;
448                  transfer8x8_copy = transfer8x8_copy_3dne;                  transfer8x8_copy = transfer8x8_copy_3dne;
# Line 480  Line 470 
470          }          }
471    
472    
473          if ((cpu_flags & XVID_CPU_SSE2) > 0) {          if ((cpu_flags & XVID_CPU_SSE2)) {
 #ifdef EXPERIMENTAL_SSE2_CODE  
474    
475                  calc_cbp = calc_cbp_sse2;                  calc_cbp = calc_cbp_sse2;
476    
# Line 491  Line 480 
480                  quant_inter   = quant_inter_sse2;                  quant_inter   = quant_inter_sse2;
481                  dequant_inter = dequant_inter_sse2;                  dequant_inter = dequant_inter_sse2;
482    
483                  /* ME */  #if defined(EXPERIMENTAL_SSE2_CODE)
484                    /* ME; slower than xmm */
485                  sad16    = sad16_sse2;                  sad16    = sad16_sse2;
486                  dev16    = dev16_sse2;                  dev16    = dev16_sse2;
487    #endif
488                  /* Forward and Inverse DCT */                  /* Forward and Inverse DCT */
489                  idct  = idct_sse2;                  idct  = idct_sse2;
490                  fdct = fdct_sse2;                  fdct = fdct_sse2;
 #endif  
491          }          }
   
492  #endif  #endif
493    
494  #ifdef ARCH_IA64  #if defined(ARCH_IS_IA64)
495          if ((cpu_flags & XVID_CPU_IA64) > 0) { //use assembler routines?          if ((cpu_flags & XVID_CPU_ASM)) { /* use assembler routines? */
496            idct_ia64_init();            idct_ia64_init();
497            fdct = fdct_ia64;            fdct = fdct_ia64;
498            idct = idct_ia64;   //not yet working, crashes            idct = idct_ia64;   /*not yet working, crashes */
499            interpolate8x8_halfpel_h = interpolate8x8_halfpel_h_ia64;            interpolate8x8_halfpel_h = interpolate8x8_halfpel_h_ia64;
500            interpolate8x8_halfpel_v = interpolate8x8_halfpel_v_ia64;            interpolate8x8_halfpel_v = interpolate8x8_halfpel_v_ia64;
501            interpolate8x8_halfpel_hv = interpolate8x8_halfpel_hv_ia64;            interpolate8x8_halfpel_hv = interpolate8x8_halfpel_hv_ia64;
# Line 515  Line 503 
503            sad16bi = sad16bi_ia64;            sad16bi = sad16bi_ia64;
504            sad8 = sad8_ia64;            sad8 = sad8_ia64;
505            dev16 = dev16_ia64;            dev16 = dev16_ia64;
506  //        Halfpel8_Refine = Halfpel8_Refine_ia64;  /*        Halfpel8_Refine = Halfpel8_Refine_ia64; */
507            quant_intra = quant_intra_ia64;            quant_intra = quant_intra_ia64;
508            dequant_intra = dequant_intra_ia64;            dequant_intra = dequant_intra_ia64;
509            quant_inter = quant_inter_ia64;            quant_inter = quant_inter_ia64;
# Line 526  Line 514 
514            transfer_8to16sub2 = transfer_8to16sub2_ia64;            transfer_8to16sub2 = transfer_8to16sub2_ia64;
515            transfer_16to8add = transfer_16to8add_ia64;            transfer_16to8add = transfer_16to8add_ia64;
516            transfer8x8_copy = transfer8x8_copy_ia64;            transfer8x8_copy = transfer8x8_copy_ia64;
517            DEBUG("Using IA-64 assembler routines.\n");            DPRINTF(DPRINTF_DEBUG, "Using IA-64 assembler routines.");
518          }          }
519  #endif  #endif
520    
521  #ifdef ARCH_PPC  #if defined(ARCH_IS_PPC)
522  #ifdef ARCH_PPC_ALTIVEC          if ((cpu_flags & XVID_CPU_ASM))
523            {
524                    calc_cbp = calc_cbp_ppc;
525            }
526    
527            if ((cpu_flags & XVID_CPU_ALTIVEC))
528            {
529          calc_cbp = calc_cbp_altivec;          calc_cbp = calc_cbp_altivec;
530          fdct = fdct_altivec;          fdct = fdct_altivec;
531          idct = idct_altivec;          idct = idct_altivec;
# Line 539  Line 533 
533          sad16 = sad16_altivec;          sad16 = sad16_altivec;
534          sad8 = sad8_altivec;          sad8 = sad8_altivec;
535          dev16 = dev16_altivec;          dev16 = dev16_altivec;
536  #else          }
         calc_cbp = calc_cbp_ppc;  
537  #endif  #endif
538    
539    #if defined(_DEBUG)
540        xvid_debug = init->debug;
541  #endif  #endif
542    
543          return XVID_ERR_OK;      return 0;
544  }  }
545    
546    
547    static int
548    xvid_gbl_info(xvid_gbl_info_t * info)
549    {
550            if (XVID_MAJOR(info->version) != 1) /* v1.x.x */
551                    return XVID_ERR_VERSION;
552    
553            info->actual_version = XVID_VERSION;
554            info->build = "dev-api-4";
555            info->cpu_flags = detect_cpu_flags();
556    
557    #if defined(_SMP) && defined(WIN32)
558            info->num_threads = pthread_num_processors_np();;
559    #else
560            info->num_threads = 0;
561    #endif
562    
563            return 0;
564    }
565    
566    
567  static int  static int
568  xvid_init_convert(XVID_INIT_CONVERTINFO* convert)  xvid_gbl_convert(xvid_gbl_convert_t* convert)
569  {  {
570          // const int flip1 = (convert->input.colorspace & XVID_CSP_VFLIP) ^ (convert->output.colorspace & XVID_CSP_VFLIP);          int width;
571          const int width = convert->width;          int height;
572          const int height = convert->height;          int width2;
573          const int width2 = convert->width/2;          int height2;
         const int height2 = convert->height/2;  
574          IMAGE img;          IMAGE img;
575    
576          switch (convert->input.colorspace & ~XVID_CSP_VFLIP)          if (XVID_MAJOR(convert->version) != 1)   /* v1.x.x */
577                  return XVID_ERR_VERSION;
578    
579            // const int flip1 = (convert->input.colorspace & XVID_CSP_VFLIP) ^ (convert->output.colorspace & XVID_CSP_VFLIP);
580            width = convert->width;
581            height = convert->height;
582            width2 = convert->width/2;
583            height2 = convert->height/2;
584    
585            switch (convert->input.csp & ~XVID_CSP_VFLIP)
586          {          {
587                  case XVID_CSP_YV12 :                  case XVID_CSP_YV12 :
588                          img.y = convert->input.y;                          img.y = convert->input.plane[0];
589                          img.v = (uint8_t*)convert->input.y + width*height;                          img.v = (uint8_t*)convert->input.plane[0] + convert->input.stride[0]*height;
590                          img.u = (uint8_t*)convert->input.y + width*height + width2*height2;                          img.u = (uint8_t*)convert->input.plane[0] + convert->input.stride[0]*height + (convert->input.stride[0]/2)*height2;
591                          image_output(&img, width, height, width,                          image_output(&img, width, height, width,
592                                                  convert->output.y, convert->output.y_stride,                                                  (uint8_t**)convert->output.plane, convert->output.stride,
593                                                  convert->output.colorspace, convert->interlacing);                                                  convert->output.csp, convert->interlacing);
594                          break;                          break;
595    
596                  default :                  default :
# Line 576  Line 599 
599    
600    
601          emms();          emms();
602          return XVID_ERR_OK;          return 0;
603  }  }
604    
605    
# Line 621  Line 644 
644          return 0;          return 0;
645  }  }
646    
647    int diff16(const int16_t * blockA, const int16_t * blockB, int size)
648    {
649            int i, diff = 0;
650            for (i = 0; i < size; i++)
651                    diff += abs(blockA[i]-blockB[i]);
652            return diff;
653    }
654    
655    
656    #define XVID_TEST_RANDOM        0x00000001      /* random input data */
657    #define XVID_TEST_VERBOSE       0x00000002      /* verbose error output */
658    
659  int test_h263_intra(quanth263_intraFunc * funcA, quanth263_intraFunc * funcB,  
660                                                  const char * nameA, const char * nameB,  #define TEST_FORWARD    0x00000001      /* intra */
661                                                  int min, int max)  #define TEST_FDCT  (TEST_FORWARD)
662    #define TEST_IDCT  (0)
663    
664    static int test_transform(void * funcA, void * funcB, const char * nameB,
665                                       int test, int flags)
666  {  {
667          int q,i;          int i;
668          int64_t timeSTART;          int64_t timeSTART;
669          int64_t timeA = 0;          int64_t timeA = 0;
670          int64_t timeB = 0;          int64_t timeB = 0;
         DECLARE_ALIGNED_MATRIX(arrayX, 1, 64, int16_t, CACHE_LINE);  
671          DECLARE_ALIGNED_MATRIX(arrayA, 1, 64, int16_t, CACHE_LINE);          DECLARE_ALIGNED_MATRIX(arrayA, 1, 64, int16_t, CACHE_LINE);
672          DECLARE_ALIGNED_MATRIX(arrayB, 1, 64, int16_t, CACHE_LINE);          DECLARE_ALIGNED_MATRIX(arrayB, 1, 64, int16_t, CACHE_LINE);
673            int min, max;
674            int count = 0;
675    
676          for (q = 1; q <= 31; q++)       /* quantizer */          int tmp;
677            int min_error = 0x10000*64;
678            int max_error = 0;
679    
680    
681            if ((test & TEST_FORWARD))      /* forward */
682          {          {
683                  for (i = min; i < max; i++)     /* input coeff */                  min = -256;
684                    max = 255;
685            }else{          /* inverse */
686                    min = -2048;
687                    max = 2047;
688            }
689    
690            for (i = 0; i < 64*64; i++)
691                  {                  {
692                          fill16(arrayX, 64, i);                  if ((flags & XVID_TEST_RANDOM))
693                    {
694                            random16(arrayA, 64, min, max);
695                    }else{
696                            fill16(arrayA, 64, i);
697                    }
698                    memcpy(arrayB, arrayA, 64*sizeof(int16_t));
699    
700                    if ((test & TEST_FORWARD))
701                    {
702                          timeSTART = read_counter();                          timeSTART = read_counter();
703                          funcA(arrayA, arrayX, q, q);                          ((fdctFunc*)funcA)(arrayA);
704                          timeA += read_counter() - timeSTART;                          timeA += read_counter() - timeSTART;
705    
706                          timeSTART = read_counter();                          timeSTART = read_counter();
707                          funcB(arrayB, arrayX, q, q);                          ((fdctFunc*)funcB)(arrayB);
708                          timeB += read_counter() - timeSTART;                          timeB += read_counter() - timeSTART;
   
                         if (compare16(arrayA, arrayB, 64))  
                         {  
                                 printf("%s/%s error: q=%i, i=%i\n", nameA?nameA:"?", nameB?nameB:"?", q, i);  
                                 return 0;  
709                          }                          }
710                    else
711                    {
712                            timeSTART = read_counter();
713                            ((idctFunc*)funcA)(arrayA);
714                            timeA += read_counter() - timeSTART;
715    
716                            timeSTART = read_counter();
717                            ((idctFunc*)funcB)(arrayB);
718                            timeB += read_counter() - timeSTART;
719                  }                  }
720    
721                    tmp = diff16(arrayA, arrayB, 64) / 64;
722                    if (tmp > max_error)
723                            max_error = tmp;
724                    if (tmp < min_error)
725                            min_error = tmp;
726    
727                    count++;
728          }          }
729    
730          if (nameA) printf("%s:\t%I64i\n", nameA, timeA);          /* print the "average difference" of best/worst transforms */
731          if (nameB) printf("%s:\t%I64i\n", nameB, timeB);          printf("%s:\t%i\t(min_error:%i, max_error:%i)\n", nameB, (int)(timeB / count), min_error, max_error);
732    
733          return 0;          return 0;
734  }  }
735    
736  int test_h263_inter(quanth263_interFunc * funcA, quanth263_interFunc * funcB,  
737                                                  const char * nameA, const char * nameB,  #define TEST_QUANT      0x00000001      /* forward quantization */
738                                                  int min, int max)  #define TEST_INTRA      0x00000002      /* intra */
739    #define TEST_QUANT_INTRA        (TEST_QUANT|TEST_INTRA)
740    #define TEST_QUANT_INTER        (TEST_QUANT)
741    #define TEST_DEQUANT_INTRA      (TEST_INTRA)
742    #define TEST_DEQUANT_INTER      (0)
743    
744    static int test_quant(void * funcA, void * funcB, const char * nameB,
745                               int test, int flags)
746  {  {
747          int q,i;          int q,i;
748          int64_t timeSTART;          int64_t timeSTART;
749          int64_t timeA = 0;          int64_t timeA = 0;
750          int64_t timeB = 0;          int64_t timeB = 0;
751            int retA = 0, retB = 0;
752          DECLARE_ALIGNED_MATRIX(arrayX, 1, 64, int16_t, CACHE_LINE);          DECLARE_ALIGNED_MATRIX(arrayX, 1, 64, int16_t, CACHE_LINE);
753          DECLARE_ALIGNED_MATRIX(arrayA, 1, 64, int16_t, CACHE_LINE);          DECLARE_ALIGNED_MATRIX(arrayA, 1, 64, int16_t, CACHE_LINE);
754          DECLARE_ALIGNED_MATRIX(arrayB, 1, 64, int16_t, CACHE_LINE);          DECLARE_ALIGNED_MATRIX(arrayB, 1, 64, int16_t, CACHE_LINE);
755            int min, max;
756            int count = 0;
757            int errors = 0;
758    
759            if ((test & TEST_QUANT))        /* quant */
760            {
761                    min = -2048;
762                    max = 2047;
763            }else{          /* dequant */
764                    min = -256;
765                    max = 255;
766            }
767    
768          for (q = 1; q <= 31; q++)       /* quantizer */          for (q = 1; q <= 31; q++)       /* quantizer */
769          {          {
770                  for (i = min; i < max; i++)     /* input coeff */                  for (i = min; i < max; i++)     /* input coeff */
771                  {                  {
772                            if ((flags & XVID_TEST_RANDOM))
773                            {
774                                    random16(arrayX, 64, min, max);
775                            }else{
776                          fill16(arrayX, 64, i);                          fill16(arrayX, 64, i);
777                            }
778    
779                            if ((test & TEST_INTRA))        /* intra */
780                            {
781                                    timeSTART = read_counter();
782                                    ((quanth263_intraFunc*)funcA)(arrayA, arrayX, q, q);
783                                    timeA += read_counter() - timeSTART;
784    
785                          timeSTART = read_counter();                          timeSTART = read_counter();
786                          funcA(arrayA, arrayX, q);                                  ((quanth263_intraFunc*)funcB)(arrayB, arrayX, q, q);
787                                    timeB += read_counter() - timeSTART;
788                            }
789                            else    /* inter */
790                            {
791                                    timeSTART = read_counter();
792                                    retA = ((quanth263_interFunc*)funcA)(arrayA, arrayX, q);
793                          timeA += read_counter() - timeSTART;                          timeA += read_counter() - timeSTART;
794    
795                          timeSTART = read_counter();                          timeSTART = read_counter();
796                          funcB(arrayB, arrayX, q);                                  retB = ((quanth263_interFunc*)funcB)(arrayB, arrayX, q);
797                          timeB += read_counter() - timeSTART;                          timeB += read_counter() - timeSTART;
798                            }
799    
800                          if (compare16(arrayA, arrayB, 64))                          /* compare return value from quant_inter, and compare (de)quantiz'd arrays */
801                            if ( ((test&TEST_QUANT) && !(test&TEST_INTRA) && retA != retB ) ||
802                                    compare16(arrayA, arrayB, 64))
803                          {                          {
804                                  printf("%s/%s error: q=%i, i=%i\n", nameA?nameA:"?", nameB?nameB:"?", q, i);                                  errors++;
805                                  return 0;                                  if ((flags & XVID_TEST_VERBOSE))
806                                            printf("%s error: q=%i, i=%i\n", nameB, q, i);
807                          }                          }
808    
809                            count++;
810                  }                  }
811          }          }
812    
813          if (nameA) printf("%s:\t%I64i\n", nameA, timeA);          printf("%s:\t%i", nameB, (int)(timeB / count));
814          if (nameB) printf("%s:\t%I64i\n", nameB, timeB);          if (errors>0)
815                    printf("\t(%i errors out of %i)", errors, count);
816            printf("\n");
817    
818          return 0;          return 0;
819  }  }
820    
821    
822    
823  int xvid_init_test()  int xvid_init_test(int flags)
824  {  {
825    #if defined(ARCH_IS_IA32)
826          int cpu_flags;          int cpu_flags;
827    #endif
828    
829            printf("XviD tests\n\n");
830    
831    #if defined(ARCH_IS_IA32)
832            cpu_flags = detect_cpu_flags();
833    #endif
834    
835            idct_int32_init();
836            emms();
837    
838            srand(time(0));
839    
840            /* fDCT test */
841            printf("--- fdct ---\n");
842                    test_transform(fdct_int32, fdct_int32, "c", TEST_FDCT, flags);
843    
844    #if defined(ARCH_IS_IA32)
845            if (cpu_flags & XVID_CPU_MMX)
846                    test_transform(fdct_int32, fdct_mmx, "mmx", TEST_FDCT, flags);
847            if (cpu_flags & XVID_CPU_SSE2)
848                    test_transform(fdct_int32, fdct_sse2, "sse2", TEST_FDCT, flags);
849    #endif
850    
851          printf("xvid_init_test\n");          /* iDCT test */
852            printf("\n--- idct ---\n");
853                    test_transform(idct_int32, idct_int32, "c", TEST_IDCT, flags);
854    
855  #if defined(ARCH_X86)  #if defined(ARCH_IS_IA32)
856          cpu_flags = check_cpu_features();          if (cpu_flags & XVID_CPU_MMX)
857                    test_transform(idct_int32, idct_mmx, "mmx", TEST_IDCT, flags);
858            if (cpu_flags & XVID_CPU_MMXEXT)
859                    test_transform(idct_int32, idct_xmm, "xmm", TEST_IDCT, flags);
860            if (cpu_flags & XVID_CPU_3DNOWEXT)
861                    test_transform(idct_int32, idct_3dne, "3dne", TEST_IDCT, flags);
862            if (cpu_flags & XVID_CPU_SSE2)
863                    test_transform(idct_int32, idct_sse2, "sse2", TEST_IDCT, flags);
864    #endif
865    
866          emms_mmx();          /* Intra quantization test */
867            printf("\n--- quant intra ---\n");
868                    test_quant(quant_intra_c, quant_intra_c, "c", TEST_QUANT_INTRA, flags);
869    
870          printf("--- quant intra ---\n");  #if defined(ARCH_IS_IA32)
871          if (cpu_flags & XVID_CPU_MMX)          if (cpu_flags & XVID_CPU_MMX)
872                  test_h263_intra(quant_intra_c, quant_intra_mmx, "c", "mmx", -2048, 2047);                  test_quant(quant_intra_c, quant_intra_mmx, "mmx", TEST_QUANT_INTRA, flags);
873          if (cpu_flags & XVID_CPU_3DNOWEXT)          if (cpu_flags & XVID_CPU_3DNOWEXT)
874                  test_h263_intra(quant_intra_c, quant_intra_3dne, NULL, "3dne", -2048, 2047);                  test_quant(quant_intra_c, quant_intra_3dne, "3dne", TEST_QUANT_INTRA, flags);
875          if (cpu_flags & XVID_CPU_SSE2)          if (cpu_flags & XVID_CPU_SSE2)
876                  test_h263_intra(quant_intra_c, quant_intra_sse2, NULL, "sse2", -2048, 2047);                  test_quant(quant_intra_c, quant_intra_sse2, "sse2", TEST_QUANT_INTRA, flags);
877    #endif
878    
879            /* Inter quantization test */
880          printf("\n--- quant inter ---\n");          printf("\n--- quant inter ---\n");
881                    test_quant(quant_inter_c, quant_inter_c, "c", TEST_QUANT_INTER, flags);
882    
883    #if defined(ARCH_IS_IA32)
884          if (cpu_flags & XVID_CPU_MMX)          if (cpu_flags & XVID_CPU_MMX)
885                  test_h263_inter(quant_inter_c, quant_inter_mmx, "c", "mmx", -2048, 2047);                  test_quant(quant_inter_c, quant_inter_mmx, "mmx", TEST_QUANT_INTER, flags);
886          if (cpu_flags & XVID_CPU_3DNOWEXT)          if (cpu_flags & XVID_CPU_3DNOWEXT)
887                  test_h263_inter(quant_inter_c, quant_inter_3dne, NULL, "3dne", -2048, 2047);                  test_quant(quant_inter_c, quant_inter_3dne, "3dne", TEST_QUANT_INTER, flags);
888          if (cpu_flags & XVID_CPU_SSE2)          if (cpu_flags & XVID_CPU_SSE2)
889                  test_h263_inter(quant_inter_c, quant_inter_sse2, NULL, "sse2", -2048, 2047);                  test_quant(quant_inter_c, quant_inter_sse2, "sse2", TEST_QUANT_INTER, flags);
890    #endif
891    
892          printf("\n--- dequan intra ---\n");          /* Intra dequantization test */
893            printf("\n--- dequant intra ---\n");
894                    test_quant(dequant_intra_c, dequant_intra_c, "c", TEST_DEQUANT_INTRA, flags);
895    
896    #if defined(ARCH_IS_IA32)
897          if (cpu_flags & XVID_CPU_MMX)          if (cpu_flags & XVID_CPU_MMX)
898                  test_h263_intra(dequant_intra_c, dequant_intra_mmx, "c", "mmx", -256, 255);                  test_quant(dequant_intra_c, dequant_intra_mmx, "mmx", TEST_DEQUANT_INTRA, flags);
899          if (cpu_flags & XVID_CPU_MMXEXT)          if (cpu_flags & XVID_CPU_MMXEXT)
900                  test_h263_intra(dequant_intra_c, dequant_intra_xmm, NULL, "xmm", -256, 255);                  test_quant(dequant_intra_c, dequant_intra_xmm, "xmm", TEST_DEQUANT_INTRA, flags);
901          if (cpu_flags & XVID_CPU_3DNOWEXT)          if (cpu_flags & XVID_CPU_3DNOWEXT)
902                  test_h263_intra(dequant_intra_c, dequant_intra_3dne, NULL, "3dne", -256, 255);                  test_quant(dequant_intra_c, dequant_intra_3dne, "3dne", TEST_DEQUANT_INTRA, flags);
903          if (cpu_flags & XVID_CPU_SSE2)          if (cpu_flags & XVID_CPU_SSE2)
904                  test_h263_intra(dequant_intra_c, dequant_intra_sse2, NULL, "sse2", -256, 255);                  test_quant(dequant_intra_c, dequant_intra_sse2, "sse2", TEST_DEQUANT_INTRA, flags);
905    #endif
906    
907            /* Inter dequantization test */
908          printf("\n--- dequant inter ---\n");          printf("\n--- dequant inter ---\n");
909          if (cpu_flags & XVID_CPU_MMX)                  test_quant(dequant_inter_c, dequant_inter_c, "c", TEST_DEQUANT_INTER, flags);
                 test_h263_inter((quanth263_interFunc*)dequant_inter_c,  
                                                 (quanth263_interFunc*)dequant_inter_mmx, "c", "mmx", -256, 255);  
910    
911    #if defined(ARCH_IS_IA32)
912            if (cpu_flags & XVID_CPU_MMX)
913                    test_quant(dequant_inter_c, dequant_inter_mmx, "mmx", TEST_DEQUANT_INTER, flags);
914          if (cpu_flags & XVID_CPU_MMXEXT)          if (cpu_flags & XVID_CPU_MMXEXT)
915                  test_h263_inter((quanth263_interFunc*)dequant_inter_c,                  test_quant(dequant_inter_c, dequant_inter_xmm, "xmm", TEST_DEQUANT_INTER, flags);
                                                 (quanth263_interFunc*)dequant_inter_xmm, NULL, "xmm", -256, 255);  
916          if (cpu_flags & XVID_CPU_3DNOWEXT)          if (cpu_flags & XVID_CPU_3DNOWEXT)
917                  test_h263_inter((quanth263_interFunc*)dequant_inter_c,                  test_quant(dequant_inter_c, dequant_inter_3dne, "3dne", TEST_DEQUANT_INTER, flags);
                                                 (quanth263_interFunc*)dequant_inter_3dne, NULL, "3dne", -256, 255);  
918          if (cpu_flags & XVID_CPU_SSE2)          if (cpu_flags & XVID_CPU_SSE2)
919                  test_h263_inter((quanth263_interFunc*)dequant_inter_c,                  test_quant(dequant_inter_c, dequant_inter_sse2, "sse2", TEST_DEQUANT_INTER, flags);
920                                                  (quanth263_interFunc*)dequant_inter_sse2, NULL, "sse2", -256, 255);  #endif
921    
922            /* Intra quantization test */
923            printf("\n--- quant4 intra ---\n");
924                    test_quant(quant4_intra_c, quant4_intra_c, "c", TEST_QUANT_INTRA, flags);
925    
926          printf("\n--- quant4_intra ---\n");  #if defined(ARCH_IS_IA32)
927          if (cpu_flags & XVID_CPU_MMX)          if (cpu_flags & XVID_CPU_MMX)
928                  test_h263_intra((quanth263_intraFunc*)quant4_intra_c,                  test_quant(quant4_intra_c, quant4_intra_mmx, "mmx", TEST_QUANT_INTRA, flags);
                                                 (quanth263_intraFunc*)quant4_intra_mmx, "c", "mmx", -2048, 2047);  
929          if (cpu_flags & XVID_CPU_MMXEXT)          if (cpu_flags & XVID_CPU_MMXEXT)
930                  test_h263_intra((quanth263_intraFunc*)quant4_intra_c,                  test_quant(quant4_intra_c, quant4_intra_xmm, "xmm", TEST_QUANT_INTRA, flags);
931                                                  (quanth263_intraFunc*)quant4_intra_xmm, NULL, "xmm", -2048, 2047);  #endif
932    
933            /* Inter quantization test */
934            printf("\n--- quant4 inter ---\n");
935                    test_quant(quant4_inter_c, quant4_inter_c, "c", TEST_QUANT_INTER, flags);
936    
937          printf("\n--- quant4_inter ---\n");  #if defined(ARCH_IS_IA32)
938          if (cpu_flags & XVID_CPU_MMX)          if (cpu_flags & XVID_CPU_MMX)
939                  test_h263_inter((quanth263_interFunc*)quant4_inter_c,                  test_quant(quant4_inter_c, quant4_inter_mmx, "mmx", TEST_QUANT_INTER, flags);
                                                 (quanth263_interFunc*)quant4_inter_mmx, "c", "mmx", -2048, 2047);  
940          if (cpu_flags & XVID_CPU_MMXEXT)          if (cpu_flags & XVID_CPU_MMXEXT)
941                  test_h263_inter((quanth263_interFunc*)quant4_inter_c,                  test_quant(quant4_inter_c, quant4_inter_xmm, "xmm", TEST_QUANT_INTER, flags);
942                                                  (quanth263_interFunc*)quant4_inter_xmm, NULL, "xmm", -2048, 2047);  #endif
943    
944            /* Intra dequantization test */
945            printf("\n--- dequant4 intra ---\n");
946                    test_quant(dequant4_intra_c, dequant4_intra_c, "c", TEST_DEQUANT_INTRA, flags);
947    
948          printf("\n--- dequant4_intra ---\n");  #if defined(ARCH_IS_IA32)
949          if (cpu_flags & XVID_CPU_MMX)          if (cpu_flags & XVID_CPU_MMX)
950                  test_h263_intra((quanth263_intraFunc*)dequant4_intra_c,                  test_quant(dequant4_intra_c, dequant4_intra_mmx, "mmx", TEST_DEQUANT_INTRA, flags);
                                                 (quanth263_intraFunc*)dequant4_intra_mmx, "c", "mmx", -256, 255);  
951          if (cpu_flags & XVID_CPU_3DNOWEXT)          if (cpu_flags & XVID_CPU_3DNOWEXT)
952                  test_h263_intra((quanth263_intraFunc*)dequant4_intra_c,                  test_quant(dequant4_intra_c, dequant4_intra_3dne, "3dne", TEST_DEQUANT_INTRA, flags);
953                                                  (quanth263_intraFunc*)dequant4_intra_3dne, NULL, "sse2", -256, 255);  #endif
954    
955            /* Inter dequantization test */
956            printf("\n--- dequant4 inter ---\n");
957                    test_quant(dequant4_inter_c, dequant4_inter_c, "c", TEST_DEQUANT_INTER, flags);
958    
959          printf("\n--- dequant4_inter ---\n");  #if defined(ARCH_IS_IA32)
960          if (cpu_flags & XVID_CPU_MMX)          if (cpu_flags & XVID_CPU_MMX)
961                  test_h263_inter((quanth263_interFunc*)dequant4_inter_c,                  test_quant(dequant4_inter_c, dequant4_inter_mmx, "mmx", TEST_DEQUANT_INTER, flags);
                                                 (quanth263_interFunc*)dequant4_inter_mmx, "c", "mmx", -256, 255);  
962          if (cpu_flags & XVID_CPU_3DNOWEXT)          if (cpu_flags & XVID_CPU_3DNOWEXT)
963                  test_h263_inter((quanth263_interFunc*)dequant4_inter_c,                  test_quant(dequant4_inter_c, dequant4_inter_3dne, "3dne", TEST_DEQUANT_INTER, flags);
                                                 (quanth263_interFunc*)dequant4_inter_3dne, NULL, "sse2", -256, 255);  
   
         emms_mmx();  
   
964  #endif  #endif
965    
966          return XVID_ERR_OK;          emms();
967    
968            return 0;
969  }  }
970    
971    
972    /*****************************************************************************
973     * XviD Global Entry point
974     *
975     * Well this function initialize all internal function pointers according
976     * to the CPU features forced by the library client or autodetected (depending
977     * on the XVID_CPU_FORCE flag). It also initializes vlc coding tables and all
978     * image colorspace transformation tables.
979     *
980     ****************************************************************************/
981    
982    
983  int  int
984  xvid_init(void *handle,  xvid_global(void *handle,
985                    int opt,                    int opt,
986                    void *param1,                    void *param1,
987                    void *param2)                    void *param2)
988  {  {
989          switch(opt)          switch(opt)
990          {          {
991                  case XVID_INIT_INIT :                  case XVID_GBL_INIT :
992                          return xvid_init_init((XVID_INIT_PARAM*)param1);                          return xvid_gbl_init((xvid_gbl_init_t*)param1);
993    
994                  case XVID_INIT_CONVERT :          case XVID_GBL_INFO :
995                          return xvid_init_convert((XVID_INIT_CONVERTINFO*)param1);              return xvid_gbl_info((xvid_gbl_info_t*)param1);
996    
997                  case XVID_INIT_TEST :                  case XVID_GBL_CONVERT :
998                          return xvid_init_test();                          return xvid_gbl_convert((xvid_gbl_convert_t*)param1);
999    
1000                    case XVID_GBL_TEST :
1001                    {
1002                            ptr_t flags = (ptr_t)param1;
1003                            return xvid_init_test((int)flags);
1004                    }
1005                  default :                  default :
1006                          return XVID_ERR_FAIL;                          return XVID_ERR_FAIL;
1007          }          }
# Line 837  Line 1024 
1024                          void *param2)                          void *param2)
1025  {  {
1026          switch (opt) {          switch (opt) {
         case XVID_DEC_DECODE:  
                 return decoder_decode((DECODER *) handle, (XVID_DEC_FRAME *) param1, (XVID_DEC_STATS*) param2);  
   
1027          case XVID_DEC_CREATE:          case XVID_DEC_CREATE:
1028                  return decoder_create((XVID_DEC_PARAM *) param1);                  return decoder_create((xvid_dec_create_t *) param1);
1029    
1030          case XVID_DEC_DESTROY:          case XVID_DEC_DESTROY:
1031                  return decoder_destroy((DECODER *) handle);                  return decoder_destroy((DECODER *) handle);
1032    
1033            case XVID_DEC_DECODE:
1034                    return decoder_decode((DECODER *) handle, (xvid_dec_frame_t *) param1, (xvid_dec_stats_t*) param2);
1035    
1036          default:          default:
1037                  return XVID_ERR_FAIL;                  return XVID_ERR_FAIL;
1038          }          }
# Line 871  Line 1058 
1058          switch (opt) {          switch (opt) {
1059          case XVID_ENC_ENCODE:          case XVID_ENC_ENCODE:
1060    
1061                  if (((Encoder *) handle)->mbParam.max_bframes >= 0)                  return enc_encode((Encoder *) handle,
1062                  return encoder_encode_bframes((Encoder *) handle, (XVID_ENC_FRAME *) param1,                                                            (xvid_enc_frame_t *) param1,
1063                                                            (XVID_ENC_STATS *) param2);                                                            (xvid_enc_stats_t *) param2);
                 else  
                 return encoder_encode((Encoder *) handle, (XVID_ENC_FRAME *) param1,  
                                                           (XVID_ENC_STATS *) param2);  
1064    
1065          case XVID_ENC_CREATE:          case XVID_ENC_CREATE:
1066                  return encoder_create((XVID_ENC_PARAM *) param1);                  return enc_create((xvid_enc_create_t *) param1);
1067    
1068          case XVID_ENC_DESTROY:          case XVID_ENC_DESTROY:
1069                  return encoder_destroy((Encoder *) handle);                  return enc_destroy((Encoder *) handle);
1070    
1071          default:          default:
1072                  return XVID_ERR_FAIL;                  return XVID_ERR_FAIL;
Legend:
Removed from v.1.33.2.19  
changed lines
  Added in v.1.45.2.4
No admin address has been configured
 ViewVC Help
Powered by ViewVC 1.0.4