diff options
author | Michael Niedermayer <michael@niedermayer.cc> | 2017-07-09 04:10:51 +0300 |
---|---|---|
committer | Michael Niedermayer <michael@niedermayer.cc> | 2017-08-11 13:06:41 +0300 |
commit | 84786e928f9e54bdc3622f68db0ce35bae43fb46 (patch) | |
tree | b8f33d6bdd49bfaf625394585a49344c15a5d6a1 /libavcodec/tests | |
parent | 86cbffdc4db268bab0d798dca1e30dc46606dad0 (diff) |
avcodec/tests/dct: Add Mean square error test
based on quotes of IEEE 1180 / ISO/IEC 23002-1
Signed-off-by: Michael Niedermayer <michael@niedermayer.cc>
Diffstat (limited to 'libavcodec/tests')
-rw-r--r-- | libavcodec/tests/dct.c | 10 |
1 files changed, 7 insertions, 3 deletions
diff --git a/libavcodec/tests/dct.c b/libavcodec/tests/dct.c index cf71b96508..b44c66f427 100644 --- a/libavcodec/tests/dct.c +++ b/libavcodec/tests/dct.c @@ -182,6 +182,7 @@ static int dct_error(const struct algo *dct, int test, int is_idct, int speed, c int err_inf, v; int64_t err2, ti, ti1, it1, err_sum = 0; int64_t sysErr[64], sysErrMax = 0; + int64_t err2_matrix[64], err2_max = 0; int maxout = 0; int blockSumErrMax = 0, blockSumErr; AVLFG prng; @@ -194,7 +195,7 @@ static int dct_error(const struct algo *dct, int test, int is_idct, int speed, c err_inf = 0; err2 = 0; for (i = 0; i < 64; i++) - sysErr[i] = 0; + err2_matrix[i] = sysErr[i] = 0; for (it = 0; it < NB_ITS; it++) { init_block(block1, test, is_idct, &prng, vals); permute(block, block1, dct->perm_type); @@ -221,6 +222,7 @@ static int dct_error(const struct algo *dct, int test, int is_idct, int speed, c v = abs(err); if (v > err_inf) err_inf = v; + err2_matrix[i] += v * v; err2 += v * v; sysErr[i] += block[i] - block1[i]; blockSumErr += v; @@ -230,8 +232,10 @@ static int dct_error(const struct algo *dct, int test, int is_idct, int speed, c if (blockSumErrMax < blockSumErr) blockSumErrMax = blockSumErr; } - for (i = 0; i < 64; i++) + for (i = 0; i < 64; i++) { sysErrMax = FFMAX(sysErrMax, FFABS(sysErr[i])); + err2_max = FFMAX(err2_max , FFABS(err2_matrix[i])); + } for (i = 0; i < 64; i++) { if (i % 8 == 0) @@ -245,7 +249,7 @@ static int dct_error(const struct algo *dct, int test, int is_idct, int speed, c spec_err = is_idct && (err_inf > 1 || omse > 0.02 || fabs(ome) > 0.0015); if (test < 2) - spec_err = is_idct && ((double) sysErrMax / NB_ITS > 0.015); + spec_err = is_idct && ((double) err2_max / NB_ITS > 0.06 || (double) sysErrMax / NB_ITS > 0.015); printf("%s %s: max_err=%d omse=%0.8f ome=%0.8f syserr=%0.8f maxout=%d blockSumErr=%d\n", is_idct ? "IDCT" : "DCT", dct->name, err_inf, |