Mercurial > libavcodec.hg
annotate cook.c @ 5344:b0d49997bd92 libavcodec
separate the actual math for recovering stereo from a signal channel, this is to allow fixpoint implementation
author | mhoffman |
---|---|
date | Mon, 16 Jul 2007 11:25:56 +0000 |
parents | 1aa26b974dc3 |
children | af1c85da8982 |
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 | |
665 | |
666 /** | |
4653 | 667 * The modulated lapped transform, this takes transform coefficients |
668 * and transforms them into timedomain samples. | |
669 * Apply transform window, overlap buffers, apply gain profile | |
670 * and buffer management. | |
2956 | 671 * |
672 * @param q pointer to the COOKContext | |
4653 | 673 * @param inbuffer pointer to the mltcoefficients |
4639 | 674 * @param gains_ptr current and previous gains |
2956 | 675 * @param previous_buffer pointer to the previous buffer to be used for overlapping |
676 */ | |
677 | |
4653 | 678 static void imlt_gain(COOKContext *q, float *inbuffer, |
679 cook_gains *gains_ptr, float* previous_buffer) | |
4639 | 680 { |
681 const float fc = q->pow2tab[gains_ptr->previous[0] + 63]; | |
4653 | 682 float *buffer0 = q->mono_mdct_output; |
683 float *buffer1 = q->mono_mdct_output + q->samples_per_channel; | |
2956 | 684 int i; |
4639 | 685 |
4653 | 686 /* Inverse modified discrete cosine transform */ |
687 q->mdct_ctx.fft.imdct_calc(&q->mdct_ctx, q->mono_mdct_output, | |
688 inbuffer, q->mdct_tmp); | |
689 | |
690 /* The weird thing here, is that the two halves of the time domain | |
691 * buffer are swapped. Also, the newest data, that we save away for | |
692 * next frame, has the wrong sign. Hence the subtraction below. | |
693 * Almost sounds like a complex conjugate/reverse data/FFT effect. | |
694 */ | |
695 | |
696 /* Apply window and overlap */ | |
697 for(i = 0; i < q->samples_per_channel; i++){ | |
698 buffer1[i] = buffer1[i] * fc * q->mlt_window[i] - | |
699 previous_buffer[i] * q->mlt_window[q->samples_per_channel - 1 - i]; | |
2956 | 700 } |
701 | |
4639 | 702 /* Apply gain profile */ |
703 for (i = 0; i < 8; i++) { | |
704 if (gains_ptr->now[i] || gains_ptr->now[i + 1]) | |
4653 | 705 interpolate(q, &buffer1[q->gain_size_factor * i], |
4639 | 706 gains_ptr->now[i], gains_ptr->now[i + 1]); |
707 } | |
2956 | 708 |
709 /* Save away the current to be previous block. */ | |
4653 | 710 memcpy(previous_buffer, buffer0, sizeof(float)*q->samples_per_channel); |
2956 | 711 } |
712 | |
713 | |
714 /** | |
715 * function for getting the jointstereo coupling information | |
716 * | |
717 * @param q pointer to the COOKContext | |
718 * @param decouple_tab decoupling array | |
719 * | |
720 */ | |
721 | |
722 static void decouple_info(COOKContext *q, int* decouple_tab){ | |
723 int length, i; | |
724 | |
725 if(get_bits1(&q->gb)) { | |
726 if(cplband[q->js_subband_start] > cplband[q->subbands-1]) return; | |
727 | |
728 length = cplband[q->subbands-1] - cplband[q->js_subband_start] + 1; | |
729 for (i=0 ; i<length ; i++) { | |
730 decouple_tab[cplband[q->js_subband_start] + i] = get_vlc2(&q->gb, q->ccpl.table, q->ccpl.bits, 2); | |
731 } | |
732 return; | |
733 } | |
734 | |
735 if(cplband[q->js_subband_start] > cplband[q->subbands-1]) return; | |
736 | |
737 length = cplband[q->subbands-1] - cplband[q->js_subband_start] + 1; | |
738 for (i=0 ; i<length ; i++) { | |
739 decouple_tab[cplband[q->js_subband_start] + i] = get_bits(&q->gb, q->js_vlc_bits); | |
740 } | |
741 return; | |
742 } | |
743 | |
5344
b0d49997bd92
separate the actual math for recovering stereo from a signal channel, this is to allow fixpoint implementation
mhoffman
parents:
5343
diff
changeset
|
744 /* |
b0d49997bd92
separate the actual math for recovering stereo from a signal channel, this is to allow fixpoint implementation
mhoffman
parents:
5343
diff
changeset
|
745 * 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
|
746 * |
b0d49997bd92
separate the actual math for recovering stereo from a signal channel, this is to allow fixpoint implementation
mhoffman
parents:
5343
diff
changeset
|
747 * @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
|
748 * @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
|
749 * @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
|
750 * @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
|
751 * @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
|
752 * @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
|
753 * @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
|
754 */ |
b0d49997bd92
separate the actual math for recovering stereo from a signal channel, this is to allow fixpoint implementation
mhoffman
parents:
5343
diff
changeset
|
755 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
|
756 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
|
757 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
|
758 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
|
759 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
|
760 { |
b0d49997bd92
separate the actual math for recovering stereo from a signal channel, this is to allow fixpoint implementation
mhoffman
parents:
5343
diff
changeset
|
761 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
|
762 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
|
763 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
|
764 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
|
765 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
|
766 } |
b0d49997bd92
separate the actual math for recovering stereo from a signal channel, this is to allow fixpoint implementation
mhoffman
parents:
5343
diff
changeset
|
767 } |
2956 | 768 |
769 /** | |
770 * function for decoding joint stereo data | |
771 * | |
772 * @param q pointer to the COOKContext | |
773 * @param mlt_buffer1 pointer to left channel mlt coefficients | |
774 * @param mlt_buffer2 pointer to right channel mlt coefficients | |
775 */ | |
776 | |
777 static void joint_decode(COOKContext *q, float* mlt_buffer1, | |
778 float* mlt_buffer2) { | |
779 int i,j; | |
780 int decouple_tab[SUBBAND_SIZE]; | |
5342
887509b42d3f
moving automatic allocation of joint_decode/decode_buffer
mhoffman
parents:
5341
diff
changeset
|
781 float *decode_buffer = q->decode_buffer_0; |
2956 | 782 int idx, cpl_tmp,tmp_idx; |
783 float f1,f2; | |
784 float* cplscale; | |
785 | |
786 memset(decouple_tab, 0, sizeof(decouple_tab)); | |
787 memset(decode_buffer, 0, sizeof(decode_buffer)); | |
788 | |
789 /* Make sure the buffers are zeroed out. */ | |
790 memset(mlt_buffer1,0, 1024*sizeof(float)); | |
791 memset(mlt_buffer2,0, 1024*sizeof(float)); | |
792 decouple_info(q, decouple_tab); | |
793 mono_decode(q, decode_buffer); | |
794 | |
795 /* The two channels are stored interleaved in decode_buffer. */ | |
796 for (i=0 ; i<q->js_subband_start ; i++) { | |
797 for (j=0 ; j<SUBBAND_SIZE ; j++) { | |
798 mlt_buffer1[i*20+j] = decode_buffer[i*40+j]; | |
799 mlt_buffer2[i*20+j] = decode_buffer[i*40+20+j]; | |
800 } | |
801 } | |
802 | |
803 /* When we reach js_subband_start (the higher frequencies) | |
804 the coefficients are stored in a coupling scheme. */ | |
805 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
|
806 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
|
807 cpl_tmp = cplband[i]; |
f5898b9b8a8a
Fix an out of array access and some minor cleanup of the code.
diego
parents:
2959
diff
changeset
|
808 idx -=decouple_tab[cpl_tmp]; |
f5898b9b8a8a
Fix an out of array access and some minor cleanup of the code.
diego
parents:
2959
diff
changeset
|
809 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
|
810 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
|
811 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
|
812 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
|
813 idx = (1 << q->js_vlc_bits) - 1; |
2956 | 814 } |
815 } | |
816 | |
817 /** | |
4424
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
818 * First part of subpacket decoding: |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
819 * decode raw stream bytes and read gain info. |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
820 * |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
821 * @param q pointer to the COOKContext |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
822 * @param inbuffer pointer to raw stream data |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
823 * @param gain_ptr array of current/prev gain pointers |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
824 */ |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
825 |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
826 static inline void |
4428
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
827 decode_bytes_and_gain(COOKContext *q, uint8_t *inbuffer, |
4639 | 828 cook_gains *gains_ptr) |
4424
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
829 { |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
830 int offset; |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
831 |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
832 offset = decode_bytes(inbuffer, q->decoded_bytes_buffer, |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
833 q->bits_per_subpacket/8); |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
834 init_get_bits(&q->gb, q->decoded_bytes_buffer + offset, |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
835 q->bits_per_subpacket); |
4639 | 836 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
|
837 |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
838 /* Swap current and previous gains */ |
4639 | 839 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
|
840 } |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
841 |
5343
1aa26b974dc3
separating saturation codes so that we can support other data formats
mhoffman
parents:
5342
diff
changeset
|
842 /** |
1aa26b974dc3
separating saturation codes so that we can support other data formats
mhoffman
parents:
5342
diff
changeset
|
843 * 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
|
844 * |
1aa26b974dc3
separating saturation codes so that we can support other data formats
mhoffman
parents:
5342
diff
changeset
|
845 * @param q pointer to the COOKContext |
1aa26b974dc3
separating saturation codes so that we can support other data formats
mhoffman
parents:
5342
diff
changeset
|
846 * @param chan channel to saturate |
1aa26b974dc3
separating saturation codes so that we can support other data formats
mhoffman
parents:
5342
diff
changeset
|
847 * @param out pointer to the output vector |
1aa26b974dc3
separating saturation codes so that we can support other data formats
mhoffman
parents:
5342
diff
changeset
|
848 */ |
1aa26b974dc3
separating saturation codes so that we can support other data formats
mhoffman
parents:
5342
diff
changeset
|
849 static void |
1aa26b974dc3
separating saturation codes so that we can support other data formats
mhoffman
parents:
5342
diff
changeset
|
850 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
|
851 { |
1aa26b974dc3
separating saturation codes so that we can support other data formats
mhoffman
parents:
5342
diff
changeset
|
852 int j; |
1aa26b974dc3
separating saturation codes so that we can support other data formats
mhoffman
parents:
5342
diff
changeset
|
853 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
|
854 /* Clip and convert floats to 16 bits. |
1aa26b974dc3
separating saturation codes so that we can support other data formats
mhoffman
parents:
5342
diff
changeset
|
855 */ |
1aa26b974dc3
separating saturation codes so that we can support other data formats
mhoffman
parents:
5342
diff
changeset
|
856 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
|
857 out[chan + q->nb_channels * j] = |
1aa26b974dc3
separating saturation codes so that we can support other data formats
mhoffman
parents:
5342
diff
changeset
|
858 av_clip(lrintf(output[j]), -32768, 32767); |
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 } |
1aa26b974dc3
separating saturation codes so that we can support other data formats
mhoffman
parents:
5342
diff
changeset
|
861 |
4428
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
862 /** |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
863 * Final part of subpacket decoding: |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
864 * Apply modulated lapped transform, gain compensation, |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
865 * clip and convert to integer. |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
866 * |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
867 * @param q pointer to the COOKContext |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
868 * @param decode_buffer pointer to the mlt coefficients |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
869 * @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
|
870 * @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
|
871 * @param out pointer to the output buffer |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
872 * @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
|
873 */ |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
874 |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
875 static inline void |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
876 mlt_compensate_output(COOKContext *q, float *decode_buffer, |
4639 | 877 cook_gains *gains, float *previous_buffer, |
4428
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
878 int16_t *out, int chan) |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
879 { |
4653 | 880 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
|
881 saturate_output_float (q, chan, out); |
4424
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
882 } |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
883 |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
884 |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
885 /** |
2956 | 886 * Cook subpacket decoding. This function returns one decoded subpacket, |
887 * usually 1024 samples per channel. | |
888 * | |
889 * @param q pointer to the COOKContext | |
890 * @param inbuffer pointer to the inbuffer | |
891 * @param sub_packet_size subpacket size | |
892 * @param outbuffer pointer to the outbuffer | |
893 */ | |
894 | |
895 | |
896 static int decode_subpacket(COOKContext *q, uint8_t *inbuffer, | |
897 int sub_packet_size, int16_t *outbuffer) { | |
898 /* packet dump */ | |
899 // for (i=0 ; i<sub_packet_size ; i++) { | |
900 // av_log(NULL, AV_LOG_ERROR, "%02x", inbuffer[i]); | |
901 // } | |
902 // av_log(NULL, AV_LOG_ERROR, "\n"); | |
903 | |
4639 | 904 decode_bytes_and_gain(q, inbuffer, &q->gains1); |
2956 | 905 |
4428
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
906 if (q->joint_stereo) { |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
907 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
|
908 } else { |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
909 mono_decode(q, q->decode_buffer_1); |
2956 | 910 |
4428
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
911 if (q->nb_channels == 2) { |
4639 | 912 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
|
913 mono_decode(q, q->decode_buffer_2); |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
914 } |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
915 } |
3014
959b8ad880dc
Dual mono stereo strems sound ok now, added sanity checks and removed
rtognimp
parents:
3009
diff
changeset
|
916 |
4639 | 917 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
|
918 q->mono_previous_buffer1, outbuffer, 0); |
2956 | 919 |
4428
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
920 if (q->nb_channels == 2) { |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
921 if (q->joint_stereo) { |
4639 | 922 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
|
923 q->mono_previous_buffer2, outbuffer, 1); |
eeb16216b454
decode_subpacket cleanup by Ian Braithwaite ian braithwaite dot dk.
banan
parents:
4425
diff
changeset
|
924 } else { |
4639 | 925 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
|
926 q->mono_previous_buffer2, outbuffer, 1); |
2956 | 927 } |
928 } | |
2959
24805f4d1b84
This patch adds some support for non-joint stereo streams. It also
rtognimp
parents:
2956
diff
changeset
|
929 return q->samples_per_frame * sizeof(int16_t); |
2956 | 930 } |
931 | |
932 | |
933 /** | |
934 * Cook frame decoding | |
935 * | |
936 * @param avctx pointer to the AVCodecContext | |
937 */ | |
938 | |
939 static int cook_decode_frame(AVCodecContext *avctx, | |
940 void *data, int *data_size, | |
941 uint8_t *buf, int buf_size) { | |
942 COOKContext *q = avctx->priv_data; | |
943 | |
944 if (buf_size < avctx->block_align) | |
945 return buf_size; | |
946 | |
947 *data_size = decode_subpacket(q, buf, avctx->block_align, data); | |
948 | |
4638
9f74306d4ac7
Don't output the first two frames, since they don't contain valid audio.
banan
parents:
4594
diff
changeset
|
949 /* 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
|
950 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
|
951 |
2956 | 952 return avctx->block_align; |
953 } | |
3090 | 954 |
2956 | 955 #ifdef COOKDEBUG |
4430
407e7fd9b4f4
Get rid of the COOKextradata struct. And use valid C to parse the extradata.
banan
parents:
4429
diff
changeset
|
956 static void dump_cook_context(COOKContext *q) |
2956 | 957 { |
958 //int i=0; | |
959 #define PRINT(a,b) av_log(NULL,AV_LOG_ERROR," %s = %d\n", a, b); | |
960 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
|
961 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
|
962 if (q->cookversion > STEREO) { |
407e7fd9b4f4
Get rid of the COOKextradata struct. And use valid C to parse the extradata.
banan
parents:
4429
diff
changeset
|
963 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
|
964 PRINT("js_vlc_bits",q->js_vlc_bits); |
2956 | 965 } |
966 av_log(NULL,AV_LOG_ERROR,"COOKContext\n"); | |
967 PRINT("nb_channels",q->nb_channels); | |
968 PRINT("bit_rate",q->bit_rate); | |
969 PRINT("sample_rate",q->sample_rate); | |
970 PRINT("samples_per_channel",q->samples_per_channel); | |
971 PRINT("samples_per_frame",q->samples_per_frame); | |
972 PRINT("subbands",q->subbands); | |
973 PRINT("random_state",q->random_state); | |
974 PRINT("js_subband_start",q->js_subband_start); | |
3090 | 975 PRINT("log2_numvector_size",q->log2_numvector_size); |
2956 | 976 PRINT("numvector_size",q->numvector_size); |
977 PRINT("total_subbands",q->total_subbands); | |
978 } | |
979 #endif | |
3090 | 980 |
2956 | 981 /** |
982 * Cook initialization | |
983 * | |
984 * @param avctx pointer to the AVCodecContext | |
985 */ | |
986 | |
987 static int cook_decode_init(AVCodecContext *avctx) | |
988 { | |
989 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
|
990 uint8_t *edata_ptr = avctx->extradata; |
2956 | 991 |
992 /* Take care of the codec specific extradata. */ | |
993 if (avctx->extradata_size <= 0) { | |
4302
a0f83004d485
av_log(NULL,... -> av_log(avctx,.. where appropriate.
banan
parents:
3947
diff
changeset
|
994 av_log(avctx,AV_LOG_ERROR,"Necessary extradata missing!\n"); |
2956 | 995 return -1; |
996 } else { | |
997 /* 8 for mono, 16 for stereo, ? for multichannel | |
998 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
|
999 av_log(avctx,AV_LOG_DEBUG,"codecdata_length=%d\n",avctx->extradata_size); |
2956 | 1000 if (avctx->extradata_size >= 8){ |
4542 | 1001 q->cookversion = bytestream_get_be32(&edata_ptr); |
1002 q->samples_per_frame = bytestream_get_be16(&edata_ptr); | |
1003 q->subbands = bytestream_get_be16(&edata_ptr); | |
2956 | 1004 } |
1005 if (avctx->extradata_size >= 16){ | |
4542 | 1006 bytestream_get_be32(&edata_ptr); //Unknown unused |
1007 q->js_subband_start = bytestream_get_be16(&edata_ptr); | |
1008 q->js_vlc_bits = bytestream_get_be16(&edata_ptr); | |
2956 | 1009 } |
1010 } | |
1011 | |
1012 /* Take data from the AVCodecContext (RM container). */ | |
1013 q->sample_rate = avctx->sample_rate; | |
1014 q->nb_channels = avctx->channels; | |
1015 q->bit_rate = avctx->bit_rate; | |
1016 | |
4692
3318e3f6470f
Small simplifications to subband coefficient handling and use av_random().
banan
parents:
4674
diff
changeset
|
1017 /* Initialize RNG. */ |
3318e3f6470f
Small simplifications to subband coefficient handling and use av_random().
banan
parents:
4674
diff
changeset
|
1018 av_init_random(1, &q->random_state); |
2956 | 1019 |
1020 /* 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
|
1021 q->samples_per_channel = q->samples_per_frame / q->nb_channels; |
2956 | 1022 q->bits_per_subpacket = avctx->block_align * 8; |
1023 | |
1024 /* Initialize default data states. */ | |
3090 | 1025 q->log2_numvector_size = 5; |
2956 | 1026 q->total_subbands = q->subbands; |
1027 | |
1028 /* 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
|
1029 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
|
1030 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
|
1031 switch (q->cookversion) { |
4425 | 1032 case MONO: |
2956 | 1033 if (q->nb_channels != 1) { |
4302
a0f83004d485
av_log(NULL,... -> av_log(avctx,.. where appropriate.
banan
parents:
3947
diff
changeset
|
1034 av_log(avctx,AV_LOG_ERROR,"Container channels != 1, report sample!\n"); |
2956 | 1035 return -1; |
1036 } | |
4425 | 1037 av_log(avctx,AV_LOG_DEBUG,"MONO\n"); |
2956 | 1038 break; |
4425 | 1039 case STEREO: |
2956 | 1040 if (q->nb_channels != 1) { |
2959
24805f4d1b84
This patch adds some support for non-joint stereo streams. It also
rtognimp
parents:
2956
diff
changeset
|
1041 q->bits_per_subpacket = q->bits_per_subpacket/2; |
2956 | 1042 } |
4425 | 1043 av_log(avctx,AV_LOG_DEBUG,"STEREO\n"); |
2956 | 1044 break; |
1045 case JOINT_STEREO: | |
1046 if (q->nb_channels != 2) { | |
4302
a0f83004d485
av_log(NULL,... -> av_log(avctx,.. where appropriate.
banan
parents:
3947
diff
changeset
|
1047 av_log(avctx,AV_LOG_ERROR,"Container channels != 2, report sample!\n"); |
2956 | 1048 return -1; |
1049 } | |
4302
a0f83004d485
av_log(NULL,... -> av_log(avctx,.. where appropriate.
banan
parents:
3947
diff
changeset
|
1050 av_log(avctx,AV_LOG_DEBUG,"JOINT_STEREO\n"); |
2956 | 1051 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
|
1052 q->total_subbands = q->subbands + q->js_subband_start; |
2956 | 1053 q->joint_stereo = 1; |
1054 } | |
1055 if (q->samples_per_channel > 256) { | |
3091
0284d5b34916
Fix broken cosmetics commit and add a check for valid headers.
banan
parents:
3090
diff
changeset
|
1056 q->log2_numvector_size = 6; |
2956 | 1057 } |
1058 if (q->samples_per_channel > 512) { | |
3091
0284d5b34916
Fix broken cosmetics commit and add a check for valid headers.
banan
parents:
3090
diff
changeset
|
1059 q->log2_numvector_size = 7; |
2956 | 1060 } |
1061 break; | |
1062 case MC_COOK: | |
4302
a0f83004d485
av_log(NULL,... -> av_log(avctx,.. where appropriate.
banan
parents:
3947
diff
changeset
|
1063 av_log(avctx,AV_LOG_ERROR,"MC_COOK not supported!\n"); |
2956 | 1064 return -1; |
1065 break; | |
1066 default: | |
4302
a0f83004d485
av_log(NULL,... -> av_log(avctx,.. where appropriate.
banan
parents:
3947
diff
changeset
|
1067 av_log(avctx,AV_LOG_ERROR,"Unknown Cook version, report sample!\n"); |
2956 | 1068 return -1; |
1069 break; | |
1070 } | |
1071 | |
1072 /* Initialize variable relations */ | |
3090 | 1073 q->numvector_size = (1 << q->log2_numvector_size); |
2956 | 1074 |
1075 /* Generate tables */ | |
1076 init_rootpow2table(q); | |
1077 init_pow2table(q); | |
1078 init_gain_table(q); | |
1079 | |
1080 if (init_cook_vlc_tables(q) != 0) | |
1081 return -1; | |
1082 | |
3303
68721b62a528
sanity checks, some might have been exploitable ...
michael
parents:
3106
diff
changeset
|
1083 |
68721b62a528
sanity checks, some might have been exploitable ...
michael
parents:
3106
diff
changeset
|
1084 if(avctx->block_align >= UINT_MAX/2) |
68721b62a528
sanity checks, some might have been exploitable ...
michael
parents:
3106
diff
changeset
|
1085 return -1; |
68721b62a528
sanity checks, some might have been exploitable ...
michael
parents:
3106
diff
changeset
|
1086 |
4424
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
1087 /* Pad the databuffer with: |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
1088 DECODE_BYTES_PAD1 or DECODE_BYTES_PAD2 for decode_bytes(), |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
1089 FF_INPUT_BUFFER_PADDING_SIZE, for the bitstreamreader. */ |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
1090 if (q->nb_channels==2 && q->joint_stereo==0) { |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
1091 q->decoded_bytes_buffer = |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
1092 av_mallocz(avctx->block_align/2 |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
1093 + DECODE_BYTES_PAD2(avctx->block_align/2) |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
1094 + FF_INPUT_BUFFER_PADDING_SIZE); |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
1095 } else { |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
1096 q->decoded_bytes_buffer = |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
1097 av_mallocz(avctx->block_align |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
1098 + DECODE_BYTES_PAD1(avctx->block_align) |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
1099 + FF_INPUT_BUFFER_PADDING_SIZE); |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
1100 } |
8c830bde8006
Cook stereo (MONO_COOK2) bugfix, by Ian Braithwaite.
banan
parents:
4423
diff
changeset
|
1101 if (q->decoded_bytes_buffer == NULL) |
2956 | 1102 return -1; |
1103 | |
4639 | 1104 q->gains1.now = q->gain_1; |
1105 q->gains1.previous = q->gain_2; | |
1106 q->gains2.now = q->gain_3; | |
1107 q->gains2.previous = q->gain_4; | |
2956 | 1108 |
1109 /* Initialize transform. */ | |
4650
31bf54d9353d
Replace custom modified discrete cosine transform with ffmpeg's own.
banan
parents:
4639
diff
changeset
|
1110 if ( init_cook_mlt(q) != 0 ) |
2956 | 1111 return -1; |
3014
959b8ad880dc
Dual mono stereo strems sound ok now, added sanity checks and removed
rtognimp
parents:
3009
diff
changeset
|
1112 |
959b8ad880dc
Dual mono stereo strems sound ok now, added sanity checks and removed
rtognimp
parents:
3009
diff
changeset
|
1113 /* 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
|
1114 if (q->total_subbands > 53) { |
4302
a0f83004d485
av_log(NULL,... -> av_log(avctx,.. where appropriate.
banan
parents:
3947
diff
changeset
|
1115 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
|
1116 return -1; |
959b8ad880dc
Dual mono stereo strems sound ok now, added sanity checks and removed
rtognimp
parents:
3009
diff
changeset
|
1117 } |
959b8ad880dc
Dual mono stereo strems sound ok now, added sanity checks and removed
rtognimp
parents:
3009
diff
changeset
|
1118 if (q->subbands > 50) { |
4302
a0f83004d485
av_log(NULL,... -> av_log(avctx,.. where appropriate.
banan
parents:
3947
diff
changeset
|
1119 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
|
1120 return -1; |
959b8ad880dc
Dual mono stereo strems sound ok now, added sanity checks and removed
rtognimp
parents:
3009
diff
changeset
|
1121 } |
3091
0284d5b34916
Fix broken cosmetics commit and add a check for valid headers.
banan
parents:
3090
diff
changeset
|
1122 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
|
1123 } else { |
4302
a0f83004d485
av_log(NULL,... -> av_log(avctx,.. where appropriate.
banan
parents:
3947
diff
changeset
|
1124 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
|
1125 return -1; |
0284d5b34916
Fix broken cosmetics commit and add a check for valid headers.
banan
parents:
3090
diff
changeset
|
1126 } |
4431
85ac154efd99
Check that js_vlc_bits from the extradata is in a valid range.
banan
parents:
4430
diff
changeset
|
1127 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
|
1128 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
|
1129 return -1; |
85ac154efd99
Check that js_vlc_bits from the extradata is in a valid range.
banan
parents:
4430
diff
changeset
|
1130 } |
3014
959b8ad880dc
Dual mono stereo strems sound ok now, added sanity checks and removed
rtognimp
parents:
3009
diff
changeset
|
1131 |
3009
f5898b9b8a8a
Fix an out of array access and some minor cleanup of the code.
diego
parents:
2959
diff
changeset
|
1132 #ifdef COOKDEBUG |
4430
407e7fd9b4f4
Get rid of the COOKextradata struct. And use valid C to parse the extradata.
banan
parents:
4429
diff
changeset
|
1133 dump_cook_context(q); |
3009
f5898b9b8a8a
Fix an out of array access and some minor cleanup of the code.
diego
parents:
2959
diff
changeset
|
1134 #endif |
2956 | 1135 return 0; |
1136 } | |
1137 | |
1138 | |
1139 AVCodec cook_decoder = | |
1140 { | |
1141 .name = "cook", | |
1142 .type = CODEC_TYPE_AUDIO, | |
1143 .id = CODEC_ID_COOK, | |
1144 .priv_data_size = sizeof(COOKContext), | |
1145 .init = cook_decode_init, | |
1146 .close = cook_decode_close, | |
1147 .decode = cook_decode_frame, | |
1148 }; |