[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.23, Wed Feb 12 11:46:18 2003 UTC revision 1.36, Wed Sep 4 22:08:45 2002 UTC
# Line 3  Line 3 
3   *  XVID MPEG-4 VIDEO CODEC   *  XVID MPEG-4 VIDEO CODEC
4   *  - Native API implementation  -   *  - Native API implementation  -
5   *   *
6     *  Copyright(C) 2001-2002 Peter Ross <pross@cs.rmit.edu.au>
7     *
8   *  This program is an implementation of a part of one or more MPEG-4   *  This program is an implementation of a part of one or more MPEG-4
9   *  Video tools as specified in ISO/IEC 14496-2 standard.  Those intending   *  Video tools as specified in ISO/IEC 14496-2 standard.  Those intending
10   *  to use this software module in hardware or software products are   *  to use this software module in hardware or software products are
# Line 26  Line 28 
28   *  along with this program ; if not, write to the Free Software   *  along with this program ; if not, write to the Free Software
29   *  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA   *  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
30   *   *
  ****************************************************************************/  
   
 /*****************************************************************************  
  *  
  *  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>  
  *  
31   *  $Id$   *  $Id$
32   *   *
33   ****************************************************************************/   ****************************************************************************/
34    
 #include <stdio.h>  
 #include <stdlib.h>  
 #include <string.h>  
 #include <time.h>  
   
35  #include "xvid.h"  #include "xvid.h"
36  #include "decoder.h"  #include "decoder.h"
37  #include "encoder.h"  #include "encoder.h"
# Line 54  Line 40 
40  #include "dct/fdct.h"  #include "dct/fdct.h"
41  #include "image/colorspace.h"  #include "image/colorspace.h"
42  #include "image/interpolate8x8.h"  #include "image/interpolate8x8.h"
 #include "image/reduced.h"  
43  #include "utils/mem_transfer.h"  #include "utils/mem_transfer.h"
 #include "utils/mbfunctions.h"  
44  #include "quant/quant_h263.h"  #include "quant/quant_h263.h"
45  #include "quant/quant_mpeg4.h"  #include "quant/quant_mpeg4.h"
46  #include "motion/motion.h"  #include "motion/motion.h"
# Line 65  Line 49 
49  #include "utils/timer.h"  #include "utils/timer.h"
50  #include "bitstream/mbcoding.h"  #include "bitstream/mbcoding.h"
51    
52  #if defined(ARCH_X86)  #if defined(ARCH_X86) && defined(EXPERIMENTAL_SSE2_CODE)
53    
54  #if defined(_MSC_VER)  #ifdef WIN32
55  #       include <windows.h>  #       include <windows.h>
56  #else  #else
57  #       include <signal.h>  #       include <signal.h>
58  #       include <setjmp.h>  #       include <setjmp.h>
59    #endif
60    
61    
62    #ifndef WIN32
63    
64          static jmp_buf mark;          static jmp_buf mark;
65    
# Line 84  Line 72 
72    
73    
74  /*  /*
75  calls the funcptr, and returns whether SIGILL (illegal instruction) was signalled   * Calls the funcptr, and returns whether SIGILL (illegal instruction) was signalled
76  return values:   * Return values:
77  -1 : could not determine   * -1 : could not determine
78  0  : SIGILL was *not* signalled   * 0  : SIGILL was *not* signalled
79  1  : SIGILL was signalled   * 1  : SIGILL was signalled
80  */  */
81    
82  int  int
83  sigill_check(void (*func)())  sigill_check(void (*func)())
84  {  {
85  #if defined(_MSC_VER)  #ifdef WIN32
86          _try {          _try {
87                  func();                  func();
88          }          }
# Line 128  Line 116 
116  }  }
117  #endif  #endif
118    
   
 /* detect cpu flags  */  
 static unsigned int  
 detect_cpu_flags()  
 {  
         /* enable native assembly optimizations by default */  
         unsigned int cpu_flags = XVID_CPU_ASM;  
   
 #if defined(ARCH_X86)  
         cpu_flags |= check_cpu_features();  
         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 defined(ARCH_PPC)  
 #if defined(ARCH_PPC_ALTIVEC)  
         cpu_flags |= XVID_CPU_ALTIVEC;  
 #endif  
 #endif  
   
         return cpu_flags;  
 }  
   
   
119  /*****************************************************************************  /*****************************************************************************
120   * XviD Init Entry point   * XviD Init Entry point
121   *   *
# Line 169  Line 130 
130   *   *
131   ****************************************************************************/   ****************************************************************************/
132    
133    int
134  static  xvid_init(void *handle,
135  int xvid_init_init(XVID_INIT_PARAM * init_param)                    int opt,
136                      void *param1,
137                      void *param2)
138  {  {
139          int cpu_flags;          int cpu_flags;
140            XVID_INIT_PARAM *init_param;
141    
142            init_param = (XVID_INIT_PARAM *) param1;
143    
144          /* Inform the client the API version */          /* Inform the client the API version */
145          init_param->api_version = API_VERSION;          init_param->api_version = API_VERSION;
# Line 188  Line 154 
154    
155          } else {          } else {
156    
157                  cpu_flags = detect_cpu_flags();                  cpu_flags = check_cpu_features();
158    
159    #if defined(ARCH_X86) && defined(EXPERIMENTAL_SSE2_CODE)
160                    if ((cpu_flags & XVID_CPU_SSE) && sigill_check(sse_os_trigger))
161                            cpu_flags &= ~XVID_CPU_SSE;
162    
163                    if ((cpu_flags & XVID_CPU_SSE2) && sigill_check(sse2_os_trigger))
164                            cpu_flags &= ~XVID_CPU_SSE2;
165    #endif
166          }          }
167    
168          if ((init_param->cpu_flags & XVID_CPU_CHKONLY))          if ((init_param->cpu_flags & XVID_CPU_CHKONLY))
# Line 229  Line 203 
203          transfer_8to16copy = transfer_8to16copy_c;          transfer_8to16copy = transfer_8to16copy_c;
204          transfer_16to8copy = transfer_16to8copy_c;          transfer_16to8copy = transfer_16to8copy_c;
205          transfer_8to16sub  = transfer_8to16sub_c;          transfer_8to16sub  = transfer_8to16sub_c;
         transfer_8to16subro  = transfer_8to16subro_c;  
206          transfer_8to16sub2 = transfer_8to16sub2_c;          transfer_8to16sub2 = transfer_8to16sub2_c;
207          transfer_16to8add  = transfer_16to8add_c;          transfer_16to8add  = transfer_16to8add_c;
208          transfer8x8_copy   = transfer8x8_copy_c;          transfer8x8_copy   = transfer8x8_copy_c;
209    
         /* Interlacing functions */  
         MBFieldTest = MBFieldTest_c;  
   
210          /* Image interpolation related functions */          /* Image interpolation related functions */
211          interpolate8x8_halfpel_h  = interpolate8x8_halfpel_h_c;          interpolate8x8_halfpel_h  = interpolate8x8_halfpel_h_c;
212          interpolate8x8_halfpel_v  = interpolate8x8_halfpel_v_c;          interpolate8x8_halfpel_v  = interpolate8x8_halfpel_v_c;
213          interpolate8x8_halfpel_hv = interpolate8x8_halfpel_hv_c;          interpolate8x8_halfpel_hv = interpolate8x8_halfpel_hv_c;
214    
         interpolate16x16_lowpass_h = interpolate16x16_lowpass_h_c;  
         interpolate16x16_lowpass_v = interpolate16x16_lowpass_v_c;  
         interpolate16x16_lowpass_hv = interpolate16x16_lowpass_hv_c;  
   
         interpolate8x8_lowpass_h = interpolate8x8_lowpass_h_c;  
         interpolate8x8_lowpass_v = interpolate8x8_lowpass_v_c;  
         interpolate8x8_lowpass_hv = interpolate8x8_lowpass_hv_c;  
   
         interpolate8x8_6tap_lowpass_h = interpolate8x8_6tap_lowpass_h_c;  
         interpolate8x8_6tap_lowpass_v = interpolate8x8_6tap_lowpass_v_c;  
   
         interpolate8x8_avg2 = interpolate8x8_avg2_c;  
         interpolate8x8_avg4 = interpolate8x8_avg4_c;  
   
         /* reduced resoltuion */  
         copy_upsampled_8x8_16to8 = xvid_Copy_Upsampled_8x8_16To8_C;  
         add_upsampled_8x8_16to8 = xvid_Add_Upsampled_8x8_16To8_C;  
         vfilter_31 = xvid_VFilter_31_C;  
         hfilter_31 = xvid_HFilter_31_C;  
         filter_18x18_to_8x8 = xvid_Filter_18x18_To_8x8_C;  
         filter_diff_18x18_to_8x8 = xvid_Filter_Diff_18x18_To_8x8_C;  
   
215          /* Initialize internal colorspace transformation tables */          /* Initialize internal colorspace transformation tables */
216          colorspace_init();          colorspace_init();
217    
218          /* All colorspace transformation functions User Format->YV12 */          /* All colorspace transformation functions User Format->YV12 */
         yv12_to_yv12    = yv12_to_yv12_c;  
219          rgb555_to_yv12  = rgb555_to_yv12_c;          rgb555_to_yv12  = rgb555_to_yv12_c;
220          rgb565_to_yv12  = rgb565_to_yv12_c;          rgb565_to_yv12  = rgb565_to_yv12_c;
221          bgr_to_yv12     = bgr_to_yv12_c;          rgb24_to_yv12  = rgb24_to_yv12_c;
222          bgra_to_yv12    = bgra_to_yv12_c;          rgb32_to_yv12  = rgb32_to_yv12_c;
223          abgr_to_yv12    = abgr_to_yv12_c;          yuv_to_yv12    = yuv_to_yv12_c;
         rgba_to_yv12    = rgba_to_yv12_c;  
224          yuyv_to_yv12    = yuyv_to_yv12_c;          yuyv_to_yv12    = yuyv_to_yv12_c;
225          uyvy_to_yv12    = uyvy_to_yv12_c;          uyvy_to_yv12    = uyvy_to_yv12_c;
226    
         rgb555i_to_yv12 = rgb555i_to_yv12_c;  
         rgb565i_to_yv12 = rgb565i_to_yv12_c;  
         bgri_to_yv12    = bgri_to_yv12_c;  
         bgrai_to_yv12   = bgrai_to_yv12_c;  
         abgri_to_yv12   = abgri_to_yv12_c;  
         rgbai_to_yv12   = rgbai_to_yv12_c;  
         yuyvi_to_yv12   = yuyvi_to_yv12_c;  
         uyvyi_to_yv12   = uyvyi_to_yv12_c;  
   
   
227          /* All colorspace transformation functions YV12->User format */          /* All colorspace transformation functions YV12->User format */
228          yv12_to_rgb555  = yv12_to_rgb555_c;          yv12_to_rgb555  = yv12_to_rgb555_c;
229          yv12_to_rgb565  = yv12_to_rgb565_c;          yv12_to_rgb565  = yv12_to_rgb565_c;
230          yv12_to_bgr     = yv12_to_bgr_c;          yv12_to_rgb24  = yv12_to_rgb24_c;
231          yv12_to_bgra    = yv12_to_bgra_c;          yv12_to_rgb32  = yv12_to_rgb32_c;
232          yv12_to_abgr    = yv12_to_abgr_c;          yv12_to_yuv    = yv12_to_yuv_c;
         yv12_to_rgba    = yv12_to_rgba_c;  
233          yv12_to_yuyv    = yv12_to_yuyv_c;          yv12_to_yuyv    = yv12_to_yuyv_c;
234          yv12_to_uyvy    = yv12_to_uyvy_c;          yv12_to_uyvy    = yv12_to_uyvy_c;
235    
         yv12_to_rgb555i = yv12_to_rgb555i_c;  
         yv12_to_rgb565i = yv12_to_rgb565i_c;  
         yv12_to_bgri    = yv12_to_bgri_c;  
         yv12_to_bgrai   = yv12_to_bgrai_c;  
         yv12_to_abgri   = yv12_to_abgri_c;  
         yv12_to_rgbai   = yv12_to_rgbai_c;  
         yv12_to_yuyvi   = yv12_to_yuyvi_c;  
         yv12_to_uyvyi   = yv12_to_uyvyi_c;  
   
236          /* Functions used in motion estimation algorithms */          /* Functions used in motion estimation algorithms */
237          calc_cbp = calc_cbp_c;          calc_cbp = calc_cbp_c;
238          sad16    = sad16_c;          sad16    = sad16_c;
# Line 314  Line 240 
240          sad16bi  = sad16bi_c;          sad16bi  = sad16bi_c;
241          sad8bi   = sad8bi_c;          sad8bi   = sad8bi_c;
242          dev16    = dev16_c;          dev16    = dev16_c;
         sad16v   = sad16v_c;  
243    
244  //      Halfpel8_Refine = Halfpel8_Refine_c;          Halfpel8_Refine = Halfpel8_Refine_c;
   
 #if defined(ARCH_X86)  
   
         if ((cpu_flags & XVID_CPU_ASM))  
         {  
                 vfilter_31 = xvid_VFilter_31_x86;  
                 hfilter_31 = xvid_HFilter_31_x86;  
         }  
   
         if ((cpu_flags & XVID_CPU_MMX) || (cpu_flags & XVID_CPU_MMXEXT) ||  
                 (cpu_flags & XVID_CPU_3DNOW) || (cpu_flags & XVID_CPU_3DNOWEXT) ||  
                 (cpu_flags & XVID_CPU_SSE) || (cpu_flags & XVID_CPU_SSE2))  
         {  
                 /* Restore FPU context : emms_c is a nop functions */  
                 emms = emms_mmx;  
         }  
245    
246          if ((cpu_flags & XVID_CPU_MMX)) {  #ifdef ARCH_X86
247            if ((cpu_flags & XVID_CPU_MMX) > 0) {
248    
249                  /* Forward and Inverse Discrete Cosine Transformation functions */                  /* Forward and Inverse Discrete Cosine Transformation functions */
250                  fdct = fdct_mmx;                  fdct = fdct_mmx;
251                  idct = idct_mmx;                  idct = idct_mmx;
252    
253                    /* To restore FPU context after mmx use */
254                    emms = emms_mmx;
255    
256                  /* Quantization related functions */                  /* Quantization related functions */
257                  quant_intra   = quant_intra_mmx;                  quant_intra   = quant_intra_mmx;
258                  dequant_intra = dequant_intra_mmx;                  dequant_intra = dequant_intra_mmx;
# Line 355  Line 268 
268                  transfer_8to16copy = transfer_8to16copy_mmx;                  transfer_8to16copy = transfer_8to16copy_mmx;
269                  transfer_16to8copy = transfer_16to8copy_mmx;                  transfer_16to8copy = transfer_16to8copy_mmx;
270                  transfer_8to16sub  = transfer_8to16sub_mmx;                  transfer_8to16sub  = transfer_8to16sub_mmx;
                 transfer_8to16subro  = transfer_8to16subro_mmx;  
271                  transfer_8to16sub2 = transfer_8to16sub2_mmx;                  transfer_8to16sub2 = transfer_8to16sub2_mmx;
272                  transfer_16to8add  = transfer_16to8add_mmx;                  transfer_16to8add  = transfer_16to8add_mmx;
273                  transfer8x8_copy   = transfer8x8_copy_mmx;                  transfer8x8_copy   = transfer8x8_copy_mmx;
274    
                 /* Interlacing Functions */  
                 MBFieldTest = MBFieldTest_mmx;  
275    
276                  /* Image Interpolation related functions */                  /* Image Interpolation related functions */
277                  interpolate8x8_halfpel_h  = interpolate8x8_halfpel_h_mmx;                  interpolate8x8_halfpel_h  = interpolate8x8_halfpel_h_mmx;
278                  interpolate8x8_halfpel_v  = interpolate8x8_halfpel_v_mmx;                  interpolate8x8_halfpel_v  = interpolate8x8_halfpel_v_mmx;
279                  interpolate8x8_halfpel_hv = interpolate8x8_halfpel_hv_mmx;                  interpolate8x8_halfpel_hv = interpolate8x8_halfpel_hv_mmx;
280    
281                  interpolate8x8_6tap_lowpass_h = interpolate8x8_6tap_lowpass_h_mmx;                  /* Image RGB->YV12 related functions */
282                  interpolate8x8_6tap_lowpass_v = interpolate8x8_6tap_lowpass_v_mmx;                  rgb24_to_yv12 = rgb24_to_yv12_mmx;
283                    rgb32_to_yv12 = rgb32_to_yv12_mmx;
284                  interpolate8x8_avg2 = interpolate8x8_avg2_mmx;                  yuv_to_yv12   = yuv_to_yv12_mmx;
                 interpolate8x8_avg4 = interpolate8x8_avg4_mmx;  
   
                 /* reduced resolution */  
                 copy_upsampled_8x8_16to8 = xvid_Copy_Upsampled_8x8_16To8_mmx;  
                 add_upsampled_8x8_16to8 = xvid_Add_Upsampled_8x8_16To8_mmx;  
                 hfilter_31 = xvid_HFilter_31_mmx;  
                 filter_18x18_to_8x8 = xvid_Filter_18x18_To_8x8_mmx;  
                 filter_diff_18x18_to_8x8 = xvid_Filter_Diff_18x18_To_8x8_mmx;  
   
                 /* image input xxx_to_yv12 related functions */  
                 yv12_to_yv12  = yv12_to_yv12_mmx;  
                 bgr_to_yv12   = bgr_to_yv12_mmx;  
                 bgra_to_yv12  = bgra_to_yv12_mmx;  
285                  yuyv_to_yv12  = yuyv_to_yv12_mmx;                  yuyv_to_yv12  = yuyv_to_yv12_mmx;
286                  uyvy_to_yv12  = uyvy_to_yv12_mmx;                  uyvy_to_yv12  = uyvy_to_yv12_mmx;
287    
288                  /* image output yv12_to_xxx related functions */                  /* Image YV12->RGB related functions */
289                  yv12_to_bgr   = yv12_to_bgr_mmx;                  yv12_to_rgb24 = yv12_to_rgb24_mmx;
290                  yv12_to_bgra  = yv12_to_bgra_mmx;                  yv12_to_rgb32 = yv12_to_rgb32_mmx;
291                  yv12_to_yuyv  = yv12_to_yuyv_mmx;                  yv12_to_yuyv  = yv12_to_yuyv_mmx;
292                  yv12_to_uyvy  = yv12_to_uyvy_mmx;                  yv12_to_uyvy  = yv12_to_uyvy_mmx;
293    
                 yv12_to_yuyvi = yv12_to_yuyvi_mmx;  
                 yv12_to_uyvyi = yv12_to_uyvyi_mmx;  
   
294                  /* Motion estimation related functions */                  /* Motion estimation related functions */
295                  calc_cbp = calc_cbp_mmx;                  calc_cbp = calc_cbp_mmx;
296                  sad16    = sad16_mmx;                  sad16    = sad16_mmx;
# Line 404  Line 298 
298                  sad16bi = sad16bi_mmx;                  sad16bi = sad16bi_mmx;
299                  sad8bi  = sad8bi_mmx;                  sad8bi  = sad8bi_mmx;
300                  dev16    = dev16_mmx;                  dev16    = dev16_mmx;
301                  sad16v   = sad16v_mmx;  
302          }          }
303    
304          /* these 3dnow functions are faster than mmx, but slower than xmm. */          /* these 3dnow functions are faster than mmx, but slower than xmm. */
305          if ((cpu_flags & XVID_CPU_3DNOW)) {          if ((cpu_flags & XVID_CPU_3DNOW) > 0) {
   
                 emms = emms_3dn;  
306    
307                  /* ME functions */                  /* ME functions */
308                  sad16bi = sad16bi_3dn;                  sad16bi = sad16bi_3dn;
309                  sad8bi  = sad8bi_3dn;                  sad8bi  = sad8bi_3dn;
   
                 yuyv_to_yv12  = yuyv_to_yv12_3dn;  
                 uyvy_to_yv12  = uyvy_to_yv12_3dn;  
310          }          }
311    
312    
313          if ((cpu_flags & XVID_CPU_MMXEXT)) {          if ((cpu_flags & XVID_CPU_MMXEXT) > 0) {
314    
315                  /* Inverse DCT */                  /* Inverse DCT */
316                  idct = idct_xmm;                  idct = idct_xmm;
# Line 431  Line 320 
320                  interpolate8x8_halfpel_v  = interpolate8x8_halfpel_v_xmm;                  interpolate8x8_halfpel_v  = interpolate8x8_halfpel_v_xmm;
321                  interpolate8x8_halfpel_hv = interpolate8x8_halfpel_hv_xmm;                  interpolate8x8_halfpel_hv = interpolate8x8_halfpel_hv_xmm;
322    
                 /* reduced resolution */  
                 copy_upsampled_8x8_16to8 = xvid_Copy_Upsampled_8x8_16To8_xmm;  
                 add_upsampled_8x8_16to8 = xvid_Add_Upsampled_8x8_16To8_xmm;  
   
323                  /* Quantization */                  /* Quantization */
                 quant4_intra = quant4_intra_xmm;  
                 quant4_inter = quant4_inter_xmm;  
   
324                  dequant_intra = dequant_intra_xmm;                  dequant_intra = dequant_intra_xmm;
325                  dequant_inter = dequant_inter_xmm;                  dequant_inter = dequant_inter_xmm;
326    
# Line 446  Line 328 
328                  transfer_8to16sub2 = transfer_8to16sub2_xmm;                  transfer_8to16sub2 = transfer_8to16sub2_xmm;
329    
330                  /* Colorspace transformation */                  /* Colorspace transformation */
331                  yv12_to_yv12  = yv12_to_yv12_xmm;                  yuv_to_yv12 = yuv_to_yv12_xmm;
                 yuyv_to_yv12  = yuyv_to_yv12_xmm;  
                 uyvy_to_yv12  = uyvy_to_yv12_xmm;  
332    
333                  /* ME functions */                  /* ME functions */
334                  sad16 = sad16_xmm;                  sad16 = sad16_xmm;
# Line 456  Line 336 
336                  sad16bi = sad16bi_xmm;                  sad16bi = sad16bi_xmm;
337                  sad8bi  = sad8bi_xmm;                  sad8bi  = sad8bi_xmm;
338                  dev16 = dev16_xmm;                  dev16 = dev16_xmm;
339                  sad16v   = sad16v_xmm;  
340          }          }
341    
342          if ((cpu_flags & XVID_CPU_3DNOW)) {          if ((cpu_flags & XVID_CPU_3DNOW) > 0) {
343    
344                  /* Interpolation */                  /* Interpolation */
345                  interpolate8x8_halfpel_h = interpolate8x8_halfpel_h_3dn;                  interpolate8x8_halfpel_h = interpolate8x8_halfpel_h_3dn;
# Line 467  Line 347 
347                  interpolate8x8_halfpel_hv = interpolate8x8_halfpel_hv_3dn;                  interpolate8x8_halfpel_hv = interpolate8x8_halfpel_hv_3dn;
348          }          }
349    
350          if ((cpu_flags & XVID_CPU_3DNOWEXT)) {          if ((cpu_flags & XVID_CPU_SSE2) > 0) {
351    #ifdef EXPERIMENTAL_SSE2_CODE
                 /* Inverse DCT */  
                 idct =  idct_3dne;  
   
                 /* Buffer transfer */  
                 transfer_8to16copy =  transfer_8to16copy_3dne;  
                 transfer_16to8copy = transfer_16to8copy_3dne;  
                 transfer_8to16sub =  transfer_8to16sub_3dne;  
                 transfer_8to16subro =  transfer_8to16subro_3dne;  
                 transfer_8to16sub2 =  transfer_8to16sub2_3dne;  
                 transfer_16to8add = transfer_16to8add_3dne;  
                 transfer8x8_copy = transfer8x8_copy_3dne;  
   
                 /* Quantization */  
                 dequant4_intra = dequant4_intra_3dne;  
                 dequant4_inter = dequant4_inter_3dne;  
                 quant_intra = quant_intra_3dne;  
                 quant_inter = quant_inter_3dne;  
                 dequant_intra = dequant_intra_3dne;  
                 dequant_inter = dequant_inter_3dne;  
   
                 /* ME functions */  
                 calc_cbp = calc_cbp_3dne;  
                 sad16 = sad16_3dne;  
                 sad8 = sad8_3dne;  
                 sad16bi = sad16bi_3dne;  
                 sad8bi = sad8bi_3dne;  
                 dev16 = dev16_3dne;  
   
                 /* Interpolation */  
                 interpolate8x8_halfpel_h = interpolate8x8_halfpel_h_3dne;  
                 interpolate8x8_halfpel_v = interpolate8x8_halfpel_v_3dne;  
                 interpolate8x8_halfpel_hv = interpolate8x8_halfpel_hv_3dne;  
         }  
   
   
         if ((cpu_flags & XVID_CPU_SSE2)) {  
352    
353                  calc_cbp = calc_cbp_sse2;                  calc_cbp = calc_cbp_sse2;
354    
# Line 514  Line 358 
358                  quant_inter   = quant_inter_sse2;                  quant_inter   = quant_inter_sse2;
359                  dequant_inter = dequant_inter_sse2;                  dequant_inter = dequant_inter_sse2;
360    
361  #if defined(EXPERIMENTAL_SSE2_CODE)                  /* ME */
                 /* ME; slower than xmm */  
362                  sad16    = sad16_sse2;                  sad16    = sad16_sse2;
363                  dev16    = dev16_sse2;                  dev16    = dev16_sse2;
364  #endif  
365                  /* Forward and Inverse DCT */                  /* Forward and Inverse DCT */
366                  idct  = idct_sse2;                  idct  = idct_sse2;
367                  fdct = fdct_sse2;                  fdct = fdct_sse2;
368    #endif
369          }          }
370    
371  #endif  #endif
372    
373  #if defined(ARCH_IA64)  #ifdef ARCH_IA64
374          if ((cpu_flags & XVID_CPU_ASM)) { //use assembler routines?          if ((cpu_flags & XVID_CPU_IA64) > 0) { //use assembler routines?
375            idct_ia64_init();            idct_ia64_init();
376            fdct = fdct_ia64;            fdct = fdct_ia64;
377            idct = idct_ia64;   //not yet working, crashes            idct = idct_ia64;   //not yet working, crashes
# Line 537  Line 382 
382            sad16bi = sad16bi_ia64;            sad16bi = sad16bi_ia64;
383            sad8 = sad8_ia64;            sad8 = sad8_ia64;
384            dev16 = dev16_ia64;            dev16 = dev16_ia64;
385  //        Halfpel8_Refine = Halfpel8_Refine_ia64;            Halfpel8_Refine = Halfpel8_Refine_ia64;
386            quant_intra = quant_intra_ia64;            quant_intra = quant_intra_ia64;
387            dequant_intra = dequant_intra_ia64;            dequant_intra = dequant_intra_ia64;
388            quant_inter = quant_inter_ia64;            quant_inter = quant_inter_ia64;
# Line 552  Line 397 
397          }          }
398  #endif  #endif
399    
400  #if defined(ARCH_PPC)  #ifdef ARCH_PPC
401          if ((cpu_flags & XVID_CPU_ASM))  #ifdef ARCH_PPC_ALTIVEC
         {  
                 calc_cbp = calc_cbp_ppc;  
         }  
   
         if ((cpu_flags & XVID_CPU_ALTIVEC))  
         {  
402                  calc_cbp = calc_cbp_altivec;                  calc_cbp = calc_cbp_altivec;
403                  fdct = fdct_altivec;                  fdct = fdct_altivec;
404                  idct = idct_altivec;                  idct = idct_altivec;
# Line 567  Line 406 
406                  sad16 = sad16_altivec;                  sad16 = sad16_altivec;
407                  sad8 = sad8_altivec;                  sad8 = sad8_altivec;
408                  dev16 = dev16_altivec;                  dev16 = dev16_altivec;
409          }  #else
410            calc_cbp = calc_cbp_ppc;
411  #endif  #endif
   
         return XVID_ERR_OK;  
 }  
   
   
   
 static int  
 xvid_init_convert(XVID_INIT_CONVERTINFO* convert)  
 {  
         // const int flip1 = (convert->input.colorspace & XVID_CSP_VFLIP) ^ (convert->output.colorspace & XVID_CSP_VFLIP);  
         const int width = convert->width;  
         const int height = convert->height;  
         const int width2 = convert->width/2;  
         const int height2 = convert->height/2;  
         IMAGE img;  
   
         switch (convert->input.colorspace & ~XVID_CSP_VFLIP)  
         {  
                 case XVID_CSP_YV12 :  
                         img.y = convert->input.y;  
                         img.v = (uint8_t*)convert->input.y + width*height;  
                         img.u = (uint8_t*)convert->input.y + width*height + width2*height2;  
                         image_output(&img, width, height, width,  
                                                 convert->output.y, convert->output.y_stride,  
                                                 convert->output.colorspace, convert->interlacing);  
                         break;  
   
                 default :  
                         return XVID_ERR_FORMAT;  
         }  
   
   
         emms();  
         return XVID_ERR_OK;  
 }  
   
   
   
 void fill8(uint8_t * block, int size, int value)  
 {  
         int i;  
         for (i = 0; i < size; i++)  
                 block[i] = value;  
 }  
   
 void fill16(int16_t * block, int size, int value)  
 {  
         int i;  
         for (i = 0; i < size; i++)  
                 block[i] = value;  
 }  
   
 #define RANDOM(min,max) min + (rand() % (max-min))  
   
 void random8(uint8_t * block, int size, int min, int max)  
 {  
         int i;  
         for (i = 0; i < size; i++)  
                 block[i] = RANDOM(min,max);  
 }  
   
 void random16(int16_t * block, int size, int min, int max)  
 {  
         int i;  
         for (i = 0; i < size; i++)  
                 block[i] = RANDOM(min,max);  
 }  
   
 int compare16(const int16_t * blockA, const int16_t * blockB, int size)  
 {  
         int i;  
         for (i = 0; i < size; i++)  
                 if (blockA[i] != blockB[i])  
                         return 1;  
   
         return 0;  
 }  
   
 int diff16(const int16_t * blockA, const int16_t * blockB, int size)  
 {  
         int i, diff = 0;  
         for (i = 0; i < size; i++)  
                 diff += ABS(blockA[i]-blockB[i]);  
         return diff;  
 }  
   
   
 #define XVID_TEST_RANDOM        0x00000001      /* random input data */  
 #define XVID_TEST_VERBOSE       0x00000002      /* verbose error output */  
   
   
 #define TEST_FORWARD    0x00000001      /* intra */  
 #define TEST_FDCT  (TEST_FORWARD)  
 #define TEST_IDCT  (0)  
   
 int test_transform(void * funcA, void * funcB, const char * nameB,  
                                    int test, int flags)  
 {  
         int i;  
         int64_t timeSTART;  
         int64_t timeA = 0;  
         int64_t timeB = 0;  
         DECLARE_ALIGNED_MATRIX(arrayA, 1, 64, int16_t, CACHE_LINE);  
         DECLARE_ALIGNED_MATRIX(arrayB, 1, 64, int16_t, CACHE_LINE);  
         int min, max;  
         int count = 0;  
   
         int tmp;  
         int min_error = 0x10000*64;  
         int max_error = 0;  
   
   
         if ((test & TEST_FORWARD))      /* forward */  
         {  
                 min = -256;  
                 max = 255;  
         }else{          /* inverse */  
                 min = -2048;  
                 max = 2047;  
         }  
   
         for (i = 0; i < 64*64; i++)  
         {  
                 if ((flags & XVID_TEST_RANDOM))  
                 {  
                         random16(arrayA, 64, min, max);  
                 }else{  
                         fill16(arrayA, 64, i);  
                 }  
                 memcpy(arrayB, arrayA, 64*sizeof(int16_t));  
   
                 if ((test & TEST_FORWARD))  
                 {  
                         timeSTART = read_counter();  
                         ((fdctFunc*)funcA)(arrayA);  
                         timeA += read_counter() - timeSTART;  
   
                         timeSTART = read_counter();  
                         ((fdctFunc*)funcB)(arrayB);  
                         timeB += read_counter() - timeSTART;  
                 }  
                 else  
                 {  
                         timeSTART = read_counter();  
                         ((idctFunc*)funcA)(arrayA);  
                         timeA += read_counter() - timeSTART;  
   
                         timeSTART = read_counter();  
                         ((idctFunc*)funcB)(arrayB);  
                         timeB += read_counter() - timeSTART;  
                 }  
   
                 tmp = diff16(arrayA, arrayB, 64) / 64;  
                 if (tmp > max_error)  
                         max_error = tmp;  
                 if (tmp < min_error)  
                         min_error = tmp;  
   
                 count++;  
         }  
   
         /* print the "average difference" of best/worst transforms */  
         printf("%s:\t%i\t(min_error:%i, max_error:%i)\n", nameB, (int)(timeB / count), min_error, max_error);  
   
         return 0;  
 }  
   
   
 #define TEST_QUANT      0x00000001      /* forward quantization */  
 #define TEST_INTRA      0x00000002      /* intra */  
 #define TEST_QUANT_INTRA        (TEST_QUANT|TEST_INTRA)  
 #define TEST_QUANT_INTER        (TEST_QUANT)  
 #define TEST_DEQUANT_INTRA      (TEST_INTRA)  
 #define TEST_DEQUANT_INTER      (0)  
   
 int test_quant(void * funcA, void * funcB, const char * nameB,  
                            int test, int flags)  
 {  
         int q,i;  
         int64_t timeSTART;  
         int64_t timeA = 0;  
         int64_t timeB = 0;  
         int retA, retB;  
         DECLARE_ALIGNED_MATRIX(arrayX, 1, 64, int16_t, CACHE_LINE);  
         DECLARE_ALIGNED_MATRIX(arrayA, 1, 64, int16_t, CACHE_LINE);  
         DECLARE_ALIGNED_MATRIX(arrayB, 1, 64, int16_t, CACHE_LINE);  
         int min, max;  
         int count = 0;  
         int errors = 0;  
   
         if ((test & TEST_QUANT))        /* quant */  
         {  
                 min = -2048;  
                 max = 2047;  
         }else{          /* dequant */  
                 min = -256;  
                 max = 255;  
         }  
   
         for (q = 1; q <= 31; q++)       /* quantizer */  
         {  
                 for (i = min; i < max; i++)     /* input coeff */  
                 {  
                         if ((flags & XVID_TEST_RANDOM))  
                         {  
                                 random16(arrayX, 64, min, max);  
                         }else{  
                                 fill16(arrayX, 64, i);  
                         }  
   
                         if ((test & TEST_INTRA))        /* intra */  
                         {  
                                 timeSTART = read_counter();  
                                 ((quanth263_intraFunc*)funcA)(arrayA, arrayX, q, q);  
                                 timeA += read_counter() - timeSTART;  
   
                                 timeSTART = read_counter();  
                                 ((quanth263_intraFunc*)funcB)(arrayB, arrayX, q, q);  
                                 timeB += read_counter() - timeSTART;  
                         }  
                         else    /* inter */  
                         {  
                                 timeSTART = read_counter();  
                                 retA = ((quanth263_interFunc*)funcA)(arrayA, arrayX, q);  
                                 timeA += read_counter() - timeSTART;  
   
                                 timeSTART = read_counter();  
                                 retB = ((quanth263_interFunc*)funcB)(arrayB, arrayX, q);  
                                 timeB += read_counter() - timeSTART;  
                         }  
   
                         /* compare return value from quant_inter, and compare (de)quantiz'd arrays */  
                         if ( ((test&TEST_QUANT) && !(test&TEST_INTRA) && retA != retB ) ||  
                                 compare16(arrayA, arrayB, 64))  
                         {  
                                 errors++;  
                                 if ((flags & XVID_TEST_VERBOSE))  
                                         printf("%s error: q=%i, i=%i\n", nameB, q, i);  
                         }  
   
                         count++;  
                 }  
         }  
   
         printf("%s:\t%i", nameB, (int)(timeB / count));  
         if (errors>0)  
                 printf("\t(%i errors out of %i)", errors, count);  
         printf("\n");  
   
         return 0;  
 }  
   
   
   
 int xvid_init_test(int flags)  
 {  
         int cpu_flags;  
   
         srand(time(0));  
   
         printf("xvid_init_test\n");  
   
 #if defined(ARCH_X86)  
         cpu_flags = detect_cpu_flags();  
         idct_int32_init();  
         emms_mmx();  
   
         printf("--- fdct ---\n");  
                 test_transform(fdct_int32, fdct_int32, "c", TEST_FDCT, flags);  
         if (cpu_flags & XVID_CPU_MMX)  
                 test_transform(fdct_int32, fdct_mmx, "mmx", TEST_FDCT, flags);  
         if (cpu_flags & XVID_CPU_SSE2)  
                 test_transform(fdct_int32, fdct_sse2, "sse2", TEST_FDCT, flags);  
   
         printf("\n--- idct ---\n");  
                 test_transform(idct_int32, idct_int32, "c", TEST_IDCT, flags);  
         if (cpu_flags & XVID_CPU_MMX)  
                 test_transform(idct_int32, idct_mmx, "mmx", TEST_IDCT, flags);  
         if (cpu_flags & XVID_CPU_MMXEXT)  
                 test_transform(idct_int32, idct_xmm, "xmm", TEST_IDCT, flags);  
         if (cpu_flags & XVID_CPU_3DNOWEXT)  
                 test_transform(idct_int32, idct_3dne, "3dne", TEST_IDCT, flags);  
         if (cpu_flags & XVID_CPU_SSE2)  
                 test_transform(idct_int32, idct_sse2, "sse2", TEST_IDCT, flags);  
   
         printf("\n--- quant intra ---\n");  
                 test_quant(quant_intra_c, quant_intra_c, "c", TEST_QUANT_INTRA, flags);  
         if (cpu_flags & XVID_CPU_MMX)  
                 test_quant(quant_intra_c, quant_intra_mmx, "mmx", TEST_QUANT_INTRA, flags);  
         if (cpu_flags & XVID_CPU_3DNOWEXT)  
                 test_quant(quant_intra_c, quant_intra_3dne, "3dne", TEST_QUANT_INTRA, flags);  
         if (cpu_flags & XVID_CPU_SSE2)  
                 test_quant(quant_intra_c, quant_intra_sse2, "sse2", TEST_QUANT_INTRA, flags);  
   
         printf("\n--- quant inter ---\n");  
                 test_quant(quant_inter_c, quant_inter_c, "c", TEST_QUANT_INTER, flags);  
         if (cpu_flags & XVID_CPU_MMX)  
                 test_quant(quant_inter_c, quant_inter_mmx, "mmx", TEST_QUANT_INTER, flags);  
         if (cpu_flags & XVID_CPU_3DNOWEXT)  
                 test_quant(quant_inter_c, quant_inter_3dne, "3dne", TEST_QUANT_INTER, flags);  
         if (cpu_flags & XVID_CPU_SSE2)  
                 test_quant(quant_inter_c, quant_inter_sse2, "sse2", TEST_QUANT_INTER, flags);  
   
         printf("\n--- dequant intra ---\n");  
                 test_quant(dequant_intra_c, dequant_intra_c, "c", TEST_DEQUANT_INTRA, flags);  
         if (cpu_flags & XVID_CPU_MMX)  
                 test_quant(dequant_intra_c, dequant_intra_mmx, "mmx", TEST_DEQUANT_INTRA, flags);  
         if (cpu_flags & XVID_CPU_MMXEXT)  
                 test_quant(dequant_intra_c, dequant_intra_xmm, "xmm", TEST_DEQUANT_INTRA, flags);  
         if (cpu_flags & XVID_CPU_3DNOWEXT)  
                 test_quant(dequant_intra_c, dequant_intra_3dne, "3dne", TEST_DEQUANT_INTRA, flags);  
         if (cpu_flags & XVID_CPU_SSE2)  
                 test_quant(dequant_intra_c, dequant_intra_sse2, "sse2", TEST_DEQUANT_INTRA, flags);  
   
         printf("\n--- dequant inter ---\n");  
                 test_quant(dequant_inter_c, dequant_inter_c, "c", TEST_DEQUANT_INTER, flags);  
         if (cpu_flags & XVID_CPU_MMX)  
                 test_quant(dequant_inter_c, dequant_inter_mmx, "mmx", TEST_DEQUANT_INTER, flags);  
         if (cpu_flags & XVID_CPU_MMXEXT)  
                 test_quant(dequant_inter_c, dequant_inter_xmm, "xmm", TEST_DEQUANT_INTER, flags);  
         if (cpu_flags & XVID_CPU_3DNOWEXT)  
                 test_quant(dequant_inter_c, dequant_inter_3dne, "3dne", TEST_DEQUANT_INTER, flags);  
         if (cpu_flags & XVID_CPU_SSE2)  
                 test_quant(dequant_inter_c, dequant_inter_sse2, "sse2", TEST_DEQUANT_INTER, flags);  
   
         printf("\n--- quant4_intra ---\n");  
                 test_quant(quant4_intra_c, quant4_intra_c, "c", TEST_QUANT_INTRA, flags);  
         if (cpu_flags & XVID_CPU_MMX)  
                 test_quant(quant4_intra_c, quant4_intra_mmx, "mmx", TEST_QUANT_INTRA, flags);  
         if (cpu_flags & XVID_CPU_MMXEXT)  
                 test_quant(quant4_intra_c, quant4_intra_xmm, "xmm", TEST_QUANT_INTRA, flags);  
   
         printf("\n--- quant4_inter ---\n");  
                 test_quant(quant4_inter_c, quant4_inter_c, "c", TEST_QUANT_INTER, flags);  
         if (cpu_flags & XVID_CPU_MMX)  
                 test_quant(quant4_inter_c, quant4_inter_mmx, "mmx", TEST_QUANT_INTER, flags);  
         if (cpu_flags & XVID_CPU_MMXEXT)  
                 test_quant(quant4_inter_c, quant4_inter_xmm, "xmm", TEST_QUANT_INTER, flags);  
   
         printf("\n--- dequant4_intra ---\n");  
                 test_quant(dequant4_intra_c, dequant4_intra_c, "c", TEST_DEQUANT_INTRA, flags);  
         if (cpu_flags & XVID_CPU_MMX)  
                 test_quant(dequant4_intra_c, dequant4_intra_mmx, "mmx", TEST_DEQUANT_INTRA, flags);  
         if (cpu_flags & XVID_CPU_3DNOWEXT)  
                 test_quant(dequant4_intra_c, dequant4_intra_3dne, "3dne", TEST_DEQUANT_INTRA, flags);  
   
         printf("\n--- dequant4_inter ---\n");  
                 test_quant(dequant4_inter_c, dequant4_inter_c, "c", TEST_DEQUANT_INTER, flags);  
         if (cpu_flags & XVID_CPU_MMX)  
                 test_quant(dequant4_inter_c, dequant4_inter_mmx, "mmx", TEST_DEQUANT_INTER, flags);  
         if (cpu_flags & XVID_CPU_3DNOWEXT)  
                 test_quant(dequant4_inter_c, dequant4_inter_3dne, "3dne", TEST_DEQUANT_INTER, flags);  
   
         emms_mmx();  
   
412  #endif  #endif
413    
414          return XVID_ERR_OK;          return XVID_ERR_OK;
415  }  }
416    
   
 int  
 xvid_init(void *handle,  
                   int opt,  
                   void *param1,  
                   void *param2)  
 {  
         switch(opt)  
         {  
                 case XVID_INIT_INIT :  
                         return xvid_init_init((XVID_INIT_PARAM*)param1);  
   
                 case XVID_INIT_CONVERT :  
                         return xvid_init_convert((XVID_INIT_CONVERTINFO*)param1);  
   
                 case XVID_INIT_TEST :  
                         return xvid_init_test((int)param1);  
   
                 default :  
                         return XVID_ERR_FAIL;  
         }  
 }  
   
