yading@10: /* yading@10: * RealAudio Lossless decoder yading@10: * yading@10: * Copyright (c) 2012 Konstantin Shishkov yading@10: * yading@10: * This file is part of Libav. yading@10: * yading@10: * Libav is free software; you can redistribute it and/or yading@10: * modify it under the terms of the GNU Lesser General Public yading@10: * License as published by the Free Software Foundation; either yading@10: * version 2.1 of the License, or (at your option) any later version. yading@10: * yading@10: * Libav is distributed in the hope that it will be useful, yading@10: * but WITHOUT ANY WARRANTY; without even the implied warranty of yading@10: * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU yading@10: * Lesser General Public License for more details. yading@10: * yading@10: * You should have received a copy of the GNU Lesser General Public yading@10: * License along with Libav; if not, write to the Free Software yading@10: * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA yading@10: */ yading@10: yading@10: /** yading@10: * @file yading@10: * This is a decoder for Real Audio Lossless format. yading@10: * Dedicated to the mastermind behind it, Ralph Wiggum. yading@10: */ yading@10: yading@10: #include "libavutil/channel_layout.h" yading@10: #include "avcodec.h" yading@10: #include "get_bits.h" yading@10: #include "golomb.h" yading@10: #include "internal.h" yading@10: #include "unary.h" yading@10: #include "ralfdata.h" yading@10: yading@10: #define FILTER_NONE 0 yading@10: #define FILTER_RAW 642 yading@10: yading@10: typedef struct VLCSet { yading@10: VLC filter_params; yading@10: VLC bias; yading@10: VLC coding_mode; yading@10: VLC filter_coeffs[10][11]; yading@10: VLC short_codes[15]; yading@10: VLC long_codes[125]; yading@10: } VLCSet; yading@10: yading@10: #define RALF_MAX_PKT_SIZE 8192 yading@10: yading@10: typedef struct RALFContext { yading@10: int version; yading@10: int max_frame_size; yading@10: VLCSet sets[3]; yading@10: int32_t channel_data[2][4096]; yading@10: yading@10: int filter_params; ///< combined filter parameters for the current channel data yading@10: int filter_length; ///< length of the filter for the current channel data yading@10: int filter_bits; ///< filter precision for the current channel data yading@10: int32_t filter[64]; yading@10: yading@10: int bias[2]; ///< a constant value added to channel data after filtering yading@10: yading@10: int num_blocks; ///< number of blocks inside the frame yading@10: int sample_offset; yading@10: int block_size[1 << 12]; ///< size of the blocks yading@10: int block_pts[1 << 12]; ///< block start time (in milliseconds) yading@10: yading@10: uint8_t pkt[16384]; yading@10: int has_pkt; yading@10: } RALFContext; yading@10: yading@10: #define MAX_ELEMS 644 // no RALF table uses more than that yading@10: yading@10: static int init_ralf_vlc(VLC *vlc, const uint8_t *data, int elems) yading@10: { yading@10: uint8_t lens[MAX_ELEMS]; yading@10: uint16_t codes[MAX_ELEMS]; yading@10: int counts[17], prefixes[18]; yading@10: int i, cur_len; yading@10: int max_bits = 0; yading@10: int nb = 0; yading@10: yading@10: for (i = 0; i <= 16; i++) yading@10: counts[i] = 0; yading@10: for (i = 0; i < elems; i++) { yading@10: cur_len = (nb ? *data & 0xF : *data >> 4) + 1; yading@10: counts[cur_len]++; yading@10: max_bits = FFMAX(max_bits, cur_len); yading@10: lens[i] = cur_len; yading@10: data += nb; yading@10: nb ^= 1; yading@10: } yading@10: prefixes[1] = 0; yading@10: for (i = 1; i <= 16; i++) yading@10: prefixes[i + 1] = (prefixes[i] + counts[i]) << 1; yading@10: yading@10: for (i = 0; i < elems; i++) yading@10: codes[i] = prefixes[lens[i]]++; yading@10: yading@10: return ff_init_vlc_sparse(vlc, FFMIN(max_bits, 9), elems, yading@10: lens, 1, 1, codes, 2, 2, NULL, 0, 0, 0); yading@10: } yading@10: yading@10: static av_cold int decode_close(AVCodecContext *avctx) yading@10: { yading@10: RALFContext *ctx = avctx->priv_data; yading@10: int i, j, k; yading@10: yading@10: for (i = 0; i < 3; i++) { yading@10: ff_free_vlc(&ctx->sets[i].filter_params); yading@10: ff_free_vlc(&ctx->sets[i].bias); yading@10: ff_free_vlc(&ctx->sets[i].coding_mode); yading@10: for (j = 0; j < 10; j++) yading@10: for (k = 0; k < 11; k++) yading@10: ff_free_vlc(&ctx->sets[i].filter_coeffs[j][k]); yading@10: for (j = 0; j < 15; j++) yading@10: ff_free_vlc(&ctx->sets[i].short_codes[j]); yading@10: for (j = 0; j < 125; j++) yading@10: ff_free_vlc(&ctx->sets[i].long_codes[j]); yading@10: } yading@10: yading@10: return 0; yading@10: } yading@10: yading@10: static av_cold int decode_init(AVCodecContext *avctx) yading@10: { yading@10: RALFContext *ctx = avctx->priv_data; yading@10: int i, j, k; yading@10: int ret; yading@10: yading@10: if (avctx->extradata_size < 24 || memcmp(avctx->extradata, "LSD:", 4)) { yading@10: av_log(avctx, AV_LOG_ERROR, "Extradata is not groovy, dude\n"); yading@10: return AVERROR_INVALIDDATA; yading@10: } yading@10: yading@10: ctx->version = AV_RB16(avctx->extradata + 4); yading@10: if (ctx->version != 0x103) { yading@10: avpriv_request_sample(avctx, "Unknown version %X", ctx->version); yading@10: return AVERROR_PATCHWELCOME; yading@10: } yading@10: yading@10: avctx->channels = AV_RB16(avctx->extradata + 8); yading@10: avctx->sample_rate = AV_RB32(avctx->extradata + 12); yading@10: if (avctx->channels < 1 || avctx->channels > 2 yading@10: || avctx->sample_rate < 8000 || avctx->sample_rate > 96000) { yading@10: av_log(avctx, AV_LOG_ERROR, "Invalid coding parameters %d Hz %d ch\n", yading@10: avctx->sample_rate, avctx->channels); yading@10: return AVERROR_INVALIDDATA; yading@10: } yading@10: avctx->sample_fmt = AV_SAMPLE_FMT_S16P; yading@10: avctx->channel_layout = (avctx->channels == 2) ? AV_CH_LAYOUT_STEREO yading@10: : AV_CH_LAYOUT_MONO; yading@10: yading@10: ctx->max_frame_size = AV_RB32(avctx->extradata + 16); yading@10: if (ctx->max_frame_size > (1 << 20) || !ctx->max_frame_size) { yading@10: av_log(avctx, AV_LOG_ERROR, "invalid frame size %d\n", yading@10: ctx->max_frame_size); yading@10: } yading@10: ctx->max_frame_size = FFMAX(ctx->max_frame_size, avctx->sample_rate); yading@10: yading@10: for (i = 0; i < 3; i++) { yading@10: ret = init_ralf_vlc(&ctx->sets[i].filter_params, filter_param_def[i], yading@10: FILTERPARAM_ELEMENTS); yading@10: if (ret < 0) { yading@10: decode_close(avctx); yading@10: return ret; yading@10: } yading@10: ret = init_ralf_vlc(&ctx->sets[i].bias, bias_def[i], BIAS_ELEMENTS); yading@10: if (ret < 0) { yading@10: decode_close(avctx); yading@10: return ret; yading@10: } yading@10: ret = init_ralf_vlc(&ctx->sets[i].coding_mode, coding_mode_def[i], yading@10: CODING_MODE_ELEMENTS); yading@10: if (ret < 0) { yading@10: decode_close(avctx); yading@10: return ret; yading@10: } yading@10: for (j = 0; j < 10; j++) { yading@10: for (k = 0; k < 11; k++) { yading@10: ret = init_ralf_vlc(&ctx->sets[i].filter_coeffs[j][k], yading@10: filter_coeffs_def[i][j][k], yading@10: FILTER_COEFFS_ELEMENTS); yading@10: if (ret < 0) { yading@10: decode_close(avctx); yading@10: return ret; yading@10: } yading@10: } yading@10: } yading@10: for (j = 0; j < 15; j++) { yading@10: ret = init_ralf_vlc(&ctx->sets[i].short_codes[j], yading@10: short_codes_def[i][j], SHORT_CODES_ELEMENTS); yading@10: if (ret < 0) { yading@10: decode_close(avctx); yading@10: return ret; yading@10: } yading@10: } yading@10: for (j = 0; j < 125; j++) { yading@10: ret = init_ralf_vlc(&ctx->sets[i].long_codes[j], yading@10: long_codes_def[i][j], LONG_CODES_ELEMENTS); yading@10: if (ret < 0) { yading@10: decode_close(avctx); yading@10: return ret; yading@10: } yading@10: } yading@10: } yading@10: yading@10: return 0; yading@10: } yading@10: yading@10: static inline int extend_code(GetBitContext *gb, int val, int range, int bits) yading@10: { yading@10: if (val == 0) { yading@10: val = -range - get_ue_golomb(gb); yading@10: } else if (val == range * 2) { yading@10: val = range + get_ue_golomb(gb); yading@10: } else { yading@10: val -= range; yading@10: } yading@10: if (bits) yading@10: val = (val << bits) | get_bits(gb, bits); yading@10: return val; yading@10: } yading@10: yading@10: static int decode_channel(RALFContext *ctx, GetBitContext *gb, int ch, yading@10: int length, int mode, int bits) yading@10: { yading@10: int i, t; yading@10: int code_params; yading@10: VLCSet *set = ctx->sets + mode; yading@10: VLC *code_vlc; int range, range2, add_bits; yading@10: int *dst = ctx->channel_data[ch]; yading@10: yading@10: ctx->filter_params = get_vlc2(gb, set->filter_params.table, 9, 2); yading@10: ctx->filter_bits = (ctx->filter_params - 2) >> 6; yading@10: ctx->filter_length = ctx->filter_params - (ctx->filter_bits << 6) - 1; yading@10: yading@10: if (ctx->filter_params == FILTER_RAW) { yading@10: for (i = 0; i < length; i++) yading@10: dst[i] = get_bits(gb, bits); yading@10: ctx->bias[ch] = 0; yading@10: return 0; yading@10: } yading@10: yading@10: ctx->bias[ch] = get_vlc2(gb, set->bias.table, 9, 2); yading@10: ctx->bias[ch] = extend_code(gb, ctx->bias[ch], 127, 4); yading@10: yading@10: if (ctx->filter_params == FILTER_NONE) { yading@10: memset(dst, 0, sizeof(*dst) * length); yading@10: return 0; yading@10: } yading@10: yading@10: if (ctx->filter_params > 1) { yading@10: int cmode = 0, coeff = 0; yading@10: VLC *vlc = set->filter_coeffs[ctx->filter_bits] + 5; yading@10: yading@10: add_bits = ctx->filter_bits; yading@10: yading@10: for (i = 0; i < ctx->filter_length; i++) { yading@10: t = get_vlc2(gb, vlc[cmode].table, vlc[cmode].bits, 2); yading@10: t = extend_code(gb, t, 21, add_bits); yading@10: if (!cmode) yading@10: coeff -= 12 << add_bits; yading@10: coeff = t - coeff; yading@10: ctx->filter[i] = coeff; yading@10: yading@10: cmode = coeff >> add_bits; yading@10: if (cmode < 0) { yading@10: cmode = -1 - av_log2(-cmode); yading@10: if (cmode < -5) yading@10: cmode = -5; yading@10: } else if (cmode > 0) { yading@10: cmode = 1 + av_log2(cmode); yading@10: if (cmode > 5) yading@10: cmode = 5; yading@10: } yading@10: } yading@10: } yading@10: yading@10: code_params = get_vlc2(gb, set->coding_mode.table, set->coding_mode.bits, 2); yading@10: if (code_params >= 15) { yading@10: add_bits = av_clip((code_params / 5 - 3) / 2, 0, 10); yading@10: if (add_bits > 9 && (code_params % 5) != 2) yading@10: add_bits--; yading@10: range = 10; yading@10: range2 = 21; yading@10: code_vlc = set->long_codes + code_params - 15; yading@10: } else { yading@10: add_bits = 0; yading@10: range = 6; yading@10: range2 = 13; yading@10: code_vlc = set->short_codes + code_params; yading@10: } yading@10: yading@10: for (i = 0; i < length; i += 2) { yading@10: int code1, code2; yading@10: yading@10: t = get_vlc2(gb, code_vlc->table, code_vlc->bits, 2); yading@10: code1 = t / range2; yading@10: code2 = t % range2; yading@10: dst[i] = extend_code(gb, code1, range, 0) << add_bits; yading@10: dst[i + 1] = extend_code(gb, code2, range, 0) << add_bits; yading@10: if (add_bits) { yading@10: dst[i] |= get_bits(gb, add_bits); yading@10: dst[i + 1] |= get_bits(gb, add_bits); yading@10: } yading@10: } yading@10: yading@10: return 0; yading@10: } yading@10: yading@10: static void apply_lpc(RALFContext *ctx, int ch, int length, int bits) yading@10: { yading@10: int i, j, acc; yading@10: int *audio = ctx->channel_data[ch]; yading@10: int bias = 1 << (ctx->filter_bits - 1); yading@10: int max_clip = (1 << bits) - 1, min_clip = -max_clip - 1; yading@10: yading@10: for (i = 1; i < length; i++) { yading@10: int flen = FFMIN(ctx->filter_length, i); yading@10: yading@10: acc = 0; yading@10: for (j = 0; j < flen; j++) yading@10: acc += ctx->filter[j] * audio[i - j - 1]; yading@10: if (acc < 0) { yading@10: acc = (acc + bias - 1) >> ctx->filter_bits; yading@10: acc = FFMAX(acc, min_clip); yading@10: } else { yading@10: acc = (acc + bias) >> ctx->filter_bits; yading@10: acc = FFMIN(acc, max_clip); yading@10: } yading@10: audio[i] += acc; yading@10: } yading@10: } yading@10: yading@10: static int decode_block(AVCodecContext *avctx, GetBitContext *gb, yading@10: int16_t *dst0, int16_t *dst1) yading@10: { yading@10: RALFContext *ctx = avctx->priv_data; yading@10: int len, ch, ret; yading@10: int dmode, mode[2], bits[2]; yading@10: int *ch0, *ch1; yading@10: int i, t, t2; yading@10: yading@10: len = 12 - get_unary(gb, 0, 6); yading@10: yading@10: if (len <= 7) len ^= 1; // codes for length = 6 and 7 are swapped yading@10: len = 1 << len; yading@10: yading@10: if (ctx->sample_offset + len > ctx->max_frame_size) { yading@10: av_log(avctx, AV_LOG_ERROR, yading@10: "Decoder's stomach is crying, it ate too many samples\n"); yading@10: return AVERROR_INVALIDDATA; yading@10: } yading@10: yading@10: if (avctx->channels > 1) yading@10: dmode = get_bits(gb, 2) + 1; yading@10: else yading@10: dmode = 0; yading@10: yading@10: mode[0] = (dmode == 4) ? 1 : 0; yading@10: mode[1] = (dmode >= 2) ? 2 : 0; yading@10: bits[0] = 16; yading@10: bits[1] = (mode[1] == 2) ? 17 : 16; yading@10: yading@10: for (ch = 0; ch < avctx->channels; ch++) { yading@10: if ((ret = decode_channel(ctx, gb, ch, len, mode[ch], bits[ch])) < 0) yading@10: return ret; yading@10: if (ctx->filter_params > 1 && ctx->filter_params != FILTER_RAW) { yading@10: ctx->filter_bits += 3; yading@10: apply_lpc(ctx, ch, len, bits[ch]); yading@10: } yading@10: if (get_bits_left(gb) < 0) yading@10: return AVERROR_INVALIDDATA; yading@10: } yading@10: ch0 = ctx->channel_data[0]; yading@10: ch1 = ctx->channel_data[1]; yading@10: switch (dmode) { yading@10: case 0: yading@10: for (i = 0; i < len; i++) yading@10: dst0[i] = ch0[i] + ctx->bias[0]; yading@10: break; yading@10: case 1: yading@10: for (i = 0; i < len; i++) { yading@10: dst0[i] = ch0[i] + ctx->bias[0]; yading@10: dst1[i] = ch1[i] + ctx->bias[1]; yading@10: } yading@10: break; yading@10: case 2: yading@10: for (i = 0; i < len; i++) { yading@10: ch0[i] += ctx->bias[0]; yading@10: dst0[i] = ch0[i]; yading@10: dst1[i] = ch0[i] - (ch1[i] + ctx->bias[1]); yading@10: } yading@10: break; yading@10: case 3: yading@10: for (i = 0; i < len; i++) { yading@10: t = ch0[i] + ctx->bias[0]; yading@10: t2 = ch1[i] + ctx->bias[1]; yading@10: dst0[i] = t + t2; yading@10: dst1[i] = t; yading@10: } yading@10: break; yading@10: case 4: yading@10: for (i = 0; i < len; i++) { yading@10: t = ch1[i] + ctx->bias[1]; yading@10: t2 = ((ch0[i] + ctx->bias[0]) << 1) | (t & 1); yading@10: dst0[i] = (t2 + t) / 2; yading@10: dst1[i] = (t2 - t) / 2; yading@10: } yading@10: break; yading@10: } yading@10: yading@10: ctx->sample_offset += len; yading@10: yading@10: return 0; yading@10: } yading@10: yading@10: static int decode_frame(AVCodecContext *avctx, void *data, int *got_frame_ptr, yading@10: AVPacket *avpkt) yading@10: { yading@10: RALFContext *ctx = avctx->priv_data; yading@10: AVFrame *frame = data; yading@10: int16_t *samples0; yading@10: int16_t *samples1; yading@10: int ret; yading@10: GetBitContext gb; yading@10: int table_size, table_bytes, i; yading@10: const uint8_t *src, *block_pointer; yading@10: int src_size; yading@10: int bytes_left; yading@10: yading@10: if (ctx->has_pkt) { yading@10: ctx->has_pkt = 0; yading@10: table_bytes = (AV_RB16(avpkt->data) + 7) >> 3; yading@10: if (table_bytes + 3 > avpkt->size || avpkt->size > RALF_MAX_PKT_SIZE) { yading@10: av_log(avctx, AV_LOG_ERROR, "Wrong packet's breath smells of wrong data!\n"); yading@10: return AVERROR_INVALIDDATA; yading@10: } yading@10: if (memcmp(ctx->pkt, avpkt->data, 2 + table_bytes)) { yading@10: av_log(avctx, AV_LOG_ERROR, "Wrong packet tails are wrong!\n"); yading@10: return AVERROR_INVALIDDATA; yading@10: } yading@10: yading@10: src = ctx->pkt; yading@10: src_size = RALF_MAX_PKT_SIZE + avpkt->size; yading@10: memcpy(ctx->pkt + RALF_MAX_PKT_SIZE, avpkt->data + 2 + table_bytes, yading@10: avpkt->size - 2 - table_bytes); yading@10: } else { yading@10: if (avpkt->size == RALF_MAX_PKT_SIZE) { yading@10: memcpy(ctx->pkt, avpkt->data, avpkt->size); yading@10: ctx->has_pkt = 1; yading@10: *got_frame_ptr = 0; yading@10: yading@10: return avpkt->size; yading@10: } yading@10: src = avpkt->data; yading@10: src_size = avpkt->size; yading@10: } yading@10: yading@10: frame->nb_samples = ctx->max_frame_size; yading@10: if ((ret = ff_get_buffer(avctx, frame, 0)) < 0) yading@10: return ret; yading@10: samples0 = (int16_t *)frame->data[0]; yading@10: samples1 = (int16_t *)frame->data[1]; yading@10: yading@10: if (src_size < 5) { yading@10: av_log(avctx, AV_LOG_ERROR, "too short packets are too short!\n"); yading@10: return AVERROR_INVALIDDATA; yading@10: } yading@10: table_size = AV_RB16(src); yading@10: table_bytes = (table_size + 7) >> 3; yading@10: if (src_size < table_bytes + 3) { yading@10: av_log(avctx, AV_LOG_ERROR, "short packets are short!\n"); yading@10: return AVERROR_INVALIDDATA; yading@10: } yading@10: init_get_bits(&gb, src + 2, table_size); yading@10: ctx->num_blocks = 0; yading@10: while (get_bits_left(&gb) > 0) { yading@10: ctx->block_size[ctx->num_blocks] = get_bits(&gb, 15); yading@10: if (get_bits1(&gb)) { yading@10: ctx->block_pts[ctx->num_blocks] = get_bits(&gb, 9); yading@10: } else { yading@10: ctx->block_pts[ctx->num_blocks] = 0; yading@10: } yading@10: ctx->num_blocks++; yading@10: } yading@10: yading@10: block_pointer = src + table_bytes + 2; yading@10: bytes_left = src_size - table_bytes - 2; yading@10: ctx->sample_offset = 0; yading@10: for (i = 0; i < ctx->num_blocks; i++) { yading@10: if (bytes_left < ctx->block_size[i]) { yading@10: av_log(avctx, AV_LOG_ERROR, "I'm pedaling backwards\n"); yading@10: break; yading@10: } yading@10: init_get_bits(&gb, block_pointer, ctx->block_size[i] * 8); yading@10: if (decode_block(avctx, &gb, samples0 + ctx->sample_offset, yading@10: samples1 + ctx->sample_offset) < 0) { yading@10: av_log(avctx, AV_LOG_ERROR, "Sir, I got carsick in your office. Not decoding the rest of packet.\n"); yading@10: break; yading@10: } yading@10: block_pointer += ctx->block_size[i]; yading@10: bytes_left -= ctx->block_size[i]; yading@10: } yading@10: yading@10: frame->nb_samples = ctx->sample_offset; yading@10: *got_frame_ptr = ctx->sample_offset > 0; yading@10: yading@10: return avpkt->size; yading@10: } yading@10: yading@10: static void decode_flush(AVCodecContext *avctx) yading@10: { yading@10: RALFContext *ctx = avctx->priv_data; yading@10: yading@10: ctx->has_pkt = 0; yading@10: } yading@10: yading@10: yading@10: AVCodec ff_ralf_decoder = { yading@10: .name = "ralf", yading@10: .type = AVMEDIA_TYPE_AUDIO, yading@10: .id = AV_CODEC_ID_RALF, yading@10: .priv_data_size = sizeof(RALFContext), yading@10: .init = decode_init, yading@10: .close = decode_close, yading@10: .decode = decode_frame, yading@10: .flush = decode_flush, yading@10: .capabilities = CODEC_CAP_DR1, yading@10: .long_name = NULL_IF_CONFIG_SMALL("RealAudio Lossless"), yading@10: .sample_fmts = (const enum AVSampleFormat[]) { AV_SAMPLE_FMT_S16P, yading@10: AV_SAMPLE_FMT_NONE }, yading@10: };