mpegaudiodec.c 73.1 KB
Newer Older
Fabrice Bellard's avatar
Fabrice Bellard committed
1 2
/*
 * MPEG Audio decoder
3
 * Copyright (c) 2001, 2002 Fabrice Bellard
Fabrice Bellard's avatar
Fabrice Bellard committed
4
 *
5 6 7
 * This file is part of FFmpeg.
 *
 * FFmpeg is free software; you can redistribute it and/or
8 9
 * modify it under the terms of the GNU Lesser General Public
 * License as published by the Free Software Foundation; either
10
 * version 2.1 of the License, or (at your option) any later version.
Fabrice Bellard's avatar
Fabrice Bellard committed
11
 *
12
 * FFmpeg is distributed in the hope that it will be useful,
Fabrice Bellard's avatar
Fabrice Bellard committed
13
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
14 15
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 * Lesser General Public License for more details.
Fabrice Bellard's avatar
Fabrice Bellard committed
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
Fabrice Bellard's avatar
Fabrice Bellard committed
20
 */
Michael Niedermayer's avatar
Michael Niedermayer committed
21 22

/**
23
 * @file
Michael Niedermayer's avatar
Michael Niedermayer committed
24
 * MPEG Audio decoder.
25
 */
Michael Niedermayer's avatar
Michael Niedermayer committed
26

Fabrice Bellard's avatar
Fabrice Bellard committed
27
#include "avcodec.h"
28
#include "get_bits.h"
29
#include "dsputil.h"
30
#include "libavformat/id3v1.h"
Fabrice Bellard's avatar
Fabrice Bellard committed
31 32

/*
33 34 35
 * TODO:
 *  - in low precision mode, use more 16 bit multiplies in synth filter
 *  - test lsf / mpeg25 extensively.
Fabrice Bellard's avatar
Fabrice Bellard committed
36 37
 */

Roberto Togni's avatar
Roberto Togni committed
38
#include "mpegaudio.h"
39
#include "mpegaudiodecheader.h"
40

Luca Barbato's avatar
Luca Barbato committed
41 42
#include "mathops.h"

43
#if CONFIG_FLOAT
44
#   define SHR(a,b)       ((a)*(1.0f/(1<<(b))))
45 46
#   define compute_antialias compute_antialias_float
#   define FIXR_OLD(a)    ((int)((a) * FRAC_ONE + 0.5))
47 48
#   define FIXR(x)        ((float)(x))
#   define FIXHR(x)       ((float)(x))
49 50 51 52 53 54
#   define MULH3(x, y, s) ((s)*(y)*(x))
#   define MULLx(x, y, s) ((y)*(x))
#   define RENAME(a) a ## _float
#else
#   define SHR(a,b)       ((a)>>(b))
#   define compute_antialias compute_antialias_integer
55
/* WARNING: only correct for posititive numbers */
56 57 58 59 60 61 62
#   define FIXR_OLD(a)    ((int)((a) * FRAC_ONE + 0.5))
#   define FIXR(a)        ((int)((a) * FRAC_ONE + 0.5))
#   define FIXHR(a)       ((int)((a) * (1LL<<32) + 0.5))
#   define MULH3(x, y, s) MULH((s)*(x), y)
#   define MULLx(x, y, s) MULL(x,y,s)
#   define RENAME(a)      a
#endif
63

64 65
/****************/

Fabrice Bellard's avatar
Fabrice Bellard committed
66 67
#define HEADER_SIZE 4

68
#include "mpegaudiodata.h"
69 70
#include "mpegaudiodectab.h"

71 72 73 74 75 76
#if CONFIG_FLOAT
#    include "fft.h"
#else
#    include "dct32.c"
#endif

77
static void compute_antialias(MPADecodeContext *s, GranuleDef *g);
78 79
static void apply_window_mp3_c(MPA_INT *synth_buf, MPA_INT *window,
                               int *dither_state, OUT_INT *samples, int incr);
80

81
/* vlc structure for decoding layer 3 huffman tables */
82
static VLC huff_vlc[16];
83 84 85 86 87 88 89 90
static VLC_TYPE huff_vlc_tables[
  0+128+128+128+130+128+154+166+
  142+204+190+170+542+460+662+414
  ][2];
static const int huff_vlc_tables_sizes[16] = {
  0, 128, 128, 128, 130, 128, 154, 166,
  142, 204, 190, 170, 542, 460, 662, 414
};
91
static VLC huff_quad_vlc[2];
92 93 94 95
static VLC_TYPE huff_quad_vlc_tables[128+16][2];
static const int huff_quad_vlc_tables_sizes[2] = {
  128, 16
};
96
/* computed from band_size_long */
97
static uint16_t band_index_long[9][23];
98
#include "mpegaudio_tablegen.h"
99
/* intensity stereo coef table */
100 101
static INTFLOAT is_table[2][16];
static INTFLOAT is_table_lsf[2][2][16];
102 103
static int32_t csa_table[8][4];
static float csa_table_float[8][4];
104
static INTFLOAT mdct_win[8][36];
105

106 107 108 109 110 111 112 113
static int16_t division_tab3[1<<6 ];
static int16_t division_tab5[1<<8 ];
static int16_t division_tab9[1<<11];

static int16_t * const division_tabs[4] = {
    division_tab3, division_tab5, NULL, division_tab9
};

114
/* lower 2 bits: modulo 3, higher bits: shift */
115
static uint16_t scale_factor_modshift[64];
116
/* [i][j]:  2^(-j/3) * FRAC_ONE * 2^(i+2) / (2^(i+2) - 1) */
117
static int32_t scale_factor_mult[15][3];
118 119 120
/* mult table for layer 2 group quantization */

#define SCALE_GEN(v) \
121
{ FIXR_OLD(1.0 * (v)), FIXR_OLD(0.7937005259 * (v)), FIXR_OLD(0.6299605249 * (v)) }
122

Michael Niedermayer's avatar
Michael Niedermayer committed
123
static const int32_t scale_factor_mult2[3][3] = {
124 125 126
    SCALE_GEN(4.0 / 3.0), /* 3 steps */
    SCALE_GEN(4.0 / 5.0), /* 5 steps */
    SCALE_GEN(4.0 / 9.0), /* 9 steps */
127 128
};

129
DECLARE_ALIGNED(16, MPA_INT, RENAME(ff_mpa_synth_window))[512+256];
130

131 132 133 134
/**
 * Convert region offsets to region sizes and truncate
 * size to big_values.
 */
135
static void ff_region_offset2size(GranuleDef *g){
136 137 138 139 140 141 142 143 144
    int i, k, j=0;
    g->region_size[2] = (576 / 2);
    for(i=0;i<3;i++) {
        k = FFMIN(g->region_size[i], g->big_values);
        g->region_size[i] = k - j;
        j = k;
    }
}

145
static void ff_init_short_region(MPADecodeContext *s, GranuleDef *g){
146 147 148 149 150 151 152 153 154 155 156 157 158
    if (g->block_type == 2)
        g->region_size[0] = (36 / 2);
    else {
        if (s->sample_rate_index <= 2)
            g->region_size[0] = (36 / 2);
        else if (s->sample_rate_index != 8)
            g->region_size[0] = (54 / 2);
        else
            g->region_size[0] = (108 / 2);
    }
    g->region_size[1] = (576 / 2);
}

159
static void ff_init_long_region(MPADecodeContext *s, GranuleDef *g, int ra1, int ra2){
160 161 162 163 164 165 166 167 168
    int l;
    g->region_size[0] =
        band_index_long[s->sample_rate_index][ra1 + 1] >> 1;
    /* should not overflow */
    l = FFMIN(ra1 + ra2 + 2, 22);
    g->region_size[1] =
        band_index_long[s->sample_rate_index][l] >> 1;
}

169
static void ff_compute_band_indexes(MPADecodeContext *s, GranuleDef *g){
170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192
    if (g->block_type == 2) {
        if (g->switch_point) {
            /* if switched mode, we handle the 36 first samples as
                long blocks.  For 8000Hz, we handle the 48 first
                exponents as long blocks (XXX: check this!) */
            if (s->sample_rate_index <= 2)
                g->long_end = 8;
            else if (s->sample_rate_index != 8)
                g->long_end = 6;
            else
                g->long_end = 4; /* 8000 Hz */

            g->short_start = 2 + (s->sample_rate_index != 8);
        } else {
            g->long_end = 0;
            g->short_start = 0;
        }
    } else {
        g->short_start = 13;
        g->long_end = 22;
    }
}

193 194 195 196 197
/* layer 1 unscaling */
/* n = number of bits of the mantissa minus 1 */
static inline int l1_unscale(int n, int mant, int scale_factor)
{
    int shift, mod;
198
    int64_t val;
199 200 201 202 203 204

    shift = scale_factor_modshift[scale_factor];
    mod = shift & 3;
    shift >>= 2;
    val = MUL64(mant + (-1 << n) + 1, scale_factor_mult[n-1][mod]);
    shift += n;
205 206
    /* NOTE: at this point, 1 <= shift >= 21 + 15 */
    return (int)((val + (1LL << (shift - 1))) >> shift);
207 208 209 210 211 212 213 214 215
}

static inline int l2_unscale_group(int steps, int mant, int scale_factor)
{
    int shift, mod, val;

    shift = scale_factor_modshift[scale_factor];
    mod = shift & 3;
    shift >>= 2;
216 217 218 219 220 221

    val = (mant - (steps >> 1)) * scale_factor_mult2[steps >> 2][mod];
    /* NOTE: at this point, 0 <= shift <= 21 */
    if (shift > 0)
        val = (val + (1 << (shift - 1))) >> shift;
    return val;
222 223 224 225 226 227 228 229
}

/* compute value^(4/3) * 2^(exponent/4). It normalized to FRAC_BITS */
static inline int l3_unscale(int value, int exponent)
{
    unsigned int m;
    int e;

230 231 232 233
    e = table_4_3_exp  [4*value + (exponent&3)];
    m = table_4_3_value[4*value + (exponent&3)];
    e -= (exponent >> 2);
    assert(e>=1);
234
    if (e > 31)
235
        return 0;
236
    m = (m + (1 << (e-1))) >> e;
237

238 239 240
    return m;
}

241 242 243 244 245 246
/* all integer n^(4/3) computation code */
#define DEV_ORDER 13

#define POW_FRAC_BITS 24
#define POW_FRAC_ONE    (1 << POW_FRAC_BITS)
#define POW_FIX(a)   ((int)((a) * POW_FRAC_ONE))
247
#define POW_MULL(a,b) (((int64_t)(a) * (int64_t)(b)) >> POW_FRAC_BITS)
248 249 250

static int dev_4_3_coefs[DEV_ORDER];

251
#if 0 /* unused */
252 253 254 255 256
static int pow_mult3[3] = {
    POW_FIX(1.0),
    POW_FIX(1.25992104989487316476),
    POW_FIX(1.58740105196819947474),
};
257
#endif
258

259
static av_cold void int_pow_init(void)
260 261 262 263 264 265 266 267 268 269
{
    int i, a;

    a = POW_FIX(1.0);
    for(i=0;i<DEV_ORDER;i++) {
        a = POW_MULL(a, POW_FIX(4.0 / 3.0) - i * POW_FIX(1.0)) / (i + 1);
        dev_4_3_coefs[i] = a;
    }
}

270
#if 0 /* unused, remove? */
271 272 273 274 275
/* return the mantissa and the binary exponent */
static int int_pow(int i, int *exp_ptr)
{
    int e, er, eq, j;
    int a, a1;
276

277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302
    /* renormalize */
    a = i;
    e = POW_FRAC_BITS;
    while (a < (1 << (POW_FRAC_BITS - 1))) {
        a = a << 1;
        e--;
    }
    a -= (1 << POW_FRAC_BITS);
    a1 = 0;
    for(j = DEV_ORDER - 1; j >= 0; j--)
        a1 = POW_MULL(a, dev_4_3_coefs[j] + a1);
    a = (1 << POW_FRAC_BITS) + a1;
    /* exponent compute (exact) */
    e = e * 4;
    er = e % 3;
    eq = e / 3;
    a = POW_MULL(a, pow_mult3[er]);
    while (a >= 2 * POW_FRAC_ONE) {
        a = a >> 1;
        eq++;
    }
    /* convert to float */
    while (a < POW_FRAC_ONE) {
        a = a << 1;
        eq--;
    }
303
    /* now POW_FRAC_ONE <= a < 2 * POW_FRAC_ONE */
304
#if POW_FRAC_BITS > FRAC_BITS
305 306 307 308 309 310 311
    a = (a + (1 << (POW_FRAC_BITS - FRAC_BITS - 1))) >> (POW_FRAC_BITS - FRAC_BITS);
    /* correct overflow */
    if (a >= 2 * (1 << FRAC_BITS)) {
        a = a >> 1;
        eq++;
    }
#endif
312 313 314
    *exp_ptr = eq;
    return a;
}
315
#endif
Fabrice Bellard's avatar
Fabrice Bellard committed
316

