1 /*
2 * Copyright (c) 2015 Peter Meerwald <pmeerw@pmeerw.net>
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 "g722dsp.h"
22 #include "mathops.h"
23
24 /*
25 * quadrature mirror filter (QMF) coefficients (ITU-T G.722 Table 11) inlined
26 * in code below: 3, -11, 12, 32, -210, 951, 3876, -805, 362, -156, 53, -11
27 */
28
g722_apply_qmf(const int16_t * prev_samples,int xout[2])29 static void g722_apply_qmf(const int16_t *prev_samples, int xout[2])
30 {
31 xout[1] = MUL16(*prev_samples++, 3);
32 xout[0] = MUL16(*prev_samples++, -11);
33
34 MAC16(xout[1], *prev_samples++, -11);
35 MAC16(xout[0], *prev_samples++, 53);
36
37 MAC16(xout[1], *prev_samples++, 12);
38 MAC16(xout[0], *prev_samples++, -156);
39
40 MAC16(xout[1], *prev_samples++, 32);
41 MAC16(xout[0], *prev_samples++, 362);
42
43 MAC16(xout[1], *prev_samples++, -210);
44 MAC16(xout[0], *prev_samples++, -805);
45
46 MAC16(xout[1], *prev_samples++, 951);
47 MAC16(xout[0], *prev_samples++, 3876);
48
49 MAC16(xout[1], *prev_samples++, 3876);
50 MAC16(xout[0], *prev_samples++, 951);
51
52 MAC16(xout[1], *prev_samples++, -805);
53 MAC16(xout[0], *prev_samples++, -210);
54
55 MAC16(xout[1], *prev_samples++, 362);
56 MAC16(xout[0], *prev_samples++, 32);
57
58 MAC16(xout[1], *prev_samples++, -156);
59 MAC16(xout[0], *prev_samples++, 12);
60
61 MAC16(xout[1], *prev_samples++, 53);
62 MAC16(xout[0], *prev_samples++, -11);
63
64 MAC16(xout[1], *prev_samples++, -11);
65 MAC16(xout[0], *prev_samples++, 3);
66 }
67
ff_g722dsp_init(G722DSPContext * c)68 av_cold void ff_g722dsp_init(G722DSPContext *c)
69 {
70 c->apply_qmf = g722_apply_qmf;
71
72 if (ARCH_ARM)
73 ff_g722dsp_init_arm(c);
74 if (ARCH_X86)
75 ff_g722dsp_init_x86(c);
76 }
77