1 /*
2 * This file is part of FFmpeg.
3 *
4 * FFmpeg is free software; you can redistribute it and/or modify
5 * it under the terms of the GNU General Public License as published by
6 * the Free Software Foundation; either version 2 of the License, or
7 * (at your option) any later version.
8 *
9 * FFmpeg is distributed in the hope that it will be useful,
10 * but WITHOUT ANY WARRANTY; without even the implied warranty of
11 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12 * GNU General Public License for more details.
13 *
14 * You should have received a copy of the GNU General Public License along
15 * with FFmpeg; if not, write to the Free Software Foundation, Inc.,
16 * 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
17 */
18
19 #include "libavcodec/opusdsp.h"
20
21 #include "checkasm.h"
22
23 #define randomize_float(buf, len) \
24 do { \
25 for (int i = 0; i < len; i++) { \
26 float f = (float)rnd() / (UINT_MAX >> 5) - 16.0f; \
27 buf[i] = f; \
28 } \
29 } while (0)
30
31 #define EPS 0.005
32 #define MAX_SIZE (960)
33
34 /* period is between 15 and 1022, inclusive */
test_postfilter(int period)35 static void test_postfilter(int period)
36 {
37 LOCAL_ALIGNED(16, float, data0, [MAX_SIZE + 1024]);
38 LOCAL_ALIGNED(16, float, data1, [MAX_SIZE + 1024]);
39
40 /* This filter can explode very easily, so use a tapset from the codec.
41 * In the codec these are usually multiplied by at least 0.09375f,
42 * so its outside the largest filter value, but the filter is still stable
43 * so use it. */
44 float gains[3] = { 0.3066406250f, 0.2170410156f, 0.1296386719f };
45
46 /* The codec will always call with an offset which is aligned once
47 * (period + 2) is subtracted, but here we have to align it outselves. */
48 int offset = FFALIGN(period + 2, 4);
49
50 declare_func(void, float *data, int period, float *gains, int len);
51
52 randomize_float(data0, MAX_SIZE + 1024);
53 memcpy(data1, data0, (MAX_SIZE + 1024)*sizeof(float));
54
55 call_ref(data0 + offset, period, gains, MAX_SIZE);
56 call_new(data1 + offset, period, gains, MAX_SIZE);
57
58 if (!float_near_abs_eps_array(data0 + offset, data1 + offset, EPS, MAX_SIZE))
59 fail();
60 bench_new(data1 + offset, period, gains, MAX_SIZE);
61 }
62
test_deemphasis(void)63 static void test_deemphasis(void)
64 {
65 LOCAL_ALIGNED(16, float, src, [FFALIGN(MAX_SIZE, 4)]);
66 LOCAL_ALIGNED(16, float, dst0, [FFALIGN(MAX_SIZE, 4)]);
67 LOCAL_ALIGNED(16, float, dst1, [FFALIGN(MAX_SIZE, 4)]);
68 float coeff0 = (float)rnd() / (UINT_MAX >> 5) - 16.0f, coeff1 = coeff0;
69
70 declare_func_float(float, float *out, float *in, float coeff, int len);
71
72 randomize_float(src, MAX_SIZE);
73
74 coeff0 = call_ref(dst0, src, coeff0, MAX_SIZE);
75 coeff1 = call_new(dst1, src, coeff1, MAX_SIZE);
76
77 if (!float_near_abs_eps(coeff0, coeff1, EPS) ||
78 !float_near_abs_eps_array(dst0, dst1, EPS, MAX_SIZE))
79 fail();
80 bench_new(dst1, src, coeff1, MAX_SIZE);
81 }
82
checkasm_check_opusdsp(void)83 void checkasm_check_opusdsp(void)
84 {
85 OpusDSP ctx;
86 ff_opus_dsp_init(&ctx);
87
88 if (check_func(ctx.postfilter, "postfilter_15"))
89 test_postfilter(15);
90 report("postfilter_15");
91
92 if (check_func(ctx.postfilter, "postfilter_512"))
93 test_postfilter(512);
94 report("postfilter_512");
95
96 if (check_func(ctx.postfilter, "postfilter_1022"))
97 test_postfilter(1022);
98 report("postfilter_1022");
99
100 if (check_func(ctx.deemphasis, "deemphasis"))
101 test_deemphasis();
102 report("deemphasis");
103 }
104