aacenc.c 22.6 KB
Newer Older
K
Kostya Shishkov 已提交
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22
/*
 * AAC encoder
 * Copyright (C) 2008 Konstantin Shishkov
 *
 * This file is part of FFmpeg.
 *
 * FFmpeg is free software; you can redistribute it and/or
 * modify it under the terms of the GNU Lesser General Public
 * License as published by the Free Software Foundation; either
 * version 2.1 of the License, or (at your option) any later version.
 *
 * FFmpeg is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 * Lesser General Public License for more details.
 *
 * You should have received a copy of the GNU Lesser General Public
 * License along with FFmpeg; if not, write to the Free Software
 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
 */

/**
23
 * @file
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 152 153 154 155 156 157 158
};

/**
 * 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
    flush_put_bits(&pb);
}

static av_cold int aac_encode_init(AVCodecContext *avctx)
{
    AACEncContext *s = avctx->priv_data;
    int i;
159 160
    const uint8_t *sizes[2];
    int lengths[2];
K
Kostya Shishkov 已提交
161 162 163

    avctx->frame_size = 1024;

164 165
    for (i = 0; i < 16; i++)
        if (avctx->sample_rate == ff_mpeg4audio_sample_rates[i])
K
Kostya Shishkov 已提交
166
            break;
167
    if (i == 16) {
K
Kostya Shishkov 已提交
168 169 170
        av_log(avctx, AV_LOG_ERROR, "Unsupported sample rate %d\n", avctx->sample_rate);
        return -1;
    }
171
    if (avctx->channels > AAC_MAX_CHANNELS) {
K
Kostya Shishkov 已提交
172 173 174
        av_log(avctx, AV_LOG_ERROR, "Unsupported number of channels: %d\n", avctx->channels);
        return -1;
    }
175 176 177 178
    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;
    }
179 180 181 182
    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 已提交
183 184 185
    s->samplerate_index = i;

    dsputil_init(&s->dsp, avctx);
186 187
    ff_mdct_init(&s->mdct1024, 11, 0, 1.0);
    ff_mdct_init(&s->mdct128,   8, 0, 1.0);
188 189 190
    // 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);
191 192
    ff_init_ff_sine_windows(10);
    ff_init_ff_sine_windows(7);
K
Kostya Shishkov 已提交
193

194 195
    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 已提交
196
    avctx->extradata      = av_mallocz(2 + FF_INPUT_BUFFER_PADDING_SIZE);
K
Kostya Shishkov 已提交
197 198
    avctx->extradata_size = 2;
    put_audio_specific_config(avctx);
199

200 201
    sizes[0]   = swb_size_1024[i];
    sizes[1]   = swb_size_128[i];
202 203 204 205
    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);
206
    s->coder = &ff_aac_coders[2];
207 208

    s->lambda = avctx->global_quality ? avctx->global_quality : 120;
209 210

    ff_aac_tableinit();
211

K
Kostya Shishkov 已提交
212 213 214
    return 0;
}

215 216 217 218 219 220 221 222 223 224
static void apply_window_and_mdct(AVCodecContext *avctx, AACEncContext *s,
                                  SingleChannelElement *sce, short *audio, int channel)
{
    int i, j, k;
    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;

    if (sce->ics.window_sequence[0] != EIGHT_SHORT_SEQUENCE) {
        memcpy(s->output, sce->saved, sizeof(float)*1024);
225
        if (sce->ics.window_sequence[0] == LONG_STOP_SEQUENCE) {
226
            memset(s->output, 0, sizeof(s->output[0]) * 448);
227
            for (i = 448; i < 576; i++)
228
                s->output[i] = sce->saved[i] * pwindow[i - 448];
229
            for (i = 576; i < 704; i++)
230 231
                s->output[i] = sce->saved[i];
        }
232
        if (sce->ics.window_sequence[0] != LONG_START_SEQUENCE) {
A
Alex Converse 已提交
233
            for (i = 0, j = channel; i < 1024; i++, j += avctx->channels) {
234 235 236
                s->output[i+1024]         = audio[j] * lwindow[1024 - i - 1];
                sce->saved[i] = audio[j] * lwindow[i];
            }
237
        } else {
A
Alex Converse 已提交
238
            for (i = 0, j = channel; i < 448; i++, j += avctx->channels)
239
                s->output[i+1024]         = audio[j];
A
Alex Converse 已提交
240
            for (; i < 576; i++, j += avctx->channels)
241 242
                s->output[i+1024]         = audio[j] * swindow[576 - i - 1];
            memset(s->output+1024+576, 0, sizeof(s->output[0]) * 448);
A
Alex Converse 已提交
243
            for (i = 0, j = channel; i < 1024; i++, j += avctx->channels)
244 245 246
                sce->saved[i] = audio[j];
        }
        ff_mdct_calc(&s->mdct1024, sce->coeffs, s->output);
247
    } else {
248
        for (k = 0; k < 1024; k += 128) {
249
            for (i = 448 + k; i < 448 + k + 256; i++)
250 251 252 253 254 255 256
                s->output[i - 448 - k] = (i < 1024)
                                         ? sce->saved[i]
                                         : audio[channel + (i-1024)*avctx->channels];
            s->dsp.vector_fmul        (s->output,     k ?  swindow : pwindow, 128);
            s->dsp.vector_fmul_reverse(s->output+128, s->output+128, swindow, 128);
            ff_mdct_calc(&s->mdct128, sce->coeffs + k, s->output);
        }
A
Alex Converse 已提交
257
        for (i = 0, j = channel; i < 1024; i++, j += avctx->channels)
258 259 260 261
            sce->saved[i] = audio[j];
    }
}

K
Kostya Shishkov 已提交
262 263 264 265
/**
 * Encode ics_info element.
 * @see Table 4.6 (syntax of ics_info)
 */