417  /*****************************************************************************  /*****************************************************************************
418   * XviD Native decoder entry point   * XviD Native decoder entry point
419   *   *
# Line 970  Line 432 
432  {  {
433          switch (opt) {          switch (opt) {
434          case XVID_DEC_DECODE:          case XVID_DEC_DECODE:
435                  return decoder_decode((DECODER *) handle, (XVID_DEC_FRAME *) param1, (XVID_DEC_STATS*) param2);                  return decoder_decode((DECODER *) handle, (XVID_DEC_FRAME *) param1);
436    
437          case XVID_DEC_CREATE:          case XVID_DEC_CREATE:
438                  return decoder_create((XVID_DEC_PARAM *) param1);                  return decoder_create((XVID_DEC_PARAM *) param1);
# Line 1002  Line 464 
464  {  {
465          switch (opt) {          switch (opt) {
466          case XVID_ENC_ENCODE:          case XVID_ENC_ENCODE:
   
                 if (((Encoder *) handle)->mbParam.max_bframes >= 0)  
                 return encoder_encode_bframes((Encoder *) handle, (XVID_ENC_FRAME *) param1,  
                                                           (XVID_ENC_STATS *) param2);  
                 else  
467                  return encoder_encode((Encoder *) handle, (XVID_ENC_FRAME *) param1,                  return encoder_encode((Encoder *) handle, (XVID_ENC_FRAME *) param1,
468                                                            (XVID_ENC_STATS *) param2);                                                            (XVID_ENC_STATS *) param2);
469    
Legend:
Removed from v.1.33.2.23  
changed lines
  Added in v.1.36
No admin address has been configured
 ViewVC Help
Powered by ViewVC 1.0.4