aacenc.c 22.7 KB
Newer Older
K
Kostya Shishkov 已提交
1 2 3 4
/*
 * AAC encoder
 * Copyright (C) 2008 Konstantin Shishkov
 *
5
 * This file is part of Libav.
K
Kostya Shishkov 已提交
6
 *
7
 * Libav is free software; you can redistribute it and/or
K
Kostya Shishkov 已提交
8 9 10 11
 * 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.
 *
12
 * Libav is distributed in the hope that it will be useful,
K
Kostya Shishkov 已提交
13 14 15 16 17
 * 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
18
 * License along with Libav; if not, write to the Free Software
K
Kostya Shishkov 已提交
19 20 21 22
 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
 */

/**
23
 * @file
K
Kostya Shishkov 已提交
24 25 26 27 28
 * AAC encoder
 */

/***********************************
 *              TODOs:
29
 * add sane pulse detection
K
Kostya Shishkov 已提交
30
 * add temporal noise shaping
K
Kostya Shishkov 已提交
31 32 33
 ***********************************/

#include "avcodec.h"
34
#include "put_bits.h"
K
Kostya Shishkov 已提交
35 36 37 38 39
#include "dsputil.h"
#include "mpeg4audio.h"

#include "aac.h"
#include "aactab.h"
40 41 42
#include "aacenc.h"

#include "psymodel.h"
K
Kostya Shishkov 已提交
43

44 45
#define AAC_MAX_CHANNELS 6

K
Kostya Shishkov 已提交
46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88
static const uint8_t swb_size_1024_96[] = {
    4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 8, 8, 8, 8, 8,
    12, 12, 12, 12, 12, 16, 16, 24, 28, 36, 44,
    64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64
};

static const uint8_t swb_size_1024_64[] = {
    4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 8, 8, 8, 8,
    12, 12, 12, 16, 16, 16, 20, 24, 24, 28, 36,
    40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40
};

static const uint8_t swb_size_1024_48[] = {
    4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 8, 8, 8, 8, 8, 8, 8,
    12, 12, 12, 12, 16, 16, 20, 20, 24, 24, 28, 28,
    32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32,
    96
};

static const uint8_t swb_size_1024_32[] = {
    4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 8, 8, 8, 8, 8, 8, 8,
    12, 12, 12, 12, 16, 16, 20, 20, 24, 24, 28, 28,
    32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32
};

static const uint8_t swb_size_1024_24[] = {
    4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8,
    12, 12, 12, 12, 16, 16, 16, 20, 20, 24, 24, 28, 28,
    32, 36, 36, 40, 44, 48, 52, 52, 64, 64, 64, 64, 64
};

static const uint8_t swb_size_1024_16[] = {
    8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8,
    12, 12, 12, 12, 12, 12, 12, 12, 12, 16, 16, 16, 16, 20, 20, 20, 24, 24, 28, 28,
    32, 36, 40, 40, 44, 48, 52, 56, 60, 64, 64, 64
};

static const uint8_t swb_size_1024_8[] = {
    12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12,
    16, 16, 16, 16, 16, 16, 16, 20, 20, 20, 20, 24, 24, 24, 28, 28,
    32, 36, 36, 40, 44, 48, 52, 56, 60, 64, 80
};

89
static const uint8_t *swb_size_1024[] = {
K
Kostya Shishkov 已提交
90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115
    swb_size_1024_96, swb_size_1024_96, swb_size_1024_64,
    swb_size_1024_48, swb_size_1024_48, swb_size_1024_32,
    swb_size_1024_24, swb_size_1024_24, swb_size_1024_16,
    swb_size_1024_16, swb_size_1024_16, swb_size_1024_8
};

static const uint8_t swb_size_128_96[] = {
    4, 4, 4, 4, 4, 4, 8, 8, 8, 16, 28, 36
};

static const uint8_t swb_size_128_48[] = {
    4, 4, 4, 4, 4, 8, 8, 8, 12, 12, 12, 16, 16, 16
};