K
Kostya Shishkov 已提交
266
static void put_ics_info(AACEncContext *s, IndividualChannelStream *info)
K
Kostya Shishkov 已提交
267
{
268
    int w;
K
Kostya Shishkov 已提交
269 270 271 272

    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]);
273
    if (info->window_sequence[0] != EIGHT_SHORT_SEQUENCE) {
K
Kostya Shishkov 已提交
274 275
        put_bits(&s->pb, 6, info->max_sfb);
        put_bits(&s->pb, 1, 0);            // no prediction
276
    } else {
K
Kostya Shishkov 已提交
277
        put_bits(&s->pb, 4, info->max_sfb);
278
        for (w = 1; w < 8; w++)
279
            put_bits(&s->pb, 1, !info->group_len[w]);
K
Kostya Shishkov 已提交
280 281 282
    }
}

K
Kostya Shishkov 已提交
283
/**
284 285
 * Encode MS data.
 * @see 4.6.8.1 "Joint Coding - M/S Stereo"
K
Kostya Shishkov 已提交
286
 */
287
static void encode_ms_info(PutBitContext *pb, ChannelElement *cpe)
K
Kostya Shishkov 已提交
288 289
{
    int i, w;
290 291

    put_bits(pb, 2, cpe->ms_mode);
292 293
    if (cpe->ms_mode == 1)
        for (w = 0; w < cpe->ch[0].ics.num_windows; w += cpe->ch[0].ics.group_len[w])
294
            for (i = 0; i < cpe->ch[0].ics.max_sfb; i++)
295 296 297 298 299 300 301 302 303 304 305
                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;
    int start, sum, maxsfb, cmaxsfb;

306
    for (ch = 0; ch < chans; ch++) {
307 308 309 310
        IndividualChannelStream *ics = &cpe->ch[ch].ics;
        start = 0;
        maxsfb = 0;
        cpe->ch[ch].pulse.num_pulse = 0;
311 312
        for (w = 0; w < ics->num_windows*16; w += 16) {
            for (g = 0; g < ics->num_swb; g++) {
313 314
                sum = 0;
                //apply M/S
315 316
                if (!ch && cpe->ms_mask[w + g]) {
                    for (i = 0; i < ics->swb_sizes[g]; i++) {
317 318 319 320 321 322
                        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];
            }
323 324
            for (cmaxsfb = ics->num_swb; cmaxsfb > 0 && cpe->ch[ch].zeroes[w+cmaxsfb-1]; cmaxsfb--)
                ;
325 326 327 328 329
            maxsfb = FFMAX(maxsfb, cmaxsfb);
        }
        ics->max_sfb = maxsfb;

        //adjust zero bands for window groups
330 331
        for (w = 0; w < ics->num_windows; w += ics->group_len[w]) {
            for (g = 0; g < ics->max_sfb; g++) {
332
                i = 1;
333 334
                for (w2 = w; w2 < w + ics->group_len[w]; w2++) {
                    if (!cpe->ch[ch].zeroes[w2*16 + g]) {
335 336 337 338 339 340 341 342 343
                        i = 0;
                        break;
                    }
                }
                cpe->ch[ch].zeroes[w*16 + g] = i;
            }
        }
    }

344
    if (chans > 1 && cpe->common_window) {
345 346 347 348 349
        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;
350 351
        for (w = 0; w < ics0->num_windows*16; w += 16)
            for (i = 0; i < ics0->max_sfb; i++)
352 353
                if (cpe->ms_mask[w+i])
                    msc++;
354 355 356 357
        if (msc == 0 || ics0->max_sfb == 0)
            cpe->ms_mode = 0;
        else
            cpe->ms_mode = msc < ics0->max_sfb ? 1 : 2;
358 359 360 361 362 363 364 365 366 367
    }
}

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

368
    for (w = 0; w < sce->ics.num_windows; w += sce->ics.group_len[w])
369 370 371 372 373 374
        s->coder->encode_window_bands_info(s, sce, w, sce->ics.group_len[w], s->lambda);
}