317
static av_cold int decode_init(AVCodecContext * avctx)
Fabrice Bellard's avatar
Fabrice Bellard committed
318 319
{
    MPADecodeContext *s = avctx->priv_data;
320
    static int init=0;
321
    int i, j, k;
Fabrice Bellard's avatar
Fabrice Bellard committed
322

323
    s->avctx = avctx;
324
    s->apply_window_mp3 = apply_window_mp3_c;
325
#if HAVE_MMX && CONFIG_FLOAT
326
    ff_mpegaudiodec_init_mmx(s);
327 328 329
#endif
#if CONFIG_FLOAT
    ff_dct_init(&s->dct, 5, DCT_II);
330
#endif
331 332
    if (HAVE_ALTIVEC && CONFIG_FLOAT) ff_mpegaudiodec_init_altivec(s);

333
    avctx->sample_fmt= OUT_FMT;
334
    s->error_recognition= avctx->error_recognition;
335

336
    if (!init && !avctx->parse_only) {
337 338
        int offset;

339 340 341 342
        /* scale factors table for layer 1/2 */
        for(i=0;i<64;i++) {
            int shift, mod;
            /* 1.0 (i = 3) is normalized to 2 ^ FRAC_BITS */
343
            shift = (i / 3);
344 345 346 347 348 349 350 351
            mod = i % 3;
            scale_factor_modshift[i] = mod | (shift << 2);
        }

        /* scale factor multiply for layer 1 */
        for(i=0;i<15;i++) {
            int n, norm;
            n = i + 2;
352
            norm = ((INT64_C(1) << n) * FRAC_ONE) / ((1 << n) - 1);
353 354 355
            scale_factor_mult[i][0] = MULLx(norm, FIXR(1.0          * 2.0), FRAC_BITS);
            scale_factor_mult[i][1] = MULLx(norm, FIXR(0.7937005259 * 2.0), FRAC_BITS);
            scale_factor_mult[i][2] = MULLx(norm, FIXR(0.6299605249 * 2.0), FRAC_BITS);
356
            dprintf(avctx, "%d: norm=%x s=%x %x %x\n",
357
                    i, norm,
358 359 360 361
                    scale_factor_mult[i][0],
                    scale_factor_mult[i][1],
                    scale_factor_mult[i][2]);
        }
362

363
        RENAME(ff_mpa_synth_init)(RENAME(ff_mpa_synth_window));
364

365
        /* huffman decode tables */
366
        offset = 0;
367 368
        for(i=1;i<16;i++) {
            const HuffTable *h = &mpa_huff_tables[i];
369
            int xsize, x, y;
Michael Niedermayer's avatar
Michael Niedermayer committed
370 371
            uint8_t  tmp_bits [512];
            uint16_t tmp_codes[512];
Michael Niedermayer's avatar
Michael Niedermayer committed
372 373 374

            memset(tmp_bits , 0, sizeof(tmp_bits ));
            memset(tmp_codes, 0, sizeof(tmp_codes));
375 376

            xsize = h->xsize;
377

378 379
            j = 0;
            for(x=0;x<xsize;x++) {
Michael Niedermayer's avatar
Michael Niedermayer committed
380
                for(y=0;y<xsize;y++){
Michael Niedermayer's avatar
Michael Niedermayer committed
381 382
                    tmp_bits [(x << 5) | y | ((x&&y)<<4)]= h->bits [j  ];
                    tmp_codes[(x << 5) | y | ((x&&y)<<4)]= h->codes[j++];
Michael Niedermayer's avatar
Michael Niedermayer committed
383
                }
384
            }
Michael Niedermayer's avatar
Michael Niedermayer committed
385 386

            /* XXX: fail test */
387 388
            huff_vlc[i].table = huff_vlc_tables+offset;
            huff_vlc[i].table_allocated = huff_vlc_tables_sizes[i];
Michael Niedermayer's avatar
Michael Niedermayer committed
389
            init_vlc(&huff_vlc[i], 7, 512,
390 391 392
                     tmp_bits, 1, 1, tmp_codes, 2, 2,
                     INIT_VLC_USE_NEW_STATIC);
            offset += huff_vlc_tables_sizes[i];
393
        }
394
        assert(offset == FF_ARRAY_ELEMS(huff_vlc_tables));
395 396

        offset = 0;
397
        for(i=0;i<2;i++) {
398 399
            huff_quad_vlc[i].table = huff_quad_vlc_tables+offset;
            huff_quad_vlc[i].table_allocated = huff_quad_vlc_tables_sizes[i];
400
            init_vlc(&huff_quad_vlc[i], i == 0 ? 7 : 4, 16,
401 402 403
                     mpa_quad_bits[i], 1, 1, mpa_quad_codes[i], 1, 1,
                     INIT_VLC_USE_NEW_STATIC);
            offset += huff_quad_vlc_tables_sizes[i];
404
        }
405
        assert(offset == FF_ARRAY_ELEMS(huff_quad_vlc_tables));
406 407 408 409 410 411 412 413 414 415

        for(i=0;i<9;i++) {
            k = 0;
            for(j=0;j<22;j++) {
                band_index_long[i][j] = k;
                k += band_size_long[i][j];
            }
            band_index_long[i][22] = k;
        }

416
        /* compute n ^ (4/3) and store it in mantissa/exp format */
417

418
        int_pow_init();
419
        mpegaudio_tableinit();
420

421 422 423 424 425 426 427 428 429 430 431 432 433 434
        for (i = 0; i < 4; i++)
            if (ff_mpa_quant_bits[i] < 0)
                for (j = 0; j < (1<<(-ff_mpa_quant_bits[i]+1)); j++) {
                    int val1, val2, val3, steps;
                    int val = j;
                    steps  = ff_mpa_quant_steps[i];
                    val1 = val % steps;
                    val /= steps;
                    val2 = val % steps;
                    val3 = val / steps;
                    division_tabs[i][j] = val1 + (val2 << 4) + (val3 << 8);
                }


435 436
        for(i=0;i<7;i++) {
            float f;
437
            INTFLOAT v;
438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458 459 460
            if (i != 6) {
                f = tan((double)i * M_PI / 12.0);
                v = FIXR(f / (1.0 + f));
            } else {
                v = FIXR(1.0);
            }
            is_table[0][i] = v;
            is_table[1][6 - i] = v;
        }
        /* invalid values */
        for(i=7;i<16;i++)
            is_table[0][i] = is_table[1][i] = 0.0;

        for(i=0;i<16;i++) {
            double f;
            int e, k;

            for(j=0;j<2;j++) {
                e = -(j + 1) * ((i + 1) >> 1);
                f = pow(2.0, e / 4.0);
                k = i & 1;
                is_table_lsf[j][k ^ 1][i] = FIXR(f);
                is_table_lsf[j][k][i] = FIXR(1.0);
461
                dprintf(avctx, "is_table_lsf %d %d: %x %x\n",
462 463 464 465 466 467 468 469 470
                        i, j, is_table_lsf[j][0][i], is_table_lsf[j][1][i]);
            }
        }

        for(i=0;i<8;i++) {
            float ci, cs, ca;
            ci = ci_table[i];
            cs = 1.0 / sqrt(1.0 + ci * ci);
            ca = cs * ci;
Michael Niedermayer's avatar
Michael Niedermayer committed
471 472 473
            csa_table[i][0] = FIXHR(cs/4);
            csa_table[i][1] = FIXHR(ca/4);
            csa_table[i][2] = FIXHR(ca/4) + FIXHR(cs/4);
474
            csa_table[i][3] = FIXHR(ca/4) - FIXHR(cs/4);
475 476 477
            csa_table_float[i][0] = cs;
            csa_table_float[i][1] = ca;
            csa_table_float[i][2] = ca + cs;
478
            csa_table_float[i][3] = ca - cs;
479 480 481 482
        }

        /* compute mdct windows */
        for(i=0;i<36;i++) {
483 484
            for(j=0; j<4; j++){
                double d;
485

Michael Niedermayer's avatar
Michael Niedermayer committed
486 487
                if(j==2 && i%3 != 1)
                    continue;
488

489 490 491 492 493 494 495 496 497 498 499
                d= sin(M_PI * (i + 0.5) / 36.0);
                if(j==1){
                    if     (i>=30) d= 0;
                    else if(i>=24) d= sin(M_PI * (i - 18 + 0.5) / 12.0);
                    else if(i>=18) d= 1;
                }else if(j==3){
                    if     (i<  6) d= 0;
                    else if(i< 12) d= sin(M_PI * (i -  6 + 0.5) / 12.0);
                    else if(i< 18) d= 1;
                }
                //merge last stage of imdct into the window coefficients
Michael Niedermayer's avatar
Michael Niedermayer committed
500 501 502 503 504 505
                d*= 0.5 / cos(M_PI*(2*i + 19)/72);

                if(j==2)
                    mdct_win[j][i/3] = FIXHR((d / (1<<5)));
                else
                    mdct_win[j][i  ] = FIXHR((d / (1<<5)));
506
            }
507 508 509 510 511 512 513 514 515 516 517
        }

        /* NOTE: we do frequency inversion adter the MDCT by changing
           the sign of the right window coefs */
        for(j=0;j<4;j++) {
            for(i=0;i<36;i+=2) {
                mdct_win[j + 4][i] = mdct_win[j][i];
                mdct_win[j + 4][i + 1] = -mdct_win[j][i + 1];
            }
        }

Fabrice Bellard's avatar
Fabrice Bellard committed
518 519 520
        init = 1;
    }

Roberto Togni's avatar
Roberto Togni committed
521 522
    if (avctx->codec_id == CODEC_ID_MP3ADU)
        s->adu_mode = 1;
Fabrice Bellard's avatar
Fabrice Bellard committed
523 524 525
    return 0;
}

526

527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543
#if CONFIG_FLOAT
static inline float round_sample(float *sum)
{
    float sum1=*sum;
    *sum = 0;
    return sum1;
}

/* signed 16x16 -> 32 multiply add accumulate */
#define MACS(rt, ra, rb) rt+=(ra)*(rb)

/* signed 16x16 -> 32 multiply */
#define MULS(ra, rb) ((ra)*(rb))

#define MLSS(rt, ra, rb) rt-=(ra)*(rb)

#elif FRAC_BITS <= 15
544

545
static inline int round_sample(int *sum)
546 547
{
    int sum1;
548 549
    sum1 = (*sum) >> OUT_SHIFT;
    *sum &= (1<<OUT_SHIFT)-1;
550
    return av_clip(sum1, OUT_MIN, OUT_MAX);
551 552
}

Luca Barbato's avatar
Luca Barbato committed
553 554
/* signed 16x16 -> 32 multiply add accumulate */
#define MACS(rt, ra, rb) MAC16(rt, ra, rb)
Siarhei Siamashka's avatar
Siarhei Siamashka committed
555

Luca Barbato's avatar
Luca Barbato committed
556 557
/* signed 16x16 -> 32 multiply */
#define MULS(ra, rb) MUL16(ra, rb)
558

559 560
#define MLSS(rt, ra, rb) MLS16(rt, ra, rb)

561 562
#else

563
static inline int round_sample(int64_t *sum)
564 565
{
    int sum1;
566 567
    sum1 = (int)((*sum) >> OUT_SHIFT);
    *sum &= (1<<OUT_SHIFT)-1;
568
    return av_clip(sum1, OUT_MIN, OUT_MAX);
569 570
}

571
#   define MULS(ra, rb) MUL64(ra, rb)
572 573
#   define MACS(rt, ra, rb) MAC64(rt, ra, rb)
#   define MLSS(rt, ra, rb) MLS64(rt, ra, rb)
574 575
#endif

576 577
#define SUM8(op, sum, w, p)               \
{                                         \
578 579 580 581 582 583 584 585
    op(sum, (w)[0 * 64], (p)[0 * 64]);    \
    op(sum, (w)[1 * 64], (p)[1 * 64]);    \
    op(sum, (w)[2 * 64], (p)[2 * 64]);    \
    op(sum, (w)[3 * 64], (p)[3 * 64]);    \
    op(sum, (w)[4 * 64], (p)[4 * 64]);    \
    op(sum, (w)[5 * 64], (p)[5 * 64]);    \
    op(sum, (w)[6 * 64], (p)[6 * 64]);    \
    op(sum, (w)[7 * 64], (p)[7 * 64]);    \
586 587 588 589
}

#define SUM8P2(sum1, op1, sum2, op2, w1, w2, p) \
{                                               \
590
    INTFLOAT tmp;\
591
    tmp = p[0 * 64];\
592 593
    op1(sum1, (w1)[0 * 64], tmp);\
    op2(sum2, (w2)[0 * 64], tmp);\
594
    tmp = p[1 * 64];\
595 596
    op1(sum1, (w1)[1 * 64], tmp);\
    op2(sum2, (w2)[1 * 64], tmp);\
597
    tmp = p[2 * 64];\
598 599
    op1(sum1, (w1)[2 * 64], tmp);\
    op2(sum2, (w2)[2 * 64], tmp);\
600
    tmp = p[3 * 64];\
601 602
    op1(sum1, (w1)[3 * 64], tmp);\
    op2(sum2, (w2)[3 * 64], tmp);\
603
    tmp = p[4 * 64];\
604 605
    op1(sum1, (w1)[4 * 64], tmp);\
    op2(sum2, (w2)[4 * 64], tmp);\
606
    tmp = p[5 * 64];\
607 608
    op1(sum1, (w1)[5 * 64], tmp);\
    op2(sum2, (w2)[5 * 64], tmp);\
609
    tmp = p[6 * 64];\
610 611
    op1(sum1, (w1)[6 * 64], tmp);\
    op2(sum2, (w2)[6 * 64], tmp);\
612
    tmp = p[7 * 64];\
613 614
    op1(sum1, (w1)[7 * 64], tmp);\
    op2(sum2, (w2)[7 * 64], tmp);\
615 616
}