static const uint8_t swb_size_128_24[] = {
    4, 4, 4, 4, 4, 4, 4, 8, 8, 8, 12, 12, 16, 16, 20
};

static const uint8_t swb_size_128_16[] = {
    4, 4, 4, 4, 4, 4, 4, 4, 8, 8, 12, 12, 16, 20, 20
};

static const uint8_t swb_size_128_8[] = {
    4, 4, 4, 4, 4, 4, 4, 8, 8, 8, 8, 12, 16, 20, 20
};

116
static const uint8_t *swb_size_128[] = {
K
Kostya Shishkov 已提交
117 118 119 120 121 122 123 124 125 126
    /* the last entry on the following row is swb_size_128_64 but is a
       duplicate of swb_size_128_96 */
    swb_size_128_96, swb_size_128_96, swb_size_128_96,
    swb_size_128_48, swb_size_128_48, swb_size_128_48,
    swb_size_128_24, swb_size_128_24, swb_size_128_16,
    swb_size_128_16, swb_size_128_16, swb_size_128_8
};

/** default channel configurations */
static const uint8_t aac_chan_configs[6][5] = {
127 128 129 130 131 132
 {1, TYPE_SCE},                               // 1 channel  - single channel element
 {1, TYPE_CPE},                               // 2 channels - channel pair
 {2, TYPE_SCE, TYPE_CPE},                     // 3 channels - center + stereo
 {3, TYPE_SCE, TYPE_CPE, TYPE_SCE},           // 4 channels - front center + stereo + back center
 {3, TYPE_SCE, TYPE_CPE, TYPE_CPE},           // 5 channels - front center + stereo + back stereo
 {4, TYPE_SCE, TYPE_CPE, TYPE_CPE, TYPE_LFE}, // 6 channels - front center + stereo + back stereo + LFE
K
Kostya Shishkov 已提交
133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151
};

/**
 * Make AAC audio config object.
 * @see 1.6.2.1 "Syntax - AudioSpecificConfig"
 */
static void put_audio_specific_config(AVCodecContext *avctx)
{
    PutBitContext pb;
    AACEncContext *s = avctx->priv_data;

    init_put_bits(&pb, avctx->extradata, avctx->extradata_size*8);
    put_bits(&pb, 5, 2); //object type - AAC-LC
    put_bits(&pb, 4, s->samplerate_index); //sample rate index
    put_bits(&pb, 4, avctx->channels);
    //GASpecificConfig
    put_bits(&pb, 1, 0); //frame length - 1024 samples
    put_bits(&pb, 1, 0); //does not depend on core coder
    put_bits(&pb, 1, 0); //is not extension
A
Alex Converse 已提交
152 153

    //Explicitly Mark SBR absent
154
    put_bits(&pb, 11, 0x2b7); //sync extension
A
Alex Converse 已提交
155 156
    put_bits(&pb, 5,  AOT_SBR);
    put_bits(&pb, 1,  0);
K
Kostya Shishkov 已提交
157 158 159 160 161 162 163
    flush_put_bits(&pb);
}

