Go to most recent revision | Details | Last modification | View Log | RSS feed
| Rev | Author | Line No. | Line |
|---|---|---|---|
| 4349 | Serge | 1 | /* |
| 2 | * Copyright (c) 2004 Gildas Bazin |
||
| 3 | * Copyright (c) 2010 Mans Rullgard |
||
| 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 "libavutil/attributes.h" |
||
| 24 | #include "libavutil/intreadwrite.h" |
||
| 25 | #include "dcadsp.h" |
||
| 26 | |||
| 27 | static void dca_lfe_fir_c(float *out, const float *in, const float *coefs, |
||
| 28 | int decifactor, float scale) |
||
| 29 | { |
||
| 30 | float *out2 = out + decifactor; |
||
| 31 | const float *cf0 = coefs; |
||
| 32 | const float *cf1 = coefs + 256; |
||
| 33 | int j, k; |
||
| 34 | |||
| 35 | /* One decimated sample generates 2*decifactor interpolated ones */ |
||
| 36 | for (k = 0; k < decifactor; k++) { |
||
| 37 | float v0 = 0.0; |
||
| 38 | float v1 = 0.0; |
||
| 39 | for (j = 0; j < 256 / decifactor; j++) { |
||
| 40 | float s = in[-j]; |
||
| 41 | v0 += s * *cf0++; |
||
| 42 | v1 += s * *--cf1; |
||
| 43 | } |
||
| 44 | *out++ = v0 * scale; |
||
| 45 | *out2++ = v1 * scale; |
||
| 46 | } |
||
| 47 | } |
||
| 48 | |||
| 49 | static void dca_qmf_32_subbands(float samples_in[32][8], int sb_act, |
||
| 50 | SynthFilterContext *synth, FFTContext *imdct, |
||
| 51 | float synth_buf_ptr[512], |
||
| 52 | int *synth_buf_offset, float synth_buf2[32], |
||
| 53 | const float window[512], float *samples_out, |
||
| 54 | float raXin[32], float scale) |
||
| 55 | { |
||
| 56 | int i; |
||
| 57 | int subindex; |
||
| 58 | |||
| 59 | for (i = sb_act; i < 32; i++) |
||
| 60 | raXin[i] = 0.0; |
||
| 61 | |||
| 62 | /* Reconstructed channel sample index */ |
||
| 63 | for (subindex = 0; subindex < 8; subindex++) { |
||
| 64 | /* Load in one sample from each subband and clear inactive subbands */ |
||
| 65 | for (i = 0; i < sb_act; i++) { |
||
| 66 | unsigned sign = (i - 1) & 2; |
||
| 67 | uint32_t v = AV_RN32A(&samples_in[i][subindex]) ^ sign << 30; |
||
| 68 | AV_WN32A(&raXin[i], v); |
||
| 69 | } |
||
| 70 | |||
| 71 | synth->synth_filter_float(imdct, synth_buf_ptr, synth_buf_offset, |
||
| 72 | synth_buf2, window, samples_out, raXin, scale); |
||
| 73 | samples_out += 32; |
||
| 74 | } |
||
| 75 | } |
||
| 76 | |||
| 77 | av_cold void ff_dcadsp_init(DCADSPContext *s) |
||
| 78 | { |
||
| 79 | s->lfe_fir = dca_lfe_fir_c; |
||
| 80 | s->qmf_32_subbands = dca_qmf_32_subbands; |
||
| 81 | if (ARCH_ARM) ff_dcadsp_init_arm(s); |
||
| 82 | } |