From 568a4323fbde03665b2b23a98068d02b39121812 Mon Sep 17 00:00:00 2001 From: Janne Grunau Date: Mon, 7 Dec 2015 23:38:46 +0100 Subject: [PATCH] checkasm: add synth_filter test --- tests/checkasm/Makefile | 2 +- tests/checkasm/checkasm.c | 1 + tests/checkasm/checkasm.h | 1 + tests/checkasm/synth_filter.c | 121 ++++++++++++++++++++++++++++++++++ 4 files changed, 124 insertions(+), 1 deletion(-) create mode 100644 tests/checkasm/synth_filter.c diff --git a/tests/checkasm/Makefile b/tests/checkasm/Makefile index a7d13d5720..9bd13ac56d 100644 --- a/tests/checkasm/Makefile +++ b/tests/checkasm/Makefile @@ -1,6 +1,6 @@ # libavcodec tests AVCODECOBJS-$(CONFIG_BSWAPDSP) += bswapdsp.o -AVCODECOBJS-$(CONFIG_DCA_DECODER) += dcadsp.o +AVCODECOBJS-$(CONFIG_DCA_DECODER) += dcadsp.o synth_filter.o AVCODECOBJS-$(CONFIG_H264PRED) += h264pred.o AVCODECOBJS-$(CONFIG_H264QPEL) += h264qpel.o AVCODECOBJS-$(CONFIG_HEVC_DECODER) += hevc_mc.o diff --git a/tests/checkasm/checkasm.c b/tests/checkasm/checkasm.c index a563eaf042..c61e4d4d96 100644 --- a/tests/checkasm/checkasm.c +++ b/tests/checkasm/checkasm.c @@ -63,6 +63,7 @@ static const struct { #endif #if CONFIG_DCA_DECODER { "dcadsp", checkasm_check_dcadsp }, + { "synth_filter", checkasm_check_synth_filter }, #endif #if CONFIG_H264PRED { "h264pred", checkasm_check_h264pred }, diff --git a/tests/checkasm/checkasm.h b/tests/checkasm/checkasm.h index eb8b6dd11c..a599dba82e 100644 --- a/tests/checkasm/checkasm.h +++ b/tests/checkasm/checkasm.h @@ -35,6 +35,7 @@ void checkasm_check_dcadsp(void); void checkasm_check_h264pred(void); void checkasm_check_h264qpel(void); void checkasm_check_hevc_mc(void); +void checkasm_check_synth_filter(void); void checkasm_check_v210enc(void); void *checkasm_check_func(void *func, const char *name, ...) av_printf_format(2, 3); diff --git a/tests/checkasm/synth_filter.c b/tests/checkasm/synth_filter.c new file mode 100644 index 0000000000..157400b570 --- /dev/null +++ b/tests/checkasm/synth_filter.c @@ -0,0 +1,121 @@ +/* + * Copyright (c) 2015 Janne Grunau + * + * This file is part of Libav. + * + * Libav is free software; you can redistribute it and/or modify + * it under the terms of the GNU General Public License as published by + * the Free Software Foundation; either version 2 of the License, or + * (at your option) any later version. + * + * Libav is distributed in the hope that it will be useful, + * but WITHOUT ANY WARRANTY; without even the implied warranty of + * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the + * GNU General Public License for more details. + * + * You should have received a copy of the GNU General Public License along + * with Libav; if not, write to the Free Software Foundation, Inc., + * 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA. + */ + +#include +#include +#include +#include + +#include "libavutil/internal.h" +#include "libavutil/intfloat.h" +#include "libavcodec/dcadata.h" +#include "libavcodec/synth_filter.h" + +#include "checkasm.h" + +#define BUF_SIZE 32 + +#define randomize_input() \ + do { \ + int i; \ + for (i = 0; i < BUF_SIZE; i++) { \ + float f = (float)rnd() / (UINT_MAX >> 5) - 16.0f; \ + in[i] = f; \ + } \ + } while (0) + +void checkasm_check_synth_filter(void) +{ + FFTContext imdct; + SynthFilterContext synth; + + ff_mdct_init(&imdct, 6, 1, 1.0); + ff_synth_filter_init(&synth); + + if (check_func(synth.synth_filter_float, "synth_filter_float")) { + LOCAL_ALIGNED(32, float, out0, [BUF_SIZE]); + LOCAL_ALIGNED(32, float, out1, [BUF_SIZE]); + LOCAL_ALIGNED(32, float, out_b, [BUF_SIZE]); + LOCAL_ALIGNED(32, float, in, [BUF_SIZE]); + LOCAL_ALIGNED(32, float, buf2_0, [BUF_SIZE]); + LOCAL_ALIGNED(32, float, buf2_1, [BUF_SIZE]); + LOCAL_ALIGNED(32, float, buf2_b, [BUF_SIZE]); + LOCAL_ALIGNED(32, float, buf0, [512]); + LOCAL_ALIGNED(32, float, buf1, [512]); + LOCAL_ALIGNED(32, float, buf_b, [512]); + float scale = 1.0f; + int i, offset0 = 0, offset1 = 0, offset_b = 0; + + declare_func(void, FFTContext *, float *, int *, float[32], const float[512], + float[32], float[32], float); + + memset(buf2_0, 0, sizeof(*buf2_0) * BUF_SIZE); + memset(buf2_1, 0, sizeof(*buf2_1) * BUF_SIZE); + memset(buf2_b, 0, sizeof(*buf2_b) * BUF_SIZE); + memset(buf0, 0, sizeof(*buf2_0) * 512); + memset(buf1, 0, sizeof(*buf2_1) * 512); + memset(buf_b, 0, sizeof(*buf2_b) * 512); + + /* more than 1 synth_buf_offset wrap-around */ + for (i = 0; i < 20; i++) { + int j; + const float * window = (i & 1) ? ff_dca_fir_32bands_perfect : ff_dca_fir_32bands_nonperfect; + + memset(out0, 0, sizeof(*out0) * BUF_SIZE); + memset(out1, 0, sizeof(*out1) * BUF_SIZE); + memset(out_b, 0, sizeof(*out_b) * BUF_SIZE); + + randomize_input(); + + call_ref(&imdct, buf0, &offset0, buf2_0, window, + out0, in, scale); + call_new(&imdct, buf1, &offset1, buf2_1, window, + out1, in, scale); + + if (offset0 != offset1) { + fail(); + fprintf(stderr, "offsets do not match: %d, %d", offset0, offset1); + break; + } + + for (j = 0; j < BUF_SIZE; j++) { + if (!float_near_abs_eps_ulp(out0[j], out1[j], 7.0e-7, 16) || + !float_near_abs_eps_ulp(buf2_0[j], buf2_1[j], 7.0e-7, 16)) { + union av_intfloat32 o0, o1, b0, b1; + + fail(); + o0.f = out0[j]; o1.f = out1[j]; + b0.f = buf2_0[j], b1.f = buf2_1[j]; + fprintf(stderr, "out: %11g (0x%08x); %11g (0x%08x); abs diff %11g\n", + o0.f, o0.i, o1.f, o1.i, fabsf(o0.f - o1.f)); + fprintf(stderr, "buf2: %11g (0x%08x); %11g (0x%08x); abs diff %11g\n", + b0.f, b0.i, b1.f, b1.i, fabsf(b0.f - b1.f)); + break; + } + } + + bench_new(&imdct, buf_b, &offset_b, buf2_b, window, + out_b, in, scale); + } + } + ff_mdct_end(&imdct); + + report("synth_filter"); +}