/**
 * Encode scalefactors.
 */
375 376
static void encode_scale_factors(AVCodecContext *avctx, AACEncContext *s,
                                 SingleChannelElement *sce)
377 378 379 380
{
    int off = sce->sf_idx[0], diff;
    int i, w;

381 382 383
    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]) {
384
                diff = sce->sf_idx[w*16 + i] - off + SCALE_DIFF_ZERO;
385 386
                if (diff < 0 || diff > 120)
                    av_log(avctx, AV_LOG_ERROR, "Scalefactor difference is too big to be coded\n");
387 388 389
                off = sce->sf_idx[w*16 + i];
                put_bits(&s->pb, ff_aac_scalefactor_bits[diff], ff_aac_scalefactor_code[diff]);
            }
K
Kostya Shishkov 已提交
390 391 392 393
        }
    }
}

394 395 396
/**
 * Encode pulse data.
 */
397
static void encode_pulses(AACEncContext *s, Pulse *pulse)
398 399 400 401
{
    int i;

    put_bits(&s->pb, 1, !!pulse->num_pulse);
402 403
    if (!pulse->num_pulse)
        return;
404 405 406

    put_bits(&s->pb, 2, pulse->num_pulse - 1);
    put_bits(&s->pb, 6, pulse->start);
407
    for (i = 0; i < pulse->num_pulse; i++) {
408
        put_bits(&s->pb, 5, pulse->pos[i]);
409 410 411 412 413 414 415
        put_bits(&s->pb, 4, pulse->amp[i]);
    }
}

/**
 * Encode spectral coefficients processed by psychoacoustic model.
 */
416
static void encode_spectral_coeffs(AACEncContext *s, SingleChannelElement *sce)
417
{
418
    int start, i, w, w2;
419

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

438 439 440
/**
 * Encode one channel of audio data.
 */
441 442 443
static int encode_individual_channel(AVCodecContext *avctx, AACEncContext *s,
                                     SingleChannelElement *sce,
                                     int common_window)
444 445
{
    put_bits(&s->pb, 8, sce->sf_idx[0]);
446 447
    if (!common_window)
        put_ics_info(s, &sce->ics);
448 449 450 451 452 453 454 455 456
    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 已提交
457 458 459
/**
 * Write some auxiliary information about the created AAC file.
 */
460 461
static void put_bitstream_info(AVCodecContext *avctx, AACEncContext *s,
                               const char *name)
K
Kostya Shishkov 已提交
462 463 464 465
{
    int i, namelen, padbits;

    namelen = strlen(name) + 2;
466
    put_bits(&s->pb, 3, TYPE_FIL);
K
Kostya Shishkov 已提交
467
    put_bits(&s->pb, 4, FFMIN(namelen, 15));
468
    if (namelen >= 15)
K
Kostya Shishkov 已提交
469 470 471 472
        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);
473
    for (i = 0; i < namelen - 2; i++)
K
Kostya Shishkov 已提交
474 475 476 477
        put_bits(&s->pb, 8, name[i]);
    put_bits(&s->pb, 12 - padbits, 0);
}

478 479 480 481 482 483 484 485 486
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];
487
    FFPsyWindowInfo windows[AAC_MAX_CHANNELS];
488

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

    start_ch = 0;
