[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.4, Sun Mar 17 13:00:50 2002 UTC revision 1.33.2.22, Sat Jan 4 06:14:32 2003 UTC
# Line 1  Line 1 
1  /**************************************************************************  /*****************************************************************************
2   *   *
3   *      XVID MPEG-4 VIDEO CODEC   *      XVID MPEG-4 VIDEO CODEC
4   *      native api   *  - Native API implementation  -
5   *   *
6   *      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
7   *      Video tools as specified in ISO/IEC 14496-2 standard.  Those intending   *      Video tools as specified in ISO/IEC 14496-2 standard.  Those intending
# Line 24  Line 24 
24   *   *
25   *      You should have received a copy of the GNU General Public License   *      You should have received a copy of the GNU General Public License
26   *      along with this program; if not, write to the Free Software   *      along with this program; if not, write to the Free Software
27   *      Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.   *  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
28   *   *
29   *************************************************************************/   ****************************************************************************/
30    
31  /**************************************************************************  /*****************************************************************************
32   *   *
33   *      History:   *  History
34   *   *
35   *      17.03.2002      Added interpolate8x8_halfpel_hv_xmm   *      - 23.06.2002    added XVID_CPU_CHKONLY
36   *  22.12.2001  API change: added xvid_init() - Isibaar   *  - 17.03.2002        Added interpolate8x8_halfpel_hv_xmm
37   *      16.12.2001      inital version; (c)2001 peter ross <pross@cs.rmit.edu.au>   *  - 22.12.2001  API change: added xvid_init() - Isibaar
38     *  - 16.12.2001        inital version; (c)2001 peter ross <pross@cs.rmit.edu.au>
39   *   *
40   *************************************************************************/   *  $Id$
41     *
42     ****************************************************************************/
43    
44    #include <stdio.h>
45    #include <stdlib.h>
46    #include <string.h>
47    #include <time.h>
48    
49  #include "xvid.h"  #include "xvid.h"
50  #include "decoder.h"  #include "decoder.h"
# Line 47  Line 54 
54  #include "dct/fdct.h"  #include "dct/fdct.h"
55  #include "image/colorspace.h"  #include "image/colorspace.h"
56  #include "image/interpolate8x8.h"  #include "image/interpolate8x8.h"
57    #include "image/reduced.h"
58  #include "utils/mem_transfer.h"  #include "utils/mem_transfer.h"
59    #include "utils/mbfunctions.h"
60  #include "quant/quant_h263.h"  #include "quant/quant_h263.h"
61  #include "quant/quant_mpeg4.h"  #include "quant/quant_mpeg4.h"
62    #include "motion/motion.h"
63  #include "motion/sad.h"  #include "motion/sad.h"
64  #include "utils/emms.h"  #include "utils/emms.h"
65  #include "utils/timer.h"  #include "utils/timer.h"
66    #include "bitstream/mbcoding.h"
67    
68    #if defined(ARCH_X86)
69    
70    #if defined(_MSC_VER)
71    #       include <windows.h>
72    #else
73    #       include <signal.h>
74    #       include <setjmp.h>
75    
76  int xvid_init(void *handle, int opt, void *param1, void *param2)          static jmp_buf mark;
77    
78            static void
79            sigill_handler(int signal)
80  {  {
81          int cpu_flags;             longjmp(mark, 1);
82          XVID_INIT_PARAM *init_param;          }
83    #endif
84    
         init_param = (XVID_INIT_PARAM *) param1;  
85    
86          // force specific cpu settings?  /*
87          if((init_param->cpu_flags & XVID_CPU_FORCE) > 0)  calls the funcptr, and returns whether SIGILL (illegal instruction) was signalled
88                  cpu_flags = init_param->cpu_flags;  return values:
89          else {  -1 : could not determine
90    0  : SIGILL was *not* signalled
91    1  : SIGILL was signalled
92    */
93    
94    int
95    sigill_check(void (*func)())
96    {
97    #if defined(_MSC_VER)
98            _try {
99                    func();
100            }
101            _except(EXCEPTION_EXECUTE_HANDLER) {
102    
103  #ifdef ARCH_X86                  if (_exception_code() == STATUS_ILLEGAL_INSTRUCTION)
104                  cpu_flags = check_cpu_features();                          return 1;
105            }
106            return 0;
107  #else  #else
108                  cpu_flags = 0;      void * old_handler;
109        int jmpret;
110    
111    
112        old_handler = signal(SIGILL, sigill_handler);
113        if (old_handler == SIG_ERR)
114        {
115            return -1;
116        }
117    
118        jmpret = setjmp(mark);
119        if (jmpret == 0)
120        {
121            func();
122        }
123    
124        signal(SIGILL, old_handler);
125    
126        return jmpret;
127    #endif
128    }
129    #endif
130    
131    
132    /* detect cpu flags  */
133    static unsigned int
134    detect_cpu_flags()
135    {
136            /* enable native assembly optimizations by default */
137            unsigned int cpu_flags = XVID_CPU_ASM;
138    
139    #if defined(ARCH_X86)
140            cpu_flags |= check_cpu_features();
141            if ((cpu_flags & XVID_CPU_SSE) && sigill_check(sse_os_trigger))
142                    cpu_flags &= ~XVID_CPU_SSE;
143    
144            if ((cpu_flags & XVID_CPU_SSE2) && sigill_check(sse2_os_trigger))
145                    cpu_flags &= ~XVID_CPU_SSE2;
146  #endif  #endif
147    
148    #if defined(ARCH_PPC)
149    #if defined(ARCH_PPC_ALTIVEC)
150            cpu_flags |= XVID_CPU_ALTIVEC;
151    #endif
152    #endif
153    
154            return cpu_flags;
155    }
156    
157    
158    /*****************************************************************************
159     * XviD Init Entry point
160     *
161     * Well this function initialize all internal function pointers according
162     * to the CPU features forced by the library client or autodetected (depending
163     * on the XVID_CPU_FORCE flag). It also initializes vlc coding tables and all
164     * image colorspace transformation tables.
165     *
166     * Returned value : XVID_ERR_OK
167     *                  + API_VERSION in the input XVID_INIT_PARAM structure
168     *                  + core build  "   "    "       "               "
169     *
170     ****************************************************************************/
171    
172    
173    static
174    int xvid_init_init(XVID_INIT_PARAM * init_param)
175    {
176            int cpu_flags;
177    
178            /* Inform the client the API version */
179            init_param->api_version = API_VERSION;
180    
181            /* Inform the client the core build - unused because we're still alpha */
182            init_param->core_build = 1000;
183    
184            /* Do we have to force CPU features  ? */
185            if ((init_param->cpu_flags & XVID_CPU_FORCE)) {
186    
187                    cpu_flags = init_param->cpu_flags;
188    
189            } else {
190    
191                    cpu_flags = detect_cpu_flags();
192            }
193    
194            if ((init_param->cpu_flags & XVID_CPU_CHKONLY))
195            {
196                  init_param->cpu_flags = cpu_flags;                  init_param->cpu_flags = cpu_flags;
197                    return XVID_ERR_OK;
198          }          }
199    
200          // initialize the function pointers          init_param->cpu_flags = cpu_flags;
201    
202    
203            /* Initialize the function pointers */
204          idct_int32_init();          idct_int32_init();
205            init_vlc_tables();
206    
207            /* Fixed Point Forward/Inverse DCT transformations */
208          fdct = fdct_int32;          fdct = fdct_int32;
209          idct = idct_int32;          idct = idct_int32;
210    
211            /* Only needed on PPC Altivec archs */
212            sadInit = 0;
213    
214            /* Restore FPU context : emms_c is a nop functions */
215          emms = emms_c;          emms = emms_c;
216    
217            /* Quantization functions */
218          quant_intra = quant_intra_c;          quant_intra = quant_intra_c;
219          dequant_intra = dequant_intra_c;          dequant_intra = dequant_intra_c;
220          quant_inter = quant_inter_c;          quant_inter = quant_inter_c;
# Line 92  Line 225 
225          quant4_inter = quant4_inter_c;          quant4_inter = quant4_inter_c;
226          dequant4_inter = dequant4_inter_c;          dequant4_inter = dequant4_inter_c;
227    
228            /* Block transfer related functions */
229          transfer_8to16copy = transfer_8to16copy_c;          transfer_8to16copy = transfer_8to16copy_c;
230          transfer_16to8copy = transfer_16to8copy_c;          transfer_16to8copy = transfer_16to8copy_c;
231          transfer_8to16sub = transfer_8to16sub_c;          transfer_8to16sub = transfer_8to16sub_c;
232            transfer_8to16sub2 = transfer_8to16sub2_c;
233          transfer_16to8add = transfer_16to8add_c;          transfer_16to8add = transfer_16to8add_c;
234          transfer8x8_copy = transfer8x8_copy_c;          transfer8x8_copy = transfer8x8_copy_c;
235    
236            /* Interlacing functions */
237            MBFieldTest = MBFieldTest_c;
238    
239            /* Image interpolation related functions */
240          interpolate8x8_halfpel_h = interpolate8x8_halfpel_h_c;          interpolate8x8_halfpel_h = interpolate8x8_halfpel_h_c;
241          interpolate8x8_halfpel_v = interpolate8x8_halfpel_v_c;          interpolate8x8_halfpel_v = interpolate8x8_halfpel_v_c;
242          interpolate8x8_halfpel_hv = interpolate8x8_halfpel_hv_c;          interpolate8x8_halfpel_hv = interpolate8x8_halfpel_hv_c;
243    
244            interpolate16x16_lowpass_h = interpolate16x16_lowpass_h_c;
245            interpolate16x16_lowpass_v = interpolate16x16_lowpass_v_c;
246            interpolate16x16_lowpass_hv = interpolate16x16_lowpass_hv_c;
247    
248            interpolate8x8_lowpass_h = interpolate8x8_lowpass_h_c;
249            interpolate8x8_lowpass_v = interpolate8x8_lowpass_v_c;
250            interpolate8x8_lowpass_hv = interpolate8x8_lowpass_hv_c;
251    
252            interpolate8x8_6tap_lowpass_h = interpolate8x8_6tap_lowpass_h_c;
253            interpolate8x8_6tap_lowpass_v = interpolate8x8_6tap_lowpass_v_c;
254    
255            interpolate8x8_avg2 = interpolate8x8_avg2_c;
256            interpolate8x8_avg4 = interpolate8x8_avg4_c;
257    
258            /* reduced resoltuion */
259            copy_upsampled_8x8_16to8 = xvid_Copy_Upsampled_8x8_16To8_C;
260            add_upsampled_8x8_16to8 = xvid_Add_Upsampled_8x8_16To8_C;
261            vfilter_31 = xvid_VFilter_31_C;
262            hfilter_31 = xvid_HFilter_31_C;
263            filter_18x18_to_8x8 = xvid_Filter_18x18_To_8x8_C;
264            filter_diff_18x18_to_8x8 = xvid_Filter_Diff_18x18_To_8x8_C;
265    
266            /* Initialize internal colorspace transformation tables */
267          colorspace_init();          colorspace_init();
268    
269            /* All colorspace transformation functions User Format->YV12 */
270            yv12_to_yv12    = yv12_to_yv12_c;
271          rgb555_to_yv12 = rgb555_to_yv12_c;          rgb555_to_yv12 = rgb555_to_yv12_c;
272          rgb565_to_yv12 = rgb565_to_yv12_c;          rgb565_to_yv12 = rgb565_to_yv12_c;
273          rgb24_to_yv12 = rgb24_to_yv12_c;          bgr_to_yv12     = bgr_to_yv12_c;
274          rgb32_to_yv12 = rgb32_to_yv12_c;          bgra_to_yv12    = bgra_to_yv12_c;
275          yuv_to_yv12 = yuv_to_yv12_c;          abgr_to_yv12    = abgr_to_yv12_c;
276            rgba_to_yv12    = rgba_to_yv12_c;
277          yuyv_to_yv12 = yuyv_to_yv12_c;          yuyv_to_yv12 = yuyv_to_yv12_c;
278          uyvy_to_yv12 = uyvy_to_yv12_c;          uyvy_to_yv12 = uyvy_to_yv12_c;
279    
280            rgb555i_to_yv12 = rgb555i_to_yv12_c;
281            rgb565i_to_yv12 = rgb565i_to_yv12_c;
282            bgri_to_yv12    = bgri_to_yv12_c;
283            bgrai_to_yv12   = bgrai_to_yv12_c;
284            abgri_to_yv12   = abgri_to_yv12_c;
285            rgbai_to_yv12   = rgbai_to_yv12_c;
286            yuyvi_to_yv12   = yuyvi_to_yv12_c;
287            uyvyi_to_yv12   = uyvyi_to_yv12_c;
288    
289    
290            /* All colorspace transformation functions YV12->User format */
291          yv12_to_rgb555 = yv12_to_rgb555_c;          yv12_to_rgb555 = yv12_to_rgb555_c;
292          yv12_to_rgb565 = yv12_to_rgb565_c;          yv12_to_rgb565 = yv12_to_rgb565_c;
293          yv12_to_rgb24 = yv12_to_rgb24_c;          yv12_to_bgr     = yv12_to_bgr_c;
294          yv12_to_rgb32 = yv12_to_rgb32_c;          yv12_to_bgra    = yv12_to_bgra_c;
295          yv12_to_yuv = yv12_to_yuv_c;          yv12_to_abgr    = yv12_to_abgr_c;
296            yv12_to_rgba    = yv12_to_rgba_c;
297          yv12_to_yuyv = yv12_to_yuyv_c;          yv12_to_yuyv = yv12_to_yuyv_c;
298          yv12_to_uyvy = yv12_to_uyvy_c;          yv12_to_uyvy = yv12_to_uyvy_c;
299    
300            yv12_to_rgb555i = yv12_to_rgb555i_c;
301            yv12_to_rgb565i = yv12_to_rgb565i_c;
302            yv12_to_bgri    = yv12_to_bgri_c;
303            yv12_to_bgrai   = yv12_to_bgrai_c;
304            yv12_to_abgri   = yv12_to_abgri_c;
305            yv12_to_rgbai   = yv12_to_rgbai_c;
306            yv12_to_yuyvi   = yv12_to_yuyvi_c;
307            yv12_to_uyvyi   = yv12_to_uyvyi_c;
308    
309            /* Functions used in motion estimation algorithms */
310          calc_cbp = calc_cbp_c;          calc_cbp = calc_cbp_c;
311          sad16 = sad16_c;          sad16 = sad16_c;
312          sad8 = sad8_c;          sad8 = sad8_c;
313            sad16bi  = sad16bi_c;
314            sad8bi   = sad8bi_c;
315          dev16 = dev16_c;          dev16 = dev16_c;
316            sad16v   = sad16v_c;
317    
318  #ifdef ARCH_X86  //      Halfpel8_Refine = Halfpel8_Refine_c;
319          if((cpu_flags & XVID_CPU_MMX) > 0) {  
320                  fdct = fdct_mmx;  #if defined(ARCH_X86)
                 idct = idct_mmx;  
321    
322            if ((cpu_flags & XVID_CPU_ASM))
323            {
324                    vfilter_31 = xvid_VFilter_31_x86;
325                    hfilter_31 = xvid_HFilter_31_x86;
326            }
327    
328            if ((cpu_flags & XVID_CPU_MMX) || (cpu_flags & XVID_CPU_MMXEXT) ||
329                    (cpu_flags & XVID_CPU_3DNOW) || (cpu_flags & XVID_CPU_3DNOWEXT) ||
330                    (cpu_flags & XVID_CPU_SSE) || (cpu_flags & XVID_CPU_SSE2))
331            {
332                    /* Restore FPU context : emms_c is a nop functions */
333                  emms = emms_mmx;                  emms = emms_mmx;
334            }
335    
336            if ((cpu_flags & XVID_CPU_MMX)) {
337    
338                    /* Forward and Inverse Discrete Cosine Transformation functions */
339                    fdct = fdct_mmx;
340                    idct = idct_mmx;
341    
342                    /* Quantization related functions */
343                  quant_intra = quant_intra_mmx;                  quant_intra = quant_intra_mmx;
344                  dequant_intra = dequant_intra_mmx;                  dequant_intra = dequant_intra_mmx;
345                  quant_inter = quant_inter_mmx;                  quant_inter = quant_inter_mmx;
# Line 142  Line 350 
350                  quant4_inter = quant4_inter_mmx;                  quant4_inter = quant4_inter_mmx;
351                  dequant4_inter = dequant4_inter_mmx;                  dequant4_inter = dequant4_inter_mmx;
352    
353                    /* Block related functions */
354                  transfer_8to16copy = transfer_8to16copy_mmx;                  transfer_8to16copy = transfer_8to16copy_mmx;
355                  transfer_16to8copy = transfer_16to8copy_mmx;                  transfer_16to8copy = transfer_16to8copy_mmx;
356                  transfer_8to16sub = transfer_8to16sub_mmx;                  transfer_8to16sub = transfer_8to16sub_mmx;
357                    transfer_8to16sub2 = transfer_8to16sub2_mmx;
358                  transfer_16to8add = transfer_16to8add_mmx;                  transfer_16to8add = transfer_16to8add_mmx;
359                  transfer8x8_copy = transfer8x8_copy_mmx;                  transfer8x8_copy = transfer8x8_copy_mmx;
360    
361                    /* Interlacing Functions */
362                    MBFieldTest = MBFieldTest_mmx;
363    
364                    /* Image Interpolation related functions */
365                  interpolate8x8_halfpel_h = interpolate8x8_halfpel_h_mmx;                  interpolate8x8_halfpel_h = interpolate8x8_halfpel_h_mmx;
366                  interpolate8x8_halfpel_v = interpolate8x8_halfpel_v_mmx;                  interpolate8x8_halfpel_v = interpolate8x8_halfpel_v_mmx;
367                  interpolate8x8_halfpel_hv = interpolate8x8_halfpel_hv_mmx;                  interpolate8x8_halfpel_hv = interpolate8x8_halfpel_hv_mmx;
368    
369                  rgb24_to_yv12 = rgb24_to_yv12_mmx;                  interpolate8x8_6tap_lowpass_h = interpolate8x8_6tap_lowpass_h_mmx;
370                  rgb32_to_yv12 = rgb32_to_yv12_mmx;                  interpolate8x8_6tap_lowpass_v = interpolate8x8_6tap_lowpass_v_mmx;
371                  yuv_to_yv12 = yuv_to_yv12_mmx;  
372                    interpolate8x8_avg2 = interpolate8x8_avg2_mmx;
373                    interpolate8x8_avg4 = interpolate8x8_avg4_mmx;
374    
375                    /* reduced resolution */
376                    copy_upsampled_8x8_16to8 = xvid_Copy_Upsampled_8x8_16To8_mmx;
377                    add_upsampled_8x8_16to8 = xvid_Add_Upsampled_8x8_16To8_mmx;
378                    hfilter_31 = xvid_HFilter_31_mmx;
379                    filter_18x18_to_8x8 = xvid_Filter_18x18_To_8x8_mmx;
380                    filter_diff_18x18_to_8x8 = xvid_Filter_Diff_18x18_To_8x8_mmx;
381    
382                    /* image input xxx_to_yv12 related functions */
383                    yv12_to_yv12  = yv12_to_yv12_mmx;
384                    bgr_to_yv12   = bgr_to_yv12_mmx;
385                    bgra_to_yv12  = bgra_to_yv12_mmx;
386                  yuyv_to_yv12 = yuyv_to_yv12_mmx;                  yuyv_to_yv12 = yuyv_to_yv12_mmx;
387                  uyvy_to_yv12 = uyvy_to_yv12_mmx;                  uyvy_to_yv12 = uyvy_to_yv12_mmx;
388    
389                  yv12_to_rgb24 = yv12_to_rgb24_mmx;                  /* image output yv12_to_xxx related functions */
390                  yv12_to_rgb32 = yv12_to_rgb32_mmx;                  yv12_to_bgr   = yv12_to_bgr_mmx;
391                    yv12_to_bgra  = yv12_to_bgra_mmx;
392                  yv12_to_yuyv = yv12_to_yuyv_mmx;                  yv12_to_yuyv = yv12_to_yuyv_mmx;
393                  yv12_to_uyvy = yv12_to_uyvy_mmx;                  yv12_to_uyvy = yv12_to_uyvy_mmx;
394    
395                    yv12_to_yuyvi = yv12_to_yuyvi_mmx;
396                    yv12_to_uyvyi = yv12_to_uyvyi_mmx;
397    
398                    /* Motion estimation related functions */
399                  calc_cbp = calc_cbp_mmx;                  calc_cbp = calc_cbp_mmx;
400                  sad16 = sad16_mmx;                  sad16 = sad16_mmx;
401                  sad8 = sad8_mmx;                  sad8 = sad8_mmx;
402                    sad16bi = sad16bi_mmx;
403                    sad8bi  = sad8bi_mmx;
404                  dev16 = dev16_mmx;                  dev16 = dev16_mmx;
405                    sad16v   = sad16v_mmx;
406            }
407    
408            /* these 3dnow functions are faster than mmx, but slower than xmm. */
409            if ((cpu_flags & XVID_CPU_3DNOW)) {
410    
411                    emms = emms_3dn;
412    
413                    /* ME functions */
414                    sad16bi = sad16bi_3dn;
415                    sad8bi  = sad8bi_3dn;
416    
417                    yuyv_to_yv12  = yuyv_to_yv12_3dn;
418                    uyvy_to_yv12  = uyvy_to_yv12_3dn;
419          }          }
420    
421          if((cpu_flags & XVID_CPU_MMXEXT) > 0) {  
422            if ((cpu_flags & XVID_CPU_MMXEXT)) {
423    
424                    /* Inverse DCT */
425                  idct = idct_xmm;                  idct = idct_xmm;
426    
427                    /* Interpolation */
428                  interpolate8x8_halfpel_h = interpolate8x8_halfpel_h_xmm;                  interpolate8x8_halfpel_h = interpolate8x8_halfpel_h_xmm;
429                  interpolate8x8_halfpel_v = interpolate8x8_halfpel_v_xmm;                  interpolate8x8_halfpel_v = interpolate8x8_halfpel_v_xmm;
430                  interpolate8x8_halfpel_hv = interpolate8x8_halfpel_hv_xmm;                  interpolate8x8_halfpel_hv = interpolate8x8_halfpel_hv_xmm;
                 yuv_to_yv12 = yuv_to_yv12_xmm;  
431    
432                    /* reduced resolution */
433                    copy_upsampled_8x8_16to8 = xvid_Copy_Upsampled_8x8_16To8_xmm;
434                    add_upsampled_8x8_16to8 = xvid_Add_Upsampled_8x8_16To8_xmm;
435    
436                    /* Quantization */
437                    quant4_intra = quant4_intra_xmm;
438                    quant4_inter = quant4_inter_xmm;
439    
440                    dequant_intra = dequant_intra_xmm;
441                    dequant_inter = dequant_inter_xmm;
442    
443                    /* Buffer transfer */
444                    transfer_8to16sub2 = transfer_8to16sub2_xmm;
445    
446                    /* Colorspace transformation */
447                    yv12_to_yv12  = yv12_to_yv12_xmm;
448                    yuyv_to_yv12  = yuyv_to_yv12_xmm;
449                    uyvy_to_yv12  = uyvy_to_yv12_xmm;
450    
451                    /* ME functions */
452                  sad16 = sad16_xmm;                  sad16 = sad16_xmm;
453                  sad8 = sad8_xmm;                  sad8 = sad8_xmm;
454                    sad16bi = sad16bi_xmm;
455                    sad8bi  = sad8bi_xmm;
456                  dev16 = dev16_xmm;                  dev16 = dev16_xmm;
457                    sad16v   = sad16v_xmm;
458          }          }
459    
460          if((cpu_flags & XVID_CPU_3DNOW) > 0) {          if ((cpu_flags & XVID_CPU_3DNOW)) {
461    
462                    /* Interpolation */
463                  interpolate8x8_halfpel_h = interpolate8x8_halfpel_h_3dn;                  interpolate8x8_halfpel_h = interpolate8x8_halfpel_h_3dn;
464                  interpolate8x8_halfpel_v = interpolate8x8_halfpel_v_3dn;                  interpolate8x8_halfpel_v = interpolate8x8_halfpel_v_3dn;
465                  interpolate8x8_halfpel_hv = interpolate8x8_halfpel_hv_3dn;                  interpolate8x8_halfpel_hv = interpolate8x8_halfpel_hv_3dn;
466          }          }
467    
468            if ((cpu_flags & XVID_CPU_3DNOWEXT)) {
469    
470                    /* Inverse DCT */
471                    idct =  idct_3dne;
472    
473                    /* Buffer transfer */
474                    transfer_8to16copy =  transfer_8to16copy_3dne;
475                    transfer_16to8copy = transfer_16to8copy_3dne;
476                    transfer_8to16sub =  transfer_8to16sub_3dne;
477                    transfer_8to16sub2 =  transfer_8to16sub2_3dne;
478                    transfer_16to8add = transfer_16to8add_3dne;
479                    transfer8x8_copy = transfer8x8_copy_3dne;
480    
481                    /* Quantization */
482                    dequant4_intra = dequant4_intra_3dne;
483                    dequant4_inter = dequant4_inter_3dne;
484                    quant_intra = quant_intra_3dne;
485                    quant_inter = quant_inter_3dne;
486                    dequant_intra = dequant_intra_3dne;
487                    dequant_inter = dequant_inter_3dne;
488    
489                    /* ME functions */
490                    calc_cbp = calc_cbp_3dne;
491                    sad16 = sad16_3dne;
492                    sad8 = sad8_3dne;
493                    sad16bi = sad16bi_3dne;
494                    sad8bi = sad8bi_3dne;
495                    dev16 = dev16_3dne;
496    
497                    /* Interpolation */
498                    interpolate8x8_halfpel_h = interpolate8x8_halfpel_h_3dne;
499                    interpolate8x8_halfpel_v = interpolate8x8_halfpel_v_3dne;
500                    interpolate8x8_halfpel_hv = interpolate8x8_halfpel_hv_3dne;
501            }
502    
503    
504            if ((cpu_flags & XVID_CPU_SSE2)) {
505    
506                    calc_cbp = calc_cbp_sse2;
507    
508                    /* Quantization */
509                    quant_intra   = quant_intra_sse2;
510                    dequant_intra = dequant_intra_sse2;
511                    quant_inter   = quant_inter_sse2;
512                    dequant_inter = dequant_inter_sse2;
513    
514    #if defined(EXPERIMENTAL_SSE2_CODE)
515                    /* ME; slower than xmm */
516                    sad16    = sad16_sse2;
517                    dev16    = dev16_sse2;
518    #endif
519                    /* Forward and Inverse DCT */
520                    idct  = idct_sse2;
521                    fdct = fdct_sse2;
522            }
523  #endif  #endif
524    
525          // API version  #if defined(ARCH_IA64)
526          init_param->api_version = API_VERSION;          if ((cpu_flags & XVID_CPU_ASM)) { //use assembler routines?
527              idct_ia64_init();
528              fdct = fdct_ia64;
529              idct = idct_ia64;   //not yet working, crashes
530              interpolate8x8_halfpel_h = interpolate8x8_halfpel_h_ia64;
531              interpolate8x8_halfpel_v = interpolate8x8_halfpel_v_ia64;
532              interpolate8x8_halfpel_hv = interpolate8x8_halfpel_hv_ia64;
533              sad16 = sad16_ia64;
534              sad16bi = sad16bi_ia64;
535              sad8 = sad8_ia64;
536              dev16 = dev16_ia64;
537    //        Halfpel8_Refine = Halfpel8_Refine_ia64;
538              quant_intra = quant_intra_ia64;
539              dequant_intra = dequant_intra_ia64;
540              quant_inter = quant_inter_ia64;
541              dequant_inter = dequant_inter_ia64;
542              transfer_8to16copy = transfer_8to16copy_ia64;
543              transfer_16to8copy = transfer_16to8copy_ia64;
544              transfer_8to16sub = transfer_8to16sub_ia64;
545              transfer_8to16sub2 = transfer_8to16sub2_ia64;
546              transfer_16to8add = transfer_16to8add_ia64;
547              transfer8x8_copy = transfer8x8_copy_ia64;
548              DEBUG("Using IA-64 assembler routines.\n");
549            }
550    #endif
551    
552          // something clever has to be done for this  #if defined(ARCH_PPC)
553          init_param->core_build = 1000;          if ((cpu_flags & XVID_CPU_ASM))
554            {
555                    calc_cbp = calc_cbp_ppc;
556            }
557    
558            if ((cpu_flags & XVID_CPU_ALTIVEC))
559            {
560                    calc_cbp = calc_cbp_altivec;
561                    fdct = fdct_altivec;
562                    idct = idct_altivec;
563                    sadInit = sadInit_altivec;
564                    sad16 = sad16_altivec;
565                    sad8 = sad8_altivec;
566                    dev16 = dev16_altivec;
567            }
568    #endif
569    
570            return XVID_ERR_OK;
571    }
572    
573    
574    
575    static int
576    xvid_init_convert(XVID_INIT_CONVERTINFO* convert)
577    {
578            // const int flip1 = (convert->input.colorspace & XVID_CSP_VFLIP) ^ (convert->output.colorspace & XVID_CSP_VFLIP);
579            const int width = convert->width;
580            const int height = convert->height;
581            const int width2 = convert->width/2;
582            const int height2 = convert->height/2;
583            IMAGE img;
584    
585            switch (convert->input.colorspace & ~XVID_CSP_VFLIP)
586            {
587                    case XVID_CSP_YV12 :
588                            img.y = convert->input.y;
589                            img.v = (uint8_t*)convert->input.y + width*height;
590                            img.u = (uint8_t*)convert->input.y + width*height + width2*height2;
591                            image_output(&img, width, height, width,
592                                                    convert->output.y, convert->output.y_stride,
593                                                    convert->output.colorspace, convert->interlacing);
594                            break;
595    
596                    default :
597                            return XVID_ERR_FORMAT;
598            }
599    
600    
601            emms();
602          return XVID_ERR_OK;          return XVID_ERR_OK;
603  }  }
604    
605  int xvid_decore(void * handle, int opt, void * param1, void * param2)  
606    
607    void fill8(uint8_t * block, int size, int value)
608    {
609            int i;
610            for (i = 0; i < size; i++)
611                    block[i] = value;
612    }
613    
614    void fill16(int16_t * block, int size, int value)
615    {
616            int i;
617            for (i = 0; i < size; i++)
618                    block[i] = value;
619    }
620    
621    #define RANDOM(min,max) min + (rand() % (max-min))
622    
623    void random8(uint8_t * block, int size, int min, int max)
624    {
625            int i;
626            for (i = 0; i < size; i++)
627                    block[i] = RANDOM(min,max);
628    }
629    
630    void random16(int16_t * block, int size, int min, int max)
631    {
632            int i;
633            for (i = 0; i < size; i++)
634                    block[i] = RANDOM(min,max);
635    }
636    
637    int compare16(const int16_t * blockA, const int16_t * blockB, int size)
638    {
639            int i;
640            for (i = 0; i < size; i++)
641                    if (blockA[i] != blockB[i])
642                            return 1;
643    
644            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    
660    #define TEST_FORWARD    0x00000001      /* intra */
661    #define TEST_FDCT  (TEST_FORWARD)
662    #define TEST_IDCT  (0)
663    
664    int test_transform(void * funcA, void * funcB, const char * nameB,
665                                       int test, int flags)
666    {
667            int i;
668            int64_t timeSTART;
669            int64_t timeA = 0;
670            int64_t timeB = 0;
671            DECLARE_ALIGNED_MATRIX(arrayA, 1, 64, int16_t, CACHE_LINE);
672            DECLARE_ALIGNED_MATRIX(arrayB, 1, 64, int16_t, CACHE_LINE);
673            int min, max;
674            int count = 0;
675    
676            int tmp;
677            int min_error = 0x10000*64;
678            int max_error = 0;
679    
680    
681            if ((test & TEST_FORWARD))      /* forward */
682            {
683                    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                    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();
703                            ((fdctFunc*)funcA)(arrayA);
704                            timeA += read_counter() - timeSTART;
705    
706                            timeSTART = read_counter();
707                            ((fdctFunc*)funcB)(arrayB);
708                            timeB += read_counter() - timeSTART;
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            /* print the "average difference" of best/worst transforms */
731            printf("%s:\t%i\t(min_error:%i, max_error:%i)\n", nameB, (int)(timeB / count), min_error, max_error);
732    
733            return 0;
734    }
735    
736    
737    #define TEST_QUANT      0x00000001      /* forward quantization */
738    #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    int test_quant(void * funcA, void * funcB, const char * nameB,
745                               int test, int flags)
746    {
747            int q,i;
748            int64_t timeSTART;
749            int64_t timeA = 0;
750            int64_t timeB = 0;
751            int retA, retB;
752            DECLARE_ALIGNED_MATRIX(arrayX, 1, 64, int16_t, CACHE_LINE);
753            DECLARE_ALIGNED_MATRIX(arrayA, 1, 64, int16_t, CACHE_LINE);
754            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 */
769            {
770                    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);
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();
786                                    ((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;
794    
795                                    timeSTART = read_counter();
796                                    retB = ((quanth263_interFunc*)funcB)(arrayB, arrayX, q);
797                                    timeB += read_counter() - timeSTART;
798                            }
799    
800                            /* 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                                    errors++;
805                                    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            printf("%s:\t%i", nameB, (int)(timeB / count));
814            if (errors>0)
815                    printf("\t(%i errors out of %i)", errors, count);
816            printf("\n");
817    
818            return 0;
819    }
820    
821    
822    
823    int xvid_init_test(int flags)
824    {
825            int cpu_flags;
826    
827            srand(time(0));
828    
829            printf("xvid_init_test\n");
830    
831    #if defined(ARCH_X86)
832            cpu_flags = detect_cpu_flags();
833            idct_int32_init();
834            emms_mmx();
835    
836            printf("--- fdct ---\n");
837                    test_transform(fdct_int32, fdct_int32, "c", TEST_FDCT, flags);
838            if (cpu_flags & XVID_CPU_MMX)
839                    test_transform(fdct_int32, fdct_mmx, "mmx", TEST_FDCT, flags);
840            if (cpu_flags & XVID_CPU_SSE2)
841                    test_transform(fdct_int32, fdct_sse2, "sse2", TEST_FDCT, flags);
842    
843            printf("\n--- idct ---\n");
844                    test_transform(idct_int32, idct_int32, "c", TEST_IDCT, flags);
845            if (cpu_flags & XVID_CPU_MMX)
846                    test_transform(idct_int32, idct_mmx, "mmx", TEST_IDCT, flags);
847            if (cpu_flags & XVID_CPU_MMXEXT)
848                    test_transform(idct_int32, idct_xmm, "xmm", TEST_IDCT, flags);
849            if (cpu_flags & XVID_CPU_3DNOWEXT)
850                    test_transform(idct_int32, idct_3dne, "3dne", TEST_IDCT, flags);
851            if (cpu_flags & XVID_CPU_SSE2)
852                    test_transform(idct_int32, idct_sse2, "sse2", TEST_IDCT, flags);
853    
854            printf("\n--- quant intra ---\n");
855                    test_quant(quant_intra_c, quant_intra_c, "c", TEST_QUANT_INTRA, flags);
856            if (cpu_flags & XVID_CPU_MMX)
857                    test_quant(quant_intra_c, quant_intra_mmx, "mmx", TEST_QUANT_INTRA, flags);
858            if (cpu_flags & XVID_CPU_3DNOWEXT)
859                    test_quant(quant_intra_c, quant_intra_3dne, "3dne", TEST_QUANT_INTRA, flags);
860            if (cpu_flags & XVID_CPU_SSE2)
861                    test_quant(quant_intra_c, quant_intra_sse2, "sse2", TEST_QUANT_INTRA, flags);
862    
863            printf("\n--- quant inter ---\n");
864                    test_quant(quant_inter_c, quant_inter_c, "c", TEST_QUANT_INTER, flags);
865            if (cpu_flags & XVID_CPU_MMX)
866                    test_quant(quant_inter_c, quant_inter_mmx, "mmx", TEST_QUANT_INTER, flags);
867            if (cpu_flags & XVID_CPU_3DNOWEXT)
868                    test_quant(quant_inter_c, quant_inter_3dne, "3dne", TEST_QUANT_INTER, flags);
869            if (cpu_flags & XVID_CPU_SSE2)
870                    test_quant(quant_inter_c, quant_inter_sse2, "sse2", TEST_QUANT_INTER, flags);
871    
872            printf("\n--- dequant intra ---\n");
873                    test_quant(dequant_intra_c, dequant_intra_c, "c", TEST_DEQUANT_INTRA, flags);
874            if (cpu_flags & XVID_CPU_MMX)
875                    test_quant(dequant_intra_c, dequant_intra_mmx, "mmx", TEST_DEQUANT_INTRA, flags);
876            if (cpu_flags & XVID_CPU_MMXEXT)
877                    test_quant(dequant_intra_c, dequant_intra_xmm, "xmm", TEST_DEQUANT_INTRA, flags);
878            if (cpu_flags & XVID_CPU_3DNOWEXT)
879                    test_quant(dequant_intra_c, dequant_intra_3dne, "3dne", TEST_DEQUANT_INTRA, flags);
880            if (cpu_flags & XVID_CPU_SSE2)
881                    test_quant(dequant_intra_c, dequant_intra_sse2, "sse2", TEST_DEQUANT_INTRA, flags);
882    
883            printf("\n--- dequant inter ---\n");
884                    test_quant(dequant_inter_c, dequant_inter_c, "c", TEST_DEQUANT_INTER, flags);
885            if (cpu_flags & XVID_CPU_MMX)
886                    test_quant(dequant_inter_c, dequant_inter_mmx, "mmx", TEST_DEQUANT_INTER, flags);
887            if (cpu_flags & XVID_CPU_MMXEXT)
888                    test_quant(dequant_inter_c, dequant_inter_xmm, "xmm", TEST_DEQUANT_INTER, flags);
889            if (cpu_flags & XVID_CPU_3DNOWEXT)
890                    test_quant(dequant_inter_c, dequant_inter_3dne, "3dne", TEST_DEQUANT_INTER, flags);
891            if (cpu_flags & XVID_CPU_SSE2)
892                    test_quant(dequant_inter_c, dequant_inter_sse2, "sse2", TEST_DEQUANT_INTER, flags);
893    
894            printf("\n--- quant4_intra ---\n");
895                    test_quant(quant4_intra_c, quant4_intra_c, "c", TEST_QUANT_INTRA, flags);
896            if (cpu_flags & XVID_CPU_MMX)
897                    test_quant(quant4_intra_c, quant4_intra_mmx, "mmx", TEST_QUANT_INTRA, flags);
898            if (cpu_flags & XVID_CPU_MMXEXT)
899                    test_quant(quant4_intra_c, quant4_intra_xmm, "xmm", TEST_QUANT_INTRA, flags);
900    
901            printf("\n--- quant4_inter ---\n");
902                    test_quant(quant4_inter_c, quant4_inter_c, "c", TEST_QUANT_INTER, flags);
903            if (cpu_flags & XVID_CPU_MMX)
904                    test_quant(quant4_inter_c, quant4_inter_mmx, "mmx", TEST_QUANT_INTER, flags);
905            if (cpu_flags & XVID_CPU_MMXEXT)
906                    test_quant(quant4_inter_c, quant4_inter_xmm, "xmm", TEST_QUANT_INTER, flags);
907    
908            printf("\n--- dequant4_intra ---\n");
909                    test_quant(dequant4_intra_c, dequant4_intra_c, "c", TEST_DEQUANT_INTRA, flags);
910            if (cpu_flags & XVID_CPU_MMX)
911                    test_quant(dequant4_intra_c, dequant4_intra_mmx, "mmx", TEST_DEQUANT_INTRA, flags);
912            if (cpu_flags & XVID_CPU_3DNOWEXT)
913                    test_quant(dequant4_intra_c, dequant4_intra_3dne, "3dne", TEST_DEQUANT_INTRA, flags);
914    
915            printf("\n--- dequant4_inter ---\n");
916                    test_quant(dequant4_inter_c, dequant4_inter_c, "c", TEST_DEQUANT_INTER, flags);
917            if (cpu_flags & XVID_CPU_MMX)
918                    test_quant(dequant4_inter_c, dequant4_inter_mmx, "mmx", TEST_DEQUANT_INTER, flags);
919            if (cpu_flags & XVID_CPU_3DNOWEXT)
920                    test_quant(dequant4_inter_c, dequant4_inter_3dne, "3dne", TEST_DEQUANT_INTER, flags);
921    
922            emms_mmx();
923    
924    #endif
925    
926            return XVID_ERR_OK;
927    }
928    
929    
930    int
931    xvid_init(void *handle,
932                      int opt,
933                      void *param1,
934                      void *param2)
935  {  {
936          switch (opt)          switch (opt)
937          {          {
938                    case XVID_INIT_INIT :
939                            return xvid_init_init((XVID_INIT_PARAM*)param1);
940    
941                    case XVID_INIT_CONVERT :
942                            return xvid_init_convert((XVID_INIT_CONVERTINFO*)param1);
943    
944                    case XVID_INIT_TEST :
945                            return xvid_init_test((int)param1);
946    
947                    default :
948                            return XVID_ERR_FAIL;
949            }
950    }
951    
952    /*****************************************************************************
953     * XviD Native decoder entry point
954     *
955     * This function is just a wrapper to all the option cases.
956     *
957     * Returned values : XVID_ERR_FAIL when opt is invalid
958     *                   else returns the wrapped function result
959     *
960     ****************************************************************************/
961    
962    int
963    xvid_decore(void *handle,
964                            int opt,
965                            void *param1,
966                            void *param2)
967    {
968            switch (opt) {
969          case XVID_DEC_DECODE :          case XVID_DEC_DECODE :
970          return decoder_decode((DECODER *) handle, (XVID_DEC_FRAME *) param1);                  return decoder_decode((DECODER *) handle, (XVID_DEC_FRAME *) param1, (XVID_DEC_STATS*) param2);
971    
972          case XVID_DEC_CREATE :          case XVID_DEC_CREATE :
973          return decoder_create((XVID_DEC_PARAM *) param1);          return decoder_create((XVID_DEC_PARAM *) param1);
# Line 219  Line 981 
981  }  }
982    
983    
984  int xvid_encore(void * handle, int opt, void * param1, void * param2)  /*****************************************************************************
985  {   * XviD Native encoder entry point
986          switch (opt)   *
987     * This function is just a wrapper to all the option cases.
988     *
989     * Returned values : XVID_ERR_FAIL when opt is invalid
990     *                   else returns the wrapped function result
991     *
992     ****************************************************************************/
993    
994    int
995    xvid_encore(void *handle,
996                            int opt,
997                            void *param1,
998                            void *param2)
999          {          {
1000            switch (opt) {
1001          case XVID_ENC_ENCODE :          case XVID_ENC_ENCODE :
1002          return encoder_encode((Encoder *) handle, (XVID_ENC_FRAME *) param1, (XVID_ENC_STATS *) param2);  
1003                    if (((Encoder *) handle)->mbParam.max_bframes >= 0)
1004                    return encoder_encode_bframes((Encoder *) handle, (XVID_ENC_FRAME *) param1,
1005                                                              (XVID_ENC_STATS *) param2);
1006                    else
1007                    return encoder_encode((Encoder *) handle, (XVID_ENC_FRAME *) param1,
1008                                                              (XVID_ENC_STATS *) param2);
1009    
1010          case XVID_ENC_CREATE :          case XVID_ENC_CREATE :
1011          return encoder_create((XVID_ENC_PARAM *) param1);          return encoder_create((XVID_ENC_PARAM *) param1);

Legend:
Removed from v.1.4  
changed lines
  Added in v.1.33.2.22

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