X-Git-Url: https://git.sesse.net/?a=blobdiff_plain;f=libavcodec%2Fjpeg2000dwt.c;h=188cc261a415a3c023f3e16745610671bbc75b97;hb=e826ceeb5594c5e63860ec5e5050f86307456898;hp=a46c93a9b20ea98984cfab89ca2770e1945cc76c;hpb=d62984d8a0816e05fb2d9dbe9c899a5563065074;p=ffmpeg diff --git a/libavcodec/jpeg2000dwt.c b/libavcodec/jpeg2000dwt.c index a46c93a9b20..188cc261a41 100644 --- a/libavcodec/jpeg2000dwt.c +++ b/libavcodec/jpeg2000dwt.c @@ -622,125 +622,3 @@ void ff_dwt_destroy(DWTContext *s) av_freep(&s->f_linebuf); av_freep(&s->i_linebuf); } - -#ifdef TEST - -#include "libavutil/lfg.h" - -#define MAX_W 256 - -static int test_dwt(int *array, int *ref, int border[2][2], int decomp_levels, int type, int max_diff) { - int ret, j; - DWTContext s1={{{0}}}, *s= &s1; - int64_t err2 = 0; - - ret = ff_jpeg2000_dwt_init(s, border, decomp_levels, type); - if (ret < 0) { - fprintf(stderr, "ff_jpeg2000_dwt_init failed\n"); - return 1; - } - ret = ff_dwt_encode(s, array); - if (ret < 0) { - fprintf(stderr, "ff_dwt_encode failed\n"); - return 1; - } - ret = ff_dwt_decode(s, array); - if (ret < 0) { - fprintf(stderr, "ff_dwt_encode failed\n"); - return 1; - } - for (j = 0; j max_diff) { - fprintf(stderr, "missmatch at %d (%d != %d) decomp:%d border %d %d %d %d\n", - j, array[j], ref[j],decomp_levels, border[0][0], border[0][1], border[1][0], border[1][1]); - return 2; - } - err2 += (array[j] - ref[j]) * (array[j] - ref[j]); - array[j] = ref[j]; - } - ff_dwt_destroy(s); - - printf("%s, decomp:%2d border %3d %3d %3d %3d milli-err2:%9"PRId64"\n", - type == FF_DWT53 ? "5/3i" : "9/7i", - decomp_levels, border[0][0], border[0][1], border[1][0], border[1][1], - 1000*err2 / ((border[0][1] - border[0][0])*(border[1][1] - border[1][0]))); - - return 0; -} - -static int test_dwtf(float *array, float *ref, int border[2][2], int decomp_levels, float max_diff) { - int ret, j; - DWTContext s1={{{0}}}, *s= &s1; - double err2 = 0; - - ret = ff_jpeg2000_dwt_init(s, border, decomp_levels, FF_DWT97); - if (ret < 0) { - fprintf(stderr, "ff_jpeg2000_dwt_init failed\n"); - return 1; - } - ret = ff_dwt_encode(s, array); - if (ret < 0) { - fprintf(stderr, "ff_dwt_encode failed\n"); - return 1; - } - ret = ff_dwt_decode(s, array); - if (ret < 0) { - fprintf(stderr, "ff_dwt_encode failed\n"); - return 1; - } - for (j = 0; j max_diff) { - fprintf(stderr, "missmatch at %d (%f != %f) decomp:%d border %d %d %d %d\n", - j, array[j], ref[j],decomp_levels, border[0][0], border[0][1], border[1][0], border[1][1]); - return 2; - } - err2 += (array[j] - ref[j]) * (array[j] - ref[j]); - array[j] = ref[j]; - } - ff_dwt_destroy(s); - - printf("9/7f, decomp:%2d border %3d %3d %3d %3d err2:%20.3f\n", - decomp_levels, border[0][0], border[0][1], border[1][0], border[1][1], - err2 / ((border[0][1] - border[0][0])*(border[1][1] - border[1][0]))); - - return 0; -} - -static int array[MAX_W * MAX_W]; -static int ref [MAX_W * MAX_W]; -static float arrayf[MAX_W * MAX_W]; -static float reff [MAX_W * MAX_W]; - -int main(void) { - AVLFG prng; - int i,j; - int border[2][2]; - int ret, decomp_levels; - - av_lfg_init(&prng, 1); - - for (i = 0; i>1][j&1] = av_lfg_get(&prng) % MAX_W; - if (border[0][0] >= border[0][1] || border[1][0] >= border[1][1]) - continue; - decomp_levels = av_lfg_get(&prng) % FF_DWT_MAX_DECLVLS; - - ret = test_dwt(array, ref, border, decomp_levels, FF_DWT53, 0); - if (ret) - return ret; - ret = test_dwt(array, ref, border, decomp_levels, FF_DWT97_INT, FFMIN(7+5*decomp_levels, 15+3*decomp_levels)); - if (ret) - return ret; - ret = test_dwtf(arrayf, reff, border, decomp_levels, 0.05); - if (ret) - return ret; - } - - return 0; -} - -#endif