mpc7.c 11.3 KB
Newer Older
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22
/*
 * Musepack SV7 decoder
 * Copyright (c) 2006 Konstantin Shishkov
 *
 * This file is part of FFmpeg.
 *
 * FFmpeg 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.
 *
 * FFmpeg 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 FFmpeg; if not, write to the Free Software
 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
 */

/**
23
 * @file
24 25 26 27
 * MPEG Audio Layer 1/2 -like codec with frames of 1152 samples
 * divided into 32 subbands.
 */

28
#include "libavutil/channel_layout.h"
29
#include "libavutil/internal.h"
30
#include "libavutil/lfg.h"
31
#include "avcodec.h"
32
#include "get_bits.h"
33
#include "internal.h"
34
#include "mpegaudiodsp.h"
35 36 37 38 39 40

#include "mpc.h"
#include "mpc7data.h"

static VLC scfi_vlc, dscf_vlc, hdr_vlc, quant_vlc[MPC7_QUANT_VLC_TABLES][2];

41
static const uint16_t quant_offsets[MPC7_QUANT_VLC_TABLES*2 + 1] =
42 43 44 45 46 47
{
       0, 512, 1024, 1536, 2052, 2564, 3076, 3588, 4100, 4612, 5124,
       5636, 6164, 6676, 7224
};