617
void av_cold RENAME(ff_mpa_synth_init)(MPA_INT *window)
618
{
619
    int i, j;
620 621 622

    /* max = 18760, max sum over all 16 coefs : 44736 */
    for(i=0;i<257;i++) {
623
        INTFLOAT v;
624
        v = ff_mpa_enwindow[i];
625 626 627
#if CONFIG_FLOAT
        v *= 1.0 / (1LL<<(16 + FRAC_BITS));
#elif WFRAC_BITS < 16
628 629 630 631 632 633 634
        v = (v + (1 << (16 - WFRAC_BITS - 1))) >> (16 - WFRAC_BITS);
#endif
        window[i] = v;
        if ((i & 63) != 0)
            v = -v;
        if (i != 0)
            window[512 - i] = v;
635
    }
636 637 638 639 640 641 642 643 644

    // Needed for avoiding shuffles in ASM implementations
    for(i=0; i < 8; i++)
        for(j=0; j < 16; j++)
            window[512+16*i+j] = window[64*i+32-j];

    for(i=0; i < 8; i++)
        for(j=0; j < 16; j++)
            window[512+128+16*i+j] = window[64*i+48-j];
645
}
646

647 648
static void apply_window_mp3_c(MPA_INT *synth_buf, MPA_INT *window,
                               int *dither_state, OUT_INT *samples, int incr)
649
{
Alex Beregszaszi's avatar
Alex Beregszaszi committed
650
    register const MPA_INT *w, *w2, *p;
651
    int j;
652
    OUT_INT *samples2;
653 654 655
#if CONFIG_FLOAT
    float sum, sum2;
#elif FRAC_BITS <= 15
656
    int sum, sum2;
657
#else
658
    int64_t sum, sum2;
659
#endif
660

661
    /* copy to avoid wrap */
662
    memcpy(synth_buf + 512, synth_buf, 32 * sizeof(*synth_buf));
663

664
    samples2 = samples + 31 * incr;
665
    w = window;
666 667
    w2 = window + 31;

668
    sum = *dither_state;
669
    p = synth_buf + 16;
670
    SUM8(MACS, sum, w, p);
671
    p = synth_buf + 48;
672
    SUM8(MLSS, sum, w + 32, p);
673
    *samples = round_sample(&sum);
674
    samples += incr;
675 676
    w++;

677 678 679 680 681
    /* we calculate two samples at the same time to avoid one memory
       access per two sample */
    for(j=1;j<16;j++) {
        sum2 = 0;
        p = synth_buf + 16 + j;
682
        SUM8P2(sum, MACS, sum2, MLSS, w, w2, p);
683
        p = synth_buf + 48 - j;
684
        SUM8P2(sum, MLSS, sum2, MLSS, w + 32, w2 + 32, p);
685

686
        *samples = round_sample(&sum);
687
        samples += incr;
688 689
        sum += sum2;
        *samples2 = round_sample(&sum);
690
        samples2 -= incr;
691
        w++;
692
        w2--;
693
    }
694

695
    p = synth_buf + 32;
696
    SUM8(MLSS, sum, w + 32, p);
697
    *samples = round_sample(&sum);
698
    *dither_state= sum;
699 700 701 702 703 704
}


/* 32 sub band synthesis filter. Input: 32 sub band samples, Output:
   32 samples. */
/* XXX: optimize by avoiding ring buffer usage */
705
#if !CONFIG_FLOAT
706 707 708 709 710 711 712 713 714
void ff_mpa_synth_filter(MPA_INT *synth_buf_ptr, int *synth_buf_offset,
                         MPA_INT *window, int *dither_state,
                         OUT_INT *samples, int incr,
                         INTFLOAT sb_samples[SBLIMIT])
{
    register MPA_INT *synth_buf;
    int offset;
#if FRAC_BITS <= 15
    int32_t tmp[32];
715
    int j;
716 717 718 719 720
#endif

    offset = *synth_buf_offset;
    synth_buf = synth_buf_ptr + offset;

721
#if FRAC_BITS <= 15
722 723 724 725 726 727 728 729 730 731 732
    dct32(tmp, sb_samples);
    for(j=0;j<32;j++) {
        /* NOTE: can cause a loss in precision if very high amplitude
           sound */
        synth_buf[j] = av_clip_int16(tmp[j]);
    }
#else
    dct32(synth_buf, sb_samples);
#endif

    apply_window_mp3_c(synth_buf, window, dither_state, samples, incr);
733

734
    offset = (offset - 32) & 511;
735
    *synth_buf_offset = offset;
736
}
737
#endif
738

Michael Niedermayer's avatar
Michael Niedermayer committed
739 740 741
#define C3 FIXHR(0.86602540378443864676/2)

/* 0.5 / cos(pi*(2*i+1)/36) */
742
static const INTFLOAT icos36[9] = {
Michael Niedermayer's avatar
Michael Niedermayer committed
743 744 745 746 747 748 749 750 751 752
    FIXR(0.50190991877167369479),
    FIXR(0.51763809020504152469), //0
    FIXR(0.55168895948124587824),
    FIXR(0.61038729438072803416),
    FIXR(0.70710678118654752439), //1
    FIXR(0.87172339781054900991),
    FIXR(1.18310079157624925896),
    FIXR(1.93185165257813657349), //2
    FIXR(5.73685662283492756461),
};
753

754
/* 0.5 / cos(pi*(2*i+1)/36) */
755
static const INTFLOAT icos36h[9] = {
756 757 758 759 760 761 762 763 764 765 766
    FIXHR(0.50190991877167369479/2),
    FIXHR(0.51763809020504152469/2), //0
    FIXHR(0.55168895948124587824/2),
    FIXHR(0.61038729438072803416/2),
    FIXHR(0.70710678118654752439/2), //1
    FIXHR(0.87172339781054900991/2),
    FIXHR(1.18310079157624925896/4),
    FIXHR(1.93185165257813657349/4), //2
//    FIXHR(5.73685662283492756461),
};

767 768
/* 12 points IMDCT. We compute it "by hand" by factorizing obvious
   cases. */
769
static void imdct12(INTFLOAT *out, INTFLOAT *in)
770
{
771
    INTFLOAT in0, in1, in2, in3, in4, in5, t1, t2;
772 773 774 775 776 777 778

    in0= in[0*3];
    in1= in[1*3] + in[0*3];
    in2= in[2*3] + in[1*3];
    in3= in[3*3] + in[2*3];
    in4= in[4*3] + in[3*3];
    in5= in[5*3] + in[4*3];
Michael Niedermayer's avatar
Michael Niedermayer committed
779 780 781
    in5 += in3;
    in3 += in1;

782 783
    in2= MULH3(in2, C3, 2);
    in3= MULH3(in3, C3, 4);
784

Michael Niedermayer's avatar
Michael Niedermayer committed
785
    t1 = in0 - in4;
786
    t2 = MULH3(in1 - in5, icos36h[4], 2);
Michael Niedermayer's avatar
Michael Niedermayer committed
787

788
    out[ 7]=
Michael Niedermayer's avatar
Michael Niedermayer committed
789 790 791 792
    out[10]= t1 + t2;
    out[ 1]=
    out[ 4]= t1 - t2;

793
    in0 += SHR(in4, 1);
Michael Niedermayer's avatar
Michael Niedermayer committed
794
    in4 = in0 + in2;
795
    in5 += 2*in1;
796
    in1 = MULH3(in5 + in3, icos36h[1], 1);
797
    out[ 8]=
798
    out[ 9]= in4 + in1;
Michael Niedermayer's avatar
Michael Niedermayer committed
799
    out[ 2]=
800
    out[ 3]= in4 - in1;
801

Michael Niedermayer's avatar
Michael Niedermayer committed
802
    in0 -= in2;
803
    in5 = MULH3(in5 - in3, icos36h[7], 2);
Michael Niedermayer's avatar
Michael Niedermayer committed
804
    out[ 0]=
805
    out[ 5]= in0 - in5;
Michael Niedermayer's avatar
Michael Niedermayer committed
806
    out[ 6]=
807
    out[11]= in0 + in5;
808 809 810
}

/* cos(pi*i/18) */
811 812 813 814 815 816 817 818 819
#define C1 FIXHR(0.98480775301220805936/2)
#define C2 FIXHR(0.93969262078590838405/2)
#define C3 FIXHR(0.86602540378443864676/2)
#define C4 FIXHR(0.76604444311897803520/2)
#define C5 FIXHR(0.64278760968653932632/2)
#define C6 FIXHR(0.5/2)
#define C7 FIXHR(0.34202014332566873304/2)
#define C8 FIXHR(0.17364817766693034885/2)

820 821

/* using Lee like decomposition followed by hand coded 9 points DCT */
822
static void imdct36(INTFLOAT *out, INTFLOAT *buf, INTFLOAT *in, INTFLOAT *win)
823
{
824 825 826
    int i, j;
    INTFLOAT t0, t1, t2, t3, s0, s1, s2, s3;
    INTFLOAT tmp[18], *tmp1, *in1;
827 828 829 830 831 832 833 834 835

    for(i=17;i>=1;i--)
        in[i] += in[i-1];
    for(i=17;i>=3;i-=2)
        in[i] += in[i-2];

    for(j=0;j<2;j++) {
        tmp1 = tmp + j;
        in1 = in + j;
836

837
        t2 = in1[2*4] + in1[2*8] - in1[2*2];
838

839
        t3 = in1[2*0] + SHR(in1[2*6],1);
840
        t1 = in1[2*0] - in1[2*6];
841
        tmp1[ 6] = t1 - SHR(t2,1);
842 843
        tmp1[16] = t1 + t2;

844 845 846
        t0 = MULH3(in1[2*2] + in1[2*4] ,    C2, 2);
        t1 = MULH3(in1[2*4] - in1[2*8] , -2*C8, 1);
        t2 = MULH3(in1[2*2] + in1[2*8] ,   -C4, 2);
847

848 849 850
        tmp1[10] = t3 - t0 - t2;
        tmp1[ 2] = t3 + t0 + t1;
        tmp1[14] = t3 + t2 - t1;
851

852 853 854 855
        tmp1[ 4] = MULH3(in1[2*5] + in1[2*7] - in1[2*1], -C3, 2);
        t2 = MULH3(in1[2*1] + in1[2*5],    C1, 2);
        t3 = MULH3(in1[2*5] - in1[2*7], -2*C7, 1);
        t0 = MULH3(in1[2*3], C3, 2);
856

857
        t1 = MULH3(in1[2*1] + in1[2*7],   -C5, 2);
858 859 860 861

        tmp1[ 0] = t2 + t3 + t0;
        tmp1[12] = t2 + t1 - t0;
        tmp1[ 8] = t3 - t1 - t0;
862 863 864 865 866 867 868 869 870 871 872
    }

    i = 0;
    for(j=0;j<4;j++) {
        t0 = tmp[i];
        t1 = tmp[i + 2];
        s0 = t1 + t0;
        s2 = t1 - t0;

        t2 = tmp[i + 1];
        t3 = tmp[i + 3];
873 874
        s1 = MULH3(t3 + t2, icos36h[j], 2);
        s3 = MULLx(t3 - t2, icos36[8 - j], FRAC_BITS);
875

876 877
        t0 = s0 + s1;
        t1 = s0 - s1;
878 879 880 881
        out[(9 + j)*SBLIMIT] =  MULH3(t1, win[9 + j], 1) + buf[9 + j];
        out[(8 - j)*SBLIMIT] =  MULH3(t1, win[8 - j], 1) + buf[8 - j];
        buf[9 + j] = MULH3(t0, win[18 + 9 + j], 1);
        buf[8 - j] = MULH3(t0, win[18 + 8 - j], 1);
882

883 884
        t0 = s2 + s3;
        t1 = s2 - s3;
885 886 887 888
        out[(9 + 8 - j)*SBLIMIT] =  MULH3(t1, win[9 + 8 - j], 1) + buf[9 + 8 - j];
        out[(        j)*SBLIMIT] =  MULH3(t1, win[        j], 1) + buf[        j];
        buf[9 + 8 - j] = MULH3(t0, win[18 + 9 + 8 - j], 1);
        buf[      + j] = MULH3(t0, win[18         + j], 1);
889 890 891 892
        i += 4;
    }

    s0 = tmp[16];
893
    s1 = MULH3(tmp[17], icos36h[4], 2);
894 895
    t0 = s0 + s1;
    t1 = s0 - s1;
896 897 898 899
    out[(9 + 4)*SBLIMIT] =  MULH3(t1, win[9 + 4], 1) + buf[9 + 4];
    out[(8 - 4)*SBLIMIT] =  MULH3(t1, win[8 - 4], 1) + buf[8 - 4];
    buf[9 + 4] = MULH3(t0, win[18 + 9 + 4], 1);
    buf[8 - 4] = MULH3(t0, win[18 + 8 - 4], 1);
900 901 902 903
}