static av_cold int aac_encode_init(AVCodecContext *avctx)
{
    AACEncContext *s = avctx->priv_data;
    int i;
164 165
    const uint8_t *sizes[2];
    int lengths[2];
K
Kostya Shishkov 已提交
166 167 168

    avctx->frame_size = 1024;

169 170
    for (i = 0; i < 16; i++)
        if (avctx->sample_rate == ff_mpeg4audio_sample_rates[i])
K
Kostya Shishkov 已提交
171
            break;
172
    if (i == 16) {
K
Kostya Shishkov 已提交
173 174 175
        av_log(avctx, AV_LOG_ERROR, "Unsupported sample rate %d\n", avctx->sample_rate);
        return -1;
    }
176
    if (avctx->channels > AAC_MAX_CHANNELS) {
K
Kostya Shishkov 已提交
177 178 179
        av_log(avctx, AV_LOG_ERROR, "Unsupported number of channels: %d\n", avctx->channels);
        return -1;
    }
180 181 182 183
    if (avctx->profile != FF_PROFILE_UNKNOWN && avctx->profile != FF_PROFILE_AAC_LOW) {
        av_log(avctx, AV_LOG_ERROR, "Unsupported profile %d\n", avctx->profile);
        return -1;
    }
184 185 186 187
    if (1024.0 * avctx->bit_rate / avctx->sample_rate > 6144 * avctx->channels) {
        av_log(avctx, AV_LOG_ERROR, "Too many bits per frame requested\n");
        return -1;
    }
K
Kostya Shishkov 已提交
188 189 190
    s->samplerate_index = i;

    dsputil_init(&s->dsp, avctx);
191 192
    ff_mdct_init(&s->mdct1024, 11, 0, 1.0);
    ff_mdct_init(&s->mdct128,   8, 0, 1.0);
193 194 195
    // window init
    ff_kbd_window_init(ff_aac_kbd_long_1024, 4.0, 1024);
    ff_kbd_window_init(ff_aac_kbd_short_128, 6.0, 128);
196 197
    ff_init_ff_sine_windows(10);
    ff_init_ff_sine_windows(7);
K
Kostya Shishkov 已提交
198

199 200
    s->samples            = av_malloc(2 * 1024 * avctx->channels * sizeof(s->samples[0]));
    s->cpe                = av_mallocz(sizeof(ChannelElement) * aac_chan_configs[avctx->channels-1][0]);
A
Alex Converse 已提交
201 202
    avctx->extradata      = av_mallocz(5 + FF_INPUT_BUFFER_PADDING_SIZE);
    avctx->extradata_size = 5;
K
Kostya Shishkov 已提交
203
    put_audio_specific_config(avctx);
204

205 206
    sizes[0]   = swb_size_1024[i];
    sizes[1]   = swb_size_128[i];
207 208 209 210
    lengths[0] = ff_aac_num_swb_1024[i];
    lengths[1] = ff_aac_num_swb_128[i];
    ff_psy_init(&s->psy, avctx, 2, sizes, lengths);
    s->psypp = ff_psy_preprocess_init(avctx);
211
    s->coder = &ff_aac_coders[2];
212 213

    s->lambda = avctx->global_quality ? avctx->global_quality : 120;
214 215

    ff_aac_tableinit();
216

K
Kostya Shishkov 已提交
217 218 219
    return 0;
}

220
static void apply_window_and_mdct(AVCodecContext *avctx, AACEncContext *s,
221
                                  SingleChannelElement *sce, short *audio)
222
{
223 224
    int i, k;
    const int chans = avctx->channels;
225 226 227
    const float * lwindow = sce->ics.use_kb_window[0] ? ff_aac_kbd_long_1024 : ff_sine_1024;
    const float * swindow = sce->ics.use_kb_window[0] ? ff_aac_kbd_short_128 : ff_sine_128;
    const float * pwindow = sce->ics.use_kb_window[1] ? ff_aac_kbd_short_128 : ff_sine_128;
Y
Young Han Lee 已提交
228
    float *output = sce->ret;
229 230

    if (sce->ics.window_sequence[0] != EIGHT_SHORT_SEQUENCE) {
Y
Young Han Lee 已提交
231
        memcpy(output, sce->saved, sizeof(float)*1024);
232
        if (sce->ics.window_sequence[0] == LONG_STOP_SEQUENCE) {
Y
Young Han Lee 已提交
233
            memset(output, 0, sizeof(output[0]) * 448);
234
            for (i = 448; i < 576; i++)
Y
Young Han Lee 已提交
235
                output[i] = sce->saved[i] * pwindow[i - 448];
236
            for (i = 576; i < 704; i++)
Y
Young Han Lee 已提交
237
                output[i] = sce->saved[i];
238
        }
239
        if (sce->ics.window_sequence[0] != LONG_START_SEQUENCE) {
240
            for (i = 0; i < 1024; i++) {
Y
Young Han Lee 已提交
241
                output[i+1024]         = audio[i * chans] * lwindow[1024 - i - 1];
242
                sce->saved[i] = audio[i * chans] * lwindow[i];
243
            }
244
        } else {
245
            for (i = 0; i < 448; i++)
Y
Young Han Lee 已提交
246
                output[i+1024]         = audio[i * chans];
247
            for (; i < 576; i++)
Y
Young Han Lee 已提交
248 249
                output[i+1024]         = audio[i * chans] * swindow[576 - i - 1];
            memset(output+1024+576, 0, sizeof(output[0]) * 448);
250 251
            for (i = 0; i < 1024; i++)
                sce->saved[i] = audio[i * chans];
252
        }
253
        s->mdct1024.mdct_calc(&s->mdct1024, sce->coeffs, output);
254
    } else {
255
        for (k = 0; k < 1024; k += 128) {
256
            for (i = 448 + k; i < 448 + k + 256; i++)
Y
Young Han Lee 已提交
257
                output[i - 448 - k] = (i < 1024)
258
                                         ? sce->saved[i]
259
                                         : audio[(i-1024)*chans];
Y
Young Han Lee 已提交
260 261
            s->dsp.vector_fmul        (output,     output, k ?  swindow : pwindow, 128);
            s->dsp.vector_fmul_reverse(output+128, output+128, swindow, 128);
262
            s->mdct128.mdct_calc(&s->mdct128, sce->coeffs + k, output);
263
        }
264 265
        for (i = 0; i < 1024; i++)
            sce->saved[i] = audio[i * chans];
266 267 268
    }
}

