ra144: use scalarproduct_int16
[ffmpeg.git] / libavcodec / ra144dec.c
1 /*
2  * Real Audio 1.0 (14.4K)
3  *
4  * Copyright (c) 2008 Vitor Sessak
5  * Copyright (c) 2003 Nick Kurshev
6  *     Based on public domain decoder at http://www.honeypot.net/audio
7  *
8  * This file is part of FFmpeg.
9  *
10  * FFmpeg is free software; you can redistribute it and/or
11  * modify it under the terms of the GNU Lesser General Public
12  * License as published by the Free Software Foundation; either
13  * version 2.1 of the License, or (at your option) any later version.
14  *
15  * FFmpeg is distributed in the hope that it will be useful,
16  * but WITHOUT ANY WARRANTY; without even the implied warranty of
17  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
18  * Lesser General Public License for more details.
19  *
20  * You should have received a copy of the GNU Lesser General Public
21  * License along with FFmpeg; if not, write to the Free Software
22  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
23  */
24
25 #include "libavutil/channel_layout.h"
26 #include "avcodec.h"
27 #include "get_bits.h"
28 #include "internal.h"
29 #include "ra144.h"
30
31
32 static av_cold int ra144_decode_init(AVCodecContext * avctx)
33 {
34     RA144Context *ractx = avctx->priv_data;
35
36     ractx->avctx = avctx;
37     ff_dsputil_init(&ractx->dsp, avctx);
38
39     ractx->lpc_coef[0] = ractx->lpc_tables[0];
40     ractx->lpc_coef[1] = ractx->lpc_tables[1];
41
42     AV_ZERO128(ractx->buffer_a+BLOCKSIZE);
43
44     avctx->channels       = 1;
45     avctx->channel_layout = AV_CH_LAYOUT_MONO;
46     avctx->sample_fmt     = AV_SAMPLE_FMT_S16;
47
48     return 0;
49 }
50
51 static void do_output_subblock(RA144Context *ractx, const int16_t  *lpc_coefs,
52                                int gval, GetBitContext *gb)
53 {
54     int cba_idx = get_bits(gb, 7); // index of the adaptive CB, 0 if none
55     int gain    = get_bits(gb, 8);
56     int cb1_idx = get_bits(gb, 7);
57     int cb2_idx = get_bits(gb, 7);
58
59     ff_subblock_synthesis(ractx, lpc_coefs, cba_idx, cb1_idx, cb2_idx, gval,
60                           gain);
61 }
62
63 /** Uncompress one block (20 bytes -> 160*2 bytes). */
64 static int ra144_decode_frame(AVCodecContext * avctx, void *data,
65                               int *got_frame_ptr, AVPacket *avpkt)
66 {
67     AVFrame *frame     = data;
68     const uint8_t *buf = avpkt->data;
69     int buf_size = avpkt->size;
70     static const uint8_t sizes[LPC_ORDER] = {6, 5, 5, 4, 4, 3, 3, 3, 3, 2};
71     unsigned int refl_rms[NBLOCKS];           // RMS of the reflection coefficients
72     int16_t block_coefs[NBLOCKS][LPC_ORDER];  // LPC coefficients of each sub-block
73     unsigned int lpc_refl[LPC_ORDER];         // LPC reflection coefficients of the frame
74     int i, j;
75     int ret;
76     int16_t *samples;
77     unsigned int energy;
78
79     RA144Context *ractx = avctx->priv_data;
80     GetBitContext gb;
81
82     if (buf_size < FRAME_SIZE) {
83         av_log(avctx, AV_LOG_ERROR,
84                "Frame too small (%d bytes). Truncated file?\n", buf_size);
85         *got_frame_ptr = 0;
86         return AVERROR_INVALIDDATA;
87     }
88
89     /* get output buffer */
90     frame->nb_samples = NBLOCKS * BLOCKSIZE;
91     if ((ret = ff_get_buffer(avctx, frame, 0)) < 0)
92         return ret;
93     samples = (int16_t *)frame->data[0];
94
95     init_get_bits8(&gb, buf, FRAME_SIZE);
96
97     for (i = 0; i < LPC_ORDER; i++)
98         lpc_refl[i] = ff_lpc_refl_cb[i][get_bits(&gb, sizes[i])];
99
100     ff_eval_coefs(ractx->lpc_coef[0], lpc_refl);
101     ractx->lpc_refl_rms[0] = ff_rms(lpc_refl);
102
103     energy = ff_energy_tab[get_bits(&gb, 5)];
104
105     refl_rms[0] = ff_interp(ractx, block_coefs[0], 1, 1, ractx->old_energy);
106     refl_rms[1] = ff_interp(ractx, block_coefs[1], 2,
107                             energy <= ractx->old_energy,
108                             ff_t_sqrt(energy*ractx->old_energy) >> 12);
109     refl_rms[2] = ff_interp(ractx, block_coefs[2], 3, 0, energy);
110     refl_rms[3] = ff_rescale_rms(ractx->lpc_refl_rms[0], energy);
111
112     ff_int_to_int16(block_coefs[3], ractx->lpc_coef[0]);
113
114     for (i=0; i < NBLOCKS; i++) {
115         do_output_subblock(ractx, block_coefs[i], refl_rms[i], &gb);
116
117         for (j=0; j < BLOCKSIZE; j++)
118             *samples++ = av_clip_int16(ractx->curr_sblock[j + 10] << 2);
119     }
120
121     ractx->old_energy = energy;
122     ractx->lpc_refl_rms[1] = ractx->lpc_refl_rms[0];
123
124     FFSWAP(unsigned int *, ractx->lpc_coef[0], ractx->lpc_coef[1]);
125
126     *got_frame_ptr = 1;
127
128     return FRAME_SIZE;
129 }
130
131 AVCodec ff_ra_144_decoder = {
132     .name           = "real_144",
133     .long_name      = NULL_IF_CONFIG_SMALL("RealAudio 1.0 (14.4K)"),
134     .type           = AVMEDIA_TYPE_AUDIO,
135     .id             = AV_CODEC_ID_RA_144,
136     .priv_data_size = sizeof(RA144Context),
137     .init           = ra144_decode_init,
138     .decode         = ra144_decode_frame,
139     .capabilities   = CODEC_CAP_DR1,
140 };