48
static av_cold int mpc7_decode_init(AVCodecContext * avctx)
49 50 51 52
{
    int i, j;
    MPCContext *c = avctx->priv_data;
    GetBitContext gb;
53
    LOCAL_ALIGNED_16(uint8_t, buf, [16]);
54
    static int vlc_initialized = 0;
55

56 57 58 59 60
    static VLC_TYPE scfi_table[1 << MPC7_SCFI_BITS][2];
    static VLC_TYPE dscf_table[1 << MPC7_DSCF_BITS][2];
    static VLC_TYPE hdr_table[1 << MPC7_HDR_BITS][2];
    static VLC_TYPE quant_tables[7224][2];

61 62
    /* Musepack SV7 is always stereo */
    if (avctx->channels != 2) {
63
        avpriv_request_sample(avctx, "%d channels", avctx->channels);
64 65 66
        return AVERROR_PATCHWELCOME;
    }

67 68 69 70 71
    if(avctx->extradata_size < 16){
        av_log(avctx, AV_LOG_ERROR, "Too small extradata size (%i)!\n", avctx->extradata_size);
        return -1;
    }
    memset(c->oldDSCF, 0, sizeof(c->oldDSCF));
72
    av_lfg_init(&c->rnd, 0xDEADBEEF);
73
    ff_bswapdsp_init(&c->bdsp);
74
    ff_mpadsp_init(&c->mpadsp);
75
    c->bdsp.bswap_buf((uint32_t *) buf, (const uint32_t *) avctx->extradata, 4);
76 77 78 79 80 81 82 83 84 85
    ff_mpc_init();
    init_get_bits(&gb, buf, 128);

    c->IS = get_bits1(&gb);
    c->MSS = get_bits1(&gb);
    c->maxbands = get_bits(&gb, 6);
    if(c->maxbands >= BANDS){
        av_log(avctx, AV_LOG_ERROR, "Too many bands: %i\n", c->maxbands);
        return -1;
    }
86
    skip_bits_long(&gb, 88);
87 88 89 90 91 92
    c->gapless = get_bits1(&gb);
    c->lastframelen = get_bits(&gb, 11);
    av_log(avctx, AV_LOG_DEBUG, "IS: %d, MSS: %d, TG: %d, LFL: %d, bands: %d\n",
            c->IS, c->MSS, c->gapless, c->lastframelen, c->maxbands);
    c->frames_to_skip = 0;

93
    avctx->sample_fmt = AV_SAMPLE_FMT_S16P;
94
    avctx->channel_layout = AV_CH_LAYOUT_STEREO;
95

96
    if(vlc_initialized) return 0;
97
    av_log(avctx, AV_LOG_DEBUG, "Initing VLC\n");
98 99
    scfi_vlc.table = scfi_table;
    scfi_vlc.table_allocated = 1 << MPC7_SCFI_BITS;
100 101
    if(init_vlc(&scfi_vlc, MPC7_SCFI_BITS, MPC7_SCFI_SIZE,
                &mpc7_scfi[1], 2, 1,
102
                &mpc7_scfi[0], 2, 1, INIT_VLC_USE_NEW_STATIC)){
103 104 105
        av_log(avctx, AV_LOG_ERROR, "Cannot init SCFI VLC\n");
        return -1;
    }
106 107
    dscf_vlc.table = dscf_table;
    dscf_vlc.table_allocated = 1 << MPC7_DSCF_BITS;
108 109
    if(init_vlc(&dscf_vlc, MPC7_DSCF_BITS, MPC7_DSCF_SIZE,
                &mpc7_dscf[1], 2, 1,
110
                &mpc7_dscf[0], 2, 1, INIT_VLC_USE_NEW_STATIC)){
111 112 113
        av_log(avctx, AV_LOG_ERROR, "Cannot init DSCF VLC\n");
        return -1;
    }
114 115
    hdr_vlc.table = hdr_table;
    hdr_vlc.table_allocated = 1 << MPC7_HDR_BITS;
116 117
    if(init_vlc(&hdr_vlc, MPC7_HDR_BITS, MPC7_HDR_SIZE,
                &mpc7_hdr[1], 2, 1,
118
                &mpc7_hdr[0], 2, 1, INIT_VLC_USE_NEW_STATIC)){
119 120 121 122 123
        av_log(avctx, AV_LOG_ERROR, "Cannot init HDR VLC\n");
        return -1;
    }
    for(i = 0; i < MPC7_QUANT_VLC_TABLES; i++){
        for(j = 0; j < 2; j++){
124 125
            quant_vlc[i][j].table = &quant_tables[quant_offsets[i*2 + j]];
            quant_vlc[i][j].table_allocated = quant_offsets[i*2 + j + 1] - quant_offsets[i*2 + j];
126 127
            if(init_vlc(&quant_vlc[i][j], 9, mpc7_quant_vlc_sizes[i],
                        &mpc7_quant_vlc[i][j][1], 4, 2,
128
                        &mpc7_quant_vlc[i][j][0], 4, 2, INIT_VLC_USE_NEW_STATIC)){
129 130 131 132 133
                av_log(avctx, AV_LOG_ERROR, "Cannot init QUANT VLC %i,%i\n",i,j);
                return -1;
            }
        }
    }
134
    vlc_initialized = 1;
135

136 137 138 139 140 141 142 143 144 145 146 147
    return 0;
}

/**
 * Fill samples for given subband
 */
static inline void idx_to_quant(MPCContext *c, GetBitContext *gb, int idx, int *dst)
{
    int i, i1, t;
    switch(idx){
    case -1:
        for(i = 0; i < SAMPLES_PER_BAND; i++){
148
            *dst++ = (av_lfg_get(&c->rnd) & 0x3FC) - 510;
149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183
        }
        break;
    case 1:
        i1 = get_bits1(gb);
        for(i = 0; i < SAMPLES_PER_BAND/3; i++){
            t = get_vlc2(gb, quant_vlc[0][i1].table, 9, 2);
            *dst++ = mpc7_idx30[t];
            *dst++ = mpc7_idx31[t];
            *dst++ = mpc7_idx32[t];
        }
        break;
    case 2:
        i1 = get_bits1(gb);
        for(i = 0; i < SAMPLES_PER_BAND/2; i++){
            t = get_vlc2(gb, quant_vlc[1][i1].table, 9, 2);
            *dst++ = mpc7_idx50[t];
            *dst++ = mpc7_idx51[t];
        }
        break;
    case  3: case  4: case  5: case  6: case  7:
        i1 = get_bits1(gb);
        for(i = 0; i < SAMPLES_PER_BAND; i++)
            *dst++ = get_vlc2(gb, quant_vlc[idx-1][i1].table, 9, 2) - mpc7_quant_vlc_off[idx-1];
        break;
    case  8: case  9: case 10: case 11: case 12:
    case 13: case 14: case 15: case 16: case 17:
        t = (1 << (idx - 2)) - 1;
        for(i = 0; i < SAMPLES_PER_BAND; i++)
            *dst++ = get_bits(gb, idx - 1) - t;
        break;
    default: // case 0 and -2..-17
        return;
    }
}