K
Kostya Shishkov 已提交
269 270 271 272
/**
 * Encode ics_info element.
 * @see Table 4.6 (syntax of ics_info)
 */
K
Kostya Shishkov 已提交
273
static void put_ics_info(AACEncContext *s, IndividualChannelStream *info)
K
Kostya Shishkov 已提交
274
{
275
    int w;
K
Kostya Shishkov 已提交
276 277 278 279

    put_bits(&s->pb, 1, 0);                // ics_reserved bit
    put_bits(&s->pb, 2, info->window_sequence[0]);
    put_bits(&s->pb, 1, info->use_kb_window[0]);
280
    if (info->window_sequence[0] != EIGHT_SHORT_SEQUENCE) {
K
Kostya Shishkov 已提交
281 282
        put_bits(&s->pb, 6, info->max_sfb);
        put_bits(&s->pb, 1, 0);            // no prediction
283
    } else {
K
Kostya Shishkov 已提交
284
        put_bits(&s->pb, 4, info->max_sfb);
285
        for (w = 1; w < 8; w++)
286
            put_bits(&s->pb, 1, !info->group_len[w]);
K
Kostya Shishkov 已提交
287 288 289
    }
}

K
Kostya Shishkov 已提交
290
/**
291 292
 * Encode MS data.
 * @see 4.6.8.1 "Joint Coding - M/S Stereo"
K
Kostya Shishkov 已提交
293
 */
294
static void encode_ms_info(PutBitContext *pb, ChannelElement *cpe)
K
Kostya Shishkov 已提交
295 296
{
    int i, w;
297 298

    put_bits(pb, 2, cpe->ms_mode);
299 300
    if (cpe->ms_mode == 1)
        for (w = 0; w < cpe->ch[0].ics.num_windows; w += cpe->ch[0].ics.group_len[w])
301
            for (i = 0; i < cpe->ch[0].ics.max_sfb; i++)
302 303 304 305 306 307 308 309 310
                put_bits(pb, 1, cpe->ms_mask[w*16 + i]);
}

/**
 * Produce integer coefficients from scalefactors provided by the model.
 */
