X-Git-Url: https://git.sesse.net/?a=blobdiff_plain;f=libavcodec%2Fatrac1.c;h=0ba2cf6bddea09f3494c79c94be732cb199f65f8;hb=5f0bb0baefd506d684adfa1ad4259c65973b455e;hp=cac82b36787876eaec3fb43ebd93909a587bb0bc;hpb=9caab878292b61cb2325dc3caa8e426df0e762ea;p=ffmpeg diff --git a/libavcodec/atrac1.c b/libavcodec/atrac1.c index cac82b36787..0ba2cf6bdde 100644 --- a/libavcodec/atrac1.c +++ b/libavcodec/atrac1.c @@ -3,27 +3,27 @@ * Copyright (c) 2009 Maxim Poliakovski * Copyright (c) 2009 Benjamin Larsson * - * This file is part of FFmpeg. + * This file is part of Libav. * - * FFmpeg is free software; you can redistribute it and/or + * Libav 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, + * Libav 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 + * License along with Libav; if not, write to the Free Software * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA */ /** - * @file libavcodec/atrac1.c + * @file * Atrac 1 compatible decoder. - * This decoder handles raw ATRAC1 data. + * This decoder handles raw ATRAC1 data and probably SDDS data. */ /* Many thanks to Tim Craig for all the help! */ @@ -35,6 +35,8 @@ #include "avcodec.h" #include "get_bits.h" #include "dsputil.h" +#include "fft.h" +#include "sinewin.h" #include "atrac.h" #include "atrac1data.h" @@ -58,11 +60,11 @@ typedef struct { int log2_block_count[AT1_QMF_BANDS]; ///< log2 number of blocks in a band int num_bfus; ///< number of Block Floating Units float* spectrum[2]; - DECLARE_ALIGNED_16(float, spec1[AT1_SU_SAMPLES]); ///< mdct buffer - DECLARE_ALIGNED_16(float, spec2[AT1_SU_SAMPLES]); ///< mdct buffer - DECLARE_ALIGNED_16(float, fst_qmf_delay[46]); ///< delay line for the 1st stacked QMF filter - DECLARE_ALIGNED_16(float, snd_qmf_delay[46]); ///< delay line for the 2nd stacked QMF filter - DECLARE_ALIGNED_16(float, last_qmf_delay[256+23]); ///< delay line for the last stacked QMF filter + DECLARE_ALIGNED(32, float, spec1)[AT1_SU_SAMPLES]; ///< mdct buffer + DECLARE_ALIGNED(32, float, spec2)[AT1_SU_SAMPLES]; ///< mdct buffer + DECLARE_ALIGNED(32, float, fst_qmf_delay)[46]; ///< delay line for the 1st stacked QMF filter + DECLARE_ALIGNED(32, float, snd_qmf_delay)[46]; ///< delay line for the 2nd stacked QMF filter + DECLARE_ALIGNED(32, float, last_qmf_delay)[256+23]; ///< delay line for the last stacked QMF filter } AT1SUCtx; /** @@ -70,20 +72,18 @@ typedef struct { */ typedef struct { AT1SUCtx SUs[AT1_MAX_CHANNELS]; ///< channel sound unit - DECLARE_ALIGNED_16(float, spec[AT1_SU_SAMPLES]); ///< the mdct spectrum buffer + DECLARE_ALIGNED(32, float, spec)[AT1_SU_SAMPLES]; ///< the mdct spectrum buffer - DECLARE_ALIGNED_16(float, low[256]); - DECLARE_ALIGNED_16(float, mid[256]); - DECLARE_ALIGNED_16(float, high[512]); + DECLARE_ALIGNED(32, float, low)[256]; + DECLARE_ALIGNED(32, float, mid)[256]; + DECLARE_ALIGNED(32, float, high)[512]; float* bands[3]; - DECLARE_ALIGNED_16(float, out_samples[AT1_MAX_CHANNELS][AT1_SU_SAMPLES]); + DECLARE_ALIGNED(32, float, out_samples)[AT1_MAX_CHANNELS][AT1_SU_SAMPLES]; FFTContext mdct_ctx[3]; int channels; DSPContext dsp; } AT1Ctx; -DECLARE_ALIGNED_16(static float, short_window[32]); - /** size of the transform in samples in the long mode for each QMF band */ static const uint16_t samples_per_band[3] = {128, 128, 256}; static const uint8_t mdct_long_nbits[3] = {7, 7, 8}; @@ -100,7 +100,7 @@ static void at1_imdct(AT1Ctx *q, float *spec, float *out, int nbits, for (i = 0; i < transf_size / 2; i++) FFSWAP(float, spec[i], spec[transf_size - 1 - i]); } - ff_imdct_half(mdct_context, out, spec); + mdct_context->imdct_half(mdct_context, out, spec); } @@ -110,6 +110,9 @@ static int at1_imdct_block(AT1SUCtx* su, AT1Ctx *q) unsigned int start_pos, ref_pos = 0, pos = 0; for (band_num = 0; band_num < AT1_QMF_BANDS; band_num++) { + float *prev_buf; + int j; + band_samples = samples_per_band[band_num]; log2_block_count = su->log2_block_count[band_num]; @@ -117,42 +120,38 @@ static int at1_imdct_block(AT1SUCtx* su, AT1Ctx *q) /* 4 for short mode(low/middle bands) and 8 for short mode(high band)*/ num_blocks = 1 << log2_block_count; - /* mdct block size in samples: 128 (long mode, low & mid bands), */ - /* 256 (long mode, high band) and 32 (short mode, all bands) */ - block_size = band_samples >> log2_block_count; + if (num_blocks == 1) { + /* mdct block size in samples: 128 (long mode, low & mid bands), */ + /* 256 (long mode, high band) and 32 (short mode, all bands) */ + block_size = band_samples >> log2_block_count; - /* calc transform size in bits according to the block_size_mode */ - nbits = mdct_long_nbits[band_num] - log2_block_count; + /* calc transform size in bits according to the block_size_mode */ + nbits = mdct_long_nbits[band_num] - log2_block_count; - if (nbits != 5 && nbits != 7 && nbits != 8) - return -1; + if (nbits != 5 && nbits != 7 && nbits != 8) + return -1; + } else { + block_size = 32; + nbits = 5; + } - if (num_blocks == 1) { - /* long blocks */ - at1_imdct(q, &q->spec[pos], &su->spectrum[0][ref_pos], nbits, band_num); - pos += block_size; // move to the next mdct block in the spectrum + start_pos = 0; + prev_buf = &su->spectrum[1][ref_pos + band_samples - 16]; + for (j=0; j < num_blocks; j++) { + at1_imdct(q, &q->spec[pos], &su->spectrum[0][ref_pos + start_pos], nbits, band_num); - /* overlap and window long blocks */ - q->dsp.vector_fmul_window(q->bands[band_num], &su->spectrum[1][ref_pos + band_samples - 16], - &su->spectrum[0][ref_pos], short_window, 0, 16); - memcpy(q->bands[band_num] + 32, &su->spectrum[0][ref_pos + 16], 240 * sizeof(float)); - } else { - /* short blocks */ - float *prev_buf; - start_pos = 0; - prev_buf = &su->spectrum[1][ref_pos + band_samples - 16]; - for (; num_blocks != 0; num_blocks--) { - at1_imdct(q, &q->spec[pos], &su->spectrum[0][ref_pos + start_pos], 5, band_num); - - /* overlap and window between short blocks */ - q->dsp.vector_fmul_window(&q->bands[band_num][start_pos], prev_buf, - &su->spectrum[0][ref_pos + start_pos], short_window, 0, 16); - - prev_buf = &su->spectrum[0][ref_pos+start_pos + 16]; - start_pos += 32; // use hardcoded block_size - pos += 32; - } + /* overlap and window */ + q->dsp.vector_fmul_window(&q->bands[band_num][start_pos], prev_buf, + &su->spectrum[0][ref_pos + start_pos], ff_sine_32, 16); + + prev_buf = &su->spectrum[0][ref_pos+start_pos + 16]; + start_pos += block_size; + pos += block_size; } + + if (num_blocks == 1) + memcpy(q->bands[band_num] + 32, &su->spectrum[0][ref_pos + 16], 240 * sizeof(float)); + ref_pos += band_samples; } @@ -225,7 +224,7 @@ static int at1_unpack_dequant(GetBitContext* gb, AT1SUCtx* su, int num_specs = specs_per_bfu[bfu_num]; int word_len = !!idwls[bfu_num] + idwls[bfu_num]; - float scale_factor = sf_table[idsfs[bfu_num]]; + float scale_factor = ff_atrac_sf_table[idsfs[bfu_num]]; bits_used += word_len * num_specs; /* add number of bits consumed by current BFU */ /* check for bitstream overflow */ @@ -254,7 +253,7 @@ static int at1_unpack_dequant(GetBitContext* gb, AT1SUCtx* su, } -void at1_subband_synthesis(AT1Ctx *q, AT1SUCtx* su, float *pOut) +static void at1_subband_synthesis(AT1Ctx *q, AT1SUCtx* su, float *pOut) { float temp[256]; float iqmf_temp[512 + 46]; @@ -307,20 +306,15 @@ static int atrac1_decode_frame(AVCodecContext *avctx, void *data, at1_subband_synthesis(q, su, q->out_samples[ch]); } - /* round, convert to 16bit and interleave */ + /* interleave; FIXME, should create/use a DSP function */ if (q->channels == 1) { /* mono */ - q->dsp.vector_clipf(samples, q->out_samples[0], -32700.0 / (1 << 15), - 32700.0 / (1 << 15), AT1_SU_SAMPLES); + memcpy(samples, q->out_samples[0], AT1_SU_SAMPLES * 4); } else { /* stereo */ for (i = 0; i < AT1_SU_SAMPLES; i++) { - samples[i * 2] = av_clipf(q->out_samples[0][i], - -32700.0 / (1 << 15), - 32700.0 / (1 << 15)); - samples[i * 2 + 1] = av_clipf(q->out_samples[1][i], - -32700.0 / (1 << 15), - 32700.0 / (1 << 15)); + samples[i * 2] = q->out_samples[0][i]; + samples[i * 2 + 1] = q->out_samples[1][i]; } } @@ -333,7 +327,7 @@ static av_cold int atrac1_decode_init(AVCodecContext *avctx) { AT1Ctx *q = avctx->priv_data; - avctx->sample_fmt = SAMPLE_FMT_FLT; + avctx->sample_fmt = AV_SAMPLE_FMT_FLT; q->channels = avctx->channels; @@ -342,7 +336,7 @@ static av_cold int atrac1_decode_init(AVCodecContext *avctx) ff_mdct_init(&q->mdct_ctx[1], 8, 1, -1.0/ (1 << 15)); ff_mdct_init(&q->mdct_ctx[2], 9, 1, -1.0/ (1 << 15)); - ff_sine_window_init(short_window, 32); + ff_init_ff_sine_windows(5); atrac_generate_tables(); @@ -372,9 +366,9 @@ static av_cold int atrac1_decode_end(AVCodecContext * avctx) { } -AVCodec atrac1_decoder = { +AVCodec ff_atrac1_decoder = { .name = "atrac1", - .type = CODEC_TYPE_AUDIO, + .type = AVMEDIA_TYPE_AUDIO, .id = CODEC_ID_ATRAC1, .priv_data_size = sizeof(AT1Ctx), .init = atrac1_decode_init,