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 "libavutil/mem_internal.h"
20
21 #include "libavcodec/opusdsp.h"
22
23 #include "checkasm.h"
24
25 #define randomize_float(buf, len) \
26 do { \
27 for (int i = 0; i < len; i++) { \
28 float f = (float)rnd() / (UINT_MAX >> 5) - 16.0f; \
29 buf[i] = f; \
30 } \
31 } while (0)
32
33 #define EPS 0.005
34 #define MAX_SIZE (960)
35
36 /* period is between 15 and 1022, inclusive */
test_postfilter(int period)37 static void test_postfilter(int period)
38 {
39 LOCAL_ALIGNED(16, float, data0, [MAX_SIZE + 1024]);
40 LOCAL_ALIGNED(16, float, data1, [MAX_SIZE + 1024]);
41
42 /* This filter can explode very easily, so use a tapset from the codec.
43 * In the codec these are usually multiplied by at least 0.09375f,
44 * so its outside the largest filter value, but the filter is still stable
45 * so use it. */
46 float gains[3] = { 0.3066406250f, 0.2170410156f, 0.1296386719f };
47
48 /* The codec will always call with an offset which is aligned once
49 * (period + 2) is subtracted, but here we have to align it outselves. */
50 int offset = FFALIGN(period + 2, 4);
51
52 declare_func(void, float *data, int period, float *gains, int len);
53
54 randomize_float(data0, MAX_SIZE + 1024);
55 memcpy(data1, data0, (MAX_SIZE + 1024)*sizeof(float));
56
57 call_ref(data0 + offset, period, gains, MAX_SIZE);
58 call_new(data1 + offset, period, gains, MAX_SIZE);
59
60 if (!float_near_abs_eps_array(data0 + offset, data1 + offset, EPS, MAX_SIZE))
61 fail();
62 bench_new(data1 + offset, period, gains, MAX_SIZE);
63 }
64
test_deemphasis(void)65 static void test_deemphasis(void)
66 {
67 LOCAL_ALIGNED(16, float, src, [FFALIGN(MAX_SIZE, 4)]);
68 LOCAL_ALIGNED(16, float, dst0, [FFALIGN(MAX_SIZE, 4)]);
69 LOCAL_ALIGNED(16, float, dst1, [FFALIGN(MAX_SIZE, 4)]);
70 float coeff0 = (float)rnd() / (UINT_MAX >> 5) - 16.0f, coeff1 = coeff0;
71
72 declare_func_float(float, float *out, float *in, float coeff, int len);
73
74 randomize_float(src, MAX_SIZE);
75
76 coeff0 = call_ref(dst0, src, coeff0, MAX_SIZE);
77 coeff1 = call_new(dst1, src, coeff1, MAX_SIZE);
78
79 if (!float_near_abs_eps(coeff0, coeff1, EPS) ||
80 !float_near_abs_eps_array(dst0, dst1, EPS, MAX_SIZE))
81 fail();
82 bench_new(dst1, src, coeff1, MAX_SIZE);
83 }
84
checkasm_check_opusdsp(void)85 void checkasm_check_opusdsp(void)
86 {
87 OpusDSP ctx;
88 ff_opus_dsp_init(&ctx);
89
90 if (check_func(ctx.postfilter, "postfilter_15"))
91 test_postfilter(15);
92 report("postfilter_15");
93
94 if (check_func(ctx.postfilter, "postfilter_512"))
95 test_postfilter(512);
96 report("postfilter_512");
97
98 if (check_func(ctx.postfilter, "postfilter_1022"))
99 test_postfilter(1022);
100 report("postfilter_1022");
101
102 if (check_func(ctx.deemphasis, "deemphasis"))
103 test_deemphasis();
104 report("deemphasis");
105 }
106