/* return the number of decoded frames */
static int mp_decode_layer1(MPADecodeContext *s)
Fabrice Bellard's avatar
Fabrice Bellard committed
904
{
905
    int bound, i, v, n, ch, j, mant;
906 907
    uint8_t allocation[MPA_MAX_CHANNELS][SBLIMIT];
    uint8_t scale_factors[MPA_MAX_CHANNELS][SBLIMIT];
908

909
    if (s->mode == MPA_JSTEREO)
910 911 912 913 914 915 916 917 918 919 920 921 922 923 924 925 926 927 928 929 930 931 932 933 934 935 936
        bound = (s->mode_ext + 1) * 4;
    else
        bound = SBLIMIT;

    /* allocation bits */
    for(i=0;i<bound;i++) {
        for(ch=0;ch<s->nb_channels;ch++) {
            allocation[ch][i] = get_bits(&s->gb, 4);
        }
    }
    for(i=bound;i<SBLIMIT;i++) {
        allocation[0][i] = get_bits(&s->gb, 4);
    }

    /* scale factors */
    for(i=0;i<bound;i++) {
        for(ch=0;ch<s->nb_channels;ch++) {
            if (allocation[ch][i])
                scale_factors[ch][i] = get_bits(&s->gb, 6);
        }
    }
    for(i=bound;i<SBLIMIT;i++) {
        if (allocation[0][i]) {
            scale_factors[0][i] = get_bits(&s->gb, 6);
            scale_factors[1][i] = get_bits(&s->gb, 6);
        }
    }
937

938 939 940 941 942 943 944 945 946 947 948 949 950 951 952 953 954 955 956 957 958 959 960 961 962 963 964 965 966 967 968 969 970 971 972 973 974 975 976 977
    /* compute samples */
    for(j=0;j<12;j++) {
        for(i=0;i<bound;i++) {
            for(ch=0;ch<s->nb_channels;ch++) {
                n = allocation[ch][i];
                if (n) {
                    mant = get_bits(&s->gb, n + 1);
                    v = l1_unscale(n, mant, scale_factors[ch][i]);
                } else {
                    v = 0;
                }
                s->sb_samples[ch][j][i] = v;
            }
        }
        for(i=bound;i<SBLIMIT;i++) {
            n = allocation[0][i];
            if (n) {
                mant = get_bits(&s->gb, n + 1);
                v = l1_unscale(n, mant, scale_factors[0][i]);
                s->sb_samples[0][j][i] = v;
                v = l1_unscale(n, mant, scale_factors[1][i]);
                s->sb_samples[1][j][i] = v;
            } else {
                s->sb_samples[0][j][i] = 0;
                s->sb_samples[1][j][i] = 0;
            }
        }
    }
    return 12;
}

static int mp_decode_layer2(MPADecodeContext *s)
{
    int sblimit; /* number of used subbands */
    const unsigned char *alloc_table;
    int table, bit_alloc_bits, i, j, ch, bound, v;
    unsigned char bit_alloc[MPA_MAX_CHANNELS][SBLIMIT];
    unsigned char scale_code[MPA_MAX_CHANNELS][SBLIMIT];
    unsigned char scale_factors[MPA_MAX_CHANNELS][SBLIMIT][3], *sf;
    int scale, qindex, bits, steps, k, l, m, b;
Fabrice Bellard's avatar
Fabrice Bellard committed
978

979
    /* select decoding table */
980
    table = ff_mpa_l2_select_table(s->bit_rate / 1000, s->nb_channels,
981
                            s->sample_rate, s->lsf);
982 983
    sblimit = ff_mpa_sblimit_table[table];
    alloc_table = ff_mpa_alloc_tables[table];
984

985
    if (s->mode == MPA_JSTEREO)
986 987 988 989
        bound = (s->mode_ext + 1) * 4;
    else
        bound = sblimit;

990
    dprintf(s->avctx, "bound=%d sblimit=%d\n", bound, sblimit);
991 992 993 994

    /* sanity check */
    if( bound > sblimit ) bound = sblimit;

995 996 997 998 999 1000 1001 1002 1003 1004 1005 1006 1007 1008 1009
    /* parse bit allocation */
    j = 0;
    for(i=0;i<bound;i++) {
        bit_alloc_bits = alloc_table[j];
        for(ch=0;ch<s->nb_channels;ch++) {
            bit_alloc[ch][i] = get_bits(&s->gb, bit_alloc_bits);
        }
        j += 1 << bit_alloc_bits;
    }
    for(i=bound;i<sblimit;i++) {
        bit_alloc_bits = alloc_table[j];
        v = get_bits(&s->gb, bit_alloc_bits);
        bit_alloc[0][i] = v;
        bit_alloc[1][i] = v;
        j += 1 << bit_alloc_bits;
Fabrice Bellard's avatar
Fabrice Bellard committed
1010
    }
1011 1012 1013 1014

    /* scale codes */
    for(i=0;i<sblimit;i++) {
        for(ch=0;ch<s->nb_channels;ch++) {
1015
            if (bit_alloc[ch][i])
1016 1017 1018
                scale_code[ch][i] = get_bits(&s->gb, 2);
        }
    }
1019

1020 1021 1022 1023 1024 1025 1026 1027 1028 1029 1030 1031 1032 1033 1034 1035 1036 1037 1038 1039 1040 1041 1042 1043 1044 1045 1046 1047 1048 1049 1050 1051 1052 1053 1054 1055 1056 1057 1058 1059 1060 1061 1062
    /* scale factors */
    for(i=0;i<sblimit;i++) {
        for(ch=0;ch<s->nb_channels;ch++) {
            if (bit_alloc[ch][i]) {
                sf = scale_factors[ch][i];
                switch(scale_code[ch][i]) {
                default:
                case 0:
                    sf[0] = get_bits(&s->gb, 6);
                    sf[1] = get_bits(&s->gb, 6);
                    sf[2] = get_bits(&s->gb, 6);
                    break;
                case 2:
                    sf[0] = get_bits(&s->gb, 6);
                    sf[1] = sf[0];
                    sf[2] = sf[0];
                    break;
                case 1:
                    sf[0] = get_bits(&s->gb, 6);
                    sf[2] = get_bits(&s->gb, 6);
                    sf[1] = sf[0];
                    break;
                case 3:
                    sf[0] = get_bits(&s->gb, 6);
                    sf[2] = get_bits(&s->gb, 6);
                    sf[1] = sf[2];
                    break;
                }
            }
        }
    }

    /* samples */
    for(k=0;k<3;k++) {
        for(l=0;l<12;l+=3) {
            j = 0;
            for(i=0;i<bound;i++) {
                bit_alloc_bits = alloc_table[j];
                for(ch=0;ch<s->nb_channels;ch++) {
                    b = bit_alloc[ch][i];
                    if (b) {
                        scale = scale_factors[ch][i][k];
                        qindex = alloc_table[j+b];
1063
                        bits = ff_mpa_quant_bits[qindex];
1064
                        if (bits < 0) {
1065
                            int v2;
1066 1067
                            /* 3 values at the same time */
                            v = get_bits(&s->gb, -bits);
1068 1069 1070
                            v2 = division_tabs[qindex][v];
                            steps  = ff_mpa_quant_steps[qindex];

1071
                            s->sb_samples[ch][k * 12 + l + 0][i] =
1072
                                l2_unscale_group(steps, v2        & 15, scale);
1073
                            s->sb_samples[ch][k * 12 + l + 1][i] =
1074
                                l2_unscale_group(steps, (v2 >> 4) & 15, scale);
1075
                            s->sb_samples[ch][k * 12 + l + 2][i] =
1076
                                l2_unscale_group(steps,  v2 >> 8      , scale);
1077 1078 1079 1080 1081 1082 1083 1084 1085 1086 1087 1088 1089 1090
                        } else {
                            for(m=0;m<3;m++) {
                                v = get_bits(&s->gb, bits);
                                v = l1_unscale(bits - 1, v, scale);
                                s->sb_samples[ch][k * 12 + l + m][i] = v;
                            }
                        }
                    } else {
                        s->sb_samples[ch][k * 12 + l + 0][i] = 0;
                        s->sb_samples[ch][k * 12 + l + 1][i] = 0;
                        s->sb_samples[ch][k * 12 + l + 2][i] = 0;
                    }
                }
                /* next subband in alloc table */
1091
                j += 1 << bit_alloc_bits;
1092 1093 1094 1095 1096 1097 1098 1099 1100 1101
            }
            /* XXX: find a way to avoid this duplication of code */
            for(i=bound;i<sblimit;i++) {
                bit_alloc_bits = alloc_table[j];
                b = bit_alloc[0][i];
                if (b) {
                    int mant, scale0, scale1;
                    scale0 = scale_factors[0][i][k];
                    scale1 = scale_factors[1][i][k];
                    qindex = alloc_table[j+b];
1102
                    bits = ff_mpa_quant_bits[qindex];
1103 1104 1105
                    if (bits < 0) {
                        /* 3 values at the same time */
                        v = get_bits(&s->gb, -bits);
1106
                        steps = ff_mpa_quant_steps[qindex];
1107 1108
                        mant = v % steps;
                        v = v / steps;
1109
                        s->sb_samples[0][k * 12 + l + 0][i] =
1110
                            l2_unscale_group(steps, mant, scale0);
1111
                        s->sb_samples[1][k * 12 + l + 0][i] =
1112 1113 1114
                            l2_unscale_group(steps, mant, scale1);
                        mant = v % steps;
                        v = v / steps;
1115
                        s->sb_samples[0][k * 12 + l + 1][i] =
1116
                            l2_unscale_group(steps, mant, scale0);
1117
                        s->sb_samples[1][k * 12 + l + 1][i] =
1118
                            l2_unscale_group(steps, mant, scale1);
1119
                        s->sb_samples[0][k * 12 + l + 2][i] =
1120
                            l2_unscale_group(steps, v, scale0);
1121
                        s->sb_samples[1][k * 12 + l + 2][i] =
1122 1123 1124 1125
                            l2_unscale_group(steps, v, scale1);
                    } else {
                        for(m=0;m<3;m++) {
                            mant = get_bits(&s->gb, bits);
1126
                            s->sb_samples[0][k * 12 + l + m][i] =
1127
                                l1_unscale(bits - 1, mant, scale0);
1128
                            s->sb_samples[1][k * 12 + l + m][i] =
1129 1130 1131 1132 1133 1134 1135 1136 1137 1138 1139 1140
                                l1_unscale(bits - 1, mant, scale1);
                        }
                    }
                } else {
                    s->sb_samples[0][k * 12 + l + 0][i] = 0;
                    s->sb_samples[0][k * 12 + l + 1][i] = 0;
                    s->sb_samples[0][k * 12 + l + 2][i] = 0;
                    s->sb_samples[1][k * 12 + l + 0][i] = 0;
                    s->sb_samples[1][k * 12 + l + 1][i] = 0;
                    s->sb_samples[1][k * 12 + l + 2][i] = 0;
                }
                /* next subband in alloc table */
1141
                j += 1 << bit_alloc_bits;
1142 1143 1144 1145 1146 1147 1148 1149 1150 1151 1152 1153
            }
            /* fill remaining samples to zero */
            for(i=sblimit;i<SBLIMIT;i++) {
                for(ch=0;ch<s->nb_channels;ch++) {
                    s->sb_samples[ch][k * 12 + l + 0][i] = 0;
                    s->sb_samples[ch][k * 12 + l + 1][i] = 0;
                    s->sb_samples[ch][k * 12 + l + 2][i] = 0;
                }
            }
        }
    }
    return 3 * 12;
Fabrice Bellard's avatar
Fabrice Bellard committed
1154 1155
}

1156 1157 1158 1159 1160 1161 1162 1163 1164 1165 1166 1167 1168 1169 1170 1171 1172 1173 1174 1175 1176
#define SPLIT(dst,sf,n)\
    if(n==3){\
        int m= (sf*171)>>9;\
        dst= sf - 3*m;\
        sf=m;\
    }else if(n==4){\
        dst= sf&3;\
        sf>>=2;\
    }else if(n==5){\
        int m= (sf*205)>>10;\
        dst= sf - 5*m;\
        sf=m;\
    }else if(n==6){\
        int m= (sf*171)>>10;\
        dst= sf - 6*m;\
        sf=m;\
    }else{\
        dst=0;\
    }

static av_always_inline void lsf_sf_expand(int *slen,
1177 1178
                                 int sf, int n1, int n2, int n3)
{
1179 1180 1181
    SPLIT(slen[3], sf, n3)
    SPLIT(slen[2], sf, n2)
    SPLIT(slen[1], sf, n1)
1182 1183 1184
    slen[0] = sf;
}

1185
static void exponents_from_scale_factors(MPADecodeContext *s,
1186
                                         GranuleDef *g,
1187
                                         int16_t *exponents)
1188
{
1189
    const uint8_t *bstab, *pretab;
1190
    int len, i, j, k, l, v0, shift, gain, gains[3];
1191
    int16_t *exp_ptr;
1192 1193 1194 1195 1196 1197 1198 1199

    exp_ptr = exponents;
    gain = g->global_gain - 210;
    shift = g->scalefac_scale + 1;

    bstab = band_size_long[s->sample_rate_index];
    pretab = mpa_pretab[g->preflag];
    for(i=0;i<g->long_end;i++) {
1200
        v0 = gain - ((g->scale_factors[i] + pretab[i]) << shift) + 400;
1201 1202 1203 1204 1205 1206 1207 1208 1209 1210 1211 1212 1213 1214
        len = bstab[i];
        for(j=len;j>0;j--)
            *exp_ptr++ = v0;
    }

    if (g->short_start < 13) {
        bstab = band_size_short[s->sample_rate_index];
        gains[0] = gain - (g->subblock_gain[0] << 3);
        gains[1] = gain - (g->subblock_gain[1] << 3);
        gains[2] = gain - (g->subblock_gain[2] << 3);
        k = g->long_end;
        for(i=g->short_start;i<13;i++) {
            len = bstab[i];
            for(l=0;l<3;l++) {
1215
                v0 = gains[l] - (g->scale_factors[k++] << shift) + 400;
1216 1217 1218 1219 1220 1221 1222 1223 1224 1225 1226 1227 1228 1229 1230 1231
                for(j=len;j>0;j--)
                *exp_ptr++ = v0;
            }
        }
    }
}

