From f679711c1b516786a39f9e582622a200502fff74 Mon Sep 17 00:00:00 2001 From: Clément Bœsch Date: Sun, 6 May 2018 10:57:23 +0200 Subject: checkasm: add vf_nlmeans test for ssd_integral_image --- tests/checkasm/Makefile | 1 + tests/checkasm/checkasm.c | 3 ++ tests/checkasm/checkasm.h | 1 + tests/checkasm/vf_nlmeans.c | 113 ++++++++++++++++++++++++++++++++++++++++++++ 4 files changed, 118 insertions(+) create mode 100644 tests/checkasm/vf_nlmeans.c (limited to 'tests') diff --git a/tests/checkasm/Makefile b/tests/checkasm/Makefile index 0233d2f989..9484acbbd7 100644 --- a/tests/checkasm/Makefile +++ b/tests/checkasm/Makefile @@ -35,6 +35,7 @@ AVFILTEROBJS-$(CONFIG_BLEND_FILTER) += vf_blend.o AVFILTEROBJS-$(CONFIG_COLORSPACE_FILTER) += vf_colorspace.o AVFILTEROBJS-$(CONFIG_HFLIP_FILTER) += vf_hflip.o AVFILTEROBJS-$(CONFIG_THRESHOLD_FILTER) += vf_threshold.o +AVFILTEROBJS-$(CONFIG_NLMEANS_FILTER) += vf_nlmeans.o CHECKASMOBJS-$(CONFIG_AVFILTER) += $(AVFILTEROBJS-yes) diff --git a/tests/checkasm/checkasm.c b/tests/checkasm/checkasm.c index ba1d1d0253..721a0912fb 100644 --- a/tests/checkasm/checkasm.c +++ b/tests/checkasm/checkasm.c @@ -159,6 +159,9 @@ static const struct { #if CONFIG_HFLIP_FILTER { "vf_hflip", checkasm_check_vf_hflip }, #endif + #if CONFIG_NLMEANS_FILTER + { "vf_nlmeans", checkasm_check_nlmeans }, + #endif #if CONFIG_THRESHOLD_FILTER { "vf_threshold", checkasm_check_vf_threshold }, #endif diff --git a/tests/checkasm/checkasm.h b/tests/checkasm/checkasm.h index dcab74de06..c45cfb46f8 100644 --- a/tests/checkasm/checkasm.h +++ b/tests/checkasm/checkasm.h @@ -62,6 +62,7 @@ void checkasm_check_huffyuvdsp(void); void checkasm_check_jpeg2000dsp(void); void checkasm_check_llviddsp(void); void checkasm_check_llviddspenc(void); +void checkasm_check_nlmeans(void); void checkasm_check_pixblockdsp(void); void checkasm_check_sbrdsp(void); void checkasm_check_synth_filter(void); diff --git a/tests/checkasm/vf_nlmeans.c b/tests/checkasm/vf_nlmeans.c new file mode 100644 index 0000000000..5e2c934226 --- /dev/null +++ b/tests/checkasm/vf_nlmeans.c @@ -0,0 +1,113 @@ +/* + * Copyright (c) 2018 Clément Bœsch + * + * This file is part of FFmpeg. + * + * FFmpeg is free software; you can redistribute it and/or + * modify it under the terms of the GNU Lesser General Public + * License as published by the Free Software Foundation; either + * version 2.1 of the License, or (at your option) any later version. + * + * FFmpeg 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 + * Lesser General Public License for more details. + * + * You should have received a copy of the GNU Lesser General Public + * License along with FFmpeg; if not, write to the Free Software + * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA + */ + +#include "checkasm.h" +#include "libavfilter/vf_nlmeans.h" +#include "libavutil/avassert.h" + +#define randomize_buffer(buf, size) do { \ + int i; \ + for (i = 0; i < size / 4; i++) \ + ((uint32_t *)buf)[i] = rnd(); \ +} while (0) + +void checkasm_check_nlmeans(void) +{ + NLMeansDSPContext dsp = {0}; + + const int w = 123; // source width + const int h = 45; // source height + const int p = 3; // patch half size + const int r = 2; // research window half size + + ff_nlmeans_init(&dsp); + + /* See the filter's code for the explanations on the variables */ + if (check_func(dsp.compute_safe_ssd_integral_image, "ssd_integral_image")) { + int offx, offy; + const int e = p + r; + const int ii_w = w + e*2; + const int ii_h = h + e*2; + const int ii_lz_32 = FFALIGN(ii_w + 1, 4); + uint32_t *ii_orig_ref = av_mallocz_array(ii_h + 1, ii_lz_32 * sizeof(*ii_orig_ref)); + uint32_t *ii_ref = ii_orig_ref + ii_lz_32 + 1; + uint32_t *ii_orig_new = av_mallocz_array(ii_h + 1, ii_lz_32 * sizeof(*ii_orig_new)); + uint32_t *ii_new = ii_orig_new + ii_lz_32 + 1; + const int src_lz = FFALIGN(w, 16); + uint8_t *src = av_mallocz_array(h, src_lz); + + declare_func(void, uint32_t *dst, ptrdiff_t dst_linesize_32, + const uint8_t *s1, ptrdiff_t linesize1, + const uint8_t *s2, ptrdiff_t linesize2, + int w, int h); + + randomize_buffer(src, h * src_lz); + + for (offy = -r; offy <= r; offy++) { + for (offx = -r; offx <= r; offx++) { + if (offx || offy) { + const int s1x = e; + const int s1y = e; + const int s2x = e + offx; + const int s2y = e + offy; + const int startx_safe = FFMAX(s1x, s2x); + const int starty_safe = FFMAX(s1y, s2y); + const int u_endx_safe = FFMIN(s1x + w, s2x + w); + const int endy_safe = FFMIN(s1y + h, s2y + h); + const int safe_pw = (u_endx_safe - startx_safe) & ~0xf; + const int safe_ph = endy_safe - starty_safe; + + av_assert0(safe_pw && safe_ph); + av_assert0(startx_safe - s1x >= 0); av_assert0(startx_safe - s1x < w); + av_assert0(starty_safe - s1y >= 0); av_assert0(starty_safe - s1y < h); + av_assert0(startx_safe - s2x >= 0); av_assert0(startx_safe - s2x < w); + av_assert0(starty_safe - s2y >= 0); av_assert0(starty_safe - s2y < h); + + memset(ii_ref, 0, ii_lz_32 * ii_h * sizeof(*ii_ref)); + memset(ii_new, 0, ii_lz_32 * ii_h * sizeof(*ii_new)); + + call_ref(ii_ref + starty_safe*ii_lz_32 + startx_safe, ii_lz_32, + src + (starty_safe - s1y) * src_lz + (startx_safe - s1x), src_lz, + src + (starty_safe - s2y) * src_lz + (startx_safe - s2x), src_lz, + safe_pw, safe_ph); + call_new(ii_new + starty_safe*ii_lz_32 + startx_safe, ii_lz_32, + src + (starty_safe - s1y) * src_lz + (startx_safe - s1x), src_lz, + src + (starty_safe - s2y) * src_lz + (startx_safe - s2x), src_lz, + safe_pw, safe_ph); + + if (memcmp(ii_ref, ii_new, ii_h * ii_lz_32 * 4)) + fail(); + + memset(ii_new, 0, ii_lz_32 * ii_h * sizeof(*ii_new)); + bench_new(ii_new + starty_safe*ii_lz_32 + startx_safe, ii_lz_32, + src + (starty_safe - s1y) * src_lz + (startx_safe - s1x), src_lz, + src + (starty_safe - s2y) * src_lz + (startx_safe - s2x), src_lz, + safe_pw, safe_ph); + } + } + } + + av_freep(&ii_orig_ref); + av_freep(&ii_orig_new); + av_freep(&src); + } + + report("dsp"); +} -- cgit v1.2.3