Mercurial > libavcodec.hg
annotate cook.c @ 5345:af1c85da8982 libavcodec
fixpoint: separate windowing arithmetic imlt_window_float
author | mhoffman |
---|---|
date | Mon, 16 Jul 2007 11:52:46 +0000 |
parents | b0d49997bd92 |
children | b41036edcf2e |
rev | line source |
---|---|
2956 | 1 /* |
2 * COOK compatible decoder | |
3 * Copyright (c) 2003 Sascha Sommer | |
4 * Copyright (c) 2005 Benjamin Larsson | |
5 * | |
3947
c8c591fe26f8
Change license headers to say 'FFmpeg' instead of 'this program/this library'
diego
parents:
3619
diff
changeset
|
6 * This file is part of FFmpeg. |
c8c591fe26f8
Change license headers to say 'FFmpeg' instead of 'this program/this library'
diego
parents:
3619
diff
changeset
|
7 * |
c8c591fe26f8
Change license headers to say 'FFmpeg' instead of 'this program/this library'
diego
parents:
3619
diff
changeset
|
8 * FFmpeg is free software; you can redistribute it and/or |
2956 | 9 * modify it under the terms of the GNU Lesser General Public |
10 * License as published by the Free Software Foundation; either | |
3947
c8c591fe26f8
Change license headers to say 'FFmpeg' instead of 'this program/this library'
diego
parents:
3619
diff
changeset
|
11 * version 2.1 of the License, or (at your option) any later version. |
2956 | 12 * |
3947
c8c591fe26f8
Change license headers to say 'FFmpeg' instead of 'this program/this library'
diego
parents:
3619
diff
changeset
|
13 * FFmpeg is distributed in the hope that it will be useful, |
2956 | 14 * but WITHOUT ANY WARRANTY; without even the implied warranty of |
15 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU | |
16 * Lesser General Public License for more details. | |
17 * | |
18 * You should have received a copy of the GNU Lesser General Public | |
3947
c8c591fe26f8
Change license headers to say 'FFmpeg' instead of 'this program/this library'
diego
parents:
3619
diff
changeset
|
19 * License along with FFmpeg; if not, write to the Free Software |
3036
0b546eab515d
Update licensing information: The FSF changed postal address.
diego
parents:
3019
diff
changeset
|
20 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA |
2956 | 21 */ |
22 | |
23 /** | |
24 * @file cook.c | |
4845
78bb9129231b
As usual Real actually took something existing and rebranded it.
banan
parents:
4692
diff
changeset
|
25 * Cook compatible decoder. Bastardization of the G.722.1 standard. |
2956 | 26 * This decoder handles RealNetworks, RealAudio G2 data. |
27 * Cook is identified by the codec name cook in RM files. | |
28 * | |
29 * To use this decoder, a calling application must supply the extradata | |
30 * bytes provided from the RM container; 8+ bytes for mono streams and | |
31 * 16+ for stereo streams (maybe more). | |
32 * | |
33 * Codec technicalities (all this assume a buffer length of 1024): | |
34 * Cook works with several different techniques to achieve its compression. | |
35 * In the timedomain the buffer is divided into 8 pieces and quantized. If | |
36 * two neighboring pieces have different quantization index a smooth | |
37 * quantization curve is used to get a smooth overlap between the different | |
38 * pieces. | |
39 * To get to the transformdomain Cook uses a modulated lapped transform. | |
40 * The transform domain has 50 subbands with 20 elements each. This | |
41 * means only a maximum of 50*20=1000 coefficients are used out of the 1024 | |
42 * available. | |
43 */ | |
44 | |
45 #include <math.h> | |
46 #include <stddef.h> | |
47 #include <stdio.h> | |
48 | |
49 #include "avcodec.h" | |
50 #include "bitstream.h" | |
51 #include "dsputil.h" | |
4430
407e7fd9b4f4
Get rid of the COOKextradata struct. And use valid C to parse the extradata.
banan
parents:
4429
diff
changeset
|
52 #include "bytestream.h" |
4692
3318e3f6470f
Small simplifications to subband coefficient handling and use av_random().
banan
parents:
4674
diff
changeset
|
53 #include "random.h" |
2956 | 54 |
55 #include "cookdata.h" | |
56 | |
57 /* the different Cook versions */ | |
4425 | 58 #define MONO 0x1000001 |
59 #define STEREO 0x1000002 | |
2956 | 60 #define JOINT_STEREO 0x1000003 |
61 #define MC_COOK 0x2000000 //multichannel Cook, not supported | |
62 | |
63 #define SUBBAND_SIZE 20 | |
64 //#define COOKDEBUG | |
65 | |
66 typedef struct { | |
4639 | 67 int *now; |
68 int *previous; | |
69 } cook_gains; | |
2956 | 70 |
71 typedef struct { | |
72 GetBitContext gb; | |
73 /* stream data */ | |
74 int nb_channels; | |
75 int joint_stereo; | |
76 int bit_rate; | |
77 int sample_rate; | |
78 int samples_per_channel; | |
79 int samples_per_frame; | |
80 int subbands; | |
3090 | 81 int log2_numvector_size; |
82 int numvector_size; //1 << log2_numvector_size; | |
2956 | 83 int js_subband_start; |
84 int total_subbands; | |
85 int num_vectors; | |
86 int bits_per_subpacket; | |
4430
407e7fd9b4f4
Get rid of the COOKextradata struct. And use valid C to parse the extradata.
banan
parents:
4429
diff
changeset
|
87 int cookversion; |
2956 | 88 /* states */ |
4692
3318e3f6470f
Small simplifications to subband coefficient handling and use av_random().
banan
parents:
4674
diff
changeset
|
89 AVRandomState random_state; |
2956 | 90 |
91 /* transform data */ | |
4650
31bf54d9353d
Replace custom modified discrete cosine transform with ffmpeg's own.
banan
parents:
4639
diff
changeset
|
92 MDCTContext mdct_ctx; |
31bf54d9353d
Replace custom modified discrete cosine transform with ffmpeg's own.
banan
parents:
4639
diff
changeset
|
93 DECLARE_ALIGNED_16(FFTSample, mdct_tmp[1024]); /* temporary storage for imlt */ |
2956 | 94 float* mlt_window; |
95 | |
96 /* gain buffers */ | |
4639 | 97 cook_gains gains1; |
98 cook_gains gains2; | |
99 int gain_1[9]; | |
100 int gain_2[9]; | |
101 int gain_3[9]; | |
102 int gain_4[9]; | |
2956 | 103 |
104 /* VLC data */ | |
105 int js_vlc_bits; | |
106 VLC envelope_quant_index[13]; | |
107 VLC sqvh[7]; //scalar quantization | |
108 VLC ccpl; //channel coupling | |
109 | |
110 /* generatable tables and related variables */ | |
111 int gain_size_factor; | |
112 float gain_table[23]; | |
113 float pow2tab[127]; | |
114 float rootpow2tab[127]; | |
115 | |
116 /* data buffers */ | |
117 | |
118 uint8_t* decoded_bytes_buffer; | |
4542 | 119 DECLARE_ALIGNED_16(float,mono_mdct_output[2048]); |
2956 | 120 float mono_previous_buffer1[1024]; |
121 float mono_previous_buffer2[1024]; | |
122 float decode_buffer_1[1024]; | |
123 float decode_buffer_2[1024]; | |
5342
887509b42d3f
moving automatic allocation of joint_decode/decode_buffer
mhoffman
parents:
5341
diff
changeset
|
124 float decode_buffer_0[1060]; /* static allocation for joint decode */ |
2956 | 125 } COOKContext; |
126 | |
127 /* debug functions */ | |
128 | |
129 #ifdef COOKDEBUG | |
130 static void dump_float_table(float* table, int size, int delimiter) { | |
131 int i=0; | |
132 av_log(NULL,AV_LOG_ERROR,"\n[%d]: ",i); | |
133 for (i=0 ; i<size ; i++) { | |
134 av_log(NULL, AV_LOG_ERROR, "%5.1f, ", table[i]); | |
135 if ((i+1)%delimiter == 0) av_log(NULL,AV_LOG_ERROR,"\n[%d]: ",i+1); | |
136 } | |
137 } | |
138 | |
139 static void dump_int_table(int* table, int size, int delimiter) { | |
140 int i=0; | |
141 av_log(NULL,AV_LOG_ERROR,"\n[%d]: ",i); | |
142 for (i=0 ; i<size ; i++) { | |
143 av_log(NULL, AV_LOG_ERROR, "%d, ", table[i]); | |
144 if ((i+1)%delimiter == 0) av_log(NULL,AV_LOG_ERROR,"\n[%d]: ",i+1); | |
145 } | |
146 } | |
147 | |
148 static void dump_short_table(short* table, int size, int delimiter) { | |
149 int i=0; | |
150 av_log(NULL,AV_LOG_ERROR,"\n[%d]: ",i); | |
151 for (i=0 ; i<size ; i++) { | |
152 av_log(NULL, AV_LOG_ERROR, "%d, ", table[i]); | |
153 if ((i+1)%delimiter == 0) av_log(NULL,AV_LOG_ERROR,"\n[%d]: ",i+1); | |
154 } | |
155 } | |
156 | |
157 #endif | |
158 | |
159 /*************** init functions ***************/ | |
160 | |
161 /* table generator */ | |
162 static void init_pow2table(COOKContext *q){ | |
163 int i; | |
164 q->pow2tab[63] = 1.0; | |
165 for (i=1 ; i<64 ; i++){ | |
3106
9cbd63cca826
Don't use pow/powf functions where we just need integer arithmetic.
al
parents:
3091
diff
changeset
|
166 q->pow2tab[63+i]=(float)((uint64_t)1<<i); |
9cbd63cca826
Don't use pow/powf functions where we just need integer arithmetic.
al
parents:
3091
diff
changeset
|
167 q->pow2tab[63-i]=1.0/(float)((uint64_t)1<<i); |
2956 | 168 } |
169 } | |
170 | |
171 /* table generator */ | |
172 static void init_rootpow2table(COOKContext *q){ | |
173 int i; | |
174 q->rootpow2tab[63] = 1.0; | |
175 for (i=1 ; i<64 ; i++){ | |
3106
9cbd63cca826
Don't use pow/powf functions where we just need integer arithmetic.
al
parents:
3091
diff
changeset
|
176 q->rootpow2tab[63+i]=sqrt((float)((uint64_t)1<<i)); |
9cbd63cca826
Don't use pow/powf functions where we just need integer arithmetic.
al
parents:
3091
diff
changeset
|
177 q->rootpow2tab[63-i]=sqrt(1.0/(float)((uint64_t)1<<i)); |
2956 | 178 } |
179 } | |
180 | |
181 /* table generator */ | |
182 static void init_gain_table(COOKContext *q) { | |
183 int i; | |
184 q->gain_size_factor = q->samples_per_channel/8; | |
185 for (i=0 ; i<23 ; i++) { | |
186 q->gain_table[i] = pow((double)q->pow2tab[i+52] , | |
187 (1.0/(double)q->gain_size_factor)); | |
188 } | |
189 } | |
190 | |
191 | |
192 static int init_cook_vlc_tables(COOKContext *q) { | |
193 int i, result; | |
194 | |
195 result = 0; | |
196 for (i=0 ; i<13 ; i++) { | |
4951 | 197 result |= init_vlc (&q->envelope_quant_index[i], 9, 24, |
2956 | 198 envelope_quant_index_huffbits[i], 1, 1, |
199 envelope_quant_index_huffcodes[i], 2, 2, 0); | |
200 } | |
201 av_log(NULL,AV_LOG_DEBUG,"sqvh VLC init\n"); | |
202 for (i=0 ; i<7 ; i++) { | |
4951 | 203 result |= init_vlc (&q->sqvh[i], vhvlcsize_tab[i], vhsize_tab[i], |
2956 | 204 cvh_huffbits[i], 1, 1, |
205 cvh_huffcodes[i], 2, 2, 0); | |
206 } | |
207 | |
208 if (q->nb_channels==2 && q->joint_stereo==1){ | |
4951 | 209 result |= init_vlc (&q->ccpl, 6, (1<<q->js_vlc_bits)-1, |
2956 | 210 ccpl_huffbits[q->js_vlc_bits-2], 1, 1, |
211 ccpl_huffcodes[q->js_vlc_bits-2], 2, 2, 0); | |
212 av_log(NULL,AV_LOG_DEBUG,"Joint-stereo VLC used.\n"); | |
213 } | |
214 | |
215 av_log(NULL,AV_LOG_DEBUG,"VLC tables initialized.\n"); | |
216 return result; | |
217 } | |
218 | |
219 static int init_cook_mlt(COOKContext *q) { | |
220 int j; | |
221 float alpha; | |
4650
31bf54d9353d
Replace custom modified discrete cosine transform with ffmpeg's own.
banan
parents:
4639
diff
changeset
|
222 int mlt_size = q->samples_per_channel; |
2956 | 223 |
4650
31bf54d9353d
Replace custom modified discrete cosine transform with ffmpeg's own.
banan
parents:
4639
diff
changeset
|
224 if ((q->mlt_window = av_malloc(sizeof(float)*mlt_size)) == 0) |
31bf54d9353d
Replace custom modified discrete cosine transform with ffmpeg's own.
banan
parents:
4639
diff
changeset
|
225 return -1; |
2956 | 226 |
227 /* Initialize the MLT window: simple sine window. */ | |
4650
31bf54d9353d
Replace custom modified discrete cosine transform with ffmpeg's own.
banan
parents:
4639
diff
changeset
|
228 alpha = M_PI / (2.0 * (float)mlt_size); |
31bf54d9353d
Replace custom modified discrete cosine transform with ffmpeg's own.
banan
parents:
4639
diff
changeset
|
229 for(j=0 ; j<mlt_size ; j++) |
31bf54d9353d
Replace custom modified discrete cosine transform with ffmpeg's own.
banan
parents:
4639
diff
changeset
|
230 q->mlt_window[j] = sin((j + 0.5) * alpha) * sqrt(2.0 / q->samples_per_channel); |
2956 | 231 |
4650
31bf54d9353d
Replace custom modified discrete cosine transform with ffmpeg's own.
banan
parents:
4639
diff
changeset
|
232 /* Initialize the MDCT. */ |
31bf54d9353d
Replace custom modified discrete cosine transform with ffmpeg's own.
banan
parents:
4639
diff
changeset
|
233 if (ff_mdct_init(&q->mdct_ctx, av_log2(mlt_size)+1, 1)) { |
31bf54d9353d
Replace custom modified discrete cosine transform with ffmpeg's own.
banan
parents:
4639
diff
changeset
|
234 av_free(q->mlt_window); |
31bf54d9353d
Replace custom modified discrete cosine transform with ffmpeg's own.
banan
parents:
4639
diff
changeset
|
235 return -1; |
2956 | 236 } |
4650
31bf54d9353d
Replace custom modified discrete cosine transform with ffmpeg's own.
banan
parents:
4639
diff
changeset
|
237 av_log(NULL,AV_LOG_DEBUG,"MDCT initialized, order = %d.\n", |
31bf54d9353d
Replace custom modified discrete cosine transform with ffmpeg's own.
banan
parents:
4639
diff
changeset
|
238 av_log2(mlt_size)+1); |
2956 | 239 |
4650
31bf54d9353d
Replace custom modified discrete cosine transform with ffmpeg's own.
banan
parents:
4639
diff
changeset
|
240 return 0; |
2956 | 241 } |
242 | |
243 /*************** init functions end ***********/ | |
244 | |
245 /** | |
246 * Cook indata decoding, every 32 bits are XORed with 0x37c511f2. | |
247 * Why? No idea, some checksum/error detection method maybe. | |
4424
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
248 * |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
249 * Out buffer size: extra bytes are needed to cope with |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
250 * padding/missalignment. |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
251 * Subpackets passed to the decoder can contain two, consecutive |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
252 * half-subpackets, of identical but arbitrary size. |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
253 * 1234 1234 1234 1234 extraA extraB |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
254 * Case 1: AAAA BBBB 0 0 |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
255 * Case 2: AAAA ABBB BB-- 3 3 |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
256 * Case 3: AAAA AABB BBBB 2 2 |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
257 * Case 4: AAAA AAAB BBBB BB-- 1 5 |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
258 * |
2956 | 259 * Nice way to waste CPU cycles. |
260 * | |
4424
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
261 * @param inbuffer pointer to byte array of indata |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
262 * @param out pointer to byte array of outdata |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
263 * @param bytes number of bytes |
2956 | 264 */ |
4424
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
265 #define DECODE_BYTES_PAD1(bytes) (3 - ((bytes)+3) % 4) |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
266 #define DECODE_BYTES_PAD2(bytes) ((bytes) % 4 + DECODE_BYTES_PAD1(2 * (bytes))) |
2956 | 267 |
4424
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
268 static inline int decode_bytes(uint8_t* inbuffer, uint8_t* out, int bytes){ |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
269 int i, off; |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
270 uint32_t c; |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
271 uint32_t* buf; |
2956 | 272 uint32_t* obuf = (uint32_t*) out; |
273 /* FIXME: 64 bit platforms would be able to do 64 bits at a time. | |
274 * I'm too lazy though, should be something like | |
275 * for(i=0 ; i<bitamount/64 ; i++) | |
276 * (int64_t)out[i] = 0x37c511f237c511f2^be2me_64(int64_t)in[i]); | |
277 * Buffer alignment needs to be checked. */ | |
278 | |
4429 | 279 off = (int)((long)inbuffer & 3); |
4424
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
280 buf = (uint32_t*) (inbuffer - off); |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
281 c = be2me_32((0x37c511f2 >> (off*8)) | (0x37c511f2 << (32-(off*8)))); |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
282 bytes += 3 + off; |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
283 for (i = 0; i < bytes/4; i++) |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
284 obuf[i] = c ^ buf[i]; |
2956 | 285 |
4424
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
286 return off; |
2956 | 287 } |
288 | |
289 /** | |
290 * Cook uninit | |
291 */ | |
292 | |
293 static int cook_decode_close(AVCodecContext *avctx) | |
294 { | |
295 int i; | |
296 COOKContext *q = avctx->priv_data; | |
4302
a0f83004d485
av_log(NULL,... -> av_log(avctx,.. where appropriate.
banan
parents:
3947
diff
changeset
|
297 av_log(avctx,AV_LOG_DEBUG, "Deallocating memory.\n"); |
2956 | 298 |
299 /* Free allocated memory buffers. */ | |
300 av_free(q->mlt_window); | |
301 av_free(q->decoded_bytes_buffer); | |
302 | |
303 /* Free the transform. */ | |
4650
31bf54d9353d
Replace custom modified discrete cosine transform with ffmpeg's own.
banan
parents:
4639
diff
changeset
|
304 ff_mdct_end(&q->mdct_ctx); |
2956 | 305 |
306 /* Free the VLC tables. */ | |
307 for (i=0 ; i<13 ; i++) { | |
308 free_vlc(&q->envelope_quant_index[i]); | |
309 } | |
310 for (i=0 ; i<7 ; i++) { | |
311 free_vlc(&q->sqvh[i]); | |
312 } | |
313 if(q->nb_channels==2 && q->joint_stereo==1 ){ | |
314 free_vlc(&q->ccpl); | |
315 } | |
316 | |
317 av_log(NULL,AV_LOG_DEBUG,"Memory deallocated.\n"); | |
318 | |
319 return 0; | |
320 } | |
321 | |
322 /** | |
4639 | 323 * Fill the gain array for the timedomain quantization. |
2956 | 324 * |
325 * @param q pointer to the COOKContext | |
4639 | 326 * @param gaininfo[9] array of gain indices |
2956 | 327 */ |
328 | |
4639 | 329 static void decode_gain_info(GetBitContext *gb, int *gaininfo) |
330 { | |
331 int i, n; | |
2956 | 332 |
333 while (get_bits1(gb)) {} | |
4639 | 334 n = get_bits_count(gb) - 1; //amount of elements*2 to update |
2956 | 335 |
4639 | 336 i = 0; |
337 while (n--) { | |
338 int index = get_bits(gb, 3); | |
339 int gain = get_bits1(gb) ? get_bits(gb, 4) - 7 : -1; | |
340 | |
341 while (i <= index) gaininfo[i++] = gain; | |
2956 | 342 } |
4639 | 343 while (i <= 8) gaininfo[i++] = 0; |
2956 | 344 } |
345 | |
346 /** | |
347 * Create the quant index table needed for the envelope. | |
348 * | |
349 * @param q pointer to the COOKContext | |
350 * @param quant_index_table pointer to the array | |
351 */ | |
352 | |
353 static void decode_envelope(COOKContext *q, int* quant_index_table) { | |
354 int i,j, vlc_index; | |
355 | |
356 quant_index_table[0]= get_bits(&q->gb,6) - 6; //This is used later in categorize | |
357 | |
358 for (i=1 ; i < q->total_subbands ; i++){ | |
359 vlc_index=i; | |
360 if (i >= q->js_subband_start * 2) { | |
361 vlc_index-=q->js_subband_start; | |
362 } else { | |
363 vlc_index/=2; | |
364 if(vlc_index < 1) vlc_index = 1; | |
365 } | |
366 if (vlc_index>13) vlc_index = 13; //the VLC tables >13 are identical to No. 13 | |
367 | |
368 j = get_vlc2(&q->gb, q->envelope_quant_index[vlc_index-1].table, | |
369 q->envelope_quant_index[vlc_index-1].bits,2); | |
370 quant_index_table[i] = quant_index_table[i-1] + j - 12; //differential encoding | |
371 } | |
372 } | |
373 | |
374 /** | |
375 * Calculate the category and category_index vector. | |
376 * | |
377 * @param q pointer to the COOKContext | |
378 * @param quant_index_table pointer to the array | |
379 * @param category pointer to the category array | |
380 * @param category_index pointer to the category_index array | |
381 */ | |
382 | |
383 static void categorize(COOKContext *q, int* quant_index_table, | |
384 int* category, int* category_index){ | |
4952
1900e2eaecda
Add another tmpbias variable, as bias' value will be used later
ramiro
parents:
4951
diff
changeset
|
385 int exp_idx, bias, tmpbias1, tmpbias2, bits_left, num_bits, index, v, i, j; |
2956 | 386 int exp_index2[102]; |
387 int exp_index1[102]; | |
388 | |
4955
bbe763044678
Use 1 array with double the size instead of 2 arrays with normal size
ramiro
parents:
4954
diff
changeset
|
389 int tmp_categorize_array[128*2]; |
bbe763044678
Use 1 array with double the size instead of 2 arrays with normal size
ramiro
parents:
4954
diff
changeset
|
390 int tmp_categorize_array1_idx=q->numvector_size; |
bbe763044678
Use 1 array with double the size instead of 2 arrays with normal size
ramiro
parents:
4954
diff
changeset
|
391 int tmp_categorize_array2_idx=q->numvector_size; |
2956 | 392 |
393 bits_left = q->bits_per_subpacket - get_bits_count(&q->gb); | |
394 | |
395 if(bits_left > q->samples_per_channel) { | |
396 bits_left = q->samples_per_channel + | |
397 ((bits_left - q->samples_per_channel)*5)/8; | |
398 //av_log(NULL, AV_LOG_ERROR, "bits_left = %d\n",bits_left); | |
399 } | |
400 | |
401 memset(&exp_index1,0,102*sizeof(int)); | |
402 memset(&exp_index2,0,102*sizeof(int)); | |
4955
bbe763044678
Use 1 array with double the size instead of 2 arrays with normal size
ramiro
parents:
4954
diff
changeset
|
403 memset(&tmp_categorize_array,0,128*2*sizeof(int)); |
2956 | 404 |
405 bias=-32; | |
406 | |
407 /* Estimate bias. */ | |
408 for (i=32 ; i>0 ; i=i/2){ | |
409 num_bits = 0; | |
410 index = 0; | |
411 for (j=q->total_subbands ; j>0 ; j--){ | |
4863 | 412 exp_idx = av_clip((i - quant_index_table[index] + bias) / 2, 0, 7); |
2956 | 413 index++; |
414 num_bits+=expbits_tab[exp_idx]; | |
415 } | |
416 if(num_bits >= bits_left - 32){ | |
417 bias+=i; | |
418 } | |
419 } | |
420 | |
421 /* Calculate total number of bits. */ | |
422 num_bits=0; | |
423 for (i=0 ; i<q->total_subbands ; i++) { | |
4863 | 424 exp_idx = av_clip((bias - quant_index_table[i]) / 2, 0, 7); |
2956 | 425 num_bits += expbits_tab[exp_idx]; |
426 exp_index1[i] = exp_idx; | |
427 exp_index2[i] = exp_idx; | |
428 } | |
4952
1900e2eaecda
Add another tmpbias variable, as bias' value will be used later
ramiro
parents:
4951
diff
changeset
|
429 tmpbias1 = tmpbias2 = num_bits; |
2956 | 430 |
431 for (j = 1 ; j < q->numvector_size ; j++) { | |
4952
1900e2eaecda
Add another tmpbias variable, as bias' value will be used later
ramiro
parents:
4951
diff
changeset
|
432 if (tmpbias1 + tmpbias2 > 2*bits_left) { /* ---> */ |
2956 | 433 int max = -999999; |
434 index=-1; | |
435 for (i=0 ; i<q->total_subbands ; i++){ | |
436 if (exp_index1[i] < 7) { | |
4954
1fd90978db25
Add bias instead of -32 or 0, as is done in g.722.1
ramiro
parents:
4953
diff
changeset
|
437 v = (-2*exp_index1[i]) - quant_index_table[i] + bias; |
2956 | 438 if ( v >= max) { |
439 max = v; | |
440 index = i; | |
441 } | |
442 } | |
443 } | |
444 if(index==-1)break; | |
4955
bbe763044678
Use 1 array with double the size instead of 2 arrays with normal size
ramiro
parents:
4954
diff
changeset
|
445 tmp_categorize_array[tmp_categorize_array1_idx++] = index; |
4952
1900e2eaecda
Add another tmpbias variable, as bias' value will be used later
ramiro
parents:
4951
diff
changeset
|
446 tmpbias1 -= expbits_tab[exp_index1[index]] - |
4953 | 447 expbits_tab[exp_index1[index]+1]; |
2956 | 448 ++exp_index1[index]; |
449 } else { /* <--- */ | |
450 int min = 999999; | |
451 index=-1; | |
452 for (i=0 ; i<q->total_subbands ; i++){ | |
453 if(exp_index2[i] > 0){ | |
4954
1fd90978db25
Add bias instead of -32 or 0, as is done in g.722.1
ramiro
parents:
4953
diff
changeset
|
454 v = (-2*exp_index2[i])-quant_index_table[i]+bias; |
2956 | 455 if ( v < min) { |
456 min = v; | |
457 index = i; | |
458 } | |
459 } | |
460 } | |
461 if(index == -1)break; | |
4955
bbe763044678
Use 1 array with double the size instead of 2 arrays with normal size
ramiro
parents:
4954
diff
changeset
|
462 tmp_categorize_array[--tmp_categorize_array2_idx] = index; |
4952
1900e2eaecda
Add another tmpbias variable, as bias' value will be used later
ramiro
parents:
4951
diff
changeset
|
463 tmpbias2 -= expbits_tab[exp_index2[index]] - |
4953 | 464 expbits_tab[exp_index2[index]-1]; |
2956 | 465 --exp_index2[index]; |
466 } | |
467 } | |
468 | |
469 for(i=0 ; i<q->total_subbands ; i++) | |
470 category[i] = exp_index2[i]; | |
471 | |
4955
bbe763044678
Use 1 array with double the size instead of 2 arrays with normal size
ramiro
parents:
4954
diff
changeset
|
472 for(i=0 ; i<q->numvector_size-1 ; i++) |
bbe763044678
Use 1 array with double the size instead of 2 arrays with normal size
ramiro
parents:
4954
diff
changeset
|
473 category_index[i] = tmp_categorize_array[tmp_categorize_array2_idx++]; |
2956 | 474 |
475 } | |
476 | |
477 | |
478 /** | |
479 * Expand the category vector. | |
480 * | |
481 * @param q pointer to the COOKContext | |
482 * @param category pointer to the category array | |
483 * @param category_index pointer to the category_index array | |
484 */ | |
485 | |
4908
777f250df232
Fix multiple "¡Æinline/static¡Ç is not at beginning of declaration" warnings.
diego
parents:
4863
diff
changeset
|
486 static inline void expand_category(COOKContext *q, int* category, |
2956 | 487 int* category_index){ |
488 int i; | |
489 for(i=0 ; i<q->num_vectors ; i++){ | |
490 ++category[category_index[i]]; | |
491 } | |
492 } | |
493 | |
494 /** | |
495 * The real requantization of the mltcoefs | |
496 * | |
497 * @param q pointer to the COOKContext | |
498 * @param index index | |
4692
3318e3f6470f
Small simplifications to subband coefficient handling and use av_random().
banan
parents:
4674
diff
changeset
|
499 * @param quant_index quantisation index |
2956 | 500 * @param subband_coef_index array of indexes to quant_centroid_tab |
4674 | 501 * @param subband_coef_sign signs of coefficients |
4692
3318e3f6470f
Small simplifications to subband coefficient handling and use av_random().
banan
parents:
4674
diff
changeset
|
502 * @param mlt_p pointer into the mlt buffer |
2956 | 503 */ |
504 | |
4692
3318e3f6470f
Small simplifications to subband coefficient handling and use av_random().
banan
parents:
4674
diff
changeset
|
505 static void scalar_dequant(COOKContext *q, int index, int quant_index, |
3318e3f6470f
Small simplifications to subband coefficient handling and use av_random().
banan
parents:
4674
diff
changeset
|
506 int* subband_coef_index, int* subband_coef_sign, |
3318e3f6470f
Small simplifications to subband coefficient handling and use av_random().
banan
parents:
4674
diff
changeset
|
507 float* mlt_p){ |
2956 | 508 int i; |
509 float f1; | |
510 | |
511 for(i=0 ; i<SUBBAND_SIZE ; i++) { | |
512 if (subband_coef_index[i]) { | |
4692
3318e3f6470f
Small simplifications to subband coefficient handling and use av_random().
banan
parents:
4674
diff
changeset
|
513 f1 = quant_centroid_tab[index][subband_coef_index[i]]; |
3318e3f6470f
Small simplifications to subband coefficient handling and use av_random().
banan
parents:
4674
diff
changeset
|
514 if (subband_coef_sign[i]) f1 = -f1; |
2956 | 515 } else { |
4674 | 516 /* noise coding if subband_coef_index[i] == 0 */ |
4692
3318e3f6470f
Small simplifications to subband coefficient handling and use av_random().
banan
parents:
4674
diff
changeset
|
517 f1 = dither_tab[index]; |
3318e3f6470f
Small simplifications to subband coefficient handling and use av_random().
banan
parents:
4674
diff
changeset
|
518 if (av_random(&q->random_state) < 0x80000000) f1 = -f1; |
2956 | 519 } |
4692
3318e3f6470f
Small simplifications to subband coefficient handling and use av_random().
banan
parents:
4674
diff
changeset
|
520 mlt_p[i] = f1 * q->rootpow2tab[quant_index+63]; |
2956 | 521 } |
522 } | |
523 /** | |
4674 | 524 * Unpack the subband_coef_index and subband_coef_sign vectors. |
2956 | 525 * |
526 * @param q pointer to the COOKContext | |
527 * @param category pointer to the category array | |
528 * @param subband_coef_index array of indexes to quant_centroid_tab | |
4674 | 529 * @param subband_coef_sign signs of coefficients |
2956 | 530 */ |
531 | |
532 static int unpack_SQVH(COOKContext *q, int category, int* subband_coef_index, | |
4674 | 533 int* subband_coef_sign) { |
2956 | 534 int i,j; |
535 int vlc, vd ,tmp, result; | |
536 | |
537 vd = vd_tab[category]; | |
538 result = 0; | |
539 for(i=0 ; i<vpr_tab[category] ; i++){ | |
540 vlc = get_vlc2(&q->gb, q->sqvh[category].table, q->sqvh[category].bits, 3); | |
541 if (q->bits_per_subpacket < get_bits_count(&q->gb)){ | |
542 vlc = 0; | |
543 result = 1; | |
544 } | |
545 for(j=vd-1 ; j>=0 ; j--){ | |
546 tmp = (vlc * invradix_tab[category])/0x100000; | |
547 subband_coef_index[vd*i+j] = vlc - tmp * (kmax_tab[category]+1); | |
548 vlc = tmp; | |
549 } | |
550 for(j=0 ; j<vd ; j++){ | |
551 if (subband_coef_index[i*vd + j]) { | |
552 if(get_bits_count(&q->gb) < q->bits_per_subpacket){ | |
4674 | 553 subband_coef_sign[i*vd+j] = get_bits1(&q->gb); |
2956 | 554 } else { |
555 result=1; | |
4674 | 556 subband_coef_sign[i*vd+j]=0; |
2956 | 557 } |
558 } else { | |
4674 | 559 subband_coef_sign[i*vd+j]=0; |
2956 | 560 } |
561 } | |
562 } | |
563 return result; | |
564 } | |
565 | |
566 | |
567 /** | |
568 * Fill the mlt_buffer with mlt coefficients. | |
569 * | |
570 * @param q pointer to the COOKContext | |
571 * @param category pointer to the category array | |
4692
3318e3f6470f
Small simplifications to subband coefficient handling and use av_random().
banan
parents:
4674
diff
changeset
|
572 * @param quant_index_table pointer to the array |
2956 | 573 * @param mlt_buffer pointer to mlt coefficients |
574 */ | |
575 | |
576 | |
577 static void decode_vectors(COOKContext* q, int* category, | |
4692
3318e3f6470f
Small simplifications to subband coefficient handling and use av_random().
banan
parents:
4674
diff
changeset
|
578 int *quant_index_table, float* mlt_buffer){ |
2956 | 579 /* A zero in this table means that the subband coefficient is |
580 random noise coded. */ | |
4674 | 581 int subband_coef_index[SUBBAND_SIZE]; |
2956 | 582 /* A zero in this table means that the subband coefficient is a |
583 positive multiplicator. */ | |
4674 | 584 int subband_coef_sign[SUBBAND_SIZE]; |
2956 | 585 int band, j; |
586 int index=0; | |
587 | |
588 for(band=0 ; band<q->total_subbands ; band++){ | |
589 index = category[band]; | |
590 if(category[band] < 7){ | |
4674 | 591 if(unpack_SQVH(q, category[band], subband_coef_index, subband_coef_sign)){ |
2956 | 592 index=7; |
593 for(j=0 ; j<q->total_subbands ; j++) category[band+j]=7; | |
594 } | |
595 } | |
596 if(index==7) { | |
597 memset(subband_coef_index, 0, sizeof(subband_coef_index)); | |
4674 | 598 memset(subband_coef_sign, 0, sizeof(subband_coef_sign)); |
2956 | 599 } |
4692
3318e3f6470f
Small simplifications to subband coefficient handling and use av_random().
banan
parents:
4674
diff
changeset
|
600 scalar_dequant(q, index, quant_index_table[band], |
3318e3f6470f
Small simplifications to subband coefficient handling and use av_random().
banan
parents:
4674
diff
changeset
|
601 subband_coef_index, subband_coef_sign, |
5341 | 602 &mlt_buffer[band * SUBBAND_SIZE]); |
2956 | 603 } |
604 | |
605 if(q->total_subbands*SUBBAND_SIZE >= q->samples_per_channel){ | |
606 return; | |
4674 | 607 } /* FIXME: should this be removed, or moved into loop above? */ |
2956 | 608 } |
609 | |
610 | |
611 /** | |
612 * function for decoding mono data | |
613 * | |
614 * @param q pointer to the COOKContext | |
4860 | 615 * @param mlt_buffer pointer to mlt coefficients |
2956 | 616 */ |
617 | |
618 static void mono_decode(COOKContext *q, float* mlt_buffer) { | |
619 | |
620 int category_index[128]; | |
621 int quant_index_table[102]; | |
622 int category[128]; | |
623 | |
624 memset(&category, 0, 128*sizeof(int)); | |
625 memset(&category_index, 0, 128*sizeof(int)); | |
626 | |
627 decode_envelope(q, quant_index_table); | |
3090 | 628 q->num_vectors = get_bits(&q->gb,q->log2_numvector_size); |
2956 | 629 categorize(q, quant_index_table, category, category_index); |
630 expand_category(q, category, category_index); | |
4692
3318e3f6470f
Small simplifications to subband coefficient handling and use av_random().
banan
parents:
4674
diff
changeset
|
631 decode_vectors(q, category, quant_index_table, mlt_buffer); |
2956 | 632 } |
633 | |
634 | |
635 /** | |
636 * the actual requantization of the timedomain samples | |
637 * | |
638 * @param q pointer to the COOKContext | |
639 * @param buffer pointer to the timedomain buffer | |
640 * @param gain_index index for the block multiplier | |
641 * @param gain_index_next index for the next block multiplier | |
642 */ | |
643 | |
644 static void interpolate(COOKContext *q, float* buffer, | |
645 int gain_index, int gain_index_next){ | |
646 int i; | |
647 float fc1, fc2; | |
648 fc1 = q->pow2tab[gain_index+63]; | |
649 | |
650 if(gain_index == gain_index_next){ //static gain | |
651 for(i=0 ; i<q->gain_size_factor ; i++){ | |
652 buffer[i]*=fc1; | |
653 } | |
654 return; | |
655 } else { //smooth gain | |
656 fc2 = q->gain_table[11 + (gain_index_next-gain_index)]; | |
657 for(i=0 ; i<q->gain_size_factor ; i++){ | |
658 buffer[i]*=fc1; | |
659 fc1*=fc2; | |
660 } | |
661 return; | |
662 } | |
663 } | |
664 | |
5345
af1c85da8982
fixpoint: separate windowing arithmetic imlt_window_float
mhoffman
parents:
5344
diff
changeset
|
665 /** |
af1c85da8982
fixpoint: separate windowing arithmetic imlt_window_float
mhoffman
parents:
5344
diff
changeset
|
666 * Apply transform window, overlap buffers. |
af1c85da8982
fixpoint: separate windowing arithmetic imlt_window_float
mhoffman
parents:
5344
diff
changeset
|
667 * |
af1c85da8982
fixpoint: separate windowing arithmetic imlt_window_float
mhoffman
parents:
5344
diff
changeset
|
668 * @param q pointer to the COOKContext |
af1c85da8982
fixpoint: separate windowing arithmetic imlt_window_float
mhoffman
parents:
5344
diff
changeset
|
669 * @param inbuffer pointer to the mltcoefficients |
af1c85da8982
fixpoint: separate windowing arithmetic imlt_window_float
mhoffman
parents:
5344
diff
changeset
|
670 * @param gains_ptr current and previous gains |
af1c85da8982
fixpoint: separate windowing arithmetic imlt_window_float
mhoffman
parents:
5344
diff
changeset
|
671 * @param previous_buffer pointer to the previous buffer to be used for overlapping |
af1c85da8982
fixpoint: separate windowing arithmetic imlt_window_float
mhoffman
parents:
5344
diff
changeset
|
672 */ |
af1c85da8982
fixpoint: separate windowing arithmetic imlt_window_float
mhoffman
parents:
5344
diff
changeset
|
673 |
af1c85da8982
fixpoint: separate windowing arithmetic imlt_window_float
mhoffman
parents:
5344
diff
changeset
|
674 static void imlt_window_float (COOKContext *q, float *buffer1, |
af1c85da8982
fixpoint: separate windowing arithmetic imlt_window_float
mhoffman
parents:
5344
diff
changeset
|
675 cook_gains *gains_ptr, float *previous_buffer) |
af1c85da8982
fixpoint: separate windowing arithmetic imlt_window_float
mhoffman
parents:
5344
diff
changeset
|
676 { |
af1c85da8982
fixpoint: separate windowing arithmetic imlt_window_float
mhoffman
parents:
5344
diff
changeset
|
677 const float fc = q->pow2tab[gains_ptr->previous[0] + 63]; |
af1c85da8982
fixpoint: separate windowing arithmetic imlt_window_float
mhoffman
parents:
5344
diff
changeset
|
678 int i; |
af1c85da8982
fixpoint: separate windowing arithmetic imlt_window_float
mhoffman
parents:
5344
diff
changeset
|
679 /* The weird thing here, is that the two halves of the time domain |
af1c85da8982
fixpoint: separate windowing arithmetic imlt_window_float
mhoffman
parents:
5344
diff
changeset
|
680 * buffer are swapped. Also, the newest data, that we save away for |
af1c85da8982
fixpoint: separate windowing arithmetic imlt_window_float
mhoffman
parents:
5344
diff
changeset
|
681 * next frame, has the wrong sign. Hence the subtraction below. |
af1c85da8982
fixpoint: separate windowing arithmetic imlt_window_float
mhoffman
parents:
5344
diff
changeset
|
682 * Almost sounds like a complex conjugate/reverse data/FFT effect. |
af1c85da8982
fixpoint: separate windowing arithmetic imlt_window_float
mhoffman
parents:
5344
diff
changeset
|
683 */ |
af1c85da8982
fixpoint: separate windowing arithmetic imlt_window_float
mhoffman
parents:
5344
diff
changeset
|
684 |
af1c85da8982
fixpoint: separate windowing arithmetic imlt_window_float
mhoffman
parents:
5344
diff
changeset
|
685 /* Apply window and overlap */ |
af1c85da8982
fixpoint: separate windowing arithmetic imlt_window_float
mhoffman
parents:
5344
diff
changeset
|
686 for(i = 0; i < q->samples_per_channel; i++){ |
af1c85da8982
fixpoint: separate windowing arithmetic imlt_window_float
mhoffman
parents:
5344
diff
changeset
|
687 buffer1[i] = buffer1[i] * fc * q->mlt_window[i] - |
af1c85da8982
fixpoint: separate windowing arithmetic imlt_window_float
mhoffman
parents:
5344
diff
changeset
|
688 previous_buffer[i] * q->mlt_window[q->samples_per_channel - 1 - i]; |
af1c85da8982
fixpoint: separate windowing arithmetic imlt_window_float
mhoffman
parents:
5344
diff
changeset
|
689 } |
af1c85da8982
fixpoint: separate windowing arithmetic imlt_window_float
mhoffman
parents:
5344
diff
changeset
|
690 } |
2956 | 691 |
692 /** | |
4653 | 693 * The modulated lapped transform, this takes transform coefficients |
694 * and transforms them into timedomain samples. | |
695 * Apply transform window, overlap buffers, apply gain profile | |
696 * and buffer management. | |
2956 | 697 * |
698 * @param q pointer to the COOKContext | |
4653 | 699 * @param inbuffer pointer to the mltcoefficients |
4639 | 700 * @param gains_ptr current and previous gains |
2956 | 701 * @param previous_buffer pointer to the previous buffer to be used for overlapping |
702 */ | |
703 | |
4653 | 704 static void imlt_gain(COOKContext *q, float *inbuffer, |
705 cook_gains *gains_ptr, float* previous_buffer) | |
4639 | 706 { |
4653 | 707 float *buffer0 = q->mono_mdct_output; |
708 float *buffer1 = q->mono_mdct_output + q->samples_per_channel; | |
2956 | 709 int i; |
4639 | 710 |
4653 | 711 /* Inverse modified discrete cosine transform */ |
712 q->mdct_ctx.fft.imdct_calc(&q->mdct_ctx, q->mono_mdct_output, | |
713 inbuffer, q->mdct_tmp); | |
714 | |
5345
af1c85da8982
fixpoint: separate windowing arithmetic imlt_window_float
mhoffman
parents:
5344
diff
changeset
|
715 imlt_window_float (q, buffer1, gains_ptr, previous_buffer); |
2956 | 716 |
4639 | 717 /* Apply gain profile */ |
718 for (i = 0; i < 8; i++) { | |
719 if (gains_ptr->now[i] || gains_ptr->now[i + 1]) | |
4653 | 720 interpolate(q, &buffer1[q->gain_size_factor * i], |
4639 | 721 gains_ptr->now[i], gains_ptr->now[i + 1]); |
722 } | |
2956 | 723 |
724 /* Save away the current to be previous block. */ | |
4653 | 725 memcpy(previous_buffer, buffer0, sizeof(float)*q->samples_per_channel); |
2956 | 726 } |
727 | |
728 | |
729 /** | |
730 * function for getting the jointstereo coupling information | |
731 * | |
732 * @param q pointer to the COOKContext | |
733 * @param decouple_tab decoupling array | |
734 * | |
735 */ | |
736 | |
737 static void decouple_info(COOKContext *q, int* decouple_tab){ | |
738 int length, i; | |
739 | |
740 if(get_bits1(&q->gb)) { | |
741 if(cplband[q->js_subband_start] > cplband[q->subbands-1]) return; | |
742 | |
743 length = cplband[q->subbands-1] - cplband[q->js_subband_start] + 1; | |
744 for (i=0 ; i<length ; i++) { | |
745 decouple_tab[cplband[q->js_subband_start] + i] = get_vlc2(&q->gb, q->ccpl.table, q->ccpl.bits, 2); | |
746 } | |
747 return; | |
748 } | |
749 | |
750 if(cplband[q->js_subband_start] > cplband[q->subbands-1]) return; | |
751 | |
752 length = cplband[q->subbands-1] - cplband[q->js_subband_start] + 1; | |
753 for (i=0 ; i<length ; i++) { | |
754 decouple_tab[cplband[q->js_subband_start] + i] = get_bits(&q->gb, q->js_vlc_bits); | |
755 } | |
756 return; | |
757 } | |
758 | |
5344
b0d49997bd92
separate the actual math for recovering stereo from a signal channel, this is to allow fixpoint implementation
mhoffman
parents:
5343
diff
changeset
|
759 /* |
b0d49997bd92
separate the actual math for recovering stereo from a signal channel, this is to allow fixpoint implementation
mhoffman
parents:
5343
diff
changeset
|
760 * function decouples a pair of signals from a single signal via multiplication. |
b0d49997bd92
separate the actual math for recovering stereo from a signal channel, this is to allow fixpoint implementation
mhoffman
parents:
5343
diff
changeset
|
761 * |
b0d49997bd92
separate the actual math for recovering stereo from a signal channel, this is to allow fixpoint implementation
mhoffman
parents:
5343
diff
changeset
|
762 * @param q pointer to the COOKContext |
b0d49997bd92
separate the actual math for recovering stereo from a signal channel, this is to allow fixpoint implementation
mhoffman
parents:
5343
diff
changeset
|
763 * @param subband index of the current subband |
b0d49997bd92
separate the actual math for recovering stereo from a signal channel, this is to allow fixpoint implementation
mhoffman
parents:
5343
diff
changeset
|
764 * @param f1 multiplier for channel 1 extraction |
b0d49997bd92
separate the actual math for recovering stereo from a signal channel, this is to allow fixpoint implementation
mhoffman
parents:
5343
diff
changeset
|
765 * @param f2 multiplier for channel 2 extraction |
b0d49997bd92
separate the actual math for recovering stereo from a signal channel, this is to allow fixpoint implementation
mhoffman
parents:
5343
diff
changeset
|
766 * @param decode_buffer input buffer |
b0d49997bd92
separate the actual math for recovering stereo from a signal channel, this is to allow fixpoint implementation
mhoffman
parents:
5343
diff
changeset
|
767 * @param mlt_buffer1 pointer to left channel mlt coefficients |
b0d49997bd92
separate the actual math for recovering stereo from a signal channel, this is to allow fixpoint implementation
mhoffman
parents:
5343
diff
changeset
|
768 * @param mlt_buffer2 pointer to right channel mlt coefficients |
b0d49997bd92
separate the actual math for recovering stereo from a signal channel, this is to allow fixpoint implementation
mhoffman
parents:
5343
diff
changeset
|
769 */ |
b0d49997bd92
separate the actual math for recovering stereo from a signal channel, this is to allow fixpoint implementation
mhoffman
parents:
5343
diff
changeset
|
770 static void decouple_float (COOKContext *q, |
b0d49997bd92
separate the actual math for recovering stereo from a signal channel, this is to allow fixpoint implementation
mhoffman
parents:
5343
diff
changeset
|
771 int subband, |
b0d49997bd92
separate the actual math for recovering stereo from a signal channel, this is to allow fixpoint implementation
mhoffman
parents:
5343
diff
changeset
|
772 float f1, float f2, |
b0d49997bd92
separate the actual math for recovering stereo from a signal channel, this is to allow fixpoint implementation
mhoffman
parents:
5343
diff
changeset
|
773 float *decode_buffer, |
b0d49997bd92
separate the actual math for recovering stereo from a signal channel, this is to allow fixpoint implementation
mhoffman
parents:
5343
diff
changeset
|
774 float *mlt_buffer1, float *mlt_buffer2) |
b0d49997bd92
separate the actual math for recovering stereo from a signal channel, this is to allow fixpoint implementation
mhoffman
parents:
5343
diff
changeset
|
775 { |
b0d49997bd92
separate the actual math for recovering stereo from a signal channel, this is to allow fixpoint implementation
mhoffman
parents:
5343
diff
changeset
|
776 int j, tmp_idx; |
b0d49997bd92
separate the actual math for recovering stereo from a signal channel, this is to allow fixpoint implementation
mhoffman
parents:
5343
diff
changeset
|
777 for (j=0 ; j<SUBBAND_SIZE ; j++) { |
b0d49997bd92
separate the actual math for recovering stereo from a signal channel, this is to allow fixpoint implementation
mhoffman
parents:
5343
diff
changeset
|
778 tmp_idx = ((q->js_subband_start + subband)*SUBBAND_SIZE)+j; |
b0d49997bd92
separate the actual math for recovering stereo from a signal channel, this is to allow fixpoint implementation
mhoffman
parents:
5343
diff
changeset
|
779 mlt_buffer1[SUBBAND_SIZE*subband + j] = f1 * decode_buffer[tmp_idx]; |
b0d49997bd92
separate the actual math for recovering stereo from a signal channel, this is to allow fixpoint implementation
mhoffman
parents:
5343
diff
changeset
|
780 mlt_buffer2[SUBBAND_SIZE*subband + j] = f2 * decode_buffer[tmp_idx]; |
b0d49997bd92
separate the actual math for recovering stereo from a signal channel, this is to allow fixpoint implementation
mhoffman
parents:
5343
diff
changeset
|
781 } |
b0d49997bd92
separate the actual math for recovering stereo from a signal channel, this is to allow fixpoint implementation
mhoffman
parents:
5343
diff
changeset
|
782 } |
2956 | 783 |
784 /** | |
785 * function for decoding joint stereo data | |
786 * | |
787 * @param q pointer to the COOKContext | |
788 * @param mlt_buffer1 pointer to left channel mlt coefficients | |
789 * @param mlt_buffer2 pointer to right channel mlt coefficients | |
790 */ | |
791 | |
792 static void joint_decode(COOKContext *q, float* mlt_buffer1, | |
793 float* mlt_buffer2) { | |
794 int i,j; | |
795 int decouple_tab[SUBBAND_SIZE]; | |
5342
887509b42d3f
moving automatic allocation of joint_decode/decode_buffer
mhoffman
parents:
5341
diff
changeset
|
796 float *decode_buffer = q->decode_buffer_0; |
2956 | 797 int idx, cpl_tmp,tmp_idx; |
798 float f1,f2; | |
799 float* cplscale; | |
800 | |
801 memset(decouple_tab, 0, sizeof(decouple_tab)); | |
802 memset(decode_buffer, 0, sizeof(decode_buffer)); | |
803 | |
804 /* Make sure the buffers are zeroed out. */ | |
805 memset(mlt_buffer1,0, 1024*sizeof(float)); | |
806 memset(mlt_buffer2,0, 1024*sizeof(float)); | |
807 decouple_info(q, decouple_tab); | |
808 mono_decode(q, decode_buffer); | |
809 | |
810 /* The two channels are stored interleaved in decode_buffer. */ | |
811 for (i=0 ; i<q->js_subband_start ; i++) { | |
812 for (j=0 ; j<SUBBAND_SIZE ; j++) { | |
813 mlt_buffer1[i*20+j] = decode_buffer[i*40+j]; | |
814 mlt_buffer2[i*20+j] = decode_buffer[i*40+20+j]; | |
815 } | |
816 } | |
817 | |
818 /* When we reach js_subband_start (the higher frequencies) | |
819 the coefficients are stored in a coupling scheme. */ | |
820 idx = (1 << q->js_vlc_bits) - 1; | |
3009
f5898b9b8a8a
Fix an out of array access and some minor cleanup of the code.
diego
parents:
2959
diff
changeset
|
821 for (i=q->js_subband_start ; i<q->subbands ; i++) { |
f5898b9b8a8a
Fix an out of array access and some minor cleanup of the code.
diego
parents:
2959
diff
changeset
|
822 cpl_tmp = cplband[i]; |
f5898b9b8a8a
Fix an out of array access and some minor cleanup of the code.
diego
parents:
2959
diff
changeset
|
823 idx -=decouple_tab[cpl_tmp]; |
f5898b9b8a8a
Fix an out of array access and some minor cleanup of the code.
diego
parents:
2959
diff
changeset
|
824 cplscale = (float*)cplscales[q->js_vlc_bits-2]; //choose decoupler table |
f5898b9b8a8a
Fix an out of array access and some minor cleanup of the code.
diego
parents:
2959
diff
changeset
|
825 f1 = cplscale[decouple_tab[cpl_tmp]]; |
f5898b9b8a8a
Fix an out of array access and some minor cleanup of the code.
diego
parents:
2959
diff
changeset
|
826 f2 = cplscale[idx-1]; |
5344
b0d49997bd92
separate the actual math for recovering stereo from a signal channel, this is to allow fixpoint implementation
mhoffman
parents:
5343
diff
changeset
|
827 decouple_float (q, i, f1, f2, decode_buffer, mlt_buffer1, mlt_buffer2); |
3009
f5898b9b8a8a
Fix an out of array access and some minor cleanup of the code.
diego
parents:
2959
diff
changeset
|
828 idx = (1 << q->js_vlc_bits) - 1; |
2956 | 829 } |
830 } | |
831 | |
832 /** | |
4424
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
833 * First part of subpacket decoding: |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
834 * decode raw stream bytes and read gain info. |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
835 * |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
836 * @param q pointer to the COOKContext |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
837 * @param inbuffer pointer to raw stream data |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
838 * @param gain_ptr array of current/prev gain pointers |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
839 */ |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
840 |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
841 static inline void |
4428
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
842 decode_bytes_and_gain(COOKContext *q, uint8_t *inbuffer, |
4639 | 843 cook_gains *gains_ptr) |
4424
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
844 { |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
845 int offset; |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
846 |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
847 offset = decode_bytes(inbuffer, q->decoded_bytes_buffer, |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
848 q->bits_per_subpacket/8); |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
849 init_get_bits(&q->gb, q->decoded_bytes_buffer + offset, |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
850 q->bits_per_subpacket); |
4639 | 851 decode_gain_info(&q->gb, gains_ptr->now); |
4428
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
852 |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
853 /* Swap current and previous gains */ |
4639 | 854 FFSWAP(int *, gains_ptr->now, gains_ptr->previous); |
4428
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
855 } |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
856 |
5343
1aa26b974dc3
separating saturation codes so that we can support other data formats
mhoffman
parents:
5342
diff
changeset
|
857 /** |
1aa26b974dc3
separating saturation codes so that we can support other data formats
mhoffman
parents:
5342
diff
changeset
|
858 * Saturate the output signal to signed 16bit integers. |
1aa26b974dc3
separating saturation codes so that we can support other data formats
mhoffman
parents:
5342
diff
changeset
|
859 * |
1aa26b974dc3
separating saturation codes so that we can support other data formats
mhoffman
parents:
5342
diff
changeset
|
860 * @param q pointer to the COOKContext |
1aa26b974dc3
separating saturation codes so that we can support other data formats
mhoffman
parents:
5342
diff
changeset
|
861 * @param chan channel to saturate |
1aa26b974dc3
separating saturation codes so that we can support other data formats
mhoffman
parents:
5342
diff
changeset
|
862 * @param out pointer to the output vector |
1aa26b974dc3
separating saturation codes so that we can support other data formats
mhoffman
parents:
5342
diff
changeset
|
863 */ |
1aa26b974dc3
separating saturation codes so that we can support other data formats
mhoffman
parents:
5342
diff
changeset
|
864 static void |
1aa26b974dc3
separating saturation codes so that we can support other data formats
mhoffman
parents:
5342
diff
changeset
|
865 saturate_output_float (COOKContext *q, int chan, int16_t *out) |
1aa26b974dc3
separating saturation codes so that we can support other data formats
mhoffman
parents:
5342
diff
changeset
|
866 { |
1aa26b974dc3
separating saturation codes so that we can support other data formats
mhoffman
parents:
5342
diff
changeset
|
867 int j; |
1aa26b974dc3
separating saturation codes so that we can support other data formats
mhoffman
parents:
5342
diff
changeset
|
868 float_t *output = q->mono_mdct_output + q->samples_per_channel; |
1aa26b974dc3
separating saturation codes so that we can support other data formats
mhoffman
parents:
5342
diff
changeset
|
869 /* Clip and convert floats to 16 bits. |
1aa26b974dc3
separating saturation codes so that we can support other data formats
mhoffman
parents:
5342
diff
changeset
|
870 */ |
1aa26b974dc3
separating saturation codes so that we can support other data formats
mhoffman
parents:
5342
diff
changeset
|
871 for (j = 0; j < q->samples_per_channel; j++) { |
1aa26b974dc3
separating saturation codes so that we can support other data formats
mhoffman
parents:
5342
diff
changeset
|
872 out[chan + q->nb_channels * j] = |
1aa26b974dc3
separating saturation codes so that we can support other data formats
mhoffman
parents:
5342
diff
changeset
|
873 av_clip(lrintf(output[j]), -32768, 32767); |
1aa26b974dc3
separating saturation codes so that we can support other data formats
mhoffman
parents:
5342
diff
changeset
|
874 } |
1aa26b974dc3
separating saturation codes so that we can support other data formats
mhoffman
parents:
5342
diff
changeset
|
875 } |
1aa26b974dc3
separating saturation codes so that we can support other data formats
mhoffman
parents:
5342
diff
changeset
|
876 |
4428
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
877 /** |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
878 * Final part of subpacket decoding: |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
879 * Apply modulated lapped transform, gain compensation, |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
880 * clip and convert to integer. |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
881 * |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
882 * @param q pointer to the COOKContext |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
883 * @param decode_buffer pointer to the mlt coefficients |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
884 * @param gain_ptr array of current/prev gain pointers |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
885 * @param previous_buffer pointer to the previous buffer to be used for overlapping |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
886 * @param out pointer to the output buffer |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
887 * @param chan 0: left or single channel, 1: right channel |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
888 */ |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
889 |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
890 static inline void |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
891 mlt_compensate_output(COOKContext *q, float *decode_buffer, |
4639 | 892 cook_gains *gains, float *previous_buffer, |
4428
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
893 int16_t *out, int chan) |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
894 { |
4653 | 895 imlt_gain(q, decode_buffer, gains, previous_buffer); |
5343
1aa26b974dc3
separating saturation codes so that we can support other data formats
mhoffman
parents:
5342
diff
changeset
|
896 saturate_output_float (q, chan, out); |
4424
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
897 } |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
898 |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
899 |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
900 /** |
2956 | 901 * Cook subpacket decoding. This function returns one decoded subpacket, |
902 * usually 1024 samples per channel. | |
903 * | |
904 * @param q pointer to the COOKContext | |
905 * @param inbuffer pointer to the inbuffer | |
906 * @param sub_packet_size subpacket size | |
907 * @param outbuffer pointer to the outbuffer | |
908 */ | |
909 | |
910 | |
911 static int decode_subpacket(COOKContext *q, uint8_t *inbuffer, | |
912 int sub_packet_size, int16_t *outbuffer) { | |
913 /* packet dump */ | |
914 // for (i=0 ; i<sub_packet_size ; i++) { | |
915 // av_log(NULL, AV_LOG_ERROR, "%02x", inbuffer[i]); | |
916 // } | |
917 // av_log(NULL, AV_LOG_ERROR, "\n"); | |
918 | |
4639 | 919 decode_bytes_and_gain(q, inbuffer, &q->gains1); |
2956 | 920 |
4428
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
921 if (q->joint_stereo) { |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
922 joint_decode(q, q->decode_buffer_1, q->decode_buffer_2); |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
923 } else { |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
924 mono_decode(q, q->decode_buffer_1); |
2956 | 925 |
4428
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
926 if (q->nb_channels == 2) { |
4639 | 927 decode_bytes_and_gain(q, inbuffer + sub_packet_size/2, &q->gains2); |
4428
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
928 mono_decode(q, q->decode_buffer_2); |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
929 } |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
930 } |
3014
959b8ad880dc
Dual mono stereo strems sound ok now, added sanity checks and removed
rtognimp
parents:
3009
diff
changeset
|
931 |
4639 | 932 mlt_compensate_output(q, q->decode_buffer_1, &q->gains1, |
4428
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
933 q->mono_previous_buffer1, outbuffer, 0); |
2956 | 934 |
4428
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
935 if (q->nb_channels == 2) { |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
936 if (q->joint_stereo) { |
4639 | 937 mlt_compensate_output(q, q->decode_buffer_2, &q->gains1, |
4428
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
938 q->mono_previous_buffer2, outbuffer, 1); |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
939 } else { |
4639 | 940 mlt_compensate_output(q, q->decode_buffer_2, &q->gains2, |
4428
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
941 q->mono_previous_buffer2, outbuffer, 1); |
2956 | 942 } |
943 } | |
2959
24805f4d1b84
This patch adds some support for non-joint stereo streams. It also
rtognimp
parents:
2956
diff
changeset
|
944 return q->samples_per_frame * sizeof(int16_t); |
2956 | 945 } |
946 | |
947 | |
948 /** | |
949 * Cook frame decoding | |
950 * | |
951 * @param avctx pointer to the AVCodecContext | |
952 */ | |
953 | |
954 static int cook_decode_frame(AVCodecContext *avctx, | |
955 void *data, int *data_size, | |
956 uint8_t *buf, int buf_size) { | |
957 COOKContext *q = avctx->priv_data; | |
958 | |
959 if (buf_size < avctx->block_align) | |
960 return buf_size; | |
961 | |
962 *data_size = decode_subpacket(q, buf, avctx->block_align, data); | |
963 | |
4638
9f74306d4ac7
Don't output the first two frames, since they don't contain valid audio.
banan
parents:
4594
diff
changeset
|
964 /* Discard the first two frames: no valid audio. */ |
9f74306d4ac7
Don't output the first two frames, since they don't contain valid audio.
banan
parents:
4594
diff
changeset
|
965 if (avctx->frame_number < 2) *data_size = 0; |
9f74306d4ac7
Don't output the first two frames, since they don't contain valid audio.
banan
parents:
4594
diff
changeset
|
966 |
2956 | 967 return avctx->block_align; |
968 } | |
3090 | 969 |
2956 | 970 #ifdef COOKDEBUG |
4430
407e7fd9b4f4
Get rid of the COOKextradata struct. And use valid C to parse the extradata.
banan
parents:
4429
diff
changeset
|
971 static void dump_cook_context(COOKContext *q) |
2956 | 972 { |
973 //int i=0; | |
974 #define PRINT(a,b) av_log(NULL,AV_LOG_ERROR," %s = %d\n", a, b); | |
975 av_log(NULL,AV_LOG_ERROR,"COOKextradata\n"); | |
4430
407e7fd9b4f4
Get rid of the COOKextradata struct. And use valid C to parse the extradata.
banan
parents:
4429
diff
changeset
|
976 av_log(NULL,AV_LOG_ERROR,"cookversion=%x\n",q->cookversion); |
407e7fd9b4f4
Get rid of the COOKextradata struct. And use valid C to parse the extradata.
banan
parents:
4429
diff
changeset
|
977 if (q->cookversion > STEREO) { |
407e7fd9b4f4
Get rid of the COOKextradata struct. And use valid C to parse the extradata.
banan
parents:
4429
diff
changeset
|
978 PRINT("js_subband_start",q->js_subband_start); |
407e7fd9b4f4
Get rid of the COOKextradata struct. And use valid C to parse the extradata.
banan
parents:
4429
diff
changeset
|
979 PRINT("js_vlc_bits",q->js_vlc_bits); |
2956 | 980 } |
981 av_log(NULL,AV_LOG_ERROR,"COOKContext\n"); | |
982 PRINT("nb_channels",q->nb_channels); | |
983 PRINT("bit_rate",q->bit_rate); | |
984 PRINT("sample_rate",q->sample_rate); | |
985 PRINT("samples_per_channel",q->samples_per_channel); | |
986 PRINT("samples_per_frame",q->samples_per_frame); | |
987 PRINT("subbands",q->subbands); | |
988 PRINT("random_state",q->random_state); | |
989 PRINT("js_subband_start",q->js_subband_start); | |
3090 | 990 PRINT("log2_numvector_size",q->log2_numvector_size); |
2956 | 991 PRINT("numvector_size",q->numvector_size); |
992 PRINT("total_subbands",q->total_subbands); | |
993 } | |
994 #endif | |
3090 | 995 |
2956 | 996 /** |
997 * Cook initialization | |
998 * | |
999 * @param avctx pointer to the AVCodecContext | |
1000 */ | |
1001 | |
1002 static int cook_decode_init(AVCodecContext *avctx) | |
1003 { | |
1004 COOKContext *q = avctx->priv_data; | |
4430
407e7fd9b4f4
Get rid of the COOKextradata struct. And use valid C to parse the extradata.
banan
parents:
4429
diff
changeset
|
1005 uint8_t *edata_ptr = avctx->extradata; |
2956 | 1006 |
1007 /* Take care of the codec specific extradata. */ | |
1008 if (avctx->extradata_size <= 0) { | |
4302
a0f83004d485
av_log(NULL,... -> av_log(avctx,.. where appropriate.
banan
parents:
3947
diff
changeset
|
1009 av_log(avctx,AV_LOG_ERROR,"Necessary extradata missing!\n"); |
2956 | 1010 return -1; |
1011 } else { | |
1012 /* 8 for mono, 16 for stereo, ? for multichannel | |
1013 Swap to right endianness so we don't need to care later on. */ | |
4302
a0f83004d485
av_log(NULL,... -> av_log(avctx,.. where appropriate.
banan
parents:
3947
diff
changeset
|
1014 av_log(avctx,AV_LOG_DEBUG,"codecdata_length=%d\n",avctx->extradata_size); |
2956 | 1015 if (avctx->extradata_size >= 8){ |
4542 | 1016 q->cookversion = bytestream_get_be32(&edata_ptr); |
1017 q->samples_per_frame = bytestream_get_be16(&edata_ptr); | |
1018 q->subbands = bytestream_get_be16(&edata_ptr); | |
2956 | 1019 } |
1020 if (avctx->extradata_size >= 16){ | |
4542 | 1021 bytestream_get_be32(&edata_ptr); //Unknown unused |
1022 q->js_subband_start = bytestream_get_be16(&edata_ptr); | |
1023 q->js_vlc_bits = bytestream_get_be16(&edata_ptr); | |
2956 | 1024 } |
1025 } | |
1026 | |
1027 /* Take data from the AVCodecContext (RM container). */ | |
1028 q->sample_rate = avctx->sample_rate; | |
1029 q->nb_channels = avctx->channels; | |
1030 q->bit_rate = avctx->bit_rate; | |
1031 | |
4692
3318e3f6470f
Small simplifications to subband coefficient handling and use av_random().
banan
parents:
4674
diff
changeset
|
1032 /* Initialize RNG. */ |
3318e3f6470f
Small simplifications to subband coefficient handling and use av_random().
banan
parents:
4674
diff
changeset
|
1033 av_init_random(1, &q->random_state); |
2956 | 1034 |
1035 /* Initialize extradata related variables. */ | |
4430
407e7fd9b4f4
Get rid of the COOKextradata struct. And use valid C to parse the extradata.
banan
parents:
4429
diff
changeset
|
1036 q->samples_per_channel = q->samples_per_frame / q->nb_channels; |
2956 | 1037 q->bits_per_subpacket = avctx->block_align * 8; |
1038 | |
1039 /* Initialize default data states. */ | |
3090 | 1040 q->log2_numvector_size = 5; |
2956 | 1041 q->total_subbands = q->subbands; |
1042 | |
1043 /* Initialize version-dependent variables */ | |
4430
407e7fd9b4f4
Get rid of the COOKextradata struct. And use valid C to parse the extradata.
banan
parents:
4429
diff
changeset
|
1044 av_log(NULL,AV_LOG_DEBUG,"q->cookversion=%x\n",q->cookversion); |
4428
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
1045 q->joint_stereo = 0; |
4430
407e7fd9b4f4
Get rid of the COOKextradata struct. And use valid C to parse the extradata.
banan
parents:
4429
diff
changeset
|
1046 switch (q->cookversion) { |
4425 | 1047 case MONO: |
2956 | 1048 if (q->nb_channels != 1) { |
4302
a0f83004d485
av_log(NULL,... -> av_log(avctx,.. where appropriate.
banan
parents:
3947
diff
changeset
|
1049 av_log(avctx,AV_LOG_ERROR,"Container channels != 1, report sample!\n"); |
2956 | 1050 return -1; |
1051 } | |
4425 | 1052 av_log(avctx,AV_LOG_DEBUG,"MONO\n"); |
2956 | 1053 break; |
4425 | 1054 case STEREO: |
2956 | 1055 if (q->nb_channels != 1) { |
2959
24805f4d1b84
This patch adds some support for non-joint stereo streams. It also
rtognimp
parents:
2956
diff
changeset
|
1056 q->bits_per_subpacket = q->bits_per_subpacket/2; |
2956 | 1057 } |
4425 | 1058 av_log(avctx,AV_LOG_DEBUG,"STEREO\n"); |
2956 | 1059 break; |
1060 case JOINT_STEREO: | |
1061 if (q->nb_channels != 2) { | |
4302
a0f83004d485
av_log(NULL,... -> av_log(avctx,.. where appropriate.
banan
parents:
3947
diff
changeset
|
1062 av_log(avctx,AV_LOG_ERROR,"Container channels != 2, report sample!\n"); |
2956 | 1063 return -1; |
1064 } | |
4302
a0f83004d485
av_log(NULL,... -> av_log(avctx,.. where appropriate.
banan
parents:
3947
diff
changeset
|
1065 av_log(avctx,AV_LOG_DEBUG,"JOINT_STEREO\n"); |
2956 | 1066 if (avctx->extradata_size >= 16){ |
4430
407e7fd9b4f4
Get rid of the COOKextradata struct. And use valid C to parse the extradata.
banan
parents:
4429
diff
changeset
|
1067 q->total_subbands = q->subbands + q->js_subband_start; |
2956 | 1068 q->joint_stereo = 1; |
1069 } | |
1070 if (q->samples_per_channel > 256) { | |
3091
0284d5b34916
Fix broken cosmetics commit and add a check for valid headers.
banan
parents:
3090
diff
changeset
|
1071 q->log2_numvector_size = 6; |
2956 | 1072 } |
1073 if (q->samples_per_channel > 512) { | |
3091
0284d5b34916
Fix broken cosmetics commit and add a check for valid headers.
banan
parents:
3090
diff
changeset
|
1074 q->log2_numvector_size = 7; |
2956 | 1075 } |
1076 break; | |
1077 case MC_COOK: | |
4302
a0f83004d485
av_log(NULL,... -> av_log(avctx,.. where appropriate.
banan
parents:
3947
diff
changeset
|
1078 av_log(avctx,AV_LOG_ERROR,"MC_COOK not supported!\n"); |
2956 | 1079 return -1; |
1080 break; | |
1081 default: | |
4302
a0f83004d485
av_log(NULL,... -> av_log(avctx,.. where appropriate.
banan
parents:
3947
diff
changeset
|
1082 av_log(avctx,AV_LOG_ERROR,"Unknown Cook version, report sample!\n"); |
2956 | 1083 return -1; |
1084 break; | |
1085 } | |
1086 | |
1087 /* Initialize variable relations */ | |
3090 | 1088 q->numvector_size = (1 << q->log2_numvector_size); |
2956 | 1089 |
1090 /* Generate tables */ | |
1091 init_rootpow2table(q); | |
1092 init_pow2table(q); | |
1093 init_gain_table(q); | |
1094 | |
1095 if (init_cook_vlc_tables(q) != 0) | |
1096 return -1; | |
1097 | |
3303
68721b62a528
sanity checks, some might have been exploitable ...
michael
parents:
3106
diff
changeset
|
1098 |
68721b62a528
sanity checks, some might have been exploitable ...
michael
parents:
3106
diff
changeset
|
1099 if(avctx->block_align >= UINT_MAX/2) |
68721b62a528
sanity checks, some might have been exploitable ...
michael
parents:
3106
diff
changeset
|
1100 return -1; |
68721b62a528
sanity checks, some might have been exploitable ...
michael
parents:
3106
diff
changeset
|
1101 |
4424
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
1102 /* Pad the databuffer with: |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
1103 DECODE_BYTES_PAD1 or DECODE_BYTES_PAD2 for decode_bytes(), |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
1104 FF_INPUT_BUFFER_PADDING_SIZE, for the bitstreamreader. */ |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
1105 if (q->nb_channels==2 && q->joint_stereo==0) { |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
1106 q->decoded_bytes_buffer = |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
1107 av_mallocz(avctx->block_align/2 |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
1108 + DECODE_BYTES_PAD2(avctx->block_align/2) |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
1109 + FF_INPUT_BUFFER_PADDING_SIZE); |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
1110 } else { |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
1111 q->decoded_bytes_buffer = |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
1112 av_mallocz(avctx->block_align |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
1113 + DECODE_BYTES_PAD1(avctx->block_align) |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
1114 + FF_INPUT_BUFFER_PADDING_SIZE); |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
1115 } |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
1116 if (q->decoded_bytes_buffer == NULL) |
2956 | 1117 return -1; |
1118 | |
4639 | 1119 q->gains1.now = q->gain_1; |
1120 q->gains1.previous = q->gain_2; | |
1121 q->gains2.now = q->gain_3; | |
1122 q->gains2.previous = q->gain_4; | |
2956 | 1123 |
1124 /* Initialize transform. */ | |
4650
31bf54d9353d
Replace custom modified discrete cosine transform with ffmpeg's own.
banan
parents:
4639
diff
changeset
|
1125 if ( init_cook_mlt(q) != 0 ) |
2956 | 1126 return -1; |
3014
959b8ad880dc
Dual mono stereo strems sound ok now, added sanity checks and removed
rtognimp
parents:
3009
diff
changeset
|
1127 |
959b8ad880dc
Dual mono stereo strems sound ok now, added sanity checks and removed
rtognimp
parents:
3009
diff
changeset
|
1128 /* Try to catch some obviously faulty streams, othervise it might be exploitable */ |
959b8ad880dc
Dual mono stereo strems sound ok now, added sanity checks and removed
rtognimp
parents:
3009
diff
changeset
|
1129 if (q->total_subbands > 53) { |
4302
a0f83004d485
av_log(NULL,... -> av_log(avctx,.. where appropriate.
banan
parents:
3947
diff
changeset
|
1130 av_log(avctx,AV_LOG_ERROR,"total_subbands > 53, report sample!\n"); |
3014
959b8ad880dc
Dual mono stereo strems sound ok now, added sanity checks and removed
rtognimp
parents:
3009
diff
changeset
|
1131 return -1; |
959b8ad880dc
Dual mono stereo strems sound ok now, added sanity checks and removed
rtognimp
parents:
3009
diff
changeset
|
1132 } |
959b8ad880dc
Dual mono stereo strems sound ok now, added sanity checks and removed
rtognimp
parents:
3009
diff
changeset
|
1133 if (q->subbands > 50) { |
4302
a0f83004d485
av_log(NULL,... -> av_log(avctx,.. where appropriate.
banan
parents:
3947
diff
changeset
|
1134 av_log(avctx,AV_LOG_ERROR,"subbands > 50, report sample!\n"); |
3014
959b8ad880dc
Dual mono stereo strems sound ok now, added sanity checks and removed
rtognimp
parents:
3009
diff
changeset
|
1135 return -1; |
959b8ad880dc
Dual mono stereo strems sound ok now, added sanity checks and removed
rtognimp
parents:
3009
diff
changeset
|
1136 } |
3091
0284d5b34916
Fix broken cosmetics commit and add a check for valid headers.
banan
parents:
3090
diff
changeset
|
1137 if ((q->samples_per_channel == 256) || (q->samples_per_channel == 512) || (q->samples_per_channel == 1024)) { |
0284d5b34916
Fix broken cosmetics commit and add a check for valid headers.
banan
parents:
3090
diff
changeset
|
1138 } else { |
4302
a0f83004d485
av_log(NULL,... -> av_log(avctx,.. where appropriate.
banan
parents:
3947
diff
changeset
|
1139 av_log(avctx,AV_LOG_ERROR,"unknown amount of samples_per_channel = %d, report sample!\n",q->samples_per_channel); |
3091
0284d5b34916
Fix broken cosmetics commit and add a check for valid headers.
banan
parents:
3090
diff
changeset
|
1140 return -1; |
0284d5b34916
Fix broken cosmetics commit and add a check for valid headers.
banan
parents:
3090
diff
changeset
|
1141 } |
4431
85ac154efd99
Check that js_vlc_bits from the extradata is in a valid range.
banan
parents:
4430
diff
changeset
|
1142 if ((q->js_vlc_bits > 6) || (q->js_vlc_bits < 0)) { |
85ac154efd99
Check that js_vlc_bits from the extradata is in a valid range.
banan
parents:
4430
diff
changeset
|
1143 av_log(avctx,AV_LOG_ERROR,"q->js_vlc_bits = %d, only >= 0 and <= 6 allowed!\n",q->js_vlc_bits); |
85ac154efd99
Check that js_vlc_bits from the extradata is in a valid range.
banan
parents:
4430
diff
changeset
|
1144 return -1; |
85ac154efd99
Check that js_vlc_bits from the extradata is in a valid range.
banan
parents:
4430
diff
changeset
|
1145 } |
3014
959b8ad880dc
Dual mono stereo strems sound ok now, added sanity checks and removed
rtognimp
parents:
3009
diff
changeset
|
1146 |
3009
f5898b9b8a8a
Fix an out of array access and some minor cleanup of the code.
diego
parents:
2959
diff
changeset
|
1147 #ifdef COOKDEBUG |
4430
407e7fd9b4f4
Get rid of the COOKextradata struct. And use valid C to parse the extradata.
banan
parents:
4429
diff
changeset
|
1148 dump_cook_context(q); |
3009
f5898b9b8a8a
Fix an out of array access and some minor cleanup of the code.
diego
parents:
2959
diff
changeset
|
1149 #endif |
2956 | 1150 return 0; |
1151 } | |
1152 | |
1153 | |
1154 AVCodec cook_decoder = | |
1155 { | |
1156 .name = "cook", | |
1157 .type = CODEC_TYPE_AUDIO, | |
1158 .id = CODEC_ID_COOK, | |
1159 .priv_data_size = sizeof(COOKContext), | |
1160 .init = cook_decode_init, | |
1161 .close = cook_decode_close, | |
1162 .decode = cook_decode_frame, | |
1163 }; |