/* handle n = 0 too */
static inline int get_bitsz(GetBitContext *s, int n)
{
    if (n == 0)
        return 0;
    else
        return get_bits(s, n);
}

1232 1233 1234 1235 1236 1237 1238 1239 1240 1241 1242 1243 1244

static void switch_buffer(MPADecodeContext *s, int *pos, int *end_pos, int *end_pos2){
    if(s->in_gb.buffer && *pos >= s->gb.size_in_bits){
        s->gb= s->in_gb;
        s->in_gb.buffer=NULL;
        assert((get_bits_count(&s->gb) & 7) == 0);
        skip_bits_long(&s->gb, *pos - *end_pos);
        *end_pos2=
        *end_pos= *end_pos2 + get_bits_count(&s->gb) - *pos;
        *pos= get_bits_count(&s->gb);
    }
}

1245 1246 1247 1248 1249 1250 1251 1252 1253 1254 1255 1256 1257 1258 1259 1260
/* Following is a optimized code for
            INTFLOAT v = *src
            if(get_bits1(&s->gb))
                v = -v;
            *dst = v;
*/
#if CONFIG_FLOAT
#define READ_FLIP_SIGN(dst,src)\
            v = AV_RN32A(src) ^ (get_bits1(&s->gb)<<31);\
            AV_WN32A(dst, v);
#else
#define READ_FLIP_SIGN(dst,src)\
            v= -get_bits1(&s->gb);\
            *(dst) = (*(src) ^ v) - v;
#endif

1261
static int huffman_decode(MPADecodeContext *s, GranuleDef *g,
1262
                          int16_t *exponents, int end_pos2)
1263 1264
{
    int s_index;
1265
    int i;
1266
    int last_pos, bits_left;
1267
    VLC *vlc;
1268
    int end_pos= FFMIN(end_pos2, s->gb.size_in_bits);
1269 1270 1271 1272

    /* low frequencies (called big values) */
    s_index = 0;
    for(i=0;i<3;i++) {
1273
        int j, k, l, linbits;
1274 1275 1276 1277 1278 1279 1280 1281 1282
        j = g->region_size[i];
        if (j == 0)
            continue;
        /* select vlc table */
        k = g->table_select[i];
        l = mpa_huff_data[k][0];
        linbits = mpa_huff_data[k][1];
        vlc = &huff_vlc[l];

1283
        if(!l){
1284
            memset(&g->sb_hybrid[s_index], 0, sizeof(*g->sb_hybrid)*2*j);
1285 1286 1287 1288
            s_index += 2*j;
            continue;
        }

1289 1290
        /* read huffcode and compute each couple */
        for(;j>0;j--) {
1291
            int exponent, x, y;
1292
            int v;
1293 1294 1295 1296
            int pos= get_bits_count(&s->gb);

            if (pos >= end_pos){
//                av_log(NULL, AV_LOG_ERROR, "pos: %d %d %d %d\n", pos, end_pos, end_pos2, s_index);
1297
                switch_buffer(s, &pos, &end_pos, &end_pos2);
1298 1299 1300 1301
//                av_log(NULL, AV_LOG_ERROR, "new pos: %d %d\n", pos, end_pos);
                if(pos >= end_pos)
                    break;
            }
1302
            y = get_vlc2(&s->gb, vlc->table, 7, 3);
1303 1304 1305 1306 1307 1308 1309 1310

            if(!y){
                g->sb_hybrid[s_index  ] =
                g->sb_hybrid[s_index+1] = 0;
                s_index += 2;
                continue;
            }

1311
            exponent= exponents[s_index];
1312

1313
            dprintf(s->avctx, "region=%d n=%d x=%d y=%d exp=%d\n",
1314
                    i, g->region_size[i] - j, x, y, exponent);
Michael Niedermayer's avatar
Michael Niedermayer committed
1315 1316 1317
            if(y&16){
                x = y >> 5;
                y = y & 0x0f;
1318
                if (x < 15){
1319
                    READ_FLIP_SIGN(g->sb_hybrid+s_index, RENAME(expval_table)[ exponent ]+x)
1320 1321
                }else{
                    x += get_bitsz(&s->gb, linbits);
1322
                    v = l3_unscale(x, exponent);
1323 1324 1325
                    if (get_bits1(&s->gb))
                        v = -v;
                    g->sb_hybrid[s_index] = v;
1326 1327
                }
                if (y < 15){
1328
                    READ_FLIP_SIGN(g->sb_hybrid+s_index+1, RENAME(expval_table)[ exponent ]+y)
1329 1330
                }else{
                    y += get_bitsz(&s->gb, linbits);
1331
                    v = l3_unscale(y, exponent);
1332 1333 1334
                    if (get_bits1(&s->gb))
                        v = -v;
                    g->sb_hybrid[s_index+1] = v;
1335
                }
Michael Niedermayer's avatar
Michael Niedermayer committed
1336 1337 1338 1339 1340
            }else{
                x = y >> 5;
                y = y & 0x0f;
                x += y;
                if (x < 15){
1341
                    READ_FLIP_SIGN(g->sb_hybrid+s_index+!!y, RENAME(expval_table)[ exponent ]+x)
Michael Niedermayer's avatar
Michael Niedermayer committed
1342 1343 1344
                }else{
                    x += get_bitsz(&s->gb, linbits);
                    v = l3_unscale(x, exponent);
1345 1346 1347
                    if (get_bits1(&s->gb))
                        v = -v;
                    g->sb_hybrid[s_index+!!y] = v;
Michael Niedermayer's avatar
Michael Niedermayer committed
1348
                }
1349
                g->sb_hybrid[s_index+ !y] = 0;
1350
            }
Michael Niedermayer's avatar
Michael Niedermayer committed
1351
            s_index+=2;
1352 1353
        }
    }
1354

1355 1356
    /* high frequencies */
    vlc = &huff_quad_vlc[g->count1table_select];
1357
    last_pos=0;
1358
    while (s_index <= 572) {
1359
        int pos, code;
1360 1361
        pos = get_bits_count(&s->gb);
        if (pos >= end_pos) {
1362 1363 1364 1365 1366
            if (pos > end_pos2 && last_pos){
                /* some encoders generate an incorrect size for this
                   part. We must go back into the data */
                s_index -= 4;
                skip_bits_long(&s->gb, last_pos - pos);
1367
                av_log(s->avctx, AV_LOG_INFO, "overread, skip %d enddists: %d %d\n", last_pos - pos, end_pos-pos, end_pos2-pos);
1368
                if(s->error_recognition >= FF_ER_COMPLIANT)
1369
                    s_index=0;
1370 1371
                break;
            }
1372
//                av_log(NULL, AV_LOG_ERROR, "pos2: %d %d %d %d\n", pos, end_pos, end_pos2, s_index);
1373
            switch_buffer(s, &pos, &end_pos, &end_pos2);
1374 1375 1376
//                av_log(NULL, AV_LOG_ERROR, "new pos2: %d %d %d\n", pos, end_pos, s_index);
            if(pos >= end_pos)
                break;
1377
        }
1378
        last_pos= pos;
1379

1380
        code = get_vlc2(&s->gb, vlc->table, vlc->bits, 1);
1381
        dprintf(s->avctx, "t=%d code=%d\n", g->count1table_select, code);
1382 1383 1384 1385 1386
        g->sb_hybrid[s_index+0]=
        g->sb_hybrid[s_index+1]=
        g->sb_hybrid[s_index+2]=
        g->sb_hybrid[s_index+3]= 0;
        while(code){
1387
            static const int idxtab[16]={3,3,2,2,1,1,1,1,0,0,0,0,0,0,0,0};
1388
            int v;
1389 1390
            int pos= s_index+idxtab[code];
            code ^= 8>>idxtab[code];
1391
            READ_FLIP_SIGN(g->sb_hybrid+pos, RENAME(exp_table)+exponents[pos])
1392
        }
1393
        s_index+=4;
1394
    }
1395
    /* skip extension bits */
1396
    bits_left = end_pos2 - get_bits_count(&s->gb);
1397
//av_log(NULL, AV_LOG_ERROR, "left:%d buf:%p\n", bits_left, s->in_gb.buffer);
1398
    if (bits_left < 0 && s->error_recognition >= FF_ER_COMPLIANT) {
1399
        av_log(s->avctx, AV_LOG_ERROR, "bits_left=%d\n", bits_left);
1400
        s_index=0;
1401
    }else if(bits_left > 0 && s->error_recognition >= FF_ER_AGGRESSIVE){
1402
        av_log(s->avctx, AV_LOG_ERROR, "bits_left=%d\n", bits_left);
1403
        s_index=0;
1404
    }
1405
    memset(&g->sb_hybrid[s_index], 0, sizeof(*g->sb_hybrid)*(576 - s_index));
1406 1407
    skip_bits_long(&s->gb, bits_left);

1408
    i= get_bits_count(&s->gb);
1409
    switch_buffer(s, &i, &end_pos, &end_pos2);
1410

Fabrice Bellard's avatar
Fabrice Bellard committed
1411 1412 1413
    return 0;
}

1414 1415 1416 1417 1418
/* Reorder short blocks from bitstream order to interleaved order. It
   would be faster to do it in parsing, but the code would be far more
   complicated */
static void reorder_block(MPADecodeContext *s, GranuleDef *g)
{
1419
    int i, j, len;
1420 1421
    INTFLOAT *ptr, *dst, *ptr1;
    INTFLOAT tmp[576];
1422 1423 1424 1425 1426 1427 1428 1429 1430 1431 1432 1433 1434

    if (g->block_type != 2)
        return;

    if (g->switch_point) {
        if (s->sample_rate_index != 8) {
            ptr = g->sb_hybrid + 36;
        } else {
            ptr = g->sb_hybrid + 48;
        }
    } else {
        ptr = g->sb_hybrid;
    }
1435

1436 1437 1438
    for(i=g->short_start;i<13;i++) {
        len = band_size_short[s->sample_rate_index][i];
        ptr1 = ptr;
1439 1440 1441 1442 1443 1444
        dst = tmp;
        for(j=len;j>0;j--) {
            *dst++ = ptr[0*len];
            *dst++ = ptr[1*len];
            *dst++ = ptr[2*len];
            ptr++;
1445
        }
1446 1447
        ptr+=2*len;
        memcpy(ptr1, tmp, len * 3 * sizeof(*ptr1));
1448 1449 1450 1451 1452 1453 1454 1455 1456
    }
}

#define ISQRT2 FIXR(0.70710678118654752440)

static void compute_stereo(MPADecodeContext *s,
                           GranuleDef *g0, GranuleDef *g1)
{
    int i, j, k, l;
1457 1458
    int sf_max, sf, len, non_zero_found;
    INTFLOAT (*is_tab)[16], *tab0, *tab1, tmp0, tmp1, v1, v2;
1459 1460 1461 1462 1463 1464 1465 1466 1467 1468 1469
    int non_zero_found_short[3];

    /* intensity stereo */
    if (s->mode_ext & MODE_EXT_I_STEREO) {
        if (!s->lsf) {
            is_tab = is_table;
            sf_max = 7;
        } else {
            is_tab = is_table_lsf[g1->scalefac_compress & 1];
            sf_max = 16;
        }
1470

1471 1472 1473 1474 1475 1476 1477 1478 1479 1480 1481 1482 1483 1484 1485 1486 1487 1488 1489 1490 1491 1492 1493 1494 1495 1496 1497 1498 1499 1500 1501
        tab0 = g0->sb_hybrid + 576;
        tab1 = g1->sb_hybrid + 576;

        non_zero_found_short[0] = 0;
        non_zero_found_short[1] = 0;
        non_zero_found_short[2] = 0;
        k = (13 - g1->short_start) * 3 + g1->long_end - 3;
        for(i = 12;i >= g1->short_start;i--) {
            /* for last band, use previous scale factor */
            if (i != 11)
                k -= 3;
            len = band_size_short[s->sample_rate_index][i];
            for(l=2;l>=0;l--) {
                tab0 -= len;
                tab1 -= len;
                if (!non_zero_found_short[l]) {
                    /* test if non zero band. if so, stop doing i-stereo */
                    for(j=0;j<len;j++) {
                        if (tab1[j] != 0) {
                            non_zero_found_short[l] = 1;
                            goto found1;
                        }
                    }
                    sf = g1->scale_factors[k + l];
                    if (sf >= sf_max)
                        goto found1;

                    v1 = is_tab[0][sf];
                    v2 = is_tab[1][sf];
                    for(j=0;j<len;j++) {
                        tmp0 = tab0[j];
1502 1503
                        tab0[j] = MULLx(tmp0, v1, FRAC_BITS);
                        tab1[j] = MULLx(tmp0, v2, FRAC_BITS);
1504 1505 1506 1507 1508 1509 1510 1511 1512
                    }
                } else {
                found1:
                    if (s->mode_ext & MODE_EXT_MS_STEREO) {
                        /* lower part of the spectrum : do ms stereo
                           if enabled */
                        for(j=0;j<len;j++) {
                            tmp0 = tab0[j];
                            tmp1 = tab1[j];
1513 1514
                            tab0[j] = MULLx(tmp0 + tmp1, ISQRT2, FRAC_BITS);
                            tab1[j] = MULLx(tmp0 - tmp1, ISQRT2, FRAC_BITS);
1515 1516 1517 1518 1519 1520
                        }
                    }
                }
            }
        }

1521 1522
        non_zero_found = non_zero_found_short[0] |
            non_zero_found_short[1] |
1523 1524 1525 1526 1527 1528 1529 1530 1531 1532 1533 1534 1535 1536 1537 1538 1539 1540 1541 1542 1543 1544 1545
            non_zero_found_short[2];

        for(i = g1->long_end - 1;i >= 0;i--) {
            len = band_size_long[s->sample_rate_index][i];
            tab0 -= len;
            tab1 -= len;
            /* test if non zero band. if so, stop doing i-stereo */
            if (!non_zero_found) {
                for(j=0;j<len;j++) {
                    if (tab1[j] != 0) {
                        non_zero_found = 1;
                        goto found2;
                    }
                }
                /* for last band, use previous scale factor */
                k = (i == 21) ? 20 : i;
                sf = g1->scale_factors[k];
                if (sf >= sf_max)
                    goto found2;
                v1 = is_tab[0][sf];
                v2 = is_tab[1][sf];
                for(j=0;j<len;j++) {
                    tmp0 = tab0[j];
1546 1547
                    tab0[j] = MULLx(tmp0, v1, FRAC_BITS);
                    tab1[j] = MULLx(tmp0, v2, FRAC_BITS);
1548 1549 1550 1551 1552 1553 1554 1555 1556
                }
            } else {
            found2:
                if (s->mode_ext & MODE_EXT_MS_STEREO) {
                    /* lower part of the spectrum : do ms stereo
                       if enabled */
                    for(j=0;j<len;j++) {
                        tmp0 = tab0[j];
                        tmp1 = tab1[j];
1557 1558
                        tab0[j] = MULLx(tmp0 + tmp1, ISQRT2, FRAC_BITS);
                        tab1[j] = MULLx(tmp0 - tmp1, ISQRT2, FRAC_BITS);
1559 1560 1561 1562 1563 1564 1565 1566 1567 1568 1569 1570 1571 1572 1573 1574 1575 1576 1577
                    }
                }
            }
        }
    } else if (s->mode_ext & MODE_EXT_MS_STEREO) {
        /* ms stereo ONLY */
        /* NOTE: the 1/sqrt(2) normalization factor is included in the
           global gain */
        tab0 = g0->sb_hybrid;
        tab1 = g1->sb_hybrid;
        for(i=0;i<576;i++) {
            tmp0 = tab0[i];
            tmp1 = tab1[i];
            tab0[i] = tmp0 + tmp1;
            tab1[i] = tmp0 - tmp1;
        }
    }
}

