[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.13, Fri Nov 8 10:10:48 2002 UTC revision 1.33.2.21, Fri Jan 3 16:25:14 2003 UTC
# Line 41  Line 41 
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"
51  #include "encoder.h"  #include "encoder.h"
# Line 49  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"  #include "utils/mbfunctions.h"
60  #include "quant/quant_h263.h"  #include "quant/quant_h263.h"
# Line 234  Line 240 
240          interpolate8x8_avg2 = interpolate8x8_avg2_c;          interpolate8x8_avg2 = interpolate8x8_avg2_c;
241          interpolate8x8_avg4 = interpolate8x8_avg4_c;          interpolate8x8_avg4 = interpolate8x8_avg4_c;
242    
243            /* reduced resoltuion */
244    
245            copy_upsampled_8x8_16to8 = xvid_Copy_Upsampled_8x8_16To8_C;
246            add_upsampled_8x8_16to8 = xvid_Add_Upsampled_8x8_16To8_C;
247    #ifdef ARCH_X86
248            vfilter_31 = xvid_VFilter_31_x86;
249            hfilter_31 = xvid_HFilter_31_x86;
250    #else
251            vfilter_31 = xvid_VFilter_31_C;
252            hfilter_31 = xvid_HFilter_31_C;
253    #endif
254            filter_18x18_to_8x8 = xvid_Filter_18x18_To_8x8_C;
255            filter_diff_18x18_to_8x8 = xvid_Filter_Diff_18x18_To_8x8_C;
256    
257          /* Initialize internal colorspace transformation tables */          /* Initialize internal colorspace transformation tables */
258          colorspace_init();          colorspace_init();
259    
# Line 334  Line 354 
354                  interpolate8x8_6tap_lowpass_h = interpolate8x8_6tap_lowpass_h_mmx;                  interpolate8x8_6tap_lowpass_h = interpolate8x8_6tap_lowpass_h_mmx;
355                  interpolate8x8_6tap_lowpass_v = interpolate8x8_6tap_lowpass_v_mmx;                  interpolate8x8_6tap_lowpass_v = interpolate8x8_6tap_lowpass_v_mmx;
356    
357  //              interpolate8x8_avg2 = interpolate8x8_avg2_mmx;                  interpolate8x8_avg2 = interpolate8x8_avg2_mmx;
358                  interpolate8x8_avg4 = interpolate8x8_avg4_mmx;                  interpolate8x8_avg4 = interpolate8x8_avg4_mmx;
359    
360                    /* reduced resolution */
361                    copy_upsampled_8x8_16to8 = xvid_Copy_Upsampled_8x8_16To8_mmx;
362                    add_upsampled_8x8_16to8 = xvid_Add_Upsampled_8x8_16To8_mmx;
363                    hfilter_31 = xvid_HFilter_31_mmx;
364                    filter_18x18_to_8x8 = xvid_Filter_18x18_To_8x8_mmx;
365                    filter_diff_18x18_to_8x8 = xvid_Filter_Diff_18x18_To_8x8_mmx;
366    
367                  /* image input xxx_to_yv12 related functions */                  /* image input xxx_to_yv12 related functions */
368                  yv12_to_yv12  = yv12_to_yv12_mmx;                  yv12_to_yv12  = yv12_to_yv12_mmx;
369                  bgr_to_yv12   = bgr_to_yv12_mmx;                  bgr_to_yv12   = bgr_to_yv12_mmx;
# Line 386  Line 413 
413                  interpolate8x8_halfpel_v  = interpolate8x8_halfpel_v_xmm;                  interpolate8x8_halfpel_v  = interpolate8x8_halfpel_v_xmm;
414                  interpolate8x8_halfpel_hv = interpolate8x8_halfpel_hv_xmm;                  interpolate8x8_halfpel_hv = interpolate8x8_halfpel_hv_xmm;
415    
416                    /* reduced resolution */
417                    copy_upsampled_8x8_16to8 = xvid_Copy_Upsampled_8x8_16To8_xmm;
418                    add_upsampled_8x8_16to8 = xvid_Add_Upsampled_8x8_16To8_xmm;
419    
420                  /* Quantization */                  /* Quantization */
421                    quant4_intra = quant4_intra_xmm;
422                    quant4_inter = quant4_inter_xmm;
423    
424                  dequant_intra = dequant_intra_xmm;                  dequant_intra = dequant_intra_xmm;
425                  dequant_inter = dequant_inter_xmm;                  dequant_inter = dequant_inter_xmm;
426    
# Line 415  Line 449 
449                  interpolate8x8_halfpel_hv = interpolate8x8_halfpel_hv_3dn;                  interpolate8x8_halfpel_hv = interpolate8x8_halfpel_hv_3dn;
450          }          }
451    
452            if ((cpu_flags & XVID_CPU_3DNOWEXT) > 0) {
453    
454                    /* Inverse DCT */
455                    idct =  idct_3dne;
456    
457                    /* Buffer transfer */
458                    transfer_8to16copy =  transfer_8to16copy_3dne;
459                    transfer_16to8copy = transfer_16to8copy_3dne;
460                    transfer_8to16sub =  transfer_8to16sub_3dne;
461                    transfer_8to16sub2 =  transfer_8to16sub2_3dne;
462                    transfer_16to8add = transfer_16to8add_3dne;
463                    transfer8x8_copy = transfer8x8_copy_3dne;
464    
465                    /* Quantization */
466                    dequant4_intra = dequant4_intra_3dne;
467                    dequant4_inter = dequant4_inter_3dne;
468                    quant_intra = quant_intra_3dne;
469                    quant_inter = quant_inter_3dne;
470                    dequant_intra = dequant_intra_3dne;
471                    dequant_inter = dequant_inter_3dne;
472    
473                    /* ME functions */
474                    calc_cbp = calc_cbp_3dne;
475                    sad16 = sad16_3dne;
476                    sad8 = sad8_3dne;
477                    sad16bi = sad16bi_3dne;
478                    sad8bi = sad8bi_3dne;
479                    dev16 = dev16_3dne;
480    
481                    /* Interpolation */
482                    interpolate8x8_halfpel_h = interpolate8x8_halfpel_h_3dne;
483                    interpolate8x8_halfpel_v = interpolate8x8_halfpel_v_3dne;
484                    interpolate8x8_halfpel_hv = interpolate8x8_halfpel_hv_3dne;
485            }
486    
487    
488          if ((cpu_flags & XVID_CPU_SSE2) > 0) {          if ((cpu_flags & XVID_CPU_SSE2) > 0) {
489  #ifdef EXPERIMENTAL_SSE2_CODE  #ifdef EXPERIMENTAL_SSE2_CODE
490    
# Line 487  Line 557 
557  static int  static int
558  xvid_init_convert(XVID_INIT_CONVERTINFO* convert)  xvid_init_convert(XVID_INIT_CONVERTINFO* convert)
559  {  {
560          const int flip1 = (convert->input.colorspace & XVID_CSP_VFLIP) ^ (convert->output.colorspace & XVID_CSP_VFLIP);          // const int flip1 = (convert->input.colorspace & XVID_CSP_VFLIP) ^ (convert->output.colorspace & XVID_CSP_VFLIP);
561          const int width = convert->width;          const int width = convert->width;
562          const int height = convert->height;          const int height = convert->height;
563          const int width2 = convert->width/2;          const int width2 = convert->width/2;
# Line 515  Line 585 
585  }  }
586    
587    
588    
589    void fill8(uint8_t * block, int size, int value)
590    {
591            int i;
592            for (i = 0; i < size; i++)
593                    block[i] = value;
594    }
595    
596    void fill16(int16_t * block, int size, int value)
597    {
598            int i;
599            for (i = 0; i < size; i++)
600                    block[i] = value;
601    }
602    
603    #define RANDOM(min,max) min + (rand() % (max-min))
604    
605    void random8(uint8_t * block, int size, int min, int max)
606    {
607            int i;
608            for (i = 0; i < size; i++)
609                    block[i] = RANDOM(min,max);
610    }
611    
612    void random16(int16_t * block, int size, int min, int max)
613    {
614            int i;
615            for (i = 0; i < size; i++)
616                    block[i] = RANDOM(min,max);
617    }
618    
619    int compare16(const int16_t * blockA, const int16_t * blockB, int size)
620    {
621            int i;
622            for (i = 0; i < size; i++)
623                    if (blockA[i] != blockB[i])
624                            return 1;
625    
626            return 0;
627    }
628    
629    int diff16(const int16_t * blockA, const int16_t * blockB, int size)
630    {
631            int i, diff = 0;
632            for (i = 0; i < size; i++)
633                    diff += ABS(blockA[i]-blockB[i]);
634            return diff;
635    }
636    
637    
638    #define XVID_TEST_RANDOM        0x00000001      /* random input data */
639    #define XVID_TEST_VERBOSE       0x00000002      /* verbose error output */
640    
641    
642    #define TEST_FORWARD    0x00000001      /* intra */
643    #define TEST_FDCT  (TEST_FORWARD)
644    #define TEST_IDCT  (0)
645    
646    int test_transform(void * funcA, void * funcB, const char * nameB,
647                                       int test, int flags)
648    {
649            int i;
650            int64_t timeSTART;
651            int64_t timeA = 0;
652            int64_t timeB = 0;
653            DECLARE_ALIGNED_MATRIX(arrayA, 1, 64, int16_t, CACHE_LINE);
654            DECLARE_ALIGNED_MATRIX(arrayB, 1, 64, int16_t, CACHE_LINE);
655            int min, max;
656            int count = 0;
657    
658            int tmp;
659            int min_error = 0x10000*64;
660            int max_error = 0;
661    
662    
663            if ((test & TEST_FORWARD))      /* forward */
664            {
665                    min = -256;
666                    max = 255;
667            }else{          /* inverse */
668                    min = -2048;
669                    max = 2047;
670            }
671    
672            for (i = 0; i < 64*64; i++)
673            {
674                    if ((flags & XVID_TEST_RANDOM))
675                    {
676                            random16(arrayA, 64, min, max);
677                    }else{
678                            fill16(arrayA, 64, i);
679                    }
680                    memcpy(arrayB, arrayA, 64*sizeof(int16_t));
681    
682                    if ((test & TEST_FORWARD))
683                    {
684                            timeSTART = read_counter();
685                            ((fdctFunc*)funcA)(arrayA);
686                            timeA += read_counter() - timeSTART;
687    
688                            timeSTART = read_counter();
689                            ((fdctFunc*)funcB)(arrayB);
690                            timeB += read_counter() - timeSTART;
691                    }
692                    else
693                    {
694                            timeSTART = read_counter();
695                            ((idctFunc*)funcA)(arrayA);
696                            timeA += read_counter() - timeSTART;
697    
698                            timeSTART = read_counter();
699                            ((idctFunc*)funcB)(arrayB);
700                            timeB += read_counter() - timeSTART;
701                    }
702    
703                    tmp = diff16(arrayA, arrayB, 64) / 64;
704                    if (tmp > max_error)
705                            max_error = tmp;
706                    if (tmp < min_error)
707                            min_error = tmp;
708    
709                    count++;
710            }
711    
712            /* print the "average difference" of best/worst transforms */
713            printf("%s:\t%i\t(min_error:%i, max_error:%i)\n", nameB, (int)(timeB / count), min_error, max_error);
714    
715            return 0;
716    }
717    
718    
719    #define TEST_QUANT      0x00000001      /* forward quantization */
720    #define TEST_INTRA      0x00000002      /* intra */
721    #define TEST_QUANT_INTRA        (TEST_QUANT|TEST_INTRA)
722    #define TEST_QUANT_INTER        (TEST_QUANT)
723    #define TEST_DEQUANT_INTRA      (TEST_INTRA)
724    #define TEST_DEQUANT_INTER      (0)
725    
726    int test_quant(void * funcA, void * funcB, const char * nameB,
727                               int test, int flags)
728    {
729            int q,i;
730            int64_t timeSTART;
731            int64_t timeA = 0;
732            int64_t timeB = 0;
733            int retA, retB;
734            DECLARE_ALIGNED_MATRIX(arrayX, 1, 64, int16_t, CACHE_LINE);
735            DECLARE_ALIGNED_MATRIX(arrayA, 1, 64, int16_t, CACHE_LINE);
736            DECLARE_ALIGNED_MATRIX(arrayB, 1, 64, int16_t, CACHE_LINE);
737            int min, max;
738            int count = 0;
739            int errors = 0;
740    
741            if ((test & TEST_QUANT))        /* quant */
742            {
743                    min = -2048;
744                    max = 2047;
745            }else{          /* dequant */
746                    min = -256;
747                    max = 255;
748            }
749    
750            for (q = 1; q <= 31; q++)       /* quantizer */
751            {
752                    for (i = min; i < max; i++)     /* input coeff */
753                    {
754                            if ((flags & XVID_TEST_RANDOM))
755                            {
756                                    random16(arrayX, 64, min, max);
757                            }else{
758                                    fill16(arrayX, 64, i);
759                            }
760    
761                            if ((test & TEST_INTRA))        /* intra */
762                            {
763                                    timeSTART = read_counter();
764                                    ((quanth263_intraFunc*)funcA)(arrayA, arrayX, q, q);
765                                    timeA += read_counter() - timeSTART;
766    
767                                    timeSTART = read_counter();
768                                    ((quanth263_intraFunc*)funcB)(arrayB, arrayX, q, q);
769                                    timeB += read_counter() - timeSTART;
770                            }
771                            else    /* inter */
772                            {
773                                    timeSTART = read_counter();
774                                    retA = ((quanth263_interFunc*)funcA)(arrayA, arrayX, q);
775                                    timeA += read_counter() - timeSTART;
776    
777                                    timeSTART = read_counter();
778                                    retB = ((quanth263_interFunc*)funcB)(arrayB, arrayX, q);
779                                    timeB += read_counter() - timeSTART;
780                            }
781    
782                            /* compare return value from quant_inter, and compare (de)quantiz'd arrays */
783                            if ( ((test&TEST_QUANT) && !(test&TEST_INTRA) && retA != retB ) ||
784                                    compare16(arrayA, arrayB, 64))
785                            {
786                                    errors++;
787                                    if ((flags & XVID_TEST_VERBOSE))
788                                            printf("%s error: q=%i, i=%i\n", nameB, q, i);
789                            }
790    
791                            count++;
792                    }
793            }
794    
795            printf("%s:\t%i", nameB, (int)(timeB / count));
796            if (errors>0)
797                    printf("\t(%i errors out of %i)", errors, count);
798            printf("\n");
799    
800            return 0;
801    }
802    
803    
804    
805    int xvid_init_test(int flags)
806    {
807            int cpu_flags;
808    
809            srand(time(0));
810    
811            printf("xvid_init_test\n");
812    
813    #if defined(ARCH_X86)
814            cpu_flags = check_cpu_features();
815            idct_int32_init();
816            emms_mmx();
817    
818            printf("--- fdct ---\n");
819                    test_transform(fdct_int32, fdct_int32, "c", TEST_FDCT, flags);
820            if (cpu_flags & XVID_CPU_MMX)
821                    test_transform(fdct_int32, fdct_mmx, "mmx", TEST_FDCT, flags);
822            if (cpu_flags & XVID_CPU_SSE2)
823                    test_transform(fdct_int32, fdct_sse2, "sse2", TEST_FDCT, flags);
824    
825            printf("\n--- idct ---\n");
826                    test_transform(idct_int32, idct_int32, "c", TEST_IDCT, flags);
827            if (cpu_flags & XVID_CPU_MMX)
828                    test_transform(idct_int32, idct_mmx, "mmx", TEST_IDCT, flags);
829            if (cpu_flags & XVID_CPU_MMXEXT)
830                    test_transform(idct_int32, idct_xmm, "xmm", TEST_IDCT, flags);
831            if (cpu_flags & XVID_CPU_3DNOWEXT)
832                    test_transform(idct_int32, idct_3dne, "3dne", TEST_IDCT, flags);
833            if (cpu_flags & XVID_CPU_SSE2)
834                    test_transform(idct_int32, idct_sse2, "sse2", TEST_IDCT, flags);
835    
836            printf("\n--- quant intra ---\n");
837                    test_quant(quant_intra_c, quant_intra_c, "c", TEST_QUANT_INTRA, flags);
838            if (cpu_flags & XVID_CPU_MMX)
839                    test_quant(quant_intra_c, quant_intra_mmx, "mmx", TEST_QUANT_INTRA, flags);
840            if (cpu_flags & XVID_CPU_3DNOWEXT)
841                    test_quant(quant_intra_c, quant_intra_3dne, "3dne", TEST_QUANT_INTRA, flags);
842            if (cpu_flags & XVID_CPU_SSE2)
843                    test_quant(quant_intra_c, quant_intra_sse2, "sse2", TEST_QUANT_INTRA, flags);
844    
845            printf("\n--- quant inter ---\n");
846                    test_quant(quant_inter_c, quant_inter_c, "c", TEST_QUANT_INTER, flags);
847            if (cpu_flags & XVID_CPU_MMX)
848                    test_quant(quant_inter_c, quant_inter_mmx, "mmx", TEST_QUANT_INTER, flags);
849            if (cpu_flags & XVID_CPU_3DNOWEXT)
850                    test_quant(quant_inter_c, quant_inter_3dne, "3dne", TEST_QUANT_INTER, flags);
851            if (cpu_flags & XVID_CPU_SSE2)
852                    test_quant(quant_inter_c, quant_inter_sse2, "sse2", TEST_QUANT_INTER, flags);
853    
854            printf("\n--- dequant intra ---\n");
855                    test_quant(dequant_intra_c, dequant_intra_c, "c", TEST_DEQUANT_INTRA, flags);
856            if (cpu_flags & XVID_CPU_MMX)
857                    test_quant(dequant_intra_c, dequant_intra_mmx, "mmx", TEST_DEQUANT_INTRA, flags);
858            if (cpu_flags & XVID_CPU_MMXEXT)
859                    test_quant(dequant_intra_c, dequant_intra_xmm, "xmm", TEST_DEQUANT_INTRA, flags);
860            if (cpu_flags & XVID_CPU_3DNOWEXT)
861                    test_quant(dequant_intra_c, dequant_intra_3dne, "3dne", TEST_DEQUANT_INTRA, flags);
862            if (cpu_flags & XVID_CPU_SSE2)
863                    test_quant(dequant_intra_c, dequant_intra_sse2, "sse2", TEST_DEQUANT_INTRA, flags);
864    
865            printf("\n--- dequant inter ---\n");
866                    test_quant(dequant_inter_c, dequant_inter_c, "c", TEST_DEQUANT_INTER, flags);
867            if (cpu_flags & XVID_CPU_MMX)
868                    test_quant(dequant_inter_c, dequant_inter_mmx, "mmx", TEST_DEQUANT_INTER, flags);
869            if (cpu_flags & XVID_CPU_MMXEXT)
870                    test_quant(dequant_inter_c, dequant_inter_xmm, "xmm", TEST_DEQUANT_INTER, flags);
871            if (cpu_flags & XVID_CPU_3DNOWEXT)
872                    test_quant(dequant_inter_c, dequant_inter_3dne, "3dne", TEST_DEQUANT_INTER, flags);
873            if (cpu_flags & XVID_CPU_SSE2)
874                    test_quant(dequant_inter_c, dequant_inter_sse2, "sse2", TEST_DEQUANT_INTER, flags);
875    
876            printf("\n--- quant4_intra ---\n");
877                    test_quant(quant4_intra_c, quant4_intra_c, "c", TEST_QUANT_INTRA, flags);
878            if (cpu_flags & XVID_CPU_MMX)
879                    test_quant(quant4_intra_c, quant4_intra_mmx, "mmx", TEST_QUANT_INTRA, flags);
880            if (cpu_flags & XVID_CPU_MMXEXT)
881                    test_quant(quant4_intra_c, quant4_intra_xmm, "xmm", TEST_QUANT_INTRA, flags);
882    
883            printf("\n--- quant4_inter ---\n");
884                    test_quant(quant4_inter_c, quant4_inter_c, "c", TEST_QUANT_INTER, flags);
885            if (cpu_flags & XVID_CPU_MMX)
886                    test_quant(quant4_inter_c, quant4_inter_mmx, "mmx", TEST_QUANT_INTER, flags);
887            if (cpu_flags & XVID_CPU_MMXEXT)
888                    test_quant(quant4_inter_c, quant4_inter_xmm, "xmm", TEST_QUANT_INTER, flags);
889    
890            printf("\n--- dequant4_intra ---\n");
891                    test_quant(dequant4_intra_c, dequant4_intra_c, "c", TEST_DEQUANT_INTRA, flags);
892            if (cpu_flags & XVID_CPU_MMX)
893                    test_quant(dequant4_intra_c, dequant4_intra_mmx, "mmx", TEST_DEQUANT_INTRA, flags);
894            if (cpu_flags & XVID_CPU_3DNOWEXT)
895                    test_quant(dequant4_intra_c, dequant4_intra_3dne, "3dne", TEST_DEQUANT_INTRA, flags);
896    
897            printf("\n--- dequant4_inter ---\n");
898                    test_quant(dequant4_inter_c, dequant4_inter_c, "c", TEST_DEQUANT_INTER, flags);
899            if (cpu_flags & XVID_CPU_MMX)
900                    test_quant(dequant4_inter_c, dequant4_inter_mmx, "mmx", TEST_DEQUANT_INTER, flags);
901            if (cpu_flags & XVID_CPU_3DNOWEXT)
902                    test_quant(dequant4_inter_c, dequant4_inter_3dne, "3dne", TEST_DEQUANT_INTER, flags);
903    
904            emms_mmx();
905    
906    #endif
907    
908            return XVID_ERR_OK;
909    }
910    
911    
912  int  int
913  xvid_init(void *handle,  xvid_init(void *handle,
914                    int opt,                    int opt,
# Line 529  Line 923 
923                  case XVID_INIT_CONVERT :                  case XVID_INIT_CONVERT :
924                          return xvid_init_convert((XVID_INIT_CONVERTINFO*)param1);                          return xvid_init_convert((XVID_INIT_CONVERTINFO*)param1);
925    
926                    case XVID_INIT_TEST :
927                            return xvid_init_test((int)param1);
928    
929                  default :                  default :
930                          return XVID_ERR_FAIL;                          return XVID_ERR_FAIL;
931          }          }

Legend:
Removed from v.1.33.2.13  
changed lines
  Added in v.1.33.2.21

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