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 29static 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 68av_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#elif ARCH_X86 75 ff_g722dsp_init_x86(c); 76#endif 77} 78