514
    for (i = 0; i < chan_map[0]; i++) {
515
        FFPsyWindowInfo* wi = windows + start_ch;
516 517 518
        tag      = chan_map[i+1];
        chans    = tag == TYPE_CPE ? 2 : 1;
        cpe      = &s->cpe[i];
519
        samples2 = samples + start_ch;
520
        la       = samples2 + (448+64) * avctx->channels + start_ch;
521 522
        if (!data)
            la = NULL;
523
        for (j = 0; j < chans; j++) {
524 525
            IndividualChannelStream *ics = &cpe->ch[j].ics;
            int k;
526 527 528 529 530 531
            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 {
532
            wi[j] = ff_psy_suggest_window(&s->psy, samples2, la, start_ch + j, ics->window_sequence[0]);
533
            }
534 535 536 537 538 539
            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];
540
            ics->num_swb            = tag == TYPE_LFE ? 12 : s->psy.num_bands[ics->num_windows == 8];
541
            for (k = 0; k < ics->num_windows; k++)
542 543 544 545
                ics->group_len[k] = wi[j].grouping[k];

            s->cur_channel = start_ch + j;
            apply_window_and_mdct(avctx, s, &cpe->ch[j], samples2, j);
546 547 548
        }
        start_ch += chans;
    }
549 550
    do {
        int frame_bits;
A
Alex Converse 已提交
551 552 553 554 555 556 557 558 559 560 561
        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];
            for (j = 0; j < chans; j++) {
562
                s->cur_channel = start_ch + j;
A
Alex Converse 已提交
563
                ff_psy_set_band_info(&s->psy, s->cur_channel, cpe->ch[j].coeffs, &wi[j]);
A
Alex Converse 已提交
564 565 566 567 568 569 570 571 572 573 574 575 576
                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;
                    }
577 578
                }
            }
579
            s->cur_channel = start_ch;
A
Alex Converse 已提交
580 581 582 583 584 585 586 587 588 589 590
            if (cpe->common_window && s->coder->search_for_ms)
                s->coder->search_for_ms(s, cpe, s->lambda);
            adjust_frame_information(s, cpe, chans);
            put_bits(&s->pb, 3, tag);
            put_bits(&s->pb, 4, chan_el_counter[tag]++);
            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);
                }
591
            }
A
Alex Converse 已提交
592 593 594 595 596
            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;
597 598
        }

599 600 601 602 603 604 605 606
        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);

607 608 609 610 611
    put_bits(&s->pb, 3, TYPE_END);
    flush_put_bits(&s->pb);
    avctx->frame_bits = put_bits_count(&s->pb);

    // rate control stuff
612
    if (!(avctx->flags & CODEC_FLAG_QSCALE)) {
613 614
        float ratio = avctx->bit_rate * 1024.0f / avctx->sample_rate / avctx->frame_bits;
        s->lambda *= ratio;
615
        s->lambda = FFMIN(s->lambda, 65536.f);
616 617
    }

618
    if (!data)
619
        s->last_frame = 1;
620 621
    memcpy(s->samples, s->samples + 1024 * avctx->channels,
           1024 * avctx->channels * sizeof(s->samples[0]));
622 623 624
    return put_bits_count(&s->pb)>>3;
}

K
Kostya Shishkov 已提交
625 626 627 628 629 630
static av_cold int aac_encode_end(AVCodecContext *avctx)
{
    AACEncContext *s = avctx->priv_data;

    ff_mdct_end(&s->mdct1024);
    ff_mdct_end(&s->mdct128);
631 632
    ff_psy_end(&s->psy);
    ff_psy_preprocess_end(s->psypp);
K
Kostya Shishkov 已提交
633 634 635 636 637 638 639
    av_freep(&s->samples);
    av_freep(&s->cpe);
    return 0;
}

AVCodec aac_encoder = {
    "aac",
640
    AVMEDIA_TYPE_AUDIO,
K
Kostya Shishkov 已提交
641 642 643 644 645
    CODEC_ID_AAC,
    sizeof(AACEncContext),
    aac_encode_init,
    aac_encode_frame,
    aac_encode_end,
646
    .capabilities = CODEC_CAP_SMALL_LAST_FRAME | CODEC_CAP_DELAY | CODEC_CAP_EXPERIMENTAL,
647
    .sample_fmts = (const enum SampleFormat[]){SAMPLE_FMT_S16,SAMPLE_FMT_NONE},
K
Kostya Shishkov 已提交
648 649
    .long_name = NULL_IF_CONFIG_SMALL("Advanced Audio Coding"),
};