184 185 186 187 188 189 190 191
static int get_scale_idx(GetBitContext *gb, int ref)
{
    int t = get_vlc2(gb, dscf_vlc.table, MPC7_DSCF_BITS, 1) - 7;
    if (t == 8)
        return get_bits(gb, 6);
    return ref + t;
}

192 193
static int mpc7_decode_frame(AVCodecContext * avctx, void *data,
                             int *got_frame_ptr, AVPacket *avpkt)
194
{
195
    AVFrame *frame     = data;
196
    const uint8_t *buf = avpkt->data;
197
    int buf_size;
198 199
    MPCContext *c = avctx->priv_data;
    GetBitContext gb;
200
    int i, ch;
201 202
    int mb = -1;
    Band *bands = c->bands;
203
    int off, ret, last_frame, skip;
204 205
    int bits_used, bits_avail;

206
    memset(bands, 0, sizeof(*bands) * (c->maxbands + 1));
207 208 209 210 211 212 213 214 215 216

    buf_size = avpkt->size & ~3;
    if (buf_size <= 0) {
        av_log(avctx, AV_LOG_ERROR, "packet size is too small (%i bytes)\n",
               avpkt->size);
        return AVERROR_INVALIDDATA;
    }
    if (buf_size != avpkt->size) {
        av_log(avctx, AV_LOG_WARNING, "packet size is not a multiple of 4. "
               "extra bytes at the end will be skipped.\n");
217 218
    }

219 220 221 222 223
    skip       = buf[0];
    last_frame = buf[1];
    buf       += 4;
    buf_size  -= 4;

224
    /* get output buffer */
225
    frame->nb_samples = MPC_FRAME_SIZE;
226
    if ((ret = ff_get_buffer(avctx, frame, 0)) < 0)
227
        return ret;
228

229 230
    av_fast_padded_malloc(&c->bits, &c->buf_size, buf_size);
    if (!c->bits)
231
        return AVERROR(ENOMEM);
232 233
    c->bdsp.bswap_buf((uint32_t *) c->bits, (const uint32_t *) buf,
                      buf_size >> 2);
234
    init_get_bits(&gb, c->bits, buf_size * 8);
235
    skip_bits_long(&gb, skip);
236 237 238 239

    /* read subband indexes */
    for(i = 0; i <= c->maxbands; i++){
        for(ch = 0; ch < 2; ch++){
240
            int t = 4;
241
            if(i) t = get_vlc2(&gb, hdr_vlc.table, MPC7_HDR_BITS, 1) - 5;
242
            if(t == 4) bands[i].res[ch] = get_bits(&gb, 4);
243
            else bands[i].res[ch] = bands[i-1].res[ch] + t;
244 245 246 247
            if (bands[i].res[ch] < -1 || bands[i].res[ch] > 17) {
                av_log(avctx, AV_LOG_ERROR, "subband index invalid\n");
                return AVERROR_INVALIDDATA;
            }
248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263
        }

        if(bands[i].res[0] || bands[i].res[1]){
            mb = i;
            if(c->MSS) bands[i].msf = get_bits1(&gb);
        }
    }
    /* get scale indexes coding method */
    for(i = 0; i <= mb; i++)
        for(ch = 0; ch < 2; ch++)
            if(bands[i].res[ch]) bands[i].scfi[ch] = get_vlc2(&gb, scfi_vlc.table, MPC7_SCFI_BITS, 1);
    /* get scale indexes */
    for(i = 0; i <= mb; i++){
        for(ch = 0; ch < 2; ch++){
            if(bands[i].res[ch]){
                bands[i].scf_idx[ch][2] = c->oldDSCF[ch][i];
264
                bands[i].scf_idx[ch][0] = get_scale_idx(&gb, bands[i].scf_idx[ch][2]);
265 266
                switch(bands[i].scfi[ch]){
                case 0:
267 268
                    bands[i].scf_idx[ch][1] = get_scale_idx(&gb, bands[i].scf_idx[ch][0]);
                    bands[i].scf_idx[ch][2] = get_scale_idx(&gb, bands[i].scf_idx[ch][1]);
269 270
                    break;
                case 1:
271
                    bands[i].scf_idx[ch][1] = get_scale_idx(&gb, bands[i].scf_idx[ch][0]);
272 273 274 275
                    bands[i].scf_idx[ch][2] = bands[i].scf_idx[ch][1];
                    break;
                case 2:
                    bands[i].scf_idx[ch][1] = bands[i].scf_idx[ch][0];
276
                    bands[i].scf_idx[ch][2] = get_scale_idx(&gb, bands[i].scf_idx[ch][1]);
277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292
                    break;
                case 3:
                    bands[i].scf_idx[ch][2] = bands[i].scf_idx[ch][1] = bands[i].scf_idx[ch][0];
                    break;
                }
                c->oldDSCF[ch][i] = bands[i].scf_idx[ch][2];
            }
        }
    }
    /* get quantizers */
    memset(c->Q, 0, sizeof(c->Q));
    off = 0;
    for(i = 0; i < BANDS; i++, off += SAMPLES_PER_BAND)
        for(ch = 0; ch < 2; ch++)
            idx_to_quant(c, &gb, bands[i].res[ch], c->Q[ch] + off);

293
    ff_mpc_dequantize_and_synth(c, mb, (int16_t **)frame->extended_data, 2);
294
    if(last_frame)
295
        frame->nb_samples = c->lastframelen;
296 297

    bits_used = get_bits_count(&gb);
298 299
    bits_avail = buf_size * 8;
    if (!last_frame && ((bits_avail < bits_used) || (bits_used + 32 <= bits_avail))) {
300
        av_log(avctx, AV_LOG_ERROR, "Error decoding frame: used %i of %i bits\n", bits_used, bits_avail);
301 302 303 304
        return -1;
    }
    if(c->frames_to_skip){
        c->frames_to_skip--;
305
        *got_frame_ptr = 0;
306
        return avpkt->size;
307
    }
308

309
    *got_frame_ptr = 1;
310

311
    return avpkt->size;
312 313 314 315 316 317 318 319 320 321
}

