1/*
2 * Copyright (c) 2018 Clément Bœsch <u pkh me>
3 *
4 * This file is part of FFmpeg.
5 *
6 * FFmpeg is free software; you can redistribute it and/or
7 * modify it under the terms of the GNU Lesser General Public
8 * License as published by the Free Software Foundation; either
9 * version 2.1 of the License, or (at your option) any later version.
10 *
11 * FFmpeg is distributed in the hope that it will be useful,
12 * but WITHOUT ANY WARRANTY; without even the implied warranty of
13 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
14 * Lesser General Public License for more details.
15 *
16 * You should have received a copy of the GNU Lesser General Public
17 * License along with FFmpeg; if not, write to the Free Software
18 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
19 */
20
21#include "checkasm.h"
22#include "libavfilter/vf_nlmeans_init.h"
23#include "libavutil/avassert.h"
24
25#define randomize_buffer(buf, size) do {    \
26    int i;                                  \
27    for (i = 0; i < size / 4; i++)          \
28        ((uint32_t *)buf)[i] = rnd();       \
29} while (0)
30
31void checkasm_check_nlmeans(void)
32{
33    NLMeansDSPContext dsp = {0};
34
35    const int w = 123;  // source width
36    const int h = 45;   // source height
37    const int p = 3;    // patch half size
38    const int r = 2;    // research window half size
39
40    ff_nlmeans_init(&dsp);
41
42    /* See the filter's code for the explanations on the variables */
43    if (check_func(dsp.compute_safe_ssd_integral_image, "ssd_integral_image")) {
44        int offx, offy;
45        const int e = p + r;
46        const int ii_w = w + e*2;
47        const int ii_h = h + e*2;
48        const int ii_lz_32 = FFALIGN(ii_w + 1, 4);
49        uint32_t *ii_orig_ref = av_calloc(ii_h + 1, ii_lz_32 * sizeof(*ii_orig_ref));
50        uint32_t *ii_ref = ii_orig_ref + ii_lz_32 + 1;
51        uint32_t *ii_orig_new = av_calloc(ii_h + 1, ii_lz_32 * sizeof(*ii_orig_new));
52        uint32_t *ii_new = ii_orig_new + ii_lz_32 + 1;
53        const int src_lz = FFALIGN(w, 16);
54        uint8_t *src = av_calloc(h, src_lz);
55
56        declare_func(void, uint32_t *dst, ptrdiff_t dst_linesize_32,
57                     const uint8_t *s1, ptrdiff_t linesize1,
58                     const uint8_t *s2, ptrdiff_t linesize2,
59                     int w, int h);
60
61        randomize_buffer(src, h * src_lz);
62
63        for (offy = -r; offy <= r; offy++) {
64            for (offx = -r; offx <= r; offx++) {
65                if (offx || offy) {
66                    const int s1x = e;
67                    const int s1y = e;
68                    const int s2x = e + offx;
69                    const int s2y = e + offy;
70                    const int startx_safe = FFMAX(s1x, s2x);
71                    const int starty_safe = FFMAX(s1y, s2y);
72                    const int u_endx_safe = FFMIN(s1x + w, s2x + w);
73                    const int endy_safe   = FFMIN(s1y + h, s2y + h);
74                    const int safe_pw = (u_endx_safe - startx_safe) & ~0xf;
75                    const int safe_ph = endy_safe - starty_safe;
76
77                    av_assert0(safe_pw && safe_ph);
78                    av_assert0(startx_safe - s1x >= 0); av_assert0(startx_safe - s1x < w);
79                    av_assert0(starty_safe - s1y >= 0); av_assert0(starty_safe - s1y < h);
80                    av_assert0(startx_safe - s2x >= 0); av_assert0(startx_safe - s2x < w);
81                    av_assert0(starty_safe - s2y >= 0); av_assert0(starty_safe - s2y < h);
82
83                    memset(ii_ref, 0, (ii_lz_32 * ii_h - 1) * sizeof(*ii_ref));
84                    memset(ii_new, 0, (ii_lz_32 * ii_h - 1) * sizeof(*ii_new));
85
86                    call_ref(ii_ref + starty_safe*ii_lz_32 + startx_safe, ii_lz_32,
87                             src + (starty_safe - s1y) * src_lz + (startx_safe - s1x), src_lz,
88                             src + (starty_safe - s2y) * src_lz + (startx_safe - s2x), src_lz,
89                             safe_pw, safe_ph);
90                    call_new(ii_new + starty_safe*ii_lz_32 + startx_safe, ii_lz_32,
91                             src + (starty_safe - s1y) * src_lz + (startx_safe - s1x), src_lz,
92                             src + (starty_safe - s2y) * src_lz + (startx_safe - s2x), src_lz,
93                             safe_pw, safe_ph);
94
95                    if (memcmp(ii_ref, ii_new, (ii_lz_32 * ii_h - 1) * sizeof(*ii_ref)))
96                        fail();
97
98                    memset(ii_new, 0, (ii_lz_32 * ii_h - 1) * sizeof(*ii_new));
99                    bench_new(ii_new + starty_safe*ii_lz_32 + startx_safe, ii_lz_32,
100                             src + (starty_safe - s1y) * src_lz + (startx_safe - s1x), src_lz,
101                             src + (starty_safe - s2y) * src_lz + (startx_safe - s2x), src_lz,
102                             safe_pw, safe_ph);
103                }
104            }
105        }
106
107        av_freep(&ii_orig_ref);
108        av_freep(&ii_orig_new);
109        av_freep(&src);
110    }
111
112    report("dsp");
113}
114