annotate ffmpeg/libavcodec/ac3enc_float.c @ 13:844d341cf643 tip

Back up before ISMIR
author Yading Song <yading.song@eecs.qmul.ac.uk>
date Thu, 31 Oct 2013 13:17:06 +0000
parents 6840f77b83aa
children
rev   line source
yading@10 1 /*
yading@10 2 * The simplest AC-3 encoder
yading@10 3 * Copyright (c) 2000 Fabrice Bellard
yading@10 4 * Copyright (c) 2006-2010 Justin Ruggles <justin.ruggles@gmail.com>
yading@10 5 * Copyright (c) 2006-2010 Prakash Punnoor <prakash@punnoor.de>
yading@10 6 *
yading@10 7 * This file is part of FFmpeg.
yading@10 8 *
yading@10 9 * FFmpeg is free software; you can redistribute it and/or
yading@10 10 * modify it under the terms of the GNU Lesser General Public
yading@10 11 * License as published by the Free Software Foundation; either
yading@10 12 * version 2.1 of the License, or (at your option) any later version.
yading@10 13 *
yading@10 14 * FFmpeg is distributed in the hope that it will be useful,
yading@10 15 * but WITHOUT ANY WARRANTY; without even the implied warranty of
yading@10 16 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
yading@10 17 * Lesser General Public License for more details.
yading@10 18 *
yading@10 19 * You should have received a copy of the GNU Lesser General Public
yading@10 20 * License along with FFmpeg; if not, write to the Free Software
yading@10 21 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
yading@10 22 */
yading@10 23
yading@10 24 /**
yading@10 25 * @file
yading@10 26 * floating-point AC-3 encoder.
yading@10 27 */
yading@10 28
yading@10 29 #define CONFIG_AC3ENC_FLOAT 1
yading@10 30 #include "internal.h"
yading@10 31 #include "ac3enc.h"
yading@10 32 #include "eac3enc.h"
yading@10 33 #include "kbdwin.h"
yading@10 34
yading@10 35
yading@10 36 #if CONFIG_AC3_ENCODER
yading@10 37 #define AC3ENC_TYPE AC3ENC_TYPE_AC3
yading@10 38 #include "ac3enc_opts_template.c"
yading@10 39 static const AVClass ac3enc_class = {
yading@10 40 .class_name = "AC-3 Encoder",
yading@10 41 .item_name = av_default_item_name,
yading@10 42 .option = ac3_options,
yading@10 43 .version = LIBAVUTIL_VERSION_INT,
yading@10 44 };
yading@10 45 #endif
yading@10 46
yading@10 47 #include "ac3enc_template.c"
yading@10 48
yading@10 49
yading@10 50 /**
yading@10 51 * Finalize MDCT and free allocated memory.
yading@10 52 *
yading@10 53 * @param s AC-3 encoder private context
yading@10 54 */
yading@10 55 av_cold void ff_ac3_float_mdct_end(AC3EncodeContext *s)
yading@10 56 {
yading@10 57 ff_mdct_end(&s->mdct);
yading@10 58 av_freep(&s->mdct_window);
yading@10 59 }
yading@10 60
yading@10 61
yading@10 62 /**
yading@10 63 * Initialize MDCT tables.
yading@10 64 *
yading@10 65 * @param s AC-3 encoder private context
yading@10 66 * @return 0 on success, negative error code on failure
yading@10 67 */
yading@10 68 av_cold int ff_ac3_float_mdct_init(AC3EncodeContext *s)
yading@10 69 {
yading@10 70 float *window;
yading@10 71 int i, n, n2;
yading@10 72
yading@10 73 n = 1 << 9;
yading@10 74 n2 = n >> 1;
yading@10 75
yading@10 76 window = av_malloc(n * sizeof(*window));
yading@10 77 if (!window) {
yading@10 78 av_log(s->avctx, AV_LOG_ERROR, "Cannot allocate memory.\n");
yading@10 79 return AVERROR(ENOMEM);
yading@10 80 }
yading@10 81 ff_kbd_window_init(window, 5.0, n2);
yading@10 82 for (i = 0; i < n2; i++)
yading@10 83 window[n-1-i] = window[i];
yading@10 84 s->mdct_window = window;
yading@10 85
yading@10 86 return ff_mdct_init(&s->mdct, 9, 0, -2.0 / n);
yading@10 87 }
yading@10 88
yading@10 89
yading@10 90 /*
yading@10 91 * Apply KBD window to input samples prior to MDCT.
yading@10 92 */
yading@10 93 static void apply_window(void *dsp, float *output,
yading@10 94 const float *input, const float *window,
yading@10 95 unsigned int len)
yading@10 96 {
yading@10 97 AVFloatDSPContext *fdsp = dsp;
yading@10 98 fdsp->vector_fmul(output, input, window, len);
yading@10 99 }
yading@10 100
yading@10 101
yading@10 102 /*
yading@10 103 * Normalize the input samples.
yading@10 104 * Not needed for the floating-point encoder.
yading@10 105 */
yading@10 106 static int normalize_samples(AC3EncodeContext *s)
yading@10 107 {
yading@10 108 return 0;
yading@10 109 }
yading@10 110
yading@10 111
yading@10 112 /*
yading@10 113 * Scale MDCT coefficients from float to 24-bit fixed-point.
yading@10 114 */
yading@10 115 static void scale_coefficients(AC3EncodeContext *s)
yading@10 116 {
yading@10 117 int chan_size = AC3_MAX_COEFS * s->num_blocks;
yading@10 118 int cpl = s->cpl_on;
yading@10 119 s->ac3dsp.float_to_fixed24(s->fixed_coef_buffer + (chan_size * !cpl),
yading@10 120 s->mdct_coef_buffer + (chan_size * !cpl),
yading@10 121 chan_size * (s->channels + cpl));
yading@10 122 }
yading@10 123
yading@10 124 static void sum_square_butterfly(AC3EncodeContext *s, float sum[4],
yading@10 125 const float *coef0, const float *coef1,
yading@10 126 int len)
yading@10 127 {
yading@10 128 s->ac3dsp.sum_square_butterfly_float(sum, coef0, coef1, len);
yading@10 129 }
yading@10 130
yading@10 131 /*
yading@10 132 * Clip MDCT coefficients to allowable range.
yading@10 133 */
yading@10 134 static void clip_coefficients(DSPContext *dsp, float *coef, unsigned int len)
yading@10 135 {
yading@10 136 dsp->vector_clipf(coef, coef, COEF_MIN, COEF_MAX, len);
yading@10 137 }
yading@10 138
yading@10 139
yading@10 140 /*
yading@10 141 * Calculate a single coupling coordinate.
yading@10 142 */
yading@10 143 static CoefType calc_cpl_coord(CoefSumType energy_ch, CoefSumType energy_cpl)
yading@10 144 {
yading@10 145 float coord = 0.125;
yading@10 146 if (energy_cpl > 0)
yading@10 147 coord *= sqrtf(energy_ch / energy_cpl);
yading@10 148 return FFMIN(coord, COEF_MAX);
yading@10 149 }
yading@10 150
yading@10 151
yading@10 152 #if CONFIG_AC3_ENCODER
yading@10 153 AVCodec ff_ac3_encoder = {
yading@10 154 .name = "ac3",
yading@10 155 .type = AVMEDIA_TYPE_AUDIO,
yading@10 156 .id = AV_CODEC_ID_AC3,
yading@10 157 .priv_data_size = sizeof(AC3EncodeContext),
yading@10 158 .init = ff_ac3_encode_init,
yading@10 159 .encode2 = ff_ac3_float_encode_frame,
yading@10 160 .close = ff_ac3_encode_close,
yading@10 161 .sample_fmts = (const enum AVSampleFormat[]){ AV_SAMPLE_FMT_FLTP,
yading@10 162 AV_SAMPLE_FMT_NONE },
yading@10 163 .long_name = NULL_IF_CONFIG_SMALL("ATSC A/52A (AC-3)"),
yading@10 164 .priv_class = &ac3enc_class,
yading@10 165 .channel_layouts = ff_ac3_channel_layouts,
yading@10 166 .defaults = ac3_defaults,
yading@10 167 };
yading@10 168 #endif