static void adjust_frame_information(AACEncContext *apc, ChannelElement *cpe, int chans)
{
    int i, w, w2, g, ch;
311
    int start, maxsfb, cmaxsfb;
312

313
    for (ch = 0; ch < chans; ch++) {
314 315 316 317
        IndividualChannelStream *ics = &cpe->ch[ch].ics;
        start = 0;
        maxsfb = 0;
        cpe->ch[ch].pulse.num_pulse = 0;
318 319
        for (w = 0; w < ics->num_windows*16; w += 16) {
            for (g = 0; g < ics->num_swb; g++) {
320
                //apply M/S
321
                if (cpe->common_window && !ch && cpe->ms_mask[w + g]) {
322
                    for (i = 0; i < ics->swb_sizes[g]; i++) {
323 324 325 326 327 328
                        cpe->ch[0].coeffs[start+i] = (cpe->ch[0].coeffs[start+i] + cpe->ch[1].coeffs[start+i]) / 2.0;
                        cpe->ch[1].coeffs[start+i] =  cpe->ch[0].coeffs[start+i] - cpe->ch[1].coeffs[start+i];
                    }
                }
                start += ics->swb_sizes[g];
            }
329 330
            for (cmaxsfb = ics->num_swb; cmaxsfb > 0 && cpe->ch[ch].zeroes[w+cmaxsfb-1]; cmaxsfb--)
                ;
331 332 333 334 335
            maxsfb = FFMAX(maxsfb, cmaxsfb);
        }
        ics->max_sfb = maxsfb;

        //adjust zero bands for window groups
336 337
        for (w = 0; w < ics->num_windows; w += ics->group_len[w]) {
            for (g = 0; g < ics->max_sfb; g++) {
338
                i = 1;
339 340
                for (w2 = w; w2 < w + ics->group_len[w]; w2++) {
                    if (!cpe->ch[ch].zeroes[w2*16 + g]) {
341 342 343 344 345 346 347 348 349
                        i = 0;
                        break;
                    }
                }
                cpe->ch[ch].zeroes[w*16 + g] = i;
            }
        }
    }

350
    if (chans > 1 && cpe->common_window) {
351 352 353 354 355
        IndividualChannelStream *ics0 = &cpe->ch[0].ics;
        IndividualChannelStream *ics1 = &cpe->ch[1].ics;
        int msc = 0;
        ics0->max_sfb = FFMAX(ics0->max_sfb, ics1->max_sfb);
        ics1->max_sfb = ics0->max_sfb;
356 357
        for (w = 0; w < ics0->num_windows*16; w += 16)
            for (i = 0; i < ics0->max_sfb; i++)
358 359
                if (cpe->ms_mask[w+i])
                    msc++;
360 361 362 363
        if (msc == 0 || ics0->max_sfb == 0)
            cpe->ms_mode = 0;
        else
            cpe->ms_mode = msc < ics0->max_sfb ? 1 : 2;
364 365 366 367 368 369 370 371 372 373
    }
}

/**
 * Encode scalefactor band coding type.
 */
static void encode_band_info(AACEncContext *s, SingleChannelElement *sce)
{
    int w;

374
    for (w = 0; w < sce->ics.num_windows; w += sce->ics.group_len[w])
375 376 377 378 379 380
        s->coder->encode_window_bands_info(s, sce, w, sce->ics.group_len[w], s->lambda);
}

/**
 * Encode scalefactors.
 */
381 382
static void encode_scale_factors(AVCodecContext *avctx, AACEncContext *s,
                                 SingleChannelElement *sce)
383 384 385 386
{
    int off = sce->sf_idx[0], diff;
    int i, w;

387 388 389
    for (w = 0; w < sce->ics.num_windows; w += sce->ics.group_len[w]) {
        for (i = 0; i < sce->ics.max_sfb; i++) {
            if (!sce->zeroes[w*16 + i]) {
390
                diff = sce->sf_idx[w*16 + i] - off + SCALE_DIFF_ZERO;
391 392
                if (diff < 0 || diff > 120)
                    av_log(avctx, AV_LOG_ERROR, "Scalefactor difference is too big to be coded\n");
393 394 395
                off = sce->sf_idx[w*16 + i];
                put_bits(&s->pb, ff_aac_scalefactor_bits[diff], ff_aac_scalefactor_code[diff]);
            }
K
Kostya Shishkov 已提交
396 397 398 399
        }
    }
}

400 401 402
/**
 * Encode pulse data.
 */
