Line | Branch | Exec | Source |
---|---|---|---|
1 | /* | ||
2 | * Musepack decoder core | ||
3 | * Copyright (c) 2006 Konstantin Shishkov | ||
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 | /** | ||
23 | * @file | ||
24 | * Musepack decoder core | ||
25 | * MPEG Audio Layer 1/2 -like codec with frames of 1152 samples | ||
26 | * divided into 32 subbands. | ||
27 | */ | ||
28 | |||
29 | #include <string.h> | ||
30 | #include "libavutil/common.h" | ||
31 | #include "mpegaudiodsp.h" | ||
32 | |||
33 | #include "mpc.h" | ||
34 | #include "mpcdata.h" | ||
35 | |||
36 | /** | ||
37 | * Process decoded Musepack data and produce PCM | ||
38 | */ | ||
39 | 914 | static void mpc_synth(MPCContext *c, int16_t **out, int channels) | |
40 | { | ||
41 | 914 | int dither_state = 0; | |
42 | int i, ch; | ||
43 | |||
44 |
2/2✓ Branch 0 taken 1828 times.
✓ Branch 1 taken 914 times.
|
2742 | for(ch = 0; ch < channels; ch++){ |
45 |
2/2✓ Branch 0 taken 65808 times.
✓ Branch 1 taken 1828 times.
|
67636 | for(i = 0; i < SAMPLES_PER_BAND; i++) { |
46 | 65808 | ff_mpa_synth_filter_fixed(&c->mpadsp, | |
47 | 65808 | c->synth_buf[ch], &(c->synth_buf_offset[ch]), | |
48 | ff_mpa_synth_window_fixed, &dither_state, | ||
49 | 65808 | out[ch] + 32 * i, 1, | |
50 | 65808 | c->sb_samples[ch][i]); | |
51 | } | ||
52 | } | ||
53 | 914 | } | |
54 | |||
55 | 914 | void ff_mpc_dequantize_and_synth(MPCContext * c, int maxband, int16_t **out, | |
56 | int channels) | ||
57 | { | ||
58 | int i, j, ch; | ||
59 | 914 | Band *bands = c->bands; | |
60 | int off; | ||
61 | float mul; | ||
62 | |||
63 | /* dequantize */ | ||
64 | 914 | memset(c->sb_samples, 0, sizeof(c->sb_samples)); | |
65 | 914 | off = 0; | |
66 |
2/2✓ Branch 0 taken 20828 times.
✓ Branch 1 taken 914 times.
|
21742 | for(i = 0; i <= maxband; i++, off += SAMPLES_PER_BAND){ |
67 |
2/2✓ Branch 0 taken 41656 times.
✓ Branch 1 taken 20828 times.
|
62484 | for(ch = 0; ch < 2; ch++){ |
68 |
2/2✓ Branch 0 taken 36564 times.
✓ Branch 1 taken 5092 times.
|
41656 | if(bands[i].res[ch]){ |
69 | 36564 | j = 0; | |
70 | 36564 | mul = (mpc_CC+1)[bands[i].res[ch]] * mpc_SCF[bands[i].scf_idx[ch][0] & 0xFF]; | |
71 |
2/2✓ Branch 0 taken 438768 times.
✓ Branch 1 taken 36564 times.
|
475332 | for(; j < 12; j++) |
72 | 438768 | c->sb_samples[ch][j][i] = av_clipf(mul * c->Q[ch][j + off], INT32_MIN, INT32_MAX); | |
73 | 36564 | mul = (mpc_CC+1)[bands[i].res[ch]] * mpc_SCF[bands[i].scf_idx[ch][1] & 0xFF]; | |
74 |
2/2✓ Branch 0 taken 438768 times.
✓ Branch 1 taken 36564 times.
|
475332 | for(; j < 24; j++) |
75 | 438768 | c->sb_samples[ch][j][i] = av_clipf(mul * c->Q[ch][j + off], INT32_MIN, INT32_MAX); | |
76 | 36564 | mul = (mpc_CC+1)[bands[i].res[ch]] * mpc_SCF[bands[i].scf_idx[ch][2] & 0xFF]; | |
77 |
2/2✓ Branch 0 taken 438768 times.
✓ Branch 1 taken 36564 times.
|
475332 | for(; j < 36; j++) |
78 | 438768 | c->sb_samples[ch][j][i] = av_clipf(mul * c->Q[ch][j + off], INT32_MIN, INT32_MAX); | |
79 | } | ||
80 | } | ||
81 |
2/2✓ Branch 0 taken 12250 times.
✓ Branch 1 taken 8578 times.
|
20828 | if(bands[i].msf){ |
82 | unsigned t1, t2; | ||
83 |
2/2✓ Branch 0 taken 441000 times.
✓ Branch 1 taken 12250 times.
|
453250 | for(j = 0; j < SAMPLES_PER_BAND; j++){ |
84 | 441000 | t1 = c->sb_samples[0][j][i]; | |
85 | 441000 | t2 = c->sb_samples[1][j][i]; | |
86 | 441000 | c->sb_samples[0][j][i] = t1 + t2; | |
87 | 441000 | c->sb_samples[1][j][i] = t1 - t2; | |
88 | } | ||
89 | } | ||
90 | } | ||
91 | |||
92 | 914 | mpc_synth(c, out, channels); | |
93 | 914 | } | |
94 |