/src/ffmpeg/libavcodec/opus/dsp.c
Line | Count | Source |
1 | | /* |
2 | | * This file is part of FFmpeg. |
3 | | * |
4 | | * FFmpeg is free software; you can redistribute it and/or |
5 | | * modify it under the terms of the GNU Lesser General Public |
6 | | * License as published by the Free Software Foundation; either |
7 | | * version 2.1 of the License, or (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 GNU |
12 | | * Lesser General Public License for more details. |
13 | | * |
14 | | * You should have received a copy of the GNU Lesser General Public |
15 | | * License along with FFmpeg; if not, write to the Free Software |
16 | | * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA |
17 | | */ |
18 | | |
19 | | #include "config.h" |
20 | | #include "libavutil/attributes.h" |
21 | | #include "libavutil/mem_internal.h" |
22 | | #include "dsp.h" |
23 | | |
24 | | static void postfilter_c(float *data, int period, float *gains, int len) |
25 | 31.9k | { |
26 | 31.9k | const float g0 = gains[0]; |
27 | 31.9k | const float g1 = gains[1]; |
28 | 31.9k | const float g2 = gains[2]; |
29 | | |
30 | 31.9k | float x4 = data[-period - 2]; |
31 | 31.9k | float x3 = data[-period - 1]; |
32 | 31.9k | float x2 = data[-period + 0]; |
33 | 31.9k | float x1 = data[-period + 1]; |
34 | | |
35 | 21.8M | for (int i = 0; i < len; i++) { |
36 | 21.8M | float x0 = data[i - period + 2]; |
37 | 21.8M | data[i] += g0 * x2 + |
38 | 21.8M | g1 * (x1 + x3) + |
39 | 21.8M | g2 * (x0 + x4); |
40 | 21.8M | x4 = x3; |
41 | 21.8M | x3 = x2; |
42 | 21.8M | x2 = x1; |
43 | 21.8M | x1 = x0; |
44 | 21.8M | } |
45 | 31.9k | } |
46 | | |
47 | | static float deemphasis_c(float *y, float *x, float coeff, const float *weights, int len) |
48 | 155k | { |
49 | 155k | const float c = weights[0]; |
50 | 87.9M | for (int i = 0; i < len; i++) |
51 | 87.8M | coeff = y[i] = x[i] + coeff*c; |
52 | | |
53 | 155k | return coeff; |
54 | 155k | } |
55 | | |
56 | | av_cold void ff_opus_dsp_init(OpusDSP *ctx) |
57 | 71.0k | { |
58 | 71.0k | ctx->postfilter = postfilter_c; |
59 | 71.0k | ctx->deemphasis = deemphasis_c; |
60 | | |
61 | | #if ARCH_AARCH64 |
62 | | ff_opus_dsp_init_aarch64(ctx); |
63 | | #elif ARCH_RISCV |
64 | | ff_opus_dsp_init_riscv(ctx); |
65 | | #elif ARCH_X86 && HAVE_X86ASM |
66 | | ff_opus_dsp_init_x86(ctx); |
67 | | #endif |
68 | 71.0k | } |