403
static void encode_pulses(AACEncContext *s, Pulse *pulse)
404 405 406 407
{
    int i;

    put_bits(&s->pb, 1, !!pulse->num_pulse);
408 409
    if (!pulse->num_pulse)
        return;
410 411 412

    put_bits(&s->pb, 2, pulse->num_pulse - 1);
    put_bits(&s->pb, 6, pulse->start);
413
    for (i = 0; i < pulse->num_pulse; i++) {
414
        put_bits(&s->pb, 5, pulse->pos[i]);
415 416 417 418 419 420 421
        put_bits(&s->pb, 4, pulse->amp[i]);
    }
}

/**
 * Encode spectral coefficients processed by psychoacoustic model.
 */
422
static void encode_spectral_coeffs(AACEncContext *s, SingleChannelElement *sce)
423
{
424
    int start, i, w, w2;
425

426
    for (w = 0; w < sce->ics.num_windows; w += sce->ics.group_len[w]) {
427
        start = 0;
428 429
        for (i = 0; i < sce->ics.max_sfb; i++) {
            if (sce->zeroes[w*16 + i]) {
430
                start += sce->ics.swb_sizes[i];
431 432
                continue;
            }
433
            for (w2 = w; w2 < w + sce->ics.group_len[w]; w2++)
434
                s->coder->quantize_and_encode_band(s, &s->pb, sce->coeffs + start + w2*128,
435 436 437 438
                                                   sce->ics.swb_sizes[i],
                                                   sce->sf_idx[w*16 + i],
                                                   sce->band_type[w*16 + i],
                                                   s->lambda);
439
            start += sce->ics.swb_sizes[i];
440 441 442 443
        }
    }
}

444 445 446
/**
 * Encode one channel of audio data.
 */
447 448 449
static int encode_individual_channel(AVCodecContext *avctx, AACEncContext *s,
                                     SingleChannelElement *sce,
                                     int common_window)
450 451
{
    put_bits(&s->pb, 8, sce->sf_idx[0]);
452 453
    if (!common_window)
        put_ics_info(s, &sce->ics);
454 455 456 457 458 459 460 461 462
    encode_band_info(s, sce);
    encode_scale_factors(avctx, s, sce);
    encode_pulses(s, &sce->pulse);
    put_bits(&s->pb, 1, 0); //tns
    put_bits(&s->pb, 1, 0); //ssr
    encode_spectral_coeffs(s, sce);
    return 0;
}

K
Kostya Shishkov 已提交
463 464 465
/**
 * Write some auxiliary information about the created AAC file.
 */
466 467
static void put_bitstream_info(AVCodecContext *avctx, AACEncContext *s,
                               const char *name)
K
Kostya Shishkov 已提交
468 469 470 471
{
    int i, namelen, padbits;

    namelen = strlen(name) + 2;
472
    put_bits(&s->pb, 3, TYPE_FIL);
K
Kostya Shishkov 已提交
473
    put_bits(&s->pb, 4, FFMIN(namelen, 15));
474
    if (namelen >= 15)
K
Kostya Shishkov 已提交
475 476 477 478
        put_bits(&s->pb, 8, namelen - 16);
    put_bits(&s->pb, 4, 0); //extension type - filler
    padbits = 8 - (put_bits_count(&s->pb) & 7);
    align_put_bits(&s->pb);
479
    for (i = 0; i < namelen - 2; i++)
K
Kostya Shishkov 已提交
480 481 482 483
        put_bits(&s->pb, 8, name[i]);
    put_bits(&s->pb, 12 - padbits, 0);
}