static void mpc7_decode_flush(AVCodecContext *avctx)
{
    MPCContext *c = avctx->priv_data;

    memset(c->oldDSCF, 0, sizeof(c->oldDSCF));
    c->frames_to_skip = 32;
}

322 323 324
static av_cold int mpc7_decode_close(AVCodecContext *avctx)
{
    MPCContext *c = avctx->priv_data;
325 326
    av_freep(&c->bits);
    c->buf_size = 0;
327 328 329
    return 0;
}

330
AVCodec ff_mpc7_decoder = {
331
    .name           = "mpc7",
332
    .long_name      = NULL_IF_CONFIG_SMALL("Musepack SV7"),
333
    .type           = AVMEDIA_TYPE_AUDIO,
334
    .id             = AV_CODEC_ID_MUSEPACK7,
335 336
    .priv_data_size = sizeof(MPCContext),
    .init           = mpc7_decode_init,
337
    .close          = mpc7_decode_close,
338
    .decode         = mpc7_decode_frame,
339
    .flush          = mpc7_decode_flush,
340
    .capabilities   = AV_CODEC_CAP_DR1,
341 342
    .sample_fmts    = (const enum AVSampleFormat[]) { AV_SAMPLE_FMT_S16P,
                                                      AV_SAMPLE_FMT_NONE },
343
};