1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
|
/*
* Copyright (c) 2015 Peter Meerwald <pmeerw@pmeerw.net>
*
* This file is part of Libav.
*
* Libav is free software; you can redistribute it and/or
* modify it under the terms of the GNU Lesser General Public
* License as published by the Free Software Foundation; either
* version 2.1 of the License, or (at your option) any later version.
*
* Libav is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
* Lesser General Public License for more details.
*
* You should have received a copy of the GNU Lesser General Public
* License along with Libav; if not, write to the Free Software
* Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
*/
#include "g722dsp.h"
#include "mathops.h"
/*
* quadrature mirror filter (QMF) coefficients (ITU-T G.722 Table 11) inlined
* in code below: 3, -11, 12, 32, -210, 951, 3876, -805, 362, -156, 53, -11
*/
static void g722_apply_qmf(const int16_t *prev_samples, int xout[2])
{
xout[1] = MUL16(*prev_samples++, 3);
xout[0] = MUL16(*prev_samples++, -11);
MAC16(xout[1], *prev_samples++, -11);
MAC16(xout[0], *prev_samples++, 53);
MAC16(xout[1], *prev_samples++, 12);
MAC16(xout[0], *prev_samples++, -156);
MAC16(xout[1], *prev_samples++, 32);
MAC16(xout[0], *prev_samples++, 362);
MAC16(xout[1], *prev_samples++, -210);
MAC16(xout[0], *prev_samples++, -805);
MAC16(xout[1], *prev_samples++, 951);
MAC16(xout[0], *prev_samples++, 3876);
MAC16(xout[1], *prev_samples++, 3876);
MAC16(xout[0], *prev_samples++, 951);
MAC16(xout[1], *prev_samples++, -805);
MAC16(xout[0], *prev_samples++, -210);
MAC16(xout[1], *prev_samples++, 362);
MAC16(xout[0], *prev_samples++, 32);
MAC16(xout[1], *prev_samples++, -156);
MAC16(xout[0], *prev_samples++, 12);
MAC16(xout[1], *prev_samples++, 53);
MAC16(xout[0], *prev_samples++, -11);
MAC16(xout[1], *prev_samples++, -11);
MAC16(xout[0], *prev_samples++, 3);
}
av_cold void ff_g722dsp_init(G722DSPContext *c)
{
c->apply_qmf = g722_apply_qmf;
if (ARCH_ARM)
ff_g722dsp_init_arm(c);
}
|