1578
#if !CONFIG_FLOAT
1579
static void compute_antialias_integer(MPADecodeContext *s,
1580 1581
                              GranuleDef *g)
{
Michael Niedermayer's avatar
Michael Niedermayer committed
1582 1583
    int32_t *ptr, *csa;
    int n, i;
1584 1585 1586 1587 1588 1589 1590 1591 1592 1593

    /* we antialias only "long" bands */
    if (g->block_type == 2) {
        if (!g->switch_point)
            return;
        /* XXX: check this for 8000Hz case */
        n = 1;
    } else {
        n = SBLIMIT - 1;
    }
1594

1595 1596
    ptr = g->sb_hybrid + 18;
    for(i = n;i > 0;i--) {
Michael Niedermayer's avatar
Michael Niedermayer committed
1597 1598 1599
        int tmp0, tmp1, tmp2;
        csa = &csa_table[0][0];
#define INT_AA(j) \
1600 1601
            tmp0 = ptr[-1-j];\
            tmp1 = ptr[   j];\
Michael Niedermayer's avatar
Michael Niedermayer committed
1602
            tmp2= MULH(tmp0 + tmp1, csa[0+4*j]);\
1603 1604
            ptr[-1-j] = 4*(tmp2 - MULH(tmp1, csa[2+4*j]));\
            ptr[   j] = 4*(tmp2 + MULH(tmp0, csa[3+4*j]));
Michael Niedermayer's avatar
Michael Niedermayer committed
1605 1606 1607 1608 1609 1610 1611 1612 1613

        INT_AA(0)
        INT_AA(1)
        INT_AA(2)
        INT_AA(3)
        INT_AA(4)
        INT_AA(5)
        INT_AA(6)
        INT_AA(7)
1614 1615

        ptr += 18;
1616 1617
    }
}
1618
#endif
1619 1620

static void compute_imdct(MPADecodeContext *s,
1621
                          GranuleDef *g,
1622 1623
                          INTFLOAT *sb_samples,
                          INTFLOAT *mdct_buf)
1624
{
1625 1626 1627
    INTFLOAT *win, *win1, *out_ptr, *ptr, *buf, *ptr1;
    INTFLOAT out2[12];
    int i, j, mdct_long_end, sblimit;
1628 1629 1630 1631 1632

    /* find last non zero block */
    ptr = g->sb_hybrid + 576;
    ptr1 = g->sb_hybrid + 2 * 18;
    while (ptr >= ptr1) {
1633
        int32_t *p;
1634
        ptr -= 6;
1635 1636
        p= (int32_t*)ptr;
        if(p[0] | p[1] | p[2] | p[3] | p[4] | p[5])
1637 1638 1639 1640 1641 1642 1643 1644 1645 1646 1647 1648 1649 1650 1651 1652 1653 1654 1655 1656 1657 1658 1659 1660 1661 1662
            break;
    }
    sblimit = ((ptr - g->sb_hybrid) / 18) + 1;

    if (g->block_type == 2) {
        /* XXX: check for 8000 Hz */
        if (g->switch_point)
            mdct_long_end = 2;
        else
            mdct_long_end = 0;
    } else {
        mdct_long_end = sblimit;
    }

    buf = mdct_buf;
    ptr = g->sb_hybrid;
    for(j=0;j<mdct_long_end;j++) {
        /* apply window & overlap with previous buffer */
        out_ptr = sb_samples + j;
        /* select window */
        if (g->switch_point && j < 2)
            win1 = mdct_win[0];
        else
            win1 = mdct_win[g->block_type];
        /* select frequency inversion */
        win = win1 + ((4 * 36) & -(j & 1));
1663 1664
        imdct36(out_ptr, buf, ptr, win);
        out_ptr += 18*SBLIMIT;
1665 1666 1667 1668 1669 1670 1671
        ptr += 18;
        buf += 18;
    }
    for(j=mdct_long_end;j<sblimit;j++) {
        /* select frequency inversion */
        win = mdct_win[2] + ((4 * 36) & -(j & 1));
        out_ptr = sb_samples + j;
1672

Michael Niedermayer's avatar
Michael Niedermayer committed
1673 1674 1675 1676 1677 1678
        for(i=0; i<6; i++){
            *out_ptr = buf[i];
            out_ptr += SBLIMIT;
        }
        imdct12(out2, ptr + 0);
        for(i=0;i<6;i++) {
1679 1680
            *out_ptr     = MULH3(out2[i    ], win[i    ], 1) + buf[i + 6*1];
            buf[i + 6*2] = MULH3(out2[i + 6], win[i + 6], 1);
1681 1682
            out_ptr += SBLIMIT;
        }
Michael Niedermayer's avatar
Michael Niedermayer committed
1683 1684
        imdct12(out2, ptr + 1);
        for(i=0;i<6;i++) {
1685 1686
            *out_ptr     = MULH3(out2[i    ], win[i    ], 1) + buf[i + 6*2];
            buf[i + 6*0] = MULH3(out2[i + 6], win[i + 6], 1);
Michael Niedermayer's avatar
Michael Niedermayer committed
1687 1688 1689 1690
            out_ptr += SBLIMIT;
        }
        imdct12(out2, ptr + 2);
        for(i=0;i<6;i++) {
1691 1692
            buf[i + 6*0] = MULH3(out2[i    ], win[i    ], 1) + buf[i + 6*0];
            buf[i + 6*1] = MULH3(out2[i + 6], win[i + 6], 1);
Michael Niedermayer's avatar
Michael Niedermayer committed
1693 1694
            buf[i + 6*2] = 0;
        }
1695 1696 1697 1698 1699 1700 1701 1702 1703 1704 1705 1706 1707 1708 1709 1710 1711 1712 1713 1714
        ptr += 18;
        buf += 18;
    }
    /* zero bands */
    for(j=sblimit;j<SBLIMIT;j++) {
        /* overlap */
        out_ptr = sb_samples + j;
        for(i=0;i<18;i++) {
            *out_ptr = buf[i];
            buf[i] = 0;
            out_ptr += SBLIMIT;
        }
        buf += 18;
    }
}