484 485 486 487 488 489 490 491 492
static int aac_encode_frame(AVCodecContext *avctx,
                            uint8_t *frame, int buf_size, void *data)
{
    AACEncContext *s = avctx->priv_data;
    int16_t *samples = s->samples, *samples2, *la;
    ChannelElement *cpe;
    int i, j, chans, tag, start_ch;
    const uint8_t *chan_map = aac_chan_configs[avctx->channels-1];
    int chan_el_counter[4];
493
    FFPsyWindowInfo windows[AAC_MAX_CHANNELS];
494

495
    if (s->last_frame)
496
        return 0;
497 498
    if (data) {
        if (!s->psypp) {
499 500
            memcpy(s->samples + 1024 * avctx->channels, data,
                   1024 * avctx->channels * sizeof(s->samples[0]));
501
        } else {
502 503
            start_ch = 0;
            samples2 = s->samples + 1024 * avctx->channels;
504
            for (i = 0; i < chan_map[0]; i++) {
505 506
                tag = chan_map[i+1];
                chans = tag == TYPE_CPE ? 2 : 1;
507 508
                ff_psy_preprocess(s->psypp, (uint16_t*)data + start_ch,
                                  samples2 + start_ch, start_ch, chans);
509 510 511 512
                start_ch += chans;
            }
        }
    }
513
    if (!avctx->frame_number) {
514 515
        memcpy(s->samples, s->samples + 1024 * avctx->channels,
               1024 * avctx->channels * sizeof(s->samples[0]));
516 517 518 519
        return 0;
    }

    start_ch = 0;
520
    for (i = 0; i < chan_map[0]; i++) {
521
        FFPsyWindowInfo* wi = windows + start_ch;
522 523 524
        tag      = chan_map[i+1];
        chans    = tag == TYPE_CPE ? 2 : 1;
        cpe      = &s->cpe[i];
525
        for (j = 0; j < chans; j++) {
526 527
            IndividualChannelStream *ics = &cpe->ch[j].ics;
            int k;
528 529 530 531 532
            int cur_channel = start_ch + j;
            samples2 = samples + cur_channel;
            la       = samples2 + (448+64) * avctx->channels;
            if (!data)
                la = NULL;
533 534 535 536 537 538
            if (tag == TYPE_LFE) {
                wi[j].window_type[0] = ONLY_LONG_SEQUENCE;
                wi[j].window_shape   = 0;
                wi[j].num_windows    = 1;
                wi[j].grouping[0]    = 1;
            } else {
539
                wi[j] = ff_psy_suggest_window(&s->psy, samples2, la, cur_channel,
A
Alex Converse 已提交
540
                                              ics->window_sequence[0]);
541
            }
542 543 544 545 546 547
            ics->window_sequence[1] = ics->window_sequence[0];
            ics->window_sequence[0] = wi[j].window_type[0];
            ics->use_kb_window[1]   = ics->use_kb_window[0];
            ics->use_kb_window[0]   = wi[j].window_shape;
            ics->num_windows        = wi[j].num_windows;
            ics->swb_sizes          = s->psy.bands    [ics->num_windows == 8];
548
            ics->num_swb            = tag == TYPE_LFE ? 12 : s->psy.num_bands[ics->num_windows == 8];
549
            for (k = 0; k < ics->num_windows; k++)
550 551
                ics->group_len[k] = wi[j].grouping[k];

552
            apply_window_and_mdct(avctx, s, &cpe->ch[j], samples2);
553 554 555
        }
        start_ch += chans;
    }
556 557
    do {
        int frame_bits;
A
Alex Converse 已提交
558 559 560 561 562 563 564 565 566 567
        init_put_bits(&s->pb, frame, buf_size*8);
        if ((avctx->frame_number & 0xFF)==1 && !(avctx->flags & CODEC_FLAG_BITEXACT))
            put_bitstream_info(avctx, s, LIBAVCODEC_IDENT);
        start_ch = 0;
        memset(chan_el_counter, 0, sizeof(chan_el_counter));
        for (i = 0; i < chan_map[0]; i++) {
            FFPsyWindowInfo* wi = windows + start_ch;
            tag      = chan_map[i+1];
            chans    = tag == TYPE_CPE ? 2 : 1;
            cpe      = &s->cpe[i];
A
Alex Converse 已提交
568 569
            put_bits(&s->pb, 3, tag);
            put_bits(&s->pb, 4, chan_el_counter[tag]++);
A
Alex Converse 已提交
570
            for (j = 0; j < chans; j++) {
571
                s->cur_channel = start_ch + j;
A
Alex Converse 已提交
572
                ff_psy_set_band_info(&s->psy, s->cur_channel, cpe->ch[j].coeffs, &wi[j]);
A
Alex Converse 已提交
573 574 575 576 577 578 579 580 581 582 583 584 585
                s->coder->search_for_quantizers(avctx, s, &cpe->ch[j], s->lambda);
            }
            cpe->common_window = 0;
            if (chans > 1
                && wi[0].window_type[0] == wi[1].window_type[0]
                && wi[0].window_shape   == wi[1].window_shape) {

                cpe->common_window = 1;
                for (j = 0; j < wi[0].num_windows; j++) {
                    if (wi[0].grouping[j] != wi[1].grouping[j]) {
                        cpe->common_window = 0;
                        break;
                    }
586 587
                }
            }
588
            s->cur_channel = start_ch;
A
Alex Converse 已提交
589 590 591 592 593 594 595 596 597
            if (cpe->common_window && s->coder->search_for_ms)
                s->coder->search_for_ms(s, cpe, s->lambda);
            adjust_frame_information(s, cpe, chans);
            if (chans == 2) {
                put_bits(&s->pb, 1, cpe->common_window);
                if (cpe->common_window) {
                    put_ics_info(s, &cpe->ch[0].ics);
                    encode_ms_info(&s->pb, cpe);
                }
598
            }
A
Alex Converse 已提交
599 600 601 602 603
            for (j = 0; j < chans; j++) {
                s->cur_channel = start_ch + j;
                encode_individual_channel(avctx, s, &cpe->ch[j], cpe->common_window);
            }
            start_ch += chans;
604 605
        }

606 607 608 609 610 611 612 613
        frame_bits = put_bits_count(&s->pb);
        if (frame_bits <= 6144 * avctx->channels - 3)
            break;

        s->lambda *= avctx->bit_rate * 1024.0f / avctx->sample_rate / frame_bits;

    } while (1);

614 615 616 617 618
    put_bits(&s->pb, 3, TYPE_END);
    flush_put_bits(&s->pb);
    avctx->frame_bits = put_bits_count(&s->pb);

    // rate control stuff
619
    if (!(avctx->flags & CODEC_FLAG_QSCALE)) {
620 621
        float ratio = avctx->bit_rate * 1024.0f / avctx->sample_rate / avctx->frame_bits;
        s->lambda *= ratio;
622
        s->lambda = FFMIN(s->lambda, 65536.f);
623 624
    }

625
    if (!data)
626
        s->last_frame = 1;
627 628
    memcpy(s->samples, s->samples + 1024 * avctx->channels,
           1024 * avctx->channels * sizeof(s->samples[0]));
629 630 631
    return put_bits_count(&s->pb)>>3;
}

