Mercurial > libavcodec.hg
changeset 1562:bf452704100f libavcodec
optionally merge postscale into quantization table for the float aan dct
author | michael |
---|---|
date | Thu, 23 Oct 2003 09:11:56 +0000 |
parents | aa71c7981523 |
children | 820e06c6ca9b |
files | dct-test.c faandct.c faandct.h mpegvideo.c |
diffstat | 4 files changed, 30 insertions(+), 11 deletions(-) [+] |
line wrap: on
line diff
--- a/dct-test.c Thu Oct 23 08:31:26 2003 +0000 +++ b/dct-test.c Thu Oct 23 09:11:56 2003 +0000 @@ -161,7 +161,11 @@ fdct_func(block); emms(); /* for ff_mmx_idct */ - if (fdct_func == fdct_ifast) { + if (fdct_func == fdct_ifast +#ifndef FAAN_POSTSCALE + || fdct_func == ff_faandct +#endif + ) { for(i=0; i<64; i++) { scale = 8*(1 << (AANSCALE_BITS + 11)) / aanscales[i]; block[i] = (block[i] * scale /*+ (1<<(AANSCALE_BITS-1))*/) >> AANSCALE_BITS;
--- a/faandct.c Thu Oct 23 08:31:26 2003 +0000 +++ b/faandct.c Thu Oct 23 09:11:56 2003 +0000 @@ -32,6 +32,11 @@ #include "faandct.h" #define FLOAT float +#ifdef FAAN_POSTSCALE +# define SCALE(x) postscale[x] +#else +# define SCALE(x) 1 +#endif //numbers generated by simple c code (not as accurate as they could be) /* @@ -130,12 +135,12 @@ tmp11= tmp1 + tmp2; tmp12= tmp1 - tmp2; - data[8*0 + i]= lrint(postscale[8*0 + i] * (tmp10 + tmp11)); - data[8*4 + i]= lrint(postscale[8*4 + i] * (tmp10 - tmp11)); + data[8*0 + i]= lrint(SCALE(8*0 + i) * (tmp10 + tmp11)); + data[8*4 + i]= lrint(SCALE(8*4 + i) * (tmp10 - tmp11)); z1= (tmp12 + tmp13)* A1; - data[8*2 + i]= lrint(postscale[8*2 + i] * (tmp13 + z1)); - data[8*6 + i]= lrint(postscale[8*6 + i] * (tmp13 - z1)); + data[8*2 + i]= lrint(SCALE(8*2 + i) * (tmp13 + z1)); + data[8*6 + i]= lrint(SCALE(8*6 + i) * (tmp13 - z1)); tmp10= tmp4 + tmp5; tmp11= tmp5 + tmp6; @@ -149,9 +154,9 @@ z11= tmp7 + z3; z13= tmp7 - z3; - data[8*5 + i]= lrint(postscale[8*5 + i] * (z13 + z2)); - data[8*3 + i]= lrint(postscale[8*3 + i] * (z13 - z2)); - data[8*1 + i]= lrint(postscale[8*1 + i] * (z11 + z4)); - data[8*7 + i]= lrint(postscale[8*7 + i] * (z11 - z4)); + data[8*5 + i]= lrint(SCALE(8*5 + i) * (z13 + z2)); + data[8*3 + i]= lrint(SCALE(8*3 + i) * (z13 - z2)); + data[8*1 + i]= lrint(SCALE(8*1 + i) * (z11 + z4)); + data[8*7 + i]= lrint(SCALE(8*7 + i) * (z11 - z4)); } }
--- a/faandct.h Thu Oct 23 08:31:26 2003 +0000 +++ b/faandct.h Thu Oct 23 09:11:56 2003 +0000 @@ -25,4 +25,6 @@ * @author Michael Niedermayer <michaelni@gmx.at> */ +#define FAAN_POSTSCALE + void ff_faandct(DCTELEM * data);
--- a/mpegvideo.c Thu Oct 23 08:31:26 2003 +0000 +++ b/mpegvideo.c Thu Oct 23 09:11:56 2003 +0000 @@ -101,7 +101,11 @@ for(qscale=qmin; qscale<=qmax; qscale++){ int i; - if (s->dsp.fdct == ff_jpeg_fdct_islow || s->dsp.fdct == ff_faandct) { + if (s->dsp.fdct == ff_jpeg_fdct_islow +#ifdef FAAN_POSTSCALE + || s->dsp.fdct == ff_faandct +#endif + ) { for(i=0;i<64;i++) { const int j= s->dsp.idct_permutation[i]; /* 16 <= qscale * quant_matrix[i] <= 7905 */ @@ -112,7 +116,11 @@ qmat[qscale][i] = (int)((uint64_t_C(1) << QMAT_SHIFT) / (qscale * quant_matrix[j])); } - } else if (s->dsp.fdct == fdct_ifast) { + } else if (s->dsp.fdct == fdct_ifast +#ifndef FAAN_POSTSCALE + || s->dsp.fdct == ff_faandct +#endif + ) { for(i=0;i<64;i++) { const int j= s->dsp.idct_permutation[i]; /* 16 <= qscale * quant_matrix[i] <= 7905 */