/* main layer3 decoding function */
static int mp_decode_layer3(MPADecodeContext *s)
{
    int nb_granules, main_data_begin, private_bits;
1715
    int gr, ch, blocksplit_flag, i, j, k, n, bits_pos;
1716
    GranuleDef *g;
1717
    int16_t exponents[576]; //FIXME try INTFLOAT
1718 1719 1720 1721

    /* read side info */
    if (s->lsf) {
        main_data_begin = get_bits(&s->gb, 8);
Michael Niedermayer's avatar
Michael Niedermayer committed
1722
        private_bits = get_bits(&s->gb, s->nb_channels);
1723 1724 1725 1726 1727 1728 1729 1730 1731
        nb_granules = 1;
    } else {
        main_data_begin = get_bits(&s->gb, 9);
        if (s->nb_channels == 2)
            private_bits = get_bits(&s->gb, 3);
        else
            private_bits = get_bits(&s->gb, 5);
        nb_granules = 2;
        for(ch=0;ch<s->nb_channels;ch++) {
1732 1733
            s->granules[ch][0].scfsi = 0;/* all scale factors are transmitted */
            s->granules[ch][1].scfsi = get_bits(&s->gb, 4);
1734 1735
        }
    }
1736

1737 1738
    for(gr=0;gr<nb_granules;gr++) {
        for(ch=0;ch<s->nb_channels;ch++) {
1739
            dprintf(s->avctx, "gr=%d ch=%d: side_info\n", gr, ch);
1740
            g = &s->granules[ch][gr];
1741 1742
            g->part2_3_length = get_bits(&s->gb, 12);
            g->big_values = get_bits(&s->gb, 9);
1743
            if(g->big_values > 288){
1744
                av_log(s->avctx, AV_LOG_ERROR, "big_values too big\n");
1745 1746 1747
                return -1;
            }

1748 1749 1750
            g->global_gain = get_bits(&s->gb, 8);
            /* if MS stereo only is selected, we precompute the
               1/sqrt(2) renormalization factor */
1751
            if ((s->mode_ext & (MODE_EXT_MS_STEREO | MODE_EXT_I_STEREO)) ==
1752 1753 1754 1755 1756 1757
                MODE_EXT_MS_STEREO)
                g->global_gain -= 2;
            if (s->lsf)
                g->scalefac_compress = get_bits(&s->gb, 9);
            else
                g->scalefac_compress = get_bits(&s->gb, 4);
1758
            blocksplit_flag = get_bits1(&s->gb);
1759 1760
            if (blocksplit_flag) {
                g->block_type = get_bits(&s->gb, 2);
1761
                if (g->block_type == 0){
1762
                    av_log(s->avctx, AV_LOG_ERROR, "invalid block type\n");
1763
                    return -1;
1764
                }
1765
                g->switch_point = get_bits1(&s->gb);
1766 1767
                for(i=0;i<2;i++)
                    g->table_select[i] = get_bits(&s->gb, 5);
1768
                for(i=0;i<3;i++)
1769
                    g->subblock_gain[i] = get_bits(&s->gb, 3);
1770
                ff_init_short_region(s, g);
1771
            } else {
1772
                int region_address1, region_address2;
1773 1774 1775 1776 1777 1778 1779
                g->block_type = 0;
                g->switch_point = 0;
                for(i=0;i<3;i++)
                    g->table_select[i] = get_bits(&s->gb, 5);
                /* compute huffman coded region sizes */
                region_address1 = get_bits(&s->gb, 4);
                region_address2 = get_bits(&s->gb, 3);
1780
                dprintf(s->avctx, "region1=%d region2=%d\n",
1781
                        region_address1, region_address2);
1782
                ff_init_long_region(s, g, region_address1, region_address2);
1783
            }
1784 1785
            ff_region_offset2size(g);
            ff_compute_band_indexes(s, g);
1786

1787 1788
            g->preflag = 0;
            if (!s->lsf)
1789 1790 1791
                g->preflag = get_bits1(&s->gb);
            g->scalefac_scale = get_bits1(&s->gb);
            g->count1table_select = get_bits1(&s->gb);
1792
            dprintf(s->avctx, "block_type=%d switch_point=%d\n",
1793 1794 1795 1796
                    g->block_type, g->switch_point);
        }
    }

Roberto Togni's avatar
Roberto Togni committed
1797
  if (!s->adu_mode) {
1798
    const uint8_t *ptr = s->gb.buffer + (get_bits_count(&s->gb)>>3);
1799
    assert((get_bits_count(&s->gb) & 7) == 0);
1800
    /* now we get bits from the main_data_begin offset */
1801
    dprintf(s->avctx, "seekback: %d\n", main_data_begin);
1802 1803 1804 1805
//av_log(NULL, AV_LOG_ERROR, "backstep:%d, lastbuf:%d\n", main_data_begin, s->last_buf_size);

    memcpy(s->last_buf + s->last_buf_size, ptr, EXTRABYTES);
    s->in_gb= s->gb;
1806 1807
        init_get_bits(&s->gb, s->last_buf, s->last_buf_size*8);
        skip_bits_long(&s->gb, 8*(s->last_buf_size - main_data_begin));
Roberto Togni's avatar
Roberto Togni committed
1808
  }
1809 1810 1811

    for(gr=0;gr<nb_granules;gr++) {
        for(ch=0;ch<s->nb_channels;ch++) {
1812
            g = &s->granules[ch][gr];
1813
            if(get_bits_count(&s->gb)<0){
1814
                av_log(s->avctx, AV_LOG_DEBUG, "mdb:%d, lastbuf:%d skipping granule %d\n",
1815 1816 1817 1818 1819 1820 1821 1822 1823 1824
                                            main_data_begin, s->last_buf_size, gr);
                skip_bits_long(&s->gb, g->part2_3_length);
                memset(g->sb_hybrid, 0, sizeof(g->sb_hybrid));
                if(get_bits_count(&s->gb) >= s->gb.size_in_bits && s->in_gb.buffer){
                    skip_bits_long(&s->in_gb, get_bits_count(&s->gb) - s->gb.size_in_bits);
                    s->gb= s->in_gb;
                    s->in_gb.buffer=NULL;
                }
                continue;
            }
1825

1826
            bits_pos = get_bits_count(&s->gb);
1827

1828
            if (!s->lsf) {
1829
                uint8_t *sc;
1830 1831 1832 1833 1834
                int slen, slen1, slen2;

                /* MPEG1 scale factors */
                slen1 = slen_table[0][g->scalefac_compress];
                slen2 = slen_table[1][g->scalefac_compress];
1835
                dprintf(s->avctx, "slen1=%d slen2=%d\n", slen1, slen2);
1836 1837 1838
                if (g->block_type == 2) {
                    n = g->switch_point ? 17 : 18;
                    j = 0;
1839 1840 1841 1842 1843 1844 1845 1846 1847 1848 1849 1850 1851 1852 1853 1854
                    if(slen1){
                        for(i=0;i<n;i++)
                            g->scale_factors[j++] = get_bits(&s->gb, slen1);
                    }else{
                        for(i=0;i<n;i++)
                            g->scale_factors[j++] = 0;
                    }
                    if(slen2){
                        for(i=0;i<18;i++)
                            g->scale_factors[j++] = get_bits(&s->gb, slen2);
                        for(i=0;i<3;i++)
                            g->scale_factors[j++] = 0;
                    }else{
                        for(i=0;i<21;i++)
                            g->scale_factors[j++] = 0;
                    }
1855
                } else {
1856
                    sc = s->granules[ch][0].scale_factors;
1857 1858 1859 1860 1861
                    j = 0;
                    for(k=0;k<4;k++) {
                        n = (k == 0 ? 6 : 5);
                        if ((g->scfsi & (0x8 >> k)) == 0) {
                            slen = (k < 2) ? slen1 : slen2;
1862 1863 1864 1865 1866 1867 1868
                            if(slen){
                                for(i=0;i<n;i++)
                                    g->scale_factors[j++] = get_bits(&s->gb, slen);
                            }else{
                                for(i=0;i<n;i++)
                                    g->scale_factors[j++] = 0;
                            }
1869 1870 1871 1872 1873 1874 1875 1876 1877 1878 1879 1880 1881 1882 1883 1884 1885 1886 1887 1888 1889 1890 1891 1892 1893 1894 1895 1896 1897 1898 1899 1900 1901 1902 1903 1904 1905 1906 1907 1908 1909 1910 1911 1912 1913 1914 1915 1916 1917 1918 1919 1920
                        } else {
                            /* simply copy from last granule */
                            for(i=0;i<n;i++) {
                                g->scale_factors[j] = sc[j];
                                j++;
                            }
                        }
                    }
                    g->scale_factors[j++] = 0;
                }
            } else {
                int tindex, tindex2, slen[4], sl, sf;

                /* LSF scale factors */
                if (g->block_type == 2) {
                    tindex = g->switch_point ? 2 : 1;
                } else {
                    tindex = 0;
                }
                sf = g->scalefac_compress;
                if ((s->mode_ext & MODE_EXT_I_STEREO) && ch == 1) {
                    /* intensity stereo case */
                    sf >>= 1;
                    if (sf < 180) {
                        lsf_sf_expand(slen, sf, 6, 6, 0);
                        tindex2 = 3;
                    } else if (sf < 244) {
                        lsf_sf_expand(slen, sf - 180, 4, 4, 0);
                        tindex2 = 4;
                    } else {
                        lsf_sf_expand(slen, sf - 244, 3, 0, 0);
                        tindex2 = 5;
                    }
                } else {
                    /* normal case */
                    if (sf < 400) {
                        lsf_sf_expand(slen, sf, 5, 4, 4);
                        tindex2 = 0;
                    } else if (sf < 500) {
                        lsf_sf_expand(slen, sf - 400, 5, 4, 0);
                        tindex2 = 1;
                    } else {
                        lsf_sf_expand(slen, sf - 500, 3, 0, 0);
                        tindex2 = 2;
                        g->preflag = 1;
                    }
                }

                j = 0;
                for(k=0;k<4;k++) {
                    n = lsf_nsf_table[tindex2][tindex][k];
                    sl = slen[k];
Michael Niedermayer's avatar
Michael Niedermayer committed
1921
                    if(sl){
1922 1923 1924 1925 1926 1927
                        for(i=0;i<n;i++)
                            g->scale_factors[j++] = get_bits(&s->gb, sl);
                    }else{
                        for(i=0;i<n;i++)
                            g->scale_factors[j++] = 0;
                    }
1928 1929 1930 1931 1932 1933 1934 1935 1936
                }
                /* XXX: should compute exact size */
                for(;j<40;j++)
                    g->scale_factors[j] = 0;
            }

            exponents_from_scale_factors(s, g, exponents);

            /* read Huffman coded residue */
1937
            huffman_decode(s, g, exponents, bits_pos + g->part2_3_length);
1938 1939 1940
        } /* ch */

        if (s->nb_channels == 2)
1941
            compute_stereo(s, &s->granules[0][gr], &s->granules[1][gr]);
1942 1943

        for(ch=0;ch<s->nb_channels;ch++) {
1944
            g = &s->granules[ch][gr];
1945 1946

            reorder_block(s, g);
1947
            compute_antialias(s, g);
1948
            compute_imdct(s, g, &s->sb_samples[ch][18 * gr][0], s->mdct_buf[ch]);
1949 1950
        }
    } /* gr */
1951 1952
    if(get_bits_count(&s->gb)<0)
        skip_bits_long(&s->gb, -get_bits_count(&s->gb));
1953 1954 1955
    return nb_granules * 18;
}

1956
static int mp_decode_frame(MPADecodeContext *s,
1957
                           OUT_INT *samples, const uint8_t *buf, int buf_size)
1958 1959
{
    int i, nb_frames, ch;
1960
    OUT_INT *samples_ptr;
1961

1962
    init_get_bits(&s->gb, buf + HEADER_SIZE, (buf_size - HEADER_SIZE)*8);
1963

1964 1965
    /* skip error protection field */
    if (s->error_protection)
1966
        skip_bits(&s->gb, 16);
1967

1968
    dprintf(s->avctx, "frame %d:\n", s->frame_count);
1969 1970
    switch(s->layer) {
    case 1:
1971
        s->avctx->frame_size = 384;
1972 1973 1974
        nb_frames = mp_decode_layer1(s);
        break;
    case 2:
1975
        s->avctx->frame_size = 1152;
1976 1977 1978
        nb_frames = mp_decode_layer2(s);
        break;
    case 3:
1979
        s->avctx->frame_size = s->lsf ? 576 : 1152;
1980 1981
    default:
        nb_frames = mp_decode_layer3(s);
1982

1983 1984 1985
        s->last_buf_size=0;
        if(s->in_gb.buffer){
            align_get_bits(&s->gb);
1986
            i= get_bits_left(&s->gb)>>3;
1987
            if(i >= 0 && i <= BACKSTEP_SIZE){
1988 1989
                memmove(s->last_buf, s->gb.buffer + (get_bits_count(&s->gb)>>3), i);
                s->last_buf_size=i;
1990
            }else
1991
                av_log(s->avctx, AV_LOG_ERROR, "invalid old backstep %d\n", i);
1992
            s->gb= s->in_gb;
1993
            s->in_gb.buffer= NULL;
1994 1995
        }

1996 1997
        align_get_bits(&s->gb);
        assert((get_bits_count(&s->gb) & 7) == 0);
1998
        i= get_bits_left(&s->gb)>>3;
1999

2000
        if(i<0 || i > BACKSTEP_SIZE || nb_frames<0){
2001 2002
            if(i<0)
                av_log(s->avctx, AV_LOG_ERROR, "invalid new backstep %d\n", i);
2003 2004
            i= FFMIN(BACKSTEP_SIZE, buf_size - HEADER_SIZE);
        }
2005
        assert(i <= buf_size - HEADER_SIZE && i>= 0);
2006
        memcpy(s->last_buf + s->last_buf_size, s->gb.buffer + buf_size - HEADER_SIZE - i, i);
2007
        s->last_buf_size += i;
2008

2009 2010
        break;
    }
2011

2012 2013 2014 2015
    /* apply the synthesis filter */
    for(ch=0;ch<s->nb_channels;ch++) {
        samples_ptr = samples + ch;
        for(i=0;i<nb_frames;i++) {
2016 2017 2018 2019 2020
            RENAME(ff_mpa_synth_filter)(
#if CONFIG_FLOAT
                         s,
#endif
                         s->synth_buf[ch], &(s->synth_buf_offset[ch]),
2021
                         RENAME(ff_mpa_synth_window), &s->dither_state,
2022
                         samples_ptr, s->nb_channels,
2023 2024 2025 2026
                         s->sb_samples[ch][i]);
            samples_ptr += 32 * s->nb_channels;
        }
    }
2027

2028
    return nb_frames * 32 * sizeof(OUT_INT) * s->nb_channels;
2029 2030
}

Fabrice Bellard's avatar
Fabrice Bellard committed
2031
static int decode_frame(AVCodecContext * avctx,
2032
                        void *data, int *data_size,
2033
                        AVPacket *avpkt)
Fabrice Bellard's avatar
Fabrice Bellard committed
2034
{
2035 2036
    const uint8_t *buf = avpkt->data;
    int buf_size = avpkt->size;
Fabrice Bellard's avatar
Fabrice Bellard committed
2037
    MPADecodeContext *s = avctx->priv_data;
2038
    uint32_t header;
2039
    int out_size;
2040
    OUT_INT *out_samples = data;
Fabrice Bellard's avatar
Fabrice Bellard committed
2041

2042 2043 2044
    if(buf_size < HEADER_SIZE)
        return -1;

2045
    header = AV_RB32(buf);
2046
    if(ff_mpa_check_header(header) < 0){
2047 2048 2049 2050 2051 2052 2053

        if (buf_size == ID3v1_TAG_SIZE
            && buf[0] == 'T' && buf[1] == 'A' && buf[2] == 'G') {
            *data_size = 0;
            return ID3v1_TAG_SIZE;
        }

2054 2055
        av_log(avctx, AV_LOG_ERROR, "Header missing\n");
        return -1;
2056 2057
    }

2058
    if (ff_mpegaudio_decode_header((MPADecodeHeader *)s, header) == 1) {
2059 2060 2061 2062 2063 2064
        /* free format: prepare to compute frame size */
        s->frame_size = -1;
        return -1;
    }
    /* update codec info */
    avctx->channels = s->nb_channels;
2065 2066
    if (!avctx->bit_rate)
        avctx->bit_rate = s->bit_rate;
2067 2068
    avctx->sub_id = s->layer;

2069 2070
    if(*data_size < 1152*avctx->channels*sizeof(OUT_INT))
        return -1;
2071
    *data_size = 0;
2072

2073
    if(s->frame_size<=0 || s->frame_size > buf_size){
2074 2075
        av_log(avctx, AV_LOG_ERROR, "incomplete frame\n");
        return -1;
2076 2077
    }else if(s->frame_size < buf_size){
        av_log(avctx, AV_LOG_ERROR, "incorrect frame size\n");
2078
        buf_size= s->frame_size;
Fabrice Bellard's avatar
Fabrice Bellard committed
2079
    }
2080 2081

    out_size = mp_decode_frame(s, out_samples, buf, buf_size);
2082
    if(out_size>=0){
2083
        *data_size = out_size;
2084 2085 2086
        avctx->sample_rate = s->sample_rate;
        //FIXME maybe move the other codec info stuff from above here too
    }else
Diego Biurrun's avatar
Diego Biurrun committed
2087
        av_log(avctx, AV_LOG_DEBUG, "Error while decoding MPEG audio frame.\n"); //FIXME return -1 / but also return the number of bytes consumed
2088
    s->frame_size = 0;
2089
    return buf_size;
Fabrice Bellard's avatar
Fabrice Bellard committed
2090 2091
}

2092 2093
static void flush(AVCodecContext *avctx){
    MPADecodeContext *s = avctx->priv_data;
2094
    memset(s->synth_buf, 0, sizeof(s->synth_buf));
2095 2096 2097
    s->last_buf_size= 0;
}

2098
#if CONFIG_MP3ADU_DECODER || CONFIG_MP3ADUFLOAT_DECODER
Roberto Togni's avatar
Roberto Togni committed
2099
static int decode_frame_adu(AVCodecContext * avctx,
2100
                        void *data, int *data_size,
2101
                        AVPacket *avpkt)