K
Kostya Shishkov 已提交
632 633 634 635 636 637
static av_cold int aac_encode_end(AVCodecContext *avctx)
{
    AACEncContext *s = avctx->priv_data;

    ff_mdct_end(&s->mdct1024);
    ff_mdct_end(&s->mdct128);
638 639
    ff_psy_end(&s->psy);
    ff_psy_preprocess_end(s->psypp);
K
Kostya Shishkov 已提交
640 641 642 643 644
    av_freep(&s->samples);
    av_freep(&s->cpe);
    return 0;
}

645
AVCodec ff_aac_encoder = {
K
Kostya Shishkov 已提交
646
    "aac",
647
    AVMEDIA_TYPE_AUDIO,
K
Kostya Shishkov 已提交
648 649 650 651 652
    CODEC_ID_AAC,
    sizeof(AACEncContext),
    aac_encode_init,
    aac_encode_frame,
    aac_encode_end,
653
    .capabilities = CODEC_CAP_SMALL_LAST_FRAME | CODEC_CAP_DELAY | CODEC_CAP_EXPERIMENTAL,
654
    .sample_fmts = (const enum AVSampleFormat[]){AV_SAMPLE_FMT_S16,AV_SAMPLE_FMT_NONE},
K
Kostya Shishkov 已提交
655 656
    .long_name = NULL_IF_CONFIG_SMALL("Advanced Audio Coding"),
};