# HG changeset patch # User bellard # Date 1033674563 0 # Node ID 409bdaa0b9647daa951a6bfea2c4f1ab1747e8f5 # Parent 2b7ff6dfee35a34f3467d8b7a892072b5e730a54 added IDCT248 testing diff -r 2b7ff6dfee35 -r 409bdaa0b964 dct-test.c --- a/dct-test.c Thu Oct 03 13:41:33 2002 +0000 +++ b/dct-test.c Thu Oct 03 19:49:23 2002 +0000 @@ -24,6 +24,12 @@ extern void odivx_idct_c (short *block); +void (*add_pixels_clamped)(const DCTELEM *block/*align 16*/, + UINT8 *pixels/*align 8*/, int line_size); + +void (*put_pixels_clamped)(const DCTELEM *block/*align 16*/, + UINT8 *pixels/*align 8*/, int line_size); + #define AANSCALE_BITS 12 static const unsigned short aanscales[64] = { /* precomputed values scaled up by 14 bits */ @@ -86,7 +92,6 @@ INT64 err2, ti, ti1, it1; INT64 sysErr[64], sysErrMax=0; int maxout=0; - int max_sum=0; int blockSumErrMax=0, blockSumErr; srandom(0); @@ -135,7 +140,7 @@ fdct_func == j_rev_dct || fdct_func == ff_mmxext_idct) { for(i=0;i<64;i++) block[idct_mmx_perm[i]] = block1[i]; - } else if(fdct_func == simple_idct_mmx ) { + } else if(fdct_func == ff_simple_idct_mmx ) { for(i=0;i<64;i++) block[idct_simple_mmx_perm[i]] = block1[i]; @@ -230,7 +235,7 @@ fdct_func == j_rev_dct || fdct_func == ff_mmxext_idct) { for(i=0;i<64;i++) block[idct_mmx_perm[i]] = block1[i]; - } else if(fdct_func == simple_idct_mmx ) { + } else if(fdct_func == ff_simple_idct_mmx ) { for(i=0;i<64;i++) block[idct_simple_mmx_perm[i]] = block1[i]; } else { @@ -259,19 +264,175 @@ #endif } +static UINT8 img_dest[64] __attribute__ ((aligned (8))); +static UINT8 img_dest1[64] __attribute__ ((aligned (8))); + +void idct248_ref(UINT8 *dest, int linesize, INT16 *block) +{ + static int init; + static double c8[8][8]; + static double c4[4][4]; + double block1[64], block2[64], block3[64]; + double s, sum, v; + int i, j, k; + + if (!init) { + init = 1; + + for(i=0;i<8;i++) { + sum = 0; + for(j=0;j<8;j++) { + s = (i==0) ? sqrt(1.0/8.0) : sqrt(1.0/4.0); + c8[i][j] = s * cos(M_PI * i * (j + 0.5) / 8.0); + sum += c8[i][j] * c8[i][j]; + } + } + + for(i=0;i<4;i++) { + sum = 0; + for(j=0;j<4;j++) { + s = (i==0) ? sqrt(1.0/4.0) : sqrt(1.0/2.0); + c4[i][j] = s * cos(M_PI * i * (j + 0.5) / 4.0); + sum += c4[i][j] * c4[i][j]; + } + } + } + + /* butterfly */ + for(i=0;i<4;i++) { + for(j=0;j<8;j++) { + block1[8*(2*i)+j] = (block[8*(2*i)+j] + block[8*(2*i+1)+j]) * 0.5; + block1[8*(2*i+1)+j] = (block[8*(2*i)+j] - block[8*(2*i+1)+j]) * 0.5; + } + } + + /* idct8 on lines */ + for(i=0;i<8;i++) { + for(j=0;j<8;j++) { + sum = 0; + for(k=0;k<8;k++) + sum += c8[k][j] * block1[8*i+k]; + block2[8*i+j] = sum; + } + } + + /* idct4 */ + for(i=0;i<8;i++) { + for(j=0;j<4;j++) { + /* top */ + sum = 0; + for(k=0;k<4;k++) + sum += c4[k][j] * block2[8*(2*k)+i]; + block3[8*(2*j)+i] = sum; + + /* bottom */ + sum = 0; + for(k=0;k<4;k++) + sum += c4[k][j] * block2[8*(2*k+1)+i]; + block3[8*(2*j+1)+i] = sum; + } + } + + /* clamp and store the result */ + for(i=0;i<8;i++) { + for(j=0;j<8;j++) { + v = block3[8*i+j] + 128.0; + if (v < 0) + v = 0; + else if (v > 255) + v = 255; + dest[i * linesize + j] = (int)rint(v); + } + } +} + +void idct248_error(const char *name, + void (*idct248_put)(UINT8 *dest, int line_size, INT16 *block)) +{ + int it, i, it1, ti, ti1, err_max, v; + + srandom(0); + + /* just one test to see if code is correct (precision is less + important here) */ + err_max = 0; + for(it=0;it err_max) + err_max = v; + } + } + printf("%s %s: err_inf=%d\n", + 1 ? "IDCT248" : "DCT248", + name, err_max); + + ti = gettime(); + it1 = 0; + do { + for(it=0;it]\n" "test-number 0 -> test with random matrixes\n" " 1 -> test with random sparse matrixes\n" " 2 -> do 3. test from mpeg4 std\n" - "-i test IDCT implementations\n"); + "-i test IDCT implementations\n" + "-4 test IDCT248 implementations\n"); exit(1); } int main(int argc, char **argv) { - int test_idct = 0; + int test_idct = 0, test_248_dct = 0; int c,i; int test=1; @@ -285,13 +446,16 @@ } for(;;) { - c = getopt(argc, argv, "ih"); + c = getopt(argc, argv, "ih4"); if (c == -1) break; switch(c) { case 'i': test_idct = 1; break; + case '4': + test_248_dct = 1; + break; default : case 'h': help(); @@ -303,27 +467,31 @@ printf("ffmpeg DCT/IDCT test\n"); - if (!test_idct) { - dct_error("REF-DBL", 0, fdct, fdct, test); /* only to verify code ! */ - dct_error("IJG-AAN-INT", 0, fdct_ifast, fdct, test); - dct_error("IJG-LLM-INT", 0, ff_jpeg_fdct_islow, fdct, test); - dct_error("MMX", 0, ff_fdct_mmx, fdct, test); + if (test_248_dct) { + idct248_error("SIMPLE-C", simple_idct248_put); } else { - dct_error("REF-DBL", 1, idct, idct, test); - dct_error("INT", 1, j_rev_dct, idct, test); - dct_error("LIBMPEG2-MMX", 1, ff_mmx_idct, idct, test); - dct_error("LIBMPEG2-MMXEXT", 1, ff_mmxext_idct, idct, test); - dct_error("SIMPLE-C", 1, simple_idct, idct, test); - dct_error("SIMPLE-MMX", 1, simple_idct_mmx, idct, test); -// dct_error("ODIVX-C", 1, odivx_idct_c, idct); -//printf(" test against odivx idct\n"); -// dct_error("REF", 1, idct, odivx_idct_c); -// dct_error("INT", 1, j_rev_dct, odivx_idct_c); -// dct_error("MMX", 1, ff_mmx_idct, odivx_idct_c); -// dct_error("MMXEXT", 1, ff_mmxext_idct, odivx_idct_c); -// dct_error("SIMPLE-C", 1, simple_idct, odivx_idct_c); -// dct_error("SIMPLE-MMX", 1, simple_idct_mmx, odivx_idct_c); -// dct_error("ODIVX-C", 1, odivx_idct_c, odivx_idct_c); + if (!test_idct) { + dct_error("REF-DBL", 0, fdct, fdct, test); /* only to verify code ! */ + dct_error("IJG-AAN-INT", 0, fdct_ifast, fdct, test); + dct_error("IJG-LLM-INT", 0, ff_jpeg_fdct_islow, fdct, test); + dct_error("MMX", 0, ff_fdct_mmx, fdct, test); + } else { + dct_error("REF-DBL", 1, idct, idct, test); + dct_error("INT", 1, j_rev_dct, idct, test); + dct_error("LIBMPEG2-MMX", 1, ff_mmx_idct, idct, test); + dct_error("LIBMPEG2-MMXEXT", 1, ff_mmxext_idct, idct, test); + dct_error("SIMPLE-C", 1, simple_idct, idct, test); + dct_error("SIMPLE-MMX", 1, ff_simple_idct_mmx, idct, test); + // dct_error("ODIVX-C", 1, odivx_idct_c, idct); + //printf(" test against odivx idct\n"); + // dct_error("REF", 1, idct, odivx_idct_c); + // dct_error("INT", 1, j_rev_dct, odivx_idct_c); + // dct_error("MMX", 1, ff_mmx_idct, odivx_idct_c); + // dct_error("MMXEXT", 1, ff_mmxext_idct, odivx_idct_c); + // dct_error("SIMPLE-C", 1, simple_idct, odivx_idct_c); + // dct_error("SIMPLE-MMX", 1, ff_simple_idct_mmx, odivx_idct_c); + // dct_error("ODIVX-C", 1, odivx_idct_c, odivx_idct_c); + } } return 0; }