Go to most recent revision | Details | Last modification | View Log | RSS feed
Rev | Author | Line No. | Line |
---|---|---|---|
4349 | Serge | 1 | /** |
2 | * @file |
||
3 | * Vorbis I decoder |
||
4 | * @author Denes Balatoni ( dbalatoni programozo hu ) |
||
5 | * |
||
6 | * This file is part of FFmpeg. |
||
7 | * |
||
8 | * FFmpeg is free software; you can redistribute it and/or |
||
9 | * modify it under the terms of the GNU Lesser General Public |
||
10 | * License as published by the Free Software Foundation; either |
||
11 | * version 2.1 of the License, or (at your option) any later version. |
||
12 | * |
||
13 | * FFmpeg is distributed in the hope that it will be useful, |
||
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 |
||
19 | * License along with FFmpeg; if not, write to the Free Software |
||
20 | * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA |
||
21 | */ |
||
22 | |||
23 | /** |
||
24 | * @file |
||
25 | * Vorbis I decoder |
||
26 | * @author Denes Balatoni ( dbalatoni programozo hu ) |
||
27 | */ |
||
28 | |||
29 | #include |
||
30 | #include |
||
31 | |||
32 | #define BITSTREAM_READER_LE |
||
33 | #include "libavutil/float_dsp.h" |
||
34 | #include "libavutil/avassert.h" |
||
35 | #include "avcodec.h" |
||
36 | #include "get_bits.h" |
||
37 | #include "fft.h" |
||
38 | #include "fmtconvert.h" |
||
39 | #include "internal.h" |
||
40 | |||
41 | #include "vorbis.h" |
||
42 | #include "vorbisdsp.h" |
||
43 | #include "xiph.h" |
||
44 | |||
45 | #define V_NB_BITS 8 |
||
46 | #define V_NB_BITS2 11 |
||
47 | #define V_MAX_VLCS (1 << 16) |
||
48 | #define V_MAX_PARTITIONS (1 << 20) |
||
49 | |||
50 | typedef struct { |
||
51 | uint8_t dimensions; |
||
52 | uint8_t lookup_type; |
||
53 | uint8_t maxdepth; |
||
54 | VLC vlc; |
||
55 | float *codevectors; |
||
56 | unsigned int nb_bits; |
||
57 | } vorbis_codebook; |
||
58 | |||
59 | typedef union vorbis_floor_u vorbis_floor_data; |
||
60 | typedef struct vorbis_floor0_s vorbis_floor0; |
||
61 | typedef struct vorbis_floor1_s vorbis_floor1; |
||
62 | struct vorbis_context_s; |
||
63 | typedef |
||
64 | int (* vorbis_floor_decode_func) |
||
65 | (struct vorbis_context_s *, vorbis_floor_data *, float *); |
||
66 | typedef struct { |
||
67 | uint8_t floor_type; |
||
68 | vorbis_floor_decode_func decode; |
||
69 | union vorbis_floor_u { |
||
70 | struct vorbis_floor0_s { |
||
71 | uint8_t order; |
||
72 | uint16_t rate; |
||
73 | uint16_t bark_map_size; |
||
74 | int32_t *map[2]; |
||
75 | uint32_t map_size[2]; |
||
76 | uint8_t amplitude_bits; |
||
77 | uint8_t amplitude_offset; |
||
78 | uint8_t num_books; |
||
79 | uint8_t *book_list; |
||
80 | float *lsp; |
||
81 | } t0; |
||
82 | struct vorbis_floor1_s { |
||
83 | uint8_t partitions; |
||
84 | uint8_t partition_class[32]; |
||
85 | uint8_t class_dimensions[16]; |
||
86 | uint8_t class_subclasses[16]; |
||
87 | uint8_t class_masterbook[16]; |
||
88 | int16_t subclass_books[16][8]; |
||
89 | uint8_t multiplier; |
||
90 | uint16_t x_list_dim; |
||
91 | vorbis_floor1_entry *list; |
||
92 | } t1; |
||
93 | } data; |
||
94 | } vorbis_floor; |
||
95 | |||
96 | typedef struct { |
||
97 | uint16_t type; |
||
98 | uint32_t begin; |
||
99 | uint32_t end; |
||
100 | unsigned partition_size; |
||
101 | uint8_t classifications; |
||
102 | uint8_t classbook; |
||
103 | int16_t books[64][8]; |
||
104 | uint8_t maxpass; |
||
105 | uint16_t ptns_to_read; |
||
106 | uint8_t *classifs; |
||
107 | } vorbis_residue; |
||
108 | |||
109 | typedef struct { |
||
110 | uint8_t submaps; |
||
111 | uint16_t coupling_steps; |
||
112 | uint8_t *magnitude; |
||
113 | uint8_t *angle; |
||
114 | uint8_t *mux; |
||
115 | uint8_t submap_floor[16]; |
||
116 | uint8_t submap_residue[16]; |
||
117 | } vorbis_mapping; |
||
118 | |||
119 | typedef struct { |
||
120 | uint8_t blockflag; |
||
121 | uint16_t windowtype; |
||
122 | uint16_t transformtype; |
||
123 | uint8_t mapping; |
||
124 | } vorbis_mode; |
||
125 | |||
126 | typedef struct vorbis_context_s { |
||
127 | AVCodecContext *avctx; |
||
128 | GetBitContext gb; |
||
129 | VorbisDSPContext dsp; |
||
130 | AVFloatDSPContext fdsp; |
||
131 | FmtConvertContext fmt_conv; |
||
132 | |||
133 | FFTContext mdct[2]; |
||
134 | uint8_t first_frame; |
||
135 | uint32_t version; |
||
136 | uint8_t audio_channels; |
||
137 | uint32_t audio_samplerate; |
||
138 | uint32_t bitrate_maximum; |
||
139 | uint32_t bitrate_nominal; |
||
140 | uint32_t bitrate_minimum; |
||
141 | uint32_t blocksize[2]; |
||
142 | const float *win[2]; |
||
143 | uint16_t codebook_count; |
||
144 | vorbis_codebook *codebooks; |
||
145 | uint8_t floor_count; |
||
146 | vorbis_floor *floors; |
||
147 | uint8_t residue_count; |
||
148 | vorbis_residue *residues; |
||
149 | uint8_t mapping_count; |
||
150 | vorbis_mapping *mappings; |
||
151 | uint8_t mode_count; |
||
152 | vorbis_mode *modes; |
||
153 | uint8_t mode_number; // mode number for the current packet |
||
154 | uint8_t previous_window; |
||
155 | float *channel_residues; |
||
156 | float *saved; |
||
157 | } vorbis_context; |
||
158 | |||
159 | /* Helper functions */ |
||
160 | |||
161 | #define BARK(x) \ |
||
162 | (13.1f * atan(0.00074f * (x)) + 2.24f * atan(1.85e-8f * (x) * (x)) + 1e-4f * (x)) |
||
163 | |||
164 | static const char idx_err_str[] = "Index value %d out of range (0 - %d) for %s at %s:%i\n"; |
||
165 | #define VALIDATE_INDEX(idx, limit) \ |
||
166 | if (idx >= limit) {\ |
||
167 | av_log(vc->avctx, AV_LOG_ERROR,\ |
||
168 | idx_err_str,\ |
||
169 | (int)(idx), (int)(limit - 1), #idx, __FILE__, __LINE__);\ |
||
170 | return AVERROR_INVALIDDATA;\ |
||
171 | } |
||
172 | #define GET_VALIDATED_INDEX(idx, bits, limit) \ |
||
173 | {\ |
||
174 | idx = get_bits(gb, bits);\ |
||
175 | VALIDATE_INDEX(idx, limit)\ |
||
176 | } |
||
177 | |||
178 | static float vorbisfloat2float(unsigned val) |
||
179 | { |
||
180 | double mant = val & 0x1fffff; |
||
181 | long exp = (val & 0x7fe00000L) >> 21; |
||
182 | if (val & 0x80000000) |
||
183 | mant = -mant; |
||
184 | return ldexp(mant, exp - 20 - 768); |
||
185 | } |
||
186 | |||
187 | |||
188 | // Free all allocated memory ----------------------------------------- |
||
189 | |||
190 | static void vorbis_free(vorbis_context *vc) |
||
191 | { |
||
192 | int i; |
||
193 | |||
194 | av_freep(&vc->channel_residues); |
||
195 | av_freep(&vc->saved); |
||
196 | |||
197 | if (vc->residues) |
||
198 | for (i = 0; i < vc->residue_count; i++) |
||
199 | av_free(vc->residues[i].classifs); |
||
200 | av_freep(&vc->residues); |
||
201 | av_freep(&vc->modes); |
||
202 | |||
203 | ff_mdct_end(&vc->mdct[0]); |
||
204 | ff_mdct_end(&vc->mdct[1]); |
||
205 | |||
206 | if (vc->codebooks) |
||
207 | for (i = 0; i < vc->codebook_count; ++i) { |
||
208 | av_free(vc->codebooks[i].codevectors); |
||
209 | ff_free_vlc(&vc->codebooks[i].vlc); |
||
210 | } |
||
211 | av_freep(&vc->codebooks); |
||
212 | |||
213 | if (vc->floors) |
||
214 | for (i = 0; i < vc->floor_count; ++i) { |
||
215 | if (vc->floors[i].floor_type == 0) { |
||
216 | av_free(vc->floors[i].data.t0.map[0]); |
||
217 | av_free(vc->floors[i].data.t0.map[1]); |
||
218 | av_free(vc->floors[i].data.t0.book_list); |
||
219 | av_free(vc->floors[i].data.t0.lsp); |
||
220 | } else { |
||
221 | av_free(vc->floors[i].data.t1.list); |
||
222 | } |
||
223 | } |
||
224 | av_freep(&vc->floors); |
||
225 | |||
226 | if (vc->mappings) |
||
227 | for (i = 0; i < vc->mapping_count; ++i) { |
||
228 | av_free(vc->mappings[i].magnitude); |
||
229 | av_free(vc->mappings[i].angle); |
||
230 | av_free(vc->mappings[i].mux); |
||
231 | } |
||
232 | av_freep(&vc->mappings); |
||
233 | } |
||
234 | |||
235 | // Parse setup header ------------------------------------------------- |
||
236 | |||
237 | // Process codebooks part |
||
238 | |||
239 | static int vorbis_parse_setup_hdr_codebooks(vorbis_context *vc) |
||
240 | { |
||
241 | unsigned cb; |
||
242 | uint8_t *tmp_vlc_bits = NULL; |
||
243 | uint32_t *tmp_vlc_codes = NULL; |
||
244 | GetBitContext *gb = &vc->gb; |
||
245 | uint16_t *codebook_multiplicands = NULL; |
||
246 | int ret = 0; |
||
247 | |||
248 | vc->codebook_count = get_bits(gb, 8) + 1; |
||
249 | |||
250 | av_dlog(NULL, " Codebooks: %d \n", vc->codebook_count); |
||
251 | |||
252 | vc->codebooks = av_mallocz(vc->codebook_count * sizeof(*vc->codebooks)); |
||
253 | tmp_vlc_bits = av_mallocz(V_MAX_VLCS * sizeof(*tmp_vlc_bits)); |
||
254 | tmp_vlc_codes = av_mallocz(V_MAX_VLCS * sizeof(*tmp_vlc_codes)); |
||
255 | codebook_multiplicands = av_malloc(V_MAX_VLCS * sizeof(*codebook_multiplicands)); |
||
256 | if (!vc->codebooks || |
||
257 | !tmp_vlc_bits || !tmp_vlc_codes || !codebook_multiplicands) { |
||
258 | ret = AVERROR(ENOMEM); |
||
259 | goto error; |
||
260 | } |
||
261 | |||
262 | for (cb = 0; cb < vc->codebook_count; ++cb) { |
||
263 | vorbis_codebook *codebook_setup = &vc->codebooks[cb]; |
||
264 | unsigned ordered, t, entries, used_entries = 0; |
||
265 | |||
266 | av_dlog(NULL, " %u. Codebook\n", cb); |
||
267 | |||
268 | if (get_bits(gb, 24) != 0x564342) { |
||
269 | av_log(vc->avctx, AV_LOG_ERROR, |
||
270 | " %u. Codebook setup data corrupt.\n", cb); |
||
271 | ret = AVERROR_INVALIDDATA; |
||
272 | goto error; |
||
273 | } |
||
274 | |||
275 | codebook_setup->dimensions=get_bits(gb, 16); |
||
276 | if (codebook_setup->dimensions > 16 || codebook_setup->dimensions == 0) { |
||
277 | av_log(vc->avctx, AV_LOG_ERROR, |
||
278 | " %u. Codebook's dimension is invalid (%d).\n", |
||
279 | cb, codebook_setup->dimensions); |
||
280 | ret = AVERROR_INVALIDDATA; |
||
281 | goto error; |
||
282 | } |
||
283 | entries = get_bits(gb, 24); |
||
284 | if (entries > V_MAX_VLCS) { |
||
285 | av_log(vc->avctx, AV_LOG_ERROR, |
||
286 | " %u. Codebook has too many entries (%u).\n", |
||
287 | cb, entries); |
||
288 | ret = AVERROR_INVALIDDATA; |
||
289 | goto error; |
||
290 | } |
||
291 | |||
292 | ordered = get_bits1(gb); |
||
293 | |||
294 | av_dlog(NULL, " codebook_dimensions %d, codebook_entries %u\n", |
||
295 | codebook_setup->dimensions, entries); |
||
296 | |||
297 | if (!ordered) { |
||
298 | unsigned ce, flag; |
||
299 | unsigned sparse = get_bits1(gb); |
||
300 | |||
301 | av_dlog(NULL, " not ordered \n"); |
||
302 | |||
303 | if (sparse) { |
||
304 | av_dlog(NULL, " sparse \n"); |
||
305 | |||
306 | used_entries = 0; |
||
307 | for (ce = 0; ce < entries; ++ce) { |
||
308 | flag = get_bits1(gb); |
||
309 | if (flag) { |
||
310 | tmp_vlc_bits[ce] = get_bits(gb, 5) + 1; |
||
311 | ++used_entries; |
||
312 | } else |
||
313 | tmp_vlc_bits[ce] = 0; |
||
314 | } |
||
315 | } else { |
||
316 | av_dlog(NULL, " not sparse \n"); |
||
317 | |||
318 | used_entries = entries; |
||
319 | for (ce = 0; ce < entries; ++ce) |
||
320 | tmp_vlc_bits[ce] = get_bits(gb, 5) + 1; |
||
321 | } |
||
322 | } else { |
||
323 | unsigned current_entry = 0; |
||
324 | unsigned current_length = get_bits(gb, 5) + 1; |
||
325 | |||
326 | av_dlog(NULL, " ordered, current length: %u\n", current_length); //FIXME |
||
327 | |||
328 | used_entries = entries; |
||
329 | for (; current_entry < used_entries && current_length <= 32; ++current_length) { |
||
330 | unsigned i, number; |
||
331 | |||
332 | av_dlog(NULL, " number bits: %u ", ilog(entries - current_entry)); |
||
333 | |||
334 | number = get_bits(gb, ilog(entries - current_entry)); |
||
335 | |||
336 | av_dlog(NULL, " number: %u\n", number); |
||
337 | |||
338 | for (i = current_entry; i < number+current_entry; ++i) |
||
339 | if (i < used_entries) |
||
340 | tmp_vlc_bits[i] = current_length; |
||
341 | |||
342 | current_entry+=number; |
||
343 | } |
||
344 | if (current_entry>used_entries) { |
||
345 | av_log(vc->avctx, AV_LOG_ERROR, " More codelengths than codes in codebook. \n"); |
||
346 | ret = AVERROR_INVALIDDATA; |
||
347 | goto error; |
||
348 | } |
||
349 | } |
||
350 | |||
351 | codebook_setup->lookup_type = get_bits(gb, 4); |
||
352 | |||
353 | av_dlog(NULL, " lookup type: %d : %s \n", codebook_setup->lookup_type, |
||
354 | codebook_setup->lookup_type ? "vq" : "no lookup"); |
||
355 | |||
356 | // If the codebook is used for (inverse) VQ, calculate codevectors. |
||
357 | |||
358 | if (codebook_setup->lookup_type == 1) { |
||
359 | unsigned i, j, k; |
||
360 | unsigned codebook_lookup_values = ff_vorbis_nth_root(entries, codebook_setup->dimensions); |
||
361 | |||
362 | float codebook_minimum_value = vorbisfloat2float(get_bits_long(gb, 32)); |
||
363 | float codebook_delta_value = vorbisfloat2float(get_bits_long(gb, 32)); |
||
364 | unsigned codebook_value_bits = get_bits(gb, 4) + 1; |
||
365 | unsigned codebook_sequence_p = get_bits1(gb); |
||
366 | |||
367 | av_dlog(NULL, " We expect %d numbers for building the codevectors. \n", |
||
368 | codebook_lookup_values); |
||
369 | av_dlog(NULL, " delta %f minmum %f \n", |
||
370 | codebook_delta_value, codebook_minimum_value); |
||
371 | |||
372 | for (i = 0; i < codebook_lookup_values; ++i) { |
||
373 | codebook_multiplicands[i] = get_bits(gb, codebook_value_bits); |
||
374 | |||
375 | av_dlog(NULL, " multiplicands*delta+minmum : %e \n", |
||
376 | (float)codebook_multiplicands[i] * codebook_delta_value + codebook_minimum_value); |
||
377 | av_dlog(NULL, " multiplicand %u\n", codebook_multiplicands[i]); |
||
378 | } |
||
379 | |||
380 | // Weed out unused vlcs and build codevector vector |
||
381 | codebook_setup->codevectors = used_entries ? av_mallocz(used_entries * |
||
382 | codebook_setup->dimensions * |
||
383 | sizeof(*codebook_setup->codevectors)) |
||
384 | : NULL; |
||
385 | for (j = 0, i = 0; i < entries; ++i) { |
||
386 | unsigned dim = codebook_setup->dimensions; |
||
387 | |||
388 | if (tmp_vlc_bits[i]) { |
||
389 | float last = 0.0; |
||
390 | unsigned lookup_offset = i; |
||
391 | |||
392 | av_dlog(vc->avctx, "Lookup offset %u ,", i); |
||
393 | |||
394 | for (k = 0; k < dim; ++k) { |
||
395 | unsigned multiplicand_offset = lookup_offset % codebook_lookup_values; |
||
396 | codebook_setup->codevectors[j * dim + k] = codebook_multiplicands[multiplicand_offset] * codebook_delta_value + codebook_minimum_value + last; |
||
397 | if (codebook_sequence_p) |
||
398 | last = codebook_setup->codevectors[j * dim + k]; |
||
399 | lookup_offset/=codebook_lookup_values; |
||
400 | } |
||
401 | tmp_vlc_bits[j] = tmp_vlc_bits[i]; |
||
402 | |||
403 | av_dlog(vc->avctx, "real lookup offset %u, vector: ", j); |
||
404 | for (k = 0; k < dim; ++k) |
||
405 | av_dlog(vc->avctx, " %f ", |
||
406 | codebook_setup->codevectors[j * dim + k]); |
||
407 | av_dlog(vc->avctx, "\n"); |
||
408 | |||
409 | ++j; |
||
410 | } |
||
411 | } |
||
412 | if (j != used_entries) { |
||
413 | av_log(vc->avctx, AV_LOG_ERROR, "Bug in codevector vector building code. \n"); |
||
414 | ret = AVERROR_INVALIDDATA; |
||
415 | goto error; |
||
416 | } |
||
417 | entries = used_entries; |
||
418 | } else if (codebook_setup->lookup_type >= 2) { |
||
419 | av_log(vc->avctx, AV_LOG_ERROR, "Codebook lookup type not supported. \n"); |
||
420 | ret = AVERROR_INVALIDDATA; |
||
421 | goto error; |
||
422 | } |
||
423 | |||
424 | // Initialize VLC table |
||
425 | if (ff_vorbis_len2vlc(tmp_vlc_bits, tmp_vlc_codes, entries)) { |
||
426 | av_log(vc->avctx, AV_LOG_ERROR, " Invalid code lengths while generating vlcs. \n"); |
||
427 | ret = AVERROR_INVALIDDATA; |
||
428 | goto error; |
||
429 | } |
||
430 | codebook_setup->maxdepth = 0; |
||
431 | for (t = 0; t < entries; ++t) |
||
432 | if (tmp_vlc_bits[t] >= codebook_setup->maxdepth) |
||
433 | codebook_setup->maxdepth = tmp_vlc_bits[t]; |
||
434 | |||
435 | if (codebook_setup->maxdepth > 3 * V_NB_BITS) |
||
436 | codebook_setup->nb_bits = V_NB_BITS2; |
||
437 | else |
||
438 | codebook_setup->nb_bits = V_NB_BITS; |
||
439 | |||
440 | codebook_setup->maxdepth = (codebook_setup->maxdepth+codebook_setup->nb_bits - 1) / codebook_setup->nb_bits; |
||
441 | |||
442 | if ((ret = init_vlc(&codebook_setup->vlc, codebook_setup->nb_bits, |
||
443 | entries, tmp_vlc_bits, sizeof(*tmp_vlc_bits), |
||
444 | sizeof(*tmp_vlc_bits), tmp_vlc_codes, |
||
445 | sizeof(*tmp_vlc_codes), sizeof(*tmp_vlc_codes), |
||
446 | INIT_VLC_LE))) { |
||
447 | av_log(vc->avctx, AV_LOG_ERROR, " Error generating vlc tables. \n"); |
||
448 | goto error; |
||
449 | } |
||
450 | } |
||
451 | |||
452 | av_free(tmp_vlc_bits); |
||
453 | av_free(tmp_vlc_codes); |
||
454 | av_free(codebook_multiplicands); |
||
455 | return 0; |
||
456 | |||
457 | // Error: |
||
458 | error: |
||
459 | av_free(tmp_vlc_bits); |
||
460 | av_free(tmp_vlc_codes); |
||
461 | av_free(codebook_multiplicands); |
||
462 | return ret; |
||
463 | } |
||
464 | |||
465 | // Process time domain transforms part (unused in Vorbis I) |
||
466 | |||
467 | static int vorbis_parse_setup_hdr_tdtransforms(vorbis_context *vc) |
||
468 | { |
||
469 | GetBitContext *gb = &vc->gb; |
||
470 | unsigned i, vorbis_time_count = get_bits(gb, 6) + 1; |
||
471 | |||
472 | for (i = 0; i < vorbis_time_count; ++i) { |
||
473 | unsigned vorbis_tdtransform = get_bits(gb, 16); |
||
474 | |||
475 | av_dlog(NULL, " Vorbis time domain transform %u: %u\n", |
||
476 | vorbis_time_count, vorbis_tdtransform); |
||
477 | |||
478 | if (vorbis_tdtransform) { |
||
479 | av_log(vc->avctx, AV_LOG_ERROR, "Vorbis time domain transform data nonzero. \n"); |
||
480 | return AVERROR_INVALIDDATA; |
||
481 | } |
||
482 | } |
||
483 | return 0; |
||
484 | } |
||
485 | |||
486 | // Process floors part |
||
487 | |||
488 | static int vorbis_floor0_decode(vorbis_context *vc, |
||
489 | vorbis_floor_data *vfu, float *vec); |
||
490 | static int create_map(vorbis_context *vc, unsigned floor_number); |
||
491 | static int vorbis_floor1_decode(vorbis_context *vc, |
||
492 | vorbis_floor_data *vfu, float *vec); |
||
493 | static int vorbis_parse_setup_hdr_floors(vorbis_context *vc) |
||
494 | { |
||
495 | GetBitContext *gb = &vc->gb; |
||
496 | int i, j, k, ret; |
||
497 | |||
498 | vc->floor_count = get_bits(gb, 6) + 1; |
||
499 | |||
500 | vc->floors = av_mallocz(vc->floor_count * sizeof(*vc->floors)); |
||
501 | if (!vc->floors) |
||
502 | return AVERROR(ENOMEM); |
||
503 | |||
504 | for (i = 0; i < vc->floor_count; ++i) { |
||
505 | vorbis_floor *floor_setup = &vc->floors[i]; |
||
506 | |||
507 | floor_setup->floor_type = get_bits(gb, 16); |
||
508 | |||
509 | av_dlog(NULL, " %d. floor type %d \n", i, floor_setup->floor_type); |
||
510 | |||
511 | if (floor_setup->floor_type == 1) { |
||
512 | int maximum_class = -1; |
||
513 | unsigned rangebits, rangemax, floor1_values = 2; |
||
514 | |||
515 | floor_setup->decode = vorbis_floor1_decode; |
||
516 | |||
517 | floor_setup->data.t1.partitions = get_bits(gb, 5); |
||
518 | |||
519 | av_dlog(NULL, " %d.floor: %d partitions \n", |
||
520 | i, floor_setup->data.t1.partitions); |
||
521 | |||
522 | for (j = 0; j < floor_setup->data.t1.partitions; ++j) { |
||
523 | floor_setup->data.t1.partition_class[j] = get_bits(gb, 4); |
||
524 | if (floor_setup->data.t1.partition_class[j] > maximum_class) |
||
525 | maximum_class = floor_setup->data.t1.partition_class[j]; |
||
526 | |||
527 | av_dlog(NULL, " %d. floor %d partition class %d \n", |
||
528 | i, j, floor_setup->data.t1.partition_class[j]); |
||
529 | |||
530 | } |
||
531 | |||
532 | av_dlog(NULL, " maximum class %d \n", maximum_class); |
||
533 | |||
534 | for (j = 0; j <= maximum_class; ++j) { |
||
535 | floor_setup->data.t1.class_dimensions[j] = get_bits(gb, 3) + 1; |
||
536 | floor_setup->data.t1.class_subclasses[j] = get_bits(gb, 2); |
||
537 | |||
538 | av_dlog(NULL, " %d floor %d class dim: %d subclasses %d \n", i, j, |
||
539 | floor_setup->data.t1.class_dimensions[j], |
||
540 | floor_setup->data.t1.class_subclasses[j]); |
||
541 | |||
542 | if (floor_setup->data.t1.class_subclasses[j]) { |
||
543 | GET_VALIDATED_INDEX(floor_setup->data.t1.class_masterbook[j], 8, vc->codebook_count) |
||
544 | |||
545 | av_dlog(NULL, " masterbook: %d \n", floor_setup->data.t1.class_masterbook[j]); |
||
546 | } |
||
547 | |||
548 | for (k = 0; k < (1 << floor_setup->data.t1.class_subclasses[j]); ++k) { |
||
549 | int16_t bits = get_bits(gb, 8) - 1; |
||
550 | if (bits != -1) |
||
551 | VALIDATE_INDEX(bits, vc->codebook_count) |
||
552 | floor_setup->data.t1.subclass_books[j][k] = bits; |
||
553 | |||
554 | av_dlog(NULL, " book %d. : %d \n", k, floor_setup->data.t1.subclass_books[j][k]); |
||
555 | } |
||
556 | } |
||
557 | |||
558 | floor_setup->data.t1.multiplier = get_bits(gb, 2) + 1; |
||
559 | floor_setup->data.t1.x_list_dim = 2; |
||
560 | |||
561 | for (j = 0; j < floor_setup->data.t1.partitions; ++j) |
||
562 | floor_setup->data.t1.x_list_dim+=floor_setup->data.t1.class_dimensions[floor_setup->data.t1.partition_class[j]]; |
||
563 | |||
564 | floor_setup->data.t1.list = av_mallocz(floor_setup->data.t1.x_list_dim * |
||
565 | sizeof(*floor_setup->data.t1.list)); |
||
566 | if (!floor_setup->data.t1.list) |
||
567 | return AVERROR(ENOMEM); |
||
568 | |||
569 | rangebits = get_bits(gb, 4); |
||
570 | rangemax = (1 << rangebits); |
||
571 | if (rangemax > vc->blocksize[1] / 2) { |
||
572 | av_log(vc->avctx, AV_LOG_ERROR, |
||
573 | "Floor value is too large for blocksize: %u (%"PRIu32")\n", |
||
574 | rangemax, vc->blocksize[1] / 2); |
||
575 | return AVERROR_INVALIDDATA; |
||
576 | } |
||
577 | floor_setup->data.t1.list[0].x = 0; |
||
578 | floor_setup->data.t1.list[1].x = rangemax; |
||
579 | |||
580 | for (j = 0; j < floor_setup->data.t1.partitions; ++j) { |
||
581 | for (k = 0; k < floor_setup->data.t1.class_dimensions[floor_setup->data.t1.partition_class[j]]; ++k, ++floor1_values) { |
||
582 | floor_setup->data.t1.list[floor1_values].x = get_bits(gb, rangebits); |
||
583 | |||
584 | av_dlog(NULL, " %u. floor1 Y coord. %d\n", floor1_values, |
||
585 | floor_setup->data.t1.list[floor1_values].x); |
||
586 | } |
||
587 | } |
||
588 | |||
589 | // Precalculate order of x coordinates - needed for decode |
||
590 | if (ff_vorbis_ready_floor1_list(vc->avctx, |
||
591 | floor_setup->data.t1.list, |
||
592 | floor_setup->data.t1.x_list_dim)) { |
||
593 | return AVERROR_INVALIDDATA; |
||
594 | } |
||
595 | } else if (floor_setup->floor_type == 0) { |
||
596 | unsigned max_codebook_dim = 0; |
||
597 | |||
598 | floor_setup->decode = vorbis_floor0_decode; |
||
599 | |||
600 | floor_setup->data.t0.order = get_bits(gb, 8); |
||
601 | if (!floor_setup->data.t0.order) { |
||
602 | av_log(vc->avctx, AV_LOG_ERROR, "Floor 0 order is 0.\n"); |
||
603 | return AVERROR_INVALIDDATA; |
||
604 | } |
||
605 | floor_setup->data.t0.rate = get_bits(gb, 16); |
||
606 | if (!floor_setup->data.t0.rate) { |
||
607 | av_log(vc->avctx, AV_LOG_ERROR, "Floor 0 rate is 0.\n"); |
||
608 | return AVERROR_INVALIDDATA; |
||
609 | } |
||
610 | floor_setup->data.t0.bark_map_size = get_bits(gb, 16); |
||
611 | if (!floor_setup->data.t0.bark_map_size) { |
||
612 | av_log(vc->avctx, AV_LOG_ERROR, |
||
613 | "Floor 0 bark map size is 0.\n"); |
||
614 | return AVERROR_INVALIDDATA; |
||
615 | } |
||
616 | floor_setup->data.t0.amplitude_bits = get_bits(gb, 6); |
||
617 | floor_setup->data.t0.amplitude_offset = get_bits(gb, 8); |
||
618 | floor_setup->data.t0.num_books = get_bits(gb, 4) + 1; |
||
619 | |||
620 | /* allocate mem for booklist */ |
||
621 | floor_setup->data.t0.book_list = |
||
622 | av_malloc(floor_setup->data.t0.num_books); |
||
623 | if (!floor_setup->data.t0.book_list) |
||
624 | return AVERROR(ENOMEM); |
||
625 | /* read book indexes */ |
||
626 | { |
||
627 | int idx; |
||
628 | unsigned book_idx; |
||
629 | for (idx = 0; idx < floor_setup->data.t0.num_books; ++idx) { |
||
630 | GET_VALIDATED_INDEX(book_idx, 8, vc->codebook_count) |
||
631 | floor_setup->data.t0.book_list[idx] = book_idx; |
||
632 | if (vc->codebooks[book_idx].dimensions > max_codebook_dim) |
||
633 | max_codebook_dim = vc->codebooks[book_idx].dimensions; |
||
634 | } |
||
635 | } |
||
636 | |||
637 | if ((ret = create_map(vc, i)) < 0) |
||
638 | return ret; |
||
639 | |||
640 | /* codebook dim is for padding if codebook dim doesn't * |
||
641 | * divide order+1 then we need to read more data */ |
||
642 | floor_setup->data.t0.lsp = |
||
643 | av_malloc((floor_setup->data.t0.order + 1 + max_codebook_dim) |
||
644 | * sizeof(*floor_setup->data.t0.lsp)); |
||
645 | if (!floor_setup->data.t0.lsp) |
||
646 | return AVERROR(ENOMEM); |
||
647 | |||
648 | /* debug output parsed headers */ |
||
649 | av_dlog(NULL, "floor0 order: %u\n", floor_setup->data.t0.order); |
||
650 | av_dlog(NULL, "floor0 rate: %u\n", floor_setup->data.t0.rate); |
||
651 | av_dlog(NULL, "floor0 bark map size: %u\n", |
||
652 | floor_setup->data.t0.bark_map_size); |
||
653 | av_dlog(NULL, "floor0 amplitude bits: %u\n", |
||
654 | floor_setup->data.t0.amplitude_bits); |
||
655 | av_dlog(NULL, "floor0 amplitude offset: %u\n", |
||
656 | floor_setup->data.t0.amplitude_offset); |
||
657 | av_dlog(NULL, "floor0 number of books: %u\n", |
||
658 | floor_setup->data.t0.num_books); |
||
659 | av_dlog(NULL, "floor0 book list pointer: %p\n", |
||
660 | floor_setup->data.t0.book_list); |
||
661 | { |
||
662 | int idx; |
||
663 | for (idx = 0; idx < floor_setup->data.t0.num_books; ++idx) { |
||
664 | av_dlog(NULL, " Book %d: %u\n", idx + 1, |
||
665 | floor_setup->data.t0.book_list[idx]); |
||
666 | } |
||
667 | } |
||
668 | } else { |
||
669 | av_log(vc->avctx, AV_LOG_ERROR, "Invalid floor type!\n"); |
||
670 | return AVERROR_INVALIDDATA; |
||
671 | } |
||
672 | } |
||
673 | return 0; |
||
674 | } |
||
675 | |||
676 | // Process residues part |
||
677 | |||
678 | static int vorbis_parse_setup_hdr_residues(vorbis_context *vc) |
||
679 | { |
||
680 | GetBitContext *gb = &vc->gb; |
||
681 | unsigned i, j, k; |
||
682 | |||
683 | vc->residue_count = get_bits(gb, 6)+1; |
||
684 | vc->residues = av_mallocz(vc->residue_count * sizeof(*vc->residues)); |
||
685 | if (!vc->residues) |
||
686 | return AVERROR(ENOMEM); |
||
687 | |||
688 | av_dlog(NULL, " There are %d residues. \n", vc->residue_count); |
||
689 | |||
690 | for (i = 0; i < vc->residue_count; ++i) { |
||
691 | vorbis_residue *res_setup = &vc->residues[i]; |
||
692 | uint8_t cascade[64]; |
||
693 | unsigned high_bits, low_bits; |
||
694 | |||
695 | res_setup->type = get_bits(gb, 16); |
||
696 | |||
697 | av_dlog(NULL, " %u. residue type %d\n", i, res_setup->type); |
||
698 | |||
699 | res_setup->begin = get_bits(gb, 24); |
||
700 | res_setup->end = get_bits(gb, 24); |
||
701 | res_setup->partition_size = get_bits(gb, 24) + 1; |
||
702 | /* Validations to prevent a buffer overflow later. */ |
||
703 | if (res_setup->begin>res_setup->end || |
||
704 | res_setup->end > (res_setup->type == 2 ? vc->audio_channels : 1) * vc->blocksize[1] / 2 || |
||
705 | (res_setup->end-res_setup->begin) / res_setup->partition_size > V_MAX_PARTITIONS) { |
||
706 | av_log(vc->avctx, AV_LOG_ERROR, |
||
707 | "partition out of bounds: type, begin, end, size, blocksize: %"PRIu16", %"PRIu32", %"PRIu32", %u, %"PRIu32"\n", |
||
708 | res_setup->type, res_setup->begin, res_setup->end, |
||
709 | res_setup->partition_size, vc->blocksize[1] / 2); |
||
710 | return AVERROR_INVALIDDATA; |
||
711 | } |
||
712 | |||
713 | res_setup->classifications = get_bits(gb, 6) + 1; |
||
714 | GET_VALIDATED_INDEX(res_setup->classbook, 8, vc->codebook_count) |
||
715 | |||
716 | res_setup->ptns_to_read = |
||
717 | (res_setup->end - res_setup->begin) / res_setup->partition_size; |
||
718 | res_setup->classifs = av_malloc(res_setup->ptns_to_read * |
||
719 | vc->audio_channels * |
||
720 | sizeof(*res_setup->classifs)); |
||
721 | if (!res_setup->classifs) |
||
722 | return AVERROR(ENOMEM); |
||
723 | |||
724 | av_dlog(NULL, " begin %d end %d part.size %d classif.s %d classbook %d \n", |
||
725 | res_setup->begin, res_setup->end, res_setup->partition_size, |
||
726 | res_setup->classifications, res_setup->classbook); |
||
727 | |||
728 | for (j = 0; j < res_setup->classifications; ++j) { |
||
729 | high_bits = 0; |
||
730 | low_bits = get_bits(gb, 3); |
||
731 | if (get_bits1(gb)) |
||
732 | high_bits = get_bits(gb, 5); |
||
733 | cascade[j] = (high_bits << 3) + low_bits; |
||
734 | |||
735 | av_dlog(NULL, " %u class cascade depth: %d\n", j, ilog(cascade[j])); |
||
736 | } |
||
737 | |||
738 | res_setup->maxpass = 0; |
||
739 | for (j = 0; j < res_setup->classifications; ++j) { |
||
740 | for (k = 0; k < 8; ++k) { |
||
741 | if (cascade[j]&(1 << k)) { |
||
742 | GET_VALIDATED_INDEX(res_setup->books[j][k], 8, vc->codebook_count) |
||
743 | |||
744 | av_dlog(NULL, " %u class cascade depth %u book: %d\n", |
||
745 | j, k, res_setup->books[j][k]); |
||
746 | |||
747 | if (k>res_setup->maxpass) |
||
748 | res_setup->maxpass = k; |
||
749 | } else { |
||
750 | res_setup->books[j][k] = -1; |
||
751 | } |
||
752 | } |
||
753 | } |
||
754 | } |
||
755 | return 0; |
||
756 | } |
||
757 | |||
758 | // Process mappings part |
||
759 | |||
760 | static int vorbis_parse_setup_hdr_mappings(vorbis_context *vc) |
||
761 | { |
||
762 | GetBitContext *gb = &vc->gb; |
||
763 | unsigned i, j; |
||
764 | |||
765 | vc->mapping_count = get_bits(gb, 6)+1; |
||
766 | vc->mappings = av_mallocz(vc->mapping_count * sizeof(*vc->mappings)); |
||
767 | if (!vc->mappings) |
||
768 | return AVERROR(ENOMEM); |
||
769 | |||
770 | av_dlog(NULL, " There are %d mappings. \n", vc->mapping_count); |
||
771 | |||
772 | for (i = 0; i < vc->mapping_count; ++i) { |
||
773 | vorbis_mapping *mapping_setup = &vc->mappings[i]; |
||
774 | |||
775 | if (get_bits(gb, 16)) { |
||
776 | av_log(vc->avctx, AV_LOG_ERROR, "Other mappings than type 0 are not compliant with the Vorbis I specification. \n"); |
||
777 | return AVERROR_INVALIDDATA; |
||
778 | } |
||
779 | if (get_bits1(gb)) { |
||
780 | mapping_setup->submaps = get_bits(gb, 4) + 1; |
||
781 | } else { |
||
782 | mapping_setup->submaps = 1; |
||
783 | } |
||
784 | |||
785 | if (get_bits1(gb)) { |
||
786 | mapping_setup->coupling_steps = get_bits(gb, 8) + 1; |
||
787 | mapping_setup->magnitude = av_mallocz(mapping_setup->coupling_steps * |
||
788 | sizeof(*mapping_setup->magnitude)); |
||
789 | mapping_setup->angle = av_mallocz(mapping_setup->coupling_steps * |
||
790 | sizeof(*mapping_setup->angle)); |
||
791 | if (!mapping_setup->angle || !mapping_setup->magnitude) |
||
792 | return AVERROR(ENOMEM); |
||
793 | |||
794 | for (j = 0; j < mapping_setup->coupling_steps; ++j) { |
||
795 | GET_VALIDATED_INDEX(mapping_setup->magnitude[j], ilog(vc->audio_channels - 1), vc->audio_channels) |
||
796 | GET_VALIDATED_INDEX(mapping_setup->angle[j], ilog(vc->audio_channels - 1), vc->audio_channels) |
||
797 | } |
||
798 | } else { |
||
799 | mapping_setup->coupling_steps = 0; |
||
800 | } |
||
801 | |||
802 | av_dlog(NULL, " %u mapping coupling steps: %d\n", |
||
803 | i, mapping_setup->coupling_steps); |
||
804 | |||
805 | if (get_bits(gb, 2)) { |
||
806 | av_log(vc->avctx, AV_LOG_ERROR, "%u. mapping setup data invalid.\n", i); |
||
807 | return AVERROR_INVALIDDATA; // following spec. |
||
808 | } |
||
809 | |||
810 | if (mapping_setup->submaps>1) { |
||
811 | mapping_setup->mux = av_mallocz(vc->audio_channels * |
||
812 | sizeof(*mapping_setup->mux)); |
||
813 | if (!mapping_setup->mux) |
||
814 | return AVERROR(ENOMEM); |
||
815 | |||
816 | for (j = 0; j < vc->audio_channels; ++j) |
||
817 | mapping_setup->mux[j] = get_bits(gb, 4); |
||
818 | } |
||
819 | |||
820 | for (j = 0; j < mapping_setup->submaps; ++j) { |
||
821 | skip_bits(gb, 8); // FIXME check? |
||
822 | GET_VALIDATED_INDEX(mapping_setup->submap_floor[j], 8, vc->floor_count) |
||
823 | GET_VALIDATED_INDEX(mapping_setup->submap_residue[j], 8, vc->residue_count) |
||
824 | |||
825 | av_dlog(NULL, " %u mapping %u submap : floor %d, residue %d\n", i, j, |
||
826 | mapping_setup->submap_floor[j], |
||
827 | mapping_setup->submap_residue[j]); |
||
828 | } |
||
829 | } |
||
830 | return 0; |
||
831 | } |
||
832 | |||
833 | // Process modes part |
||
834 | |||
835 | static int create_map(vorbis_context *vc, unsigned floor_number) |
||
836 | { |
||
837 | vorbis_floor *floors = vc->floors; |
||
838 | vorbis_floor0 *vf; |
||
839 | int idx; |
||
840 | int blockflag, n; |
||
841 | int32_t *map; |
||
842 | |||
843 | for (blockflag = 0; blockflag < 2; ++blockflag) { |
||
844 | n = vc->blocksize[blockflag] / 2; |
||
845 | floors[floor_number].data.t0.map[blockflag] = |
||
846 | av_malloc((n + 1) * sizeof(int32_t)); // n + sentinel |
||
847 | if (!floors[floor_number].data.t0.map[blockflag]) |
||
848 | return AVERROR(ENOMEM); |
||
849 | |||
850 | map = floors[floor_number].data.t0.map[blockflag]; |
||
851 | vf = &floors[floor_number].data.t0; |
||
852 | |||
853 | for (idx = 0; idx < n; ++idx) { |
||
854 | map[idx] = floor(BARK((vf->rate * idx) / (2.0f * n)) * |
||
855 | (vf->bark_map_size / BARK(vf->rate / 2.0f))); |
||
856 | if (vf->bark_map_size-1 < map[idx]) |
||
857 | map[idx] = vf->bark_map_size - 1; |
||
858 | } |
||
859 | map[n] = -1; |
||
860 | vf->map_size[blockflag] = n; |
||
861 | } |
||
862 | |||
863 | for (idx = 0; idx <= n; ++idx) { |
||
864 | av_dlog(NULL, "floor0 map: map at pos %d is %d\n", idx, map[idx]); |
||
865 | } |
||
866 | |||
867 | return 0; |
||
868 | } |
||
869 | |||
870 | static int vorbis_parse_setup_hdr_modes(vorbis_context *vc) |
||
871 | { |
||
872 | GetBitContext *gb = &vc->gb; |
||
873 | unsigned i; |
||
874 | |||
875 | vc->mode_count = get_bits(gb, 6) + 1; |
||
876 | vc->modes = av_mallocz(vc->mode_count * sizeof(*vc->modes)); |
||
877 | if (!vc->modes) |
||
878 | return AVERROR(ENOMEM); |
||
879 | |||
880 | av_dlog(NULL, " There are %d modes.\n", vc->mode_count); |
||
881 | |||
882 | for (i = 0; i < vc->mode_count; ++i) { |
||
883 | vorbis_mode *mode_setup = &vc->modes[i]; |
||
884 | |||
885 | mode_setup->blockflag = get_bits1(gb); |
||
886 | mode_setup->windowtype = get_bits(gb, 16); //FIXME check |
||
887 | mode_setup->transformtype = get_bits(gb, 16); //FIXME check |
||
888 | GET_VALIDATED_INDEX(mode_setup->mapping, 8, vc->mapping_count); |
||
889 | |||
890 | av_dlog(NULL, " %u mode: blockflag %d, windowtype %d, transformtype %d, mapping %d\n", |
||
891 | i, mode_setup->blockflag, mode_setup->windowtype, |
||
892 | mode_setup->transformtype, mode_setup->mapping); |
||
893 | } |
||
894 | return 0; |
||
895 | } |
||
896 | |||
897 | // Process the whole setup header using the functions above |
||
898 | |||
899 | static int vorbis_parse_setup_hdr(vorbis_context *vc) |
||
900 | { |
||
901 | GetBitContext *gb = &vc->gb; |
||
902 | int ret; |
||
903 | |||
904 | if ((get_bits(gb, 8) != 'v') || (get_bits(gb, 8) != 'o') || |
||
905 | (get_bits(gb, 8) != 'r') || (get_bits(gb, 8) != 'b') || |
||
906 | (get_bits(gb, 8) != 'i') || (get_bits(gb, 8) != 's')) { |
||
907 | av_log(vc->avctx, AV_LOG_ERROR, " Vorbis setup header packet corrupt (no vorbis signature). \n"); |
||
908 | return AVERROR_INVALIDDATA; |
||
909 | } |
||
910 | |||
911 | if ((ret = vorbis_parse_setup_hdr_codebooks(vc))) { |
||
912 | av_log(vc->avctx, AV_LOG_ERROR, " Vorbis setup header packet corrupt (codebooks). \n"); |
||
913 | return ret; |
||
914 | } |
||
915 | if ((ret = vorbis_parse_setup_hdr_tdtransforms(vc))) { |
||
916 | av_log(vc->avctx, AV_LOG_ERROR, " Vorbis setup header packet corrupt (time domain transforms). \n"); |
||
917 | return ret; |
||
918 | } |
||
919 | if ((ret = vorbis_parse_setup_hdr_floors(vc))) { |
||
920 | av_log(vc->avctx, AV_LOG_ERROR, " Vorbis setup header packet corrupt (floors). \n"); |
||
921 | return ret; |
||
922 | } |
||
923 | if ((ret = vorbis_parse_setup_hdr_residues(vc))) { |
||
924 | av_log(vc->avctx, AV_LOG_ERROR, " Vorbis setup header packet corrupt (residues). \n"); |
||
925 | return ret; |
||
926 | } |
||
927 | if ((ret = vorbis_parse_setup_hdr_mappings(vc))) { |
||
928 | av_log(vc->avctx, AV_LOG_ERROR, " Vorbis setup header packet corrupt (mappings). \n"); |
||
929 | return ret; |
||
930 | } |
||
931 | if ((ret = vorbis_parse_setup_hdr_modes(vc))) { |
||
932 | av_log(vc->avctx, AV_LOG_ERROR, " Vorbis setup header packet corrupt (modes). \n"); |
||
933 | return ret; |
||
934 | } |
||
935 | if (!get_bits1(gb)) { |
||
936 | av_log(vc->avctx, AV_LOG_ERROR, " Vorbis setup header packet corrupt (framing flag). \n"); |
||
937 | return AVERROR_INVALIDDATA; // framing flag bit unset error |
||
938 | } |
||
939 | |||
940 | return 0; |
||
941 | } |
||
942 | |||
943 | // Process the identification header |
||
944 | |||
945 | static int vorbis_parse_id_hdr(vorbis_context *vc) |
||
946 | { |
||
947 | GetBitContext *gb = &vc->gb; |
||
948 | unsigned bl0, bl1; |
||
949 | |||
950 | if ((get_bits(gb, 8) != 'v') || (get_bits(gb, 8) != 'o') || |
||
951 | (get_bits(gb, 8) != 'r') || (get_bits(gb, 8) != 'b') || |
||
952 | (get_bits(gb, 8) != 'i') || (get_bits(gb, 8) != 's')) { |
||
953 | av_log(vc->avctx, AV_LOG_ERROR, " Vorbis id header packet corrupt (no vorbis signature). \n"); |
||
954 | return AVERROR_INVALIDDATA; |
||
955 | } |
||
956 | |||
957 | vc->version = get_bits_long(gb, 32); //FIXME check 0 |
||
958 | vc->audio_channels = get_bits(gb, 8); |
||
959 | if (vc->audio_channels <= 0) { |
||
960 | av_log(vc->avctx, AV_LOG_ERROR, "Invalid number of channels\n"); |
||
961 | return AVERROR_INVALIDDATA; |
||
962 | } |
||
963 | vc->audio_samplerate = get_bits_long(gb, 32); |
||
964 | if (vc->audio_samplerate <= 0) { |
||
965 | av_log(vc->avctx, AV_LOG_ERROR, "Invalid samplerate\n"); |
||
966 | return AVERROR_INVALIDDATA; |
||
967 | } |
||
968 | vc->bitrate_maximum = get_bits_long(gb, 32); |
||
969 | vc->bitrate_nominal = get_bits_long(gb, 32); |
||
970 | vc->bitrate_minimum = get_bits_long(gb, 32); |
||
971 | bl0 = get_bits(gb, 4); |
||
972 | bl1 = get_bits(gb, 4); |
||
973 | if (bl0 > 13 || bl0 < 6 || bl1 > 13 || bl1 < 6 || bl1 < bl0) { |
||
974 | av_log(vc->avctx, AV_LOG_ERROR, " Vorbis id header packet corrupt (illegal blocksize). \n"); |
||
975 | return AVERROR_INVALIDDATA; |
||
976 | } |
||
977 | vc->blocksize[0] = (1 << bl0); |
||
978 | vc->blocksize[1] = (1 << bl1); |
||
979 | vc->win[0] = ff_vorbis_vwin[bl0 - 6]; |
||
980 | vc->win[1] = ff_vorbis_vwin[bl1 - 6]; |
||
981 | |||
982 | if ((get_bits1(gb)) == 0) { |
||
983 | av_log(vc->avctx, AV_LOG_ERROR, " Vorbis id header packet corrupt (framing flag not set). \n"); |
||
984 | return AVERROR_INVALIDDATA; |
||
985 | } |
||
986 | |||
987 | vc->channel_residues = av_malloc((vc->blocksize[1] / 2) * vc->audio_channels * sizeof(*vc->channel_residues)); |
||
988 | vc->saved = av_mallocz((vc->blocksize[1] / 4) * vc->audio_channels * sizeof(*vc->saved)); |
||
989 | if (!vc->channel_residues || !vc->saved) |
||
990 | return AVERROR(ENOMEM); |
||
991 | |||
992 | vc->previous_window = 0; |
||
993 | |||
994 | ff_mdct_init(&vc->mdct[0], bl0, 1, -1.0); |
||
995 | ff_mdct_init(&vc->mdct[1], bl1, 1, -1.0); |
||
996 | |||
997 | av_dlog(NULL, " vorbis version %d \n audio_channels %d \n audio_samplerate %d \n bitrate_max %d \n bitrate_nom %d \n bitrate_min %d \n blk_0 %d blk_1 %d \n ", |
||
998 | vc->version, vc->audio_channels, vc->audio_samplerate, vc->bitrate_maximum, vc->bitrate_nominal, vc->bitrate_minimum, vc->blocksize[0], vc->blocksize[1]); |
||
999 | |||
1000 | /* |
||
1001 | BLK = vc->blocksize[0]; |
||
1002 | for (i = 0; i < BLK / 2; ++i) { |
||
1003 | vc->win[0][i] = sin(0.5*3.14159265358*(sin(((float)i + 0.5) / (float)BLK*3.14159265358))*(sin(((float)i + 0.5) / (float)BLK*3.14159265358))); |
||
1004 | } |
||
1005 | */ |
||
1006 | |||
1007 | return 0; |
||
1008 | } |
||
1009 | |||
1010 | // Process the extradata using the functions above (identification header, setup header) |
||
1011 | |||
1012 | static av_cold int vorbis_decode_init(AVCodecContext *avctx) |
||
1013 | { |
||
1014 | vorbis_context *vc = avctx->priv_data; |
||
1015 | uint8_t *headers = avctx->extradata; |
||
1016 | int headers_len = avctx->extradata_size; |
||
1017 | uint8_t *header_start[3]; |
||
1018 | int header_len[3]; |
||
1019 | GetBitContext *gb = &vc->gb; |
||
1020 | int hdr_type, ret; |
||
1021 | |||
1022 | vc->avctx = avctx; |
||
1023 | ff_vorbisdsp_init(&vc->dsp); |
||
1024 | avpriv_float_dsp_init(&vc->fdsp, avctx->flags & CODEC_FLAG_BITEXACT); |
||
1025 | ff_fmt_convert_init(&vc->fmt_conv, avctx); |
||
1026 | |||
1027 | avctx->sample_fmt = AV_SAMPLE_FMT_FLTP; |
||
1028 | |||
1029 | if (!headers_len) { |
||
1030 | av_log(avctx, AV_LOG_ERROR, "Extradata missing.\n"); |
||
1031 | return AVERROR_INVALIDDATA; |
||
1032 | } |
||
1033 | |||
1034 | if ((ret = avpriv_split_xiph_headers(headers, headers_len, 30, header_start, header_len)) < 0) { |
||
1035 | av_log(avctx, AV_LOG_ERROR, "Extradata corrupt.\n"); |
||
1036 | return ret; |
||
1037 | } |
||
1038 | |||
1039 | init_get_bits(gb, header_start[0], header_len[0]*8); |
||
1040 | hdr_type = get_bits(gb, 8); |
||
1041 | if (hdr_type != 1) { |
||
1042 | av_log(avctx, AV_LOG_ERROR, "First header is not the id header.\n"); |
||
1043 | return AVERROR_INVALIDDATA; |
||
1044 | } |
||
1045 | if ((ret = vorbis_parse_id_hdr(vc))) { |
||
1046 | av_log(avctx, AV_LOG_ERROR, "Id header corrupt.\n"); |
||
1047 | vorbis_free(vc); |
||
1048 | return ret; |
||
1049 | } |
||
1050 | |||
1051 | init_get_bits(gb, header_start[2], header_len[2]*8); |
||
1052 | hdr_type = get_bits(gb, 8); |
||
1053 | if (hdr_type != 5) { |
||
1054 | av_log(avctx, AV_LOG_ERROR, "Third header is not the setup header.\n"); |
||
1055 | vorbis_free(vc); |
||
1056 | return AVERROR_INVALIDDATA; |
||
1057 | } |
||
1058 | if ((ret = vorbis_parse_setup_hdr(vc))) { |
||
1059 | av_log(avctx, AV_LOG_ERROR, "Setup header corrupt.\n"); |
||
1060 | vorbis_free(vc); |
||
1061 | return ret; |
||
1062 | } |
||
1063 | |||
1064 | if (vc->audio_channels > 8) |
||
1065 | avctx->channel_layout = 0; |
||
1066 | else |
||
1067 | avctx->channel_layout = ff_vorbis_channel_layouts[vc->audio_channels - 1]; |
||
1068 | |||
1069 | avctx->channels = vc->audio_channels; |
||
1070 | avctx->sample_rate = vc->audio_samplerate; |
||
1071 | |||
1072 | return 0; |
||
1073 | } |
||
1074 | |||
1075 | // Decode audiopackets ------------------------------------------------- |
||
1076 | |||
1077 | // Read and decode floor |
||
1078 | |||
1079 | static int vorbis_floor0_decode(vorbis_context *vc, |
||
1080 | vorbis_floor_data *vfu, float *vec) |
||
1081 | { |
||
1082 | vorbis_floor0 *vf = &vfu->t0; |
||
1083 | float *lsp = vf->lsp; |
||
1084 | unsigned amplitude, book_idx; |
||
1085 | unsigned blockflag = vc->modes[vc->mode_number].blockflag; |
||
1086 | |||
1087 | if (!vf->amplitude_bits) |
||
1088 | return 1; |
||
1089 | |||
1090 | amplitude = get_bits(&vc->gb, vf->amplitude_bits); |
||
1091 | if (amplitude > 0) { |
||
1092 | float last = 0; |
||
1093 | unsigned idx, lsp_len = 0; |
||
1094 | vorbis_codebook codebook; |
||
1095 | |||
1096 | book_idx = get_bits(&vc->gb, ilog(vf->num_books)); |
||
1097 | if (book_idx >= vf->num_books) { |
||
1098 | av_log(vc->avctx, AV_LOG_ERROR, "floor0 dec: booknumber too high!\n"); |
||
1099 | book_idx = 0; |
||
1100 | } |
||
1101 | av_dlog(NULL, "floor0 dec: booknumber: %u\n", book_idx); |
||
1102 | codebook = vc->codebooks[vf->book_list[book_idx]]; |
||
1103 | /* Invalid codebook! */ |
||
1104 | if (!codebook.codevectors) |
||
1105 | return AVERROR_INVALIDDATA; |
||
1106 | |||
1107 | while (lsp_len |
||
1108 | int vec_off; |
||
1109 | |||
1110 | av_dlog(NULL, "floor0 dec: book dimension: %d\n", codebook.dimensions); |
||
1111 | av_dlog(NULL, "floor0 dec: maximum depth: %d\n", codebook.maxdepth); |
||
1112 | /* read temp vector */ |
||
1113 | vec_off = get_vlc2(&vc->gb, codebook.vlc.table, |
||
1114 | codebook.nb_bits, codebook.maxdepth) |
||
1115 | * codebook.dimensions; |
||
1116 | av_dlog(NULL, "floor0 dec: vector offset: %d\n", vec_off); |
||
1117 | /* copy each vector component and add last to it */ |
||
1118 | for (idx = 0; idx < codebook.dimensions; ++idx) |
||
1119 | lsp[lsp_len+idx] = codebook.codevectors[vec_off+idx] + last; |
||
1120 | last = lsp[lsp_len+idx-1]; /* set last to last vector component */ |
||
1121 | |||
1122 | lsp_len += codebook.dimensions; |
||
1123 | } |
||
1124 | /* DEBUG: output lsp coeffs */ |
||
1125 | { |
||
1126 | int idx; |
||
1127 | for (idx = 0; idx < lsp_len; ++idx) |
||
1128 | av_dlog(NULL, "floor0 dec: coeff at %d is %f\n", idx, lsp[idx]); |
||
1129 | } |
||
1130 | |||
1131 | /* synthesize floor output vector */ |
||
1132 | { |
||
1133 | int i; |
||
1134 | int order = vf->order; |
||
1135 | float wstep = M_PI / vf->bark_map_size; |
||
1136 | |||
1137 | for (i = 0; i < order; i++) |
||
1138 | lsp[i] = 2.0f * cos(lsp[i]); |
||
1139 | |||
1140 | av_dlog(NULL, "floor0 synth: map_size = %"PRIu32"; m = %d; wstep = %f\n", |
||
1141 | vf->map_size[blockflag], order, wstep); |
||
1142 | |||
1143 | i = 0; |
||
1144 | while (i < vf->map_size[blockflag]) { |
||
1145 | int j, iter_cond = vf->map[blockflag][i]; |
||
1146 | float p = 0.5f; |
||
1147 | float q = 0.5f; |
||
1148 | float two_cos_w = 2.0f * cos(wstep * iter_cond); // needed all times |
||
1149 | |||
1150 | /* similar part for the q and p products */ |
||
1151 | for (j = 0; j + 1 < order; j += 2) { |
||
1152 | q *= lsp[j] - two_cos_w; |
||
1153 | p *= lsp[j + 1] - two_cos_w; |
||
1154 | } |
||
1155 | if (j == order) { // even order |
||
1156 | p *= p * (2.0f - two_cos_w); |
||
1157 | q *= q * (2.0f + two_cos_w); |
||
1158 | } else { // odd order |
||
1159 | q *= two_cos_w-lsp[j]; // one more time for q |
||
1160 | |||
1161 | /* final step and square */ |
||
1162 | p *= p * (4.f - two_cos_w * two_cos_w); |
||
1163 | q *= q; |
||
1164 | } |
||
1165 | |||
1166 | /* calculate linear floor value */ |
||
1167 | q = exp((((amplitude*vf->amplitude_offset) / |
||
1168 | (((1 << vf->amplitude_bits) - 1) * sqrt(p + q))) |
||
1169 | - vf->amplitude_offset) * .11512925f); |
||
1170 | |||
1171 | /* fill vector */ |
||
1172 | do { |
||
1173 | vec[i] = q; ++i; |
||
1174 | } while (vf->map[blockflag][i] == iter_cond); |
||
1175 | } |
||
1176 | } |
||
1177 | } else { |
||
1178 | /* this channel is unused */ |
||
1179 | return 1; |
||
1180 | } |
||
1181 | |||
1182 | av_dlog(NULL, " Floor0 decoded\n"); |
||
1183 | |||
1184 | return 0; |
||
1185 | } |
||
1186 | |||
1187 | static int vorbis_floor1_decode(vorbis_context *vc, |
||
1188 | vorbis_floor_data *vfu, float *vec) |
||
1189 | { |
||
1190 | vorbis_floor1 *vf = &vfu->t1; |
||
1191 | GetBitContext *gb = &vc->gb; |
||
1192 | uint16_t range_v[4] = { 256, 128, 86, 64 }; |
||
1193 | unsigned range = range_v[vf->multiplier - 1]; |
||
1194 | uint16_t floor1_Y[258]; |
||
1195 | uint16_t floor1_Y_final[258]; |
||
1196 | int floor1_flag[258]; |
||
1197 | unsigned partition_class, cdim, cbits, csub, cval, offset, i, j; |
||
1198 | int book, adx, ady, dy, off, predicted, err; |
||
1199 | |||
1200 | |||
1201 | if (!get_bits1(gb)) // silence |
||
1202 | return 1; |
||
1203 | |||
1204 | // Read values (or differences) for the floor's points |
||
1205 | |||
1206 | floor1_Y[0] = get_bits(gb, ilog(range - 1)); |
||
1207 | floor1_Y[1] = get_bits(gb, ilog(range - 1)); |
||
1208 | |||
1209 | av_dlog(NULL, "floor 0 Y %d floor 1 Y %d \n", floor1_Y[0], floor1_Y[1]); |
||
1210 | |||
1211 | offset = 2; |
||
1212 | for (i = 0; i < vf->partitions; ++i) { |
||
1213 | partition_class = vf->partition_class[i]; |
||
1214 | cdim = vf->class_dimensions[partition_class]; |
||
1215 | cbits = vf->class_subclasses[partition_class]; |
||
1216 | csub = (1 << cbits) - 1; |
||
1217 | cval = 0; |
||
1218 | |||
1219 | av_dlog(NULL, "Cbits %u\n", cbits); |
||
1220 | |||
1221 | if (cbits) // this reads all subclasses for this partition's class |
||
1222 | cval = get_vlc2(gb, vc->codebooks[vf->class_masterbook[partition_class]].vlc.table, |
||
1223 | vc->codebooks[vf->class_masterbook[partition_class]].nb_bits, 3); |
||
1224 | |||
1225 | for (j = 0; j < cdim; ++j) { |
||
1226 | book = vf->subclass_books[partition_class][cval & csub]; |
||
1227 | |||
1228 | av_dlog(NULL, "book %d Cbits %u cval %u bits:%d\n", |
||
1229 | book, cbits, cval, get_bits_count(gb)); |
||
1230 | |||
1231 | cval = cval >> cbits; |
||
1232 | if (book > -1) { |
||
1233 | int v = get_vlc2(gb, vc->codebooks[book].vlc.table, |
||
1234 | vc->codebooks[book].nb_bits, 3); |
||
1235 | if (v < 0) |
||
1236 | return AVERROR_INVALIDDATA; |
||
1237 | floor1_Y[offset+j] = v; |
||
1238 | } else { |
||
1239 | floor1_Y[offset+j] = 0; |
||
1240 | } |
||
1241 | |||
1242 | av_dlog(NULL, " floor(%d) = %d \n", |
||
1243 | vf->list[offset+j].x, floor1_Y[offset+j]); |
||
1244 | } |
||
1245 | offset+=cdim; |
||
1246 | } |
||
1247 | |||
1248 | // Amplitude calculation from the differences |
||
1249 | |||
1250 | floor1_flag[0] = 1; |
||
1251 | floor1_flag[1] = 1; |
||
1252 | floor1_Y_final[0] = floor1_Y[0]; |
||
1253 | floor1_Y_final[1] = floor1_Y[1]; |
||
1254 | |||
1255 | for (i = 2; i < vf->x_list_dim; ++i) { |
||
1256 | unsigned val, highroom, lowroom, room, high_neigh_offs, low_neigh_offs; |
||
1257 | |||
1258 | low_neigh_offs = vf->list[i].low; |
||
1259 | high_neigh_offs = vf->list[i].high; |
||
1260 | dy = floor1_Y_final[high_neigh_offs] - floor1_Y_final[low_neigh_offs]; // render_point begin |
||
1261 | adx = vf->list[high_neigh_offs].x - vf->list[low_neigh_offs].x; |
||
1262 | ady = FFABS(dy); |
||
1263 | err = ady * (vf->list[i].x - vf->list[low_neigh_offs].x); |
||
1264 | off = err / adx; |
||
1265 | if (dy < 0) { |
||
1266 | predicted = floor1_Y_final[low_neigh_offs] - off; |
||
1267 | } else { |
||
1268 | predicted = floor1_Y_final[low_neigh_offs] + off; |
||
1269 | } // render_point end |
||
1270 | |||
1271 | val = floor1_Y[i]; |
||
1272 | highroom = range-predicted; |
||
1273 | lowroom = predicted; |
||
1274 | if (highroom < lowroom) { |
||
1275 | room = highroom * 2; |
||
1276 | } else { |
||
1277 | room = lowroom * 2; // SPEC misspelling |
||
1278 | } |
||
1279 | if (val) { |
||
1280 | floor1_flag[low_neigh_offs] = 1; |
||
1281 | floor1_flag[high_neigh_offs] = 1; |
||
1282 | floor1_flag[i] = 1; |
||
1283 | if (val >= room) { |
||
1284 | if (highroom > lowroom) { |
||
1285 | floor1_Y_final[i] = av_clip_uint16(val - lowroom + predicted); |
||
1286 | } else { |
||
1287 | floor1_Y_final[i] = av_clip_uint16(predicted - val + highroom - 1); |
||
1288 | } |
||
1289 | } else { |
||
1290 | if (val & 1) { |
||
1291 | floor1_Y_final[i] = av_clip_uint16(predicted - (val + 1) / 2); |
||
1292 | } else { |
||
1293 | floor1_Y_final[i] = av_clip_uint16(predicted + val / 2); |
||
1294 | } |
||
1295 | } |
||
1296 | } else { |
||
1297 | floor1_flag[i] = 0; |
||
1298 | floor1_Y_final[i] = av_clip_uint16(predicted); |
||
1299 | } |
||
1300 | |||
1301 | av_dlog(NULL, " Decoded floor(%d) = %u / val %u\n", |
||
1302 | vf->list[i].x, floor1_Y_final[i], val); |
||
1303 | } |
||
1304 | |||
1305 | // Curve synth - connect the calculated dots and convert from dB scale FIXME optimize ? |
||
1306 | |||
1307 | ff_vorbis_floor1_render_list(vf->list, vf->x_list_dim, floor1_Y_final, floor1_flag, vf->multiplier, vec, vf->list[1].x); |
||
1308 | |||
1309 | av_dlog(NULL, " Floor decoded\n"); |
||
1310 | |||
1311 | return 0; |
||
1312 | } |
||
1313 | |||
1314 | static av_always_inline int setup_classifs(vorbis_context *vc, |
||
1315 | vorbis_residue *vr, |
||
1316 | uint8_t *do_not_decode, |
||
1317 | unsigned ch_used, |
||
1318 | int partition_count) |
||
1319 | { |
||
1320 | int p, j, i; |
||
1321 | unsigned c_p_c = vc->codebooks[vr->classbook].dimensions; |
||
1322 | unsigned inverse_class = ff_inverse[vr->classifications]; |
||
1323 | unsigned temp, temp2; |
||
1324 | for (p = 0, j = 0; j < ch_used; ++j) { |
||
1325 | if (!do_not_decode[j]) { |
||
1326 | temp = get_vlc2(&vc->gb, vc->codebooks[vr->classbook].vlc.table, |
||
1327 | vc->codebooks[vr->classbook].nb_bits, 3); |
||
1328 | |||
1329 | av_dlog(NULL, "Classword: %u\n", temp); |
||
1330 | |||
1331 | if ((int)temp < 0) |
||
1332 | return temp; |
||
1333 | |||
1334 | av_assert0(vr->classifications > 1); //needed for inverse[] |
||
1335 | |||
1336 | if (temp <= 65536) { |
||
1337 | for (i = partition_count + c_p_c - 1; i >= partition_count; i--) { |
||
1338 | temp2 = (((uint64_t)temp) * inverse_class) >> 32; |
||
1339 | |||
1340 | if (i < vr->ptns_to_read) |
||
1341 | vr->classifs[p + i] = temp - temp2 * vr->classifications; |
||
1342 | temp = temp2; |
||
1343 | } |
||
1344 | } else { |
||
1345 | for (i = partition_count + c_p_c - 1; i >= partition_count; i--) { |
||
1346 | temp2 = temp / vr->classifications; |
||
1347 | |||
1348 | if (i < vr->ptns_to_read) |
||
1349 | vr->classifs[p + i] = temp - temp2 * vr->classifications; |
||
1350 | temp = temp2; |
||
1351 | } |
||
1352 | } |
||
1353 | } |
||
1354 | p += vr->ptns_to_read; |
||
1355 | } |
||
1356 | return 0; |
||
1357 | } |
||
1358 | // Read and decode residue |
||
1359 | |||
1360 | static av_always_inline int vorbis_residue_decode_internal(vorbis_context *vc, |
||
1361 | vorbis_residue *vr, |
||
1362 | unsigned ch, |
||
1363 | uint8_t *do_not_decode, |
||
1364 | float *vec, |
||
1365 | unsigned vlen, |
||
1366 | unsigned ch_left, |
||
1367 | int vr_type) |
||
1368 | { |
||
1369 | GetBitContext *gb = &vc->gb; |
||
1370 | unsigned c_p_c = vc->codebooks[vr->classbook].dimensions; |
||
1371 | uint8_t *classifs = vr->classifs; |
||
1372 | unsigned pass, ch_used, i, j, k, l; |
||
1373 | unsigned max_output = (ch - 1) * vlen; |
||
1374 | int ptns_to_read = vr->ptns_to_read; |
||
1375 | |||
1376 | if (vr_type == 2) { |
||
1377 | for (j = 1; j < ch; ++j) |
||
1378 | do_not_decode[0] &= do_not_decode[j]; // FIXME - clobbering input |
||
1379 | if (do_not_decode[0]) |
||
1380 | return 0; |
||
1381 | ch_used = 1; |
||
1382 | max_output += vr->end / ch; |
||
1383 | } else { |
||
1384 | ch_used = ch; |
||
1385 | max_output += vr->end; |
||
1386 | } |
||
1387 | |||
1388 | if (max_output > ch_left * vlen) { |
||
1389 | av_log(vc->avctx, AV_LOG_ERROR, "Insufficient output buffer\n"); |
||
1390 | return AVERROR_INVALIDDATA; |
||
1391 | } |
||
1392 | |||
1393 | av_dlog(NULL, " residue type 0/1/2 decode begin, ch: %d cpc %d \n", ch, c_p_c); |
||
1394 | |||
1395 | for (pass = 0; pass <= vr->maxpass; ++pass) { // FIXME OPTIMIZE? |
||
1396 | int voffset, partition_count, j_times_ptns_to_read; |
||
1397 | |||
1398 | voffset = vr->begin; |
||
1399 | for (partition_count = 0; partition_count < ptns_to_read;) { // SPEC error |
||
1400 | if (!pass) { |
||
1401 | int ret; |
||
1402 | if ((ret = setup_classifs(vc, vr, do_not_decode, ch_used, partition_count)) < 0) |
||
1403 | return ret; |
||
1404 | } |
||
1405 | for (i = 0; (i < c_p_c) && (partition_count < ptns_to_read); ++i) { |
||
1406 | for (j_times_ptns_to_read = 0, j = 0; j < ch_used; ++j) { |
||
1407 | unsigned voffs; |
||
1408 | |||
1409 | if (!do_not_decode[j]) { |
||
1410 | unsigned vqclass = classifs[j_times_ptns_to_read + partition_count]; |
||
1411 | int vqbook = vr->books[vqclass][pass]; |
||
1412 | |||
1413 | if (vqbook >= 0 && vc->codebooks[vqbook].codevectors) { |
||
1414 | unsigned coffs; |
||
1415 | unsigned dim = vc->codebooks[vqbook].dimensions; |
||
1416 | unsigned step = FASTDIV(vr->partition_size << 1, dim << 1); |
||
1417 | vorbis_codebook codebook = vc->codebooks[vqbook]; |
||
1418 | |||
1419 | if (vr_type == 0) { |
||
1420 | |||
1421 | voffs = voffset+j*vlen; |
||
1422 | for (k = 0; k < step; ++k) { |
||
1423 | coffs = get_vlc2(gb, codebook.vlc.table, codebook.nb_bits, 3) * dim; |
||
1424 | for (l = 0; l < dim; ++l) |
||
1425 | vec[voffs + k + l * step] += codebook.codevectors[coffs + l]; |
||
1426 | } |
||
1427 | } else if (vr_type == 1) { |
||
1428 | voffs = voffset + j * vlen; |
||
1429 | for (k = 0; k < step; ++k) { |
||
1430 | coffs = get_vlc2(gb, codebook.vlc.table, codebook.nb_bits, 3) * dim; |
||
1431 | for (l = 0; l < dim; ++l, ++voffs) { |
||
1432 | vec[voffs]+=codebook.codevectors[coffs+l]; |
||
1433 | |||
1434 | av_dlog(NULL, " pass %d offs: %d curr: %f change: %f cv offs.: %d \n", |
||
1435 | pass, voffs, vec[voffs], codebook.codevectors[coffs+l], coffs); |
||
1436 | } |
||
1437 | } |
||
1438 | } else if (vr_type == 2 && ch == 2 && (voffset & 1) == 0 && (dim & 1) == 0) { // most frequent case optimized |
||
1439 | voffs = voffset >> 1; |
||
1440 | |||
1441 | if (dim == 2) { |
||
1442 | for (k = 0; k < step; ++k) { |
||
1443 | coffs = get_vlc2(gb, codebook.vlc.table, codebook.nb_bits, 3) * 2; |
||
1444 | vec[voffs + k ] += codebook.codevectors[coffs ]; |
||
1445 | vec[voffs + k + vlen] += codebook.codevectors[coffs + 1]; |
||
1446 | } |
||
1447 | } else if (dim == 4) { |
||
1448 | for (k = 0; k < step; ++k, voffs += 2) { |
||
1449 | coffs = get_vlc2(gb, codebook.vlc.table, codebook.nb_bits, 3) * 4; |
||
1450 | vec[voffs ] += codebook.codevectors[coffs ]; |
||
1451 | vec[voffs + 1 ] += codebook.codevectors[coffs + 2]; |
||
1452 | vec[voffs + vlen ] += codebook.codevectors[coffs + 1]; |
||
1453 | vec[voffs + vlen + 1] += codebook.codevectors[coffs + 3]; |
||
1454 | } |
||
1455 | } else |
||
1456 | for (k = 0; k < step; ++k) { |
||
1457 | coffs = get_vlc2(gb, codebook.vlc.table, codebook.nb_bits, 3) * dim; |
||
1458 | for (l = 0; l < dim; l += 2, voffs++) { |
||
1459 | vec[voffs ] += codebook.codevectors[coffs + l ]; |
||
1460 | vec[voffs + vlen] += codebook.codevectors[coffs + l + 1]; |
||
1461 | |||
1462 | av_dlog(NULL, " pass %d offs: %d curr: %f change: %f cv offs.: %d+%d \n", |
||
1463 | pass, voffset / ch + (voffs % ch) * vlen, |
||
1464 | vec[voffset / ch + (voffs % ch) * vlen], |
||
1465 | codebook.codevectors[coffs + l], coffs, l); |
||
1466 | } |
||
1467 | } |
||
1468 | |||
1469 | } else if (vr_type == 2) { |
||
1470 | unsigned voffs_div = FASTDIV(voffset << 1, ch <<1); |
||
1471 | unsigned voffs_mod = voffset - voffs_div * ch; |
||
1472 | |||
1473 | for (k = 0; k < step; ++k) { |
||
1474 | coffs = get_vlc2(gb, codebook.vlc.table, codebook.nb_bits, 3) * dim; |
||
1475 | for (l = 0; l < dim; ++l) { |
||
1476 | vec[voffs_div + voffs_mod * vlen] += |
||
1477 | codebook.codevectors[coffs + l]; |
||
1478 | |||
1479 | av_dlog(NULL, " pass %d offs: %d curr: %f change: %f cv offs.: %d+%d \n", |
||
1480 | pass, voffs_div + voffs_mod * vlen, |
||
1481 | vec[voffs_div + voffs_mod * vlen], |
||
1482 | codebook.codevectors[coffs + l], coffs, l); |
||
1483 | |||
1484 | if (++voffs_mod == ch) { |
||
1485 | voffs_div++; |
||
1486 | voffs_mod = 0; |
||
1487 | } |
||
1488 | } |
||
1489 | } |
||
1490 | } |
||
1491 | } |
||
1492 | } |
||
1493 | j_times_ptns_to_read += ptns_to_read; |
||
1494 | } |
||
1495 | ++partition_count; |
||
1496 | voffset += vr->partition_size; |
||
1497 | } |
||
1498 | } |
||
1499 | } |
||
1500 | return 0; |
||
1501 | } |
||
1502 | |||
1503 | static inline int vorbis_residue_decode(vorbis_context *vc, vorbis_residue *vr, |
||
1504 | unsigned ch, |
||
1505 | uint8_t *do_not_decode, |
||
1506 | float *vec, unsigned vlen, |
||
1507 | unsigned ch_left) |
||
1508 | { |
||
1509 | if (vr->type == 2) |
||
1510 | return vorbis_residue_decode_internal(vc, vr, ch, do_not_decode, vec, vlen, ch_left, 2); |
||
1511 | else if (vr->type == 1) |
||
1512 | return vorbis_residue_decode_internal(vc, vr, ch, do_not_decode, vec, vlen, ch_left, 1); |
||
1513 | else if (vr->type == 0) |
||
1514 | return vorbis_residue_decode_internal(vc, vr, ch, do_not_decode, vec, vlen, ch_left, 0); |
||
1515 | else { |
||
1516 | av_log(vc->avctx, AV_LOG_ERROR, " Invalid residue type while residue decode?! \n"); |
||
1517 | return AVERROR_INVALIDDATA; |
||
1518 | } |
||
1519 | } |
||
1520 | |||
1521 | void ff_vorbis_inverse_coupling(float *mag, float *ang, intptr_t blocksize) |
||
1522 | { |
||
1523 | int i; |
||
1524 | for (i = 0; i < blocksize; i++) { |
||
1525 | if (mag[i] > 0.0) { |
||
1526 | if (ang[i] > 0.0) { |
||
1527 | ang[i] = mag[i] - ang[i]; |
||
1528 | } else { |
||
1529 | float temp = ang[i]; |
||
1530 | ang[i] = mag[i]; |
||
1531 | mag[i] += temp; |
||
1532 | } |
||
1533 | } else { |
||
1534 | if (ang[i] > 0.0) { |
||
1535 | ang[i] += mag[i]; |
||
1536 | } else { |
||
1537 | float temp = ang[i]; |
||
1538 | ang[i] = mag[i]; |
||
1539 | mag[i] -= temp; |
||
1540 | } |
||
1541 | } |
||
1542 | } |
||
1543 | } |
||
1544 | |||
1545 | // Decode the audio packet using the functions above |
||
1546 | |||
1547 | static int vorbis_parse_audio_packet(vorbis_context *vc, float **floor_ptr) |
||
1548 | { |
||
1549 | GetBitContext *gb = &vc->gb; |
||
1550 | FFTContext *mdct; |
||
1551 | unsigned previous_window = vc->previous_window; |
||
1552 | unsigned mode_number, blockflag, blocksize; |
||
1553 | int i, j; |
||
1554 | uint8_t no_residue[255]; |
||
1555 | uint8_t do_not_decode[255]; |
||
1556 | vorbis_mapping *mapping; |
||
1557 | float *ch_res_ptr = vc->channel_residues; |
||
1558 | uint8_t res_chan[255]; |
||
1559 | unsigned res_num = 0; |
||
1560 | int retlen = 0; |
||
1561 | unsigned ch_left = vc->audio_channels; |
||
1562 | unsigned vlen; |
||
1563 | |||
1564 | if (get_bits1(gb)) { |
||
1565 | av_log(vc->avctx, AV_LOG_ERROR, "Not a Vorbis I audio packet.\n"); |
||
1566 | return AVERROR_INVALIDDATA; // packet type not audio |
||
1567 | } |
||
1568 | |||
1569 | if (vc->mode_count == 1) { |
||
1570 | mode_number = 0; |
||
1571 | } else { |
||
1572 | GET_VALIDATED_INDEX(mode_number, ilog(vc->mode_count-1), vc->mode_count) |
||
1573 | } |
||
1574 | vc->mode_number = mode_number; |
||
1575 | mapping = &vc->mappings[vc->modes[mode_number].mapping]; |
||
1576 | |||
1577 | av_dlog(NULL, " Mode number: %u , mapping: %d , blocktype %d\n", mode_number, |
||
1578 | vc->modes[mode_number].mapping, vc->modes[mode_number].blockflag); |
||
1579 | |||
1580 | blockflag = vc->modes[mode_number].blockflag; |
||
1581 | blocksize = vc->blocksize[blockflag]; |
||
1582 | vlen = blocksize / 2; |
||
1583 | if (blockflag) { |
||
1584 | previous_window = get_bits(gb, 1); |
||
1585 | skip_bits1(gb); // next_window |
||
1586 | } |
||
1587 | |||
1588 | memset(ch_res_ptr, 0, sizeof(float) * vc->audio_channels * vlen); //FIXME can this be removed ? |
||
1589 | for (i = 0; i < vc->audio_channels; ++i) |
||
1590 | memset(floor_ptr[i], 0, vlen * sizeof(floor_ptr[0][0])); //FIXME can this be removed ? |
||
1591 | |||
1592 | // Decode floor |
||
1593 | |||
1594 | for (i = 0; i < vc->audio_channels; ++i) { |
||
1595 | vorbis_floor *floor; |
||
1596 | int ret; |
||
1597 | if (mapping->submaps > 1) { |
||
1598 | floor = &vc->floors[mapping->submap_floor[mapping->mux[i]]]; |
||
1599 | } else { |
||
1600 | floor = &vc->floors[mapping->submap_floor[0]]; |
||
1601 | } |
||
1602 | |||
1603 | ret = floor->decode(vc, &floor->data, floor_ptr[i]); |
||
1604 | |||
1605 | if (ret < 0) { |
||
1606 | av_log(vc->avctx, AV_LOG_ERROR, "Invalid codebook in vorbis_floor_decode.\n"); |
||
1607 | return AVERROR_INVALIDDATA; |
||
1608 | } |
||
1609 | no_residue[i] = ret; |
||
1610 | } |
||
1611 | |||
1612 | // Nonzero vector propagate |
||
1613 | |||
1614 | for (i = mapping->coupling_steps - 1; i >= 0; --i) { |
||
1615 | if (!(no_residue[mapping->magnitude[i]] & no_residue[mapping->angle[i]])) { |
||
1616 | no_residue[mapping->magnitude[i]] = 0; |
||
1617 | no_residue[mapping->angle[i]] = 0; |
||
1618 | } |
||
1619 | } |
||
1620 | |||
1621 | // Decode residue |
||
1622 | |||
1623 | for (i = 0; i < mapping->submaps; ++i) { |
||
1624 | vorbis_residue *residue; |
||
1625 | unsigned ch = 0; |
||
1626 | int ret; |
||
1627 | |||
1628 | for (j = 0; j < vc->audio_channels; ++j) { |
||
1629 | if ((mapping->submaps == 1) || (i == mapping->mux[j])) { |
||
1630 | res_chan[j] = res_num; |
||
1631 | if (no_residue[j]) { |
||
1632 | do_not_decode[ch] = 1; |
||
1633 | } else { |
||
1634 | do_not_decode[ch] = 0; |
||
1635 | } |
||
1636 | ++ch; |
||
1637 | ++res_num; |
||
1638 | } |
||
1639 | } |
||
1640 | residue = &vc->residues[mapping->submap_residue[i]]; |
||
1641 | if (ch_left < ch) { |
||
1642 | av_log(vc->avctx, AV_LOG_ERROR, "Too many channels in vorbis_floor_decode.\n"); |
||
1643 | return AVERROR_INVALIDDATA; |
||
1644 | } |
||
1645 | if (ch) { |
||
1646 | ret = vorbis_residue_decode(vc, residue, ch, do_not_decode, ch_res_ptr, vlen, ch_left); |
||
1647 | if (ret < 0) |
||
1648 | return ret; |
||
1649 | } |
||
1650 | |||
1651 | ch_res_ptr += ch * vlen; |
||
1652 | ch_left -= ch; |
||
1653 | } |
||
1654 | |||
1655 | if (ch_left > 0) |
||
1656 | return AVERROR_INVALIDDATA; |
||
1657 | |||
1658 | // Inverse coupling |
||
1659 | |||
1660 | for (i = mapping->coupling_steps - 1; i >= 0; --i) { //warning: i has to be signed |
||
1661 | float *mag, *ang; |
||
1662 | |||
1663 | mag = vc->channel_residues+res_chan[mapping->magnitude[i]] * blocksize / 2; |
||
1664 | ang = vc->channel_residues+res_chan[mapping->angle[i]] * blocksize / 2; |
||
1665 | vc->dsp.vorbis_inverse_coupling(mag, ang, blocksize / 2); |
||
1666 | } |
||
1667 | |||
1668 | // Dotproduct, MDCT |
||
1669 | |||
1670 | mdct = &vc->mdct[blockflag]; |
||
1671 | |||
1672 | for (j = vc->audio_channels-1;j >= 0; j--) { |
||
1673 | ch_res_ptr = vc->channel_residues + res_chan[j] * blocksize / 2; |
||
1674 | vc->fdsp.vector_fmul(floor_ptr[j], floor_ptr[j], ch_res_ptr, blocksize / 2); |
||
1675 | mdct->imdct_half(mdct, ch_res_ptr, floor_ptr[j]); |
||
1676 | } |
||
1677 | |||
1678 | // Overlap/add, save data for next overlapping |
||
1679 | |||
1680 | retlen = (blocksize + vc->blocksize[previous_window]) / 4; |
||
1681 | for (j = 0; j < vc->audio_channels; j++) { |
||
1682 | unsigned bs0 = vc->blocksize[0]; |
||
1683 | unsigned bs1 = vc->blocksize[1]; |
||
1684 | float *residue = vc->channel_residues + res_chan[j] * blocksize / 2; |
||
1685 | float *saved = vc->saved + j * bs1 / 4; |
||
1686 | float *ret = floor_ptr[j]; |
||
1687 | float *buf = residue; |
||
1688 | const float *win = vc->win[blockflag & previous_window]; |
||
1689 | |||
1690 | if (blockflag == previous_window) { |
||
1691 | vc->fdsp.vector_fmul_window(ret, saved, buf, win, blocksize / 4); |
||
1692 | } else if (blockflag > previous_window) { |
||
1693 | vc->fdsp.vector_fmul_window(ret, saved, buf, win, bs0 / 4); |
||
1694 | memcpy(ret+bs0/2, buf+bs0/4, ((bs1-bs0)/4) * sizeof(float)); |
||
1695 | } else { |
||
1696 | memcpy(ret, saved, ((bs1 - bs0) / 4) * sizeof(float)); |
||
1697 | vc->fdsp.vector_fmul_window(ret + (bs1 - bs0) / 4, saved + (bs1 - bs0) / 4, buf, win, bs0 / 4); |
||
1698 | } |
||
1699 | memcpy(saved, buf + blocksize / 4, blocksize / 4 * sizeof(float)); |
||
1700 | } |
||
1701 | |||
1702 | vc->previous_window = blockflag; |
||
1703 | return retlen; |
||
1704 | } |
||
1705 | |||
1706 | // Return the decoded audio packet through the standard api |
||
1707 | |||
1708 | static int vorbis_decode_frame(AVCodecContext *avctx, void *data, |
||
1709 | int *got_frame_ptr, AVPacket *avpkt) |
||
1710 | { |
||
1711 | const uint8_t *buf = avpkt->data; |
||
1712 | int buf_size = avpkt->size; |
||
1713 | vorbis_context *vc = avctx->priv_data; |
||
1714 | AVFrame *frame = data; |
||
1715 | GetBitContext *gb = &vc->gb; |
||
1716 | float *channel_ptrs[255]; |
||
1717 | int i, len, ret; |
||
1718 | |||
1719 | av_dlog(NULL, "packet length %d \n", buf_size); |
||
1720 | |||
1721 | if (*buf == 1 && buf_size > 7) { |
||
1722 | init_get_bits(gb, buf+1, buf_size*8 - 8); |
||
1723 | vorbis_free(vc); |
||
1724 | if ((ret = vorbis_parse_id_hdr(vc))) { |
||
1725 | av_log(avctx, AV_LOG_ERROR, "Id header corrupt.\n"); |
||
1726 | vorbis_free(vc); |
||
1727 | return ret; |
||
1728 | } |
||
1729 | |||
1730 | if (vc->audio_channels > 8) |
||
1731 | avctx->channel_layout = 0; |
||
1732 | else |
||
1733 | avctx->channel_layout = ff_vorbis_channel_layouts[vc->audio_channels - 1]; |
||
1734 | |||
1735 | avctx->channels = vc->audio_channels; |
||
1736 | avctx->sample_rate = vc->audio_samplerate; |
||
1737 | return buf_size; |
||
1738 | } |
||
1739 | |||
1740 | if (*buf == 3 && buf_size > 7) { |
||
1741 | av_log(avctx, AV_LOG_DEBUG, "Ignoring comment header\n"); |
||
1742 | return buf_size; |
||
1743 | } |
||
1744 | |||
1745 | if (*buf == 5 && buf_size > 7 && vc->channel_residues && !vc->modes) { |
||
1746 | init_get_bits(gb, buf+1, buf_size*8 - 8); |
||
1747 | if ((ret = vorbis_parse_setup_hdr(vc))) { |
||
1748 | av_log(avctx, AV_LOG_ERROR, "Setup header corrupt.\n"); |
||
1749 | vorbis_free(vc); |
||
1750 | return ret; |
||
1751 | } |
||
1752 | return buf_size; |
||
1753 | } |
||
1754 | |||
1755 | if (!vc->channel_residues || !vc->modes) { |
||
1756 | av_log(avctx, AV_LOG_ERROR, "Data packet before valid headers\n"); |
||
1757 | return AVERROR_INVALIDDATA; |
||
1758 | } |
||
1759 | |||
1760 | /* get output buffer */ |
||
1761 | frame->nb_samples = vc->blocksize[1] / 2; |
||
1762 | if ((ret = ff_get_buffer(avctx, frame, 0)) < 0) |
||
1763 | return ret; |
||
1764 | |||
1765 | if (vc->audio_channels > 8) { |
||
1766 | for (i = 0; i < vc->audio_channels; i++) |
||
1767 | channel_ptrs[i] = (float *)frame->extended_data[i]; |
||
1768 | } else { |
||
1769 | for (i = 0; i < vc->audio_channels; i++) { |
||
1770 | int ch = ff_vorbis_channel_layout_offsets[vc->audio_channels - 1][i]; |
||
1771 | channel_ptrs[ch] = (float *)frame->extended_data[i]; |
||
1772 | } |
||
1773 | } |
||
1774 | |||
1775 | init_get_bits(gb, buf, buf_size*8); |
||
1776 | |||
1777 | if ((len = vorbis_parse_audio_packet(vc, channel_ptrs)) <= 0) |
||
1778 | return len; |
||
1779 | |||
1780 | if (!vc->first_frame) { |
||
1781 | vc->first_frame = 1; |
||
1782 | *got_frame_ptr = 0; |
||
1783 | av_frame_unref(frame); |
||
1784 | return buf_size; |
||
1785 | } |
||
1786 | |||
1787 | av_dlog(NULL, "parsed %d bytes %d bits, returned %d samples (*ch*bits) \n", |
||
1788 | get_bits_count(gb) / 8, get_bits_count(gb) % 8, len); |
||
1789 | |||
1790 | frame->nb_samples = len; |
||
1791 | *got_frame_ptr = 1; |
||
1792 | |||
1793 | return buf_size; |
||
1794 | } |
||
1795 | |||
1796 | // Close decoder |
||
1797 | |||
1798 | static av_cold int vorbis_decode_close(AVCodecContext *avctx) |
||
1799 | { |
||
1800 | vorbis_context *vc = avctx->priv_data; |
||
1801 | |||
1802 | vorbis_free(vc); |
||
1803 | |||
1804 | return 0; |
||
1805 | } |
||
1806 | |||
1807 | static av_cold void vorbis_decode_flush(AVCodecContext *avctx) |
||
1808 | { |
||
1809 | vorbis_context *vc = avctx->priv_data; |
||
1810 | |||
1811 | if (vc->saved) { |
||
1812 | memset(vc->saved, 0, (vc->blocksize[1] / 4) * vc->audio_channels * |
||
1813 | sizeof(*vc->saved)); |
||
1814 | } |
||
1815 | vc->previous_window = 0; |
||
1816 | } |
||
1817 | |||
1818 | AVCodec ff_vorbis_decoder = { |
||
1819 | .name = "vorbis", |
||
1820 | .long_name = NULL_IF_CONFIG_SMALL("Vorbis"), |
||
1821 | .type = AVMEDIA_TYPE_AUDIO, |
||
1822 | .id = AV_CODEC_ID_VORBIS, |
||
1823 | .priv_data_size = sizeof(vorbis_context), |
||
1824 | .init = vorbis_decode_init, |
||
1825 | .close = vorbis_decode_close, |
||
1826 | .decode = vorbis_decode_frame, |
||
1827 | .flush = vorbis_decode_flush, |
||
1828 | .capabilities = CODEC_CAP_DR1, |
||
1829 | .channel_layouts = ff_vorbis_channel_layouts, |
||
1830 | .sample_fmts = (const enum AVSampleFormat[]) { AV_SAMPLE_FMT_FLTP, |
||
1831 | AV_SAMPLE_FMT_NONE }, |
||
1832 | };=>>>>>>>>>>>>>>>1); |