1/* 2 * Copyright (c) 2004 Gildas Bazin 3 * Copyright (c) 2010 Mans Rullgard <mans@mansr.com> 4 * 5 * This file is part of FFmpeg. 6 * 7 * FFmpeg is free software; you can redistribute it and/or 8 * modify it under the terms of the GNU Lesser General Public 9 * License as published by the Free Software Foundation; either 10 * version 2.1 of the License, or (at your option) any later version. 11 * 12 * FFmpeg is distributed in the hope that it will be useful, 13 * but WITHOUT ANY WARRANTY; without even the implied warranty of 14 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU 15 * Lesser General Public License for more details. 16 * 17 * You should have received a copy of the GNU Lesser General Public 18 * License along with FFmpeg; if not, write to the Free Software 19 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA 20 */ 21 22#include "config.h" 23#include "dcadsp.h" 24 25static void dca_lfe_fir_c(float *out, const float *in, const float *coefs, 26 int decifactor, float scale, float bias) 27{ 28 float *out2 = out + decifactor; 29 const float *cf0 = coefs; 30 const float *cf1 = coefs + 256; 31 int j, k; 32 33 /* One decimated sample generates 2*decifactor interpolated ones */ 34 for (k = 0; k < decifactor; k++) { 35 float v0 = 0.0; 36 float v1 = 0.0; 37 for (j = 0; j < 256 / decifactor; j++) { 38 float s = in[-j]; 39 v0 += s * *cf0++; 40 v1 += s * *--cf1; 41 } 42 *out++ = (v0 * scale) + bias; 43 *out2++ = (v1 * scale) + bias; 44 } 45} 46 47void ff_dcadsp_init(DCADSPContext *s) 48{ 49 s->lfe_fir = dca_lfe_fir_c; 50 if (ARCH_ARM) ff_dcadsp_init_arm(s); 51} 52