Roberto Togni's avatar
Roberto Togni committed
2102
{
2103 2104
    const uint8_t *buf = avpkt->data;
    int buf_size = avpkt->size;
Roberto Togni's avatar
Roberto Togni committed
2105 2106 2107
    MPADecodeContext *s = avctx->priv_data;
    uint32_t header;
    int len, out_size;
2108
    OUT_INT *out_samples = data;
Roberto Togni's avatar
Roberto Togni committed
2109 2110 2111 2112 2113 2114 2115 2116 2117 2118 2119 2120 2121 2122

    len = buf_size;

    // Discard too short frames
    if (buf_size < HEADER_SIZE) {
        *data_size = 0;
        return buf_size;
    }


    if (len > MPA_MAX_CODED_FRAME_SIZE)
        len = MPA_MAX_CODED_FRAME_SIZE;

    // Get header and restore sync word
2123
    header = AV_RB32(buf) | 0xffe00000;
Roberto Togni's avatar
Roberto Togni committed
2124

2125
    if (ff_mpa_check_header(header) < 0) { // Bad header, discard frame
Roberto Togni's avatar
Roberto Togni committed
2126 2127 2128 2129
        *data_size = 0;
        return buf_size;
    }

2130
    ff_mpegaudio_decode_header((MPADecodeHeader *)s, header);
Roberto Togni's avatar
Roberto Togni committed
2131 2132 2133
    /* update codec info */
    avctx->sample_rate = s->sample_rate;
    avctx->channels = s->nb_channels;
2134 2135
    if (!avctx->bit_rate)
        avctx->bit_rate = s->bit_rate;
Roberto Togni's avatar
Roberto Togni committed
2136 2137
    avctx->sub_id = s->layer;

2138
    s->frame_size = len;
Roberto Togni's avatar
Roberto Togni committed
2139 2140

    if (avctx->parse_only) {
2141
        out_size = buf_size;
Roberto Togni's avatar
Roberto Togni committed
2142
    } else {
2143
        out_size = mp_decode_frame(s, out_samples, buf, buf_size);
Roberto Togni's avatar
Roberto Togni committed
2144 2145 2146 2147 2148
    }

    *data_size = out_size;
    return buf_size;
}
2149
#endif /* CONFIG_MP3ADU_DECODER || CONFIG_MP3ADUFLOAT_DECODER */
Roberto Togni's avatar
Roberto Togni committed
2150

2151
#if CONFIG_MP3ON4_DECODER || CONFIG_MP3ON4FLOAT_DECODER
2152

2153 2154 2155 2156 2157 2158
/**
 * Context for MP3On4 decoder
 */
typedef struct MP3On4DecodeContext {
    int frames;   ///< number of mp3 frames per block (number of mp3 decoder instances)
    int syncword; ///< syncword patch
2159
    const uint8_t *coff; ///< channels offsets in output buffer
2160 2161 2162
    MPADecodeContext *mp3decctx[5]; ///< MPADecodeContext for every decoder instance
} MP3On4DecodeContext;

2163 2164
#include "mpeg4audio.h"

2165
/* Next 3 arrays are indexed by channel config number (passed via codecdata) */
2166
static const uint8_t mp3Frames[8] = {0,1,1,2,3,3,4,5};   /* number of mp3 decoder instances */
2167
/* offsets into output buffer, assume output order is FL FR BL BR C LFE */
2168
static const uint8_t chan_offset[8][5] = {
2169 2170 2171 2172 2173 2174 2175 2176 2177 2178 2179 2180 2181 2182
    {0},
    {0},            // C
    {0},            // FLR
    {2,0},          // C FLR
    {2,0,3},        // C FLR BS
    {4,0,2},        // C FLR BLRS
    {4,0,2,5},      // C FLR BLRS LFE
    {4,0,2,6,5},    // C FLR BLRS BLR LFE
};


static int decode_init_mp3on4(AVCodecContext * avctx)
{
    MP3On4DecodeContext *s = avctx->priv_data;
2183
    MPEG4AudioConfig cfg;
2184 2185 2186 2187 2188 2189 2190
    int i;

    if ((avctx->extradata_size < 2) || (avctx->extradata == NULL)) {
        av_log(avctx, AV_LOG_ERROR, "Codec extradata missing or too short.\n");
        return -1;
    }

2191 2192
    ff_mpeg4audio_get_config(&cfg, avctx->extradata, avctx->extradata_size);
    if (!cfg.chan_config || cfg.chan_config > 7) {
2193 2194 2195
        av_log(avctx, AV_LOG_ERROR, "Invalid channel config number.\n");
        return -1;
    }
2196 2197 2198
    s->frames = mp3Frames[cfg.chan_config];
    s->coff = chan_offset[cfg.chan_config];
    avctx->channels = ff_mpeg4audio_channels[cfg.chan_config];
2199

2200 2201 2202 2203 2204
    if (cfg.sample_rate < 16000)
        s->syncword = 0xffe00000;
    else
        s->syncword = 0xfff00000;

2205 2206 2207
    /* Init the first mp3 decoder in standard way, so that all tables get builded
     * We replace avctx->priv_data with the context of the first decoder so that
     * decode_init() does not have to be changed.
2208
     * Other decoders will be initialized here copying data from the first context
2209 2210 2211 2212 2213 2214 2215 2216 2217 2218 2219 2220 2221 2222 2223 2224
     */
    // Allocate zeroed memory for the first decoder context
    s->mp3decctx[0] = av_mallocz(sizeof(MPADecodeContext));
    // Put decoder context in place to make init_decode() happy
    avctx->priv_data = s->mp3decctx[0];
    decode_init(avctx);
    // Restore mp3on4 context pointer
    avctx->priv_data = s;
    s->mp3decctx[0]->adu_mode = 1; // Set adu mode

    /* Create a separate codec/context for each frame (first is already ok).
     * Each frame is 1 or 2 channels - up to 5 frames allowed
     */
    for (i = 1; i < s->frames; i++) {
        s->mp3decctx[i] = av_mallocz(sizeof(MPADecodeContext));
        s->mp3decctx[i]->adu_mode = 1;
2225
        s->mp3decctx[i]->avctx = avctx;
2226 2227 2228 2229 2230 2231
    }

    return 0;
}


2232
static av_cold int decode_close_mp3on4(AVCodecContext * avctx)
2233 2234 2235 2236 2237 2238 2239 2240 2241 2242 2243 2244 2245
{
    MP3On4DecodeContext *s = avctx->priv_data;
    int i;

    for (i = 0; i < s->frames; i++)
        if (s->mp3decctx[i])
            av_free(s->mp3decctx[i]);

    return 0;
}


static int decode_frame_mp3on4(AVCodecContext * avctx,
2246
                        void *data, int *data_size,
2247
                        AVPacket *avpkt)
2248
{
2249 2250
    const uint8_t *buf = avpkt->data;
    int buf_size = avpkt->size;
2251 2252
    MP3On4DecodeContext *s = avctx->priv_data;
    MPADecodeContext *m;
2253
    int fsize, len = buf_size, out_size = 0;
2254 2255 2256 2257
    uint32_t header;
    OUT_INT *out_samples = data;
    OUT_INT decoded_buf[MPA_FRAME_SIZE * MPA_MAX_CHANNELS];
    OUT_INT *outptr, *bp;
2258
    int fr, j, n;
2259

2260 2261 2262
    if(*data_size < MPA_FRAME_SIZE * MPA_MAX_CHANNELS * s->frames * sizeof(OUT_INT))
        return -1;

2263
    *data_size = 0;
2264
    // Discard too short frames
2265 2266
    if (buf_size < HEADER_SIZE)
        return -1;
2267 2268 2269 2270

    // If only one decoder interleave is not needed
    outptr = s->frames == 1 ? out_samples : decoded_buf;

2271 2272
    avctx->bit_rate = 0;

2273
    for (fr = 0; fr < s->frames; fr++) {
Baptiste Coudurier's avatar
Baptiste Coudurier committed
2274
        fsize = AV_RB16(buf) >> 4;
2275
        fsize = FFMIN3(fsize, len, MPA_MAX_CODED_FRAME_SIZE);
2276 2277 2278
        m = s->mp3decctx[fr];
        assert (m != NULL);

2279
        header = (AV_RB32(buf) & 0x000fffff) | s->syncword; // patch header
2280

2281 2282
        if (ff_mpa_check_header(header) < 0) // Bad header, discard block
            break;
2283

2284
        ff_mpegaudio_decode_header((MPADecodeHeader *)m, header);
2285
        out_size += mp_decode_frame(m, outptr, buf, fsize);
Baptiste Coudurier's avatar
Baptiste Coudurier committed
2286 2287
        buf += fsize;
        len -= fsize;
2288 2289

        if(s->frames > 1) {
2290
            n = m->avctx->frame_size*m->nb_channels;
2291
            /* interleave output data */
2292
            bp = out_samples + s->coff[fr];
2293 2294 2295
            if(m->nb_channels == 1) {
                for(j = 0; j < n; j++) {
                    *bp = decoded_buf[j];
Baptiste Coudurier's avatar
Baptiste Coudurier committed
2296
                    bp += avctx->channels;
2297 2298 2299 2300 2301
                }
            } else {
                for(j = 0; j < n; j++) {
                    bp[0] = decoded_buf[j++];
                    bp[1] = decoded_buf[j];
Baptiste Coudurier's avatar
Baptiste Coudurier committed
2302
                    bp += avctx->channels;
2303 2304 2305
                }
            }
        }
2306
        avctx->bit_rate += m->bit_rate;
2307 2308 2309 2310 2311 2312 2313 2314
    }

    /* update codec info */
    avctx->sample_rate = s->mp3decctx[0]->sample_rate;

    *data_size = out_size;
    return buf_size;
}
2315
#endif /* CONFIG_MP3ON4_DECODER || CONFIG_MP3ON4FLOAT_DECODER */
2316

2317
#if !CONFIG_FLOAT
2318
#if CONFIG_MP1_DECODER
2319 2320 2321
AVCodec mp1_decoder =
{
    "mp1",
2322
    AVMEDIA_TYPE_AUDIO,
2323 2324 2325 2326 2327 2328 2329 2330 2331 2332 2333
    CODEC_ID_MP1,
    sizeof(MPADecodeContext),
    decode_init,
    NULL,
    NULL,
    decode_frame,
    CODEC_CAP_PARSE_ONLY,
    .flush= flush,
    .long_name= NULL_IF_CONFIG_SMALL("MP1 (MPEG audio layer 1)"),
};
#endif
2334
#if CONFIG_MP2_DECODER
2335
AVCodec mp2_decoder =
Fabrice Bellard's avatar
Fabrice Bellard committed
2336
{
2337
    "mp2",
2338
    AVMEDIA_TYPE_AUDIO,
Fabrice Bellard's avatar
Fabrice Bellard committed
2339 2340 2341 2342 2343 2344
    CODEC_ID_MP2,
    sizeof(MPADecodeContext),
    decode_init,
    NULL,
    NULL,
    decode_frame,
2345
    CODEC_CAP_PARSE_ONLY,
2346
    .flush= flush,
2347
    .long_name= NULL_IF_CONFIG_SMALL("MP2 (MPEG audio layer 2)"),
Fabrice Bellard's avatar
Fabrice Bellard committed
2348
};
2349
#endif
2350
#if CONFIG_MP3_DECODER
2351 2352 2353
AVCodec mp3_decoder =
{
    "mp3",
2354
    AVMEDIA_TYPE_AUDIO,
2355
    CODEC_ID_MP3,
2356 2357 2358 2359 2360
    sizeof(MPADecodeContext),
    decode_init,
    NULL,
    NULL,
    decode_frame,
2361
    CODEC_CAP_PARSE_ONLY,
2362
    .flush= flush,
2363
    .long_name= NULL_IF_CONFIG_SMALL("MP3 (MPEG audio layer 3)"),
2364
};
2365
#endif
2366
#if CONFIG_MP3ADU_DECODER
Roberto Togni's avatar
Roberto Togni committed
2367 2368 2369
AVCodec mp3adu_decoder =
{
    "mp3adu",
2370
    AVMEDIA_TYPE_AUDIO,
Roberto Togni's avatar
Roberto Togni committed
2371 2372 2373 2374 2375 2376 2377
    CODEC_ID_MP3ADU,
    sizeof(MPADecodeContext),
    decode_init,
    NULL,
    NULL,
    decode_frame_adu,
    CODEC_CAP_PARSE_ONLY,
2378
    .flush= flush,
2379
    .long_name= NULL_IF_CONFIG_SMALL("ADU (Application Data Unit) MP3 (MPEG audio layer 3)"),
Roberto Togni's avatar
Roberto Togni committed
2380
};
2381
#endif
2382
#if CONFIG_MP3ON4_DECODER
2383 2384 2385
AVCodec mp3on4_decoder =
{
    "mp3on4",
2386
    AVMEDIA_TYPE_AUDIO,
2387 2388 2389 2390 2391 2392
    CODEC_ID_MP3ON4,
    sizeof(MP3On4DecodeContext),
    decode_init_mp3on4,
    NULL,
    decode_close_mp3on4,
    decode_frame_mp3on4,
2393
    .flush= flush,
2394
    .long_name= NULL_IF_CONFIG_SMALL("MP3onMP4"),
2395
};
2396
#endif
2397
#endif