| 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 <inttypes.h> |
| 30 | #include <math.h> |
| 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 | int8_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 | av_freep(&vc->fdsp); |
| 197 | |
| 198 | if (vc->residues) |
| 199 | for (i = 0; i < vc->residue_count; i++) |
| 200 | av_freep(&vc->residues[i].classifs); |
| 201 | av_freep(&vc->residues); |
| 202 | av_freep(&vc->modes); |
| 203 | |
| 204 | ff_mdct_end(&vc->mdct[0]); |
| 205 | ff_mdct_end(&vc->mdct[1]); |
| 206 | |
| 207 | if (vc->codebooks) |
| 208 | for (i = 0; i < vc->codebook_count; ++i) { |
| 209 | av_freep(&vc->codebooks[i].codevectors); |
| 210 | ff_free_vlc(&vc->codebooks[i].vlc); |
| 211 | } |
| 212 | av_freep(&vc->codebooks); |
| 213 | |
| 214 | if (vc->floors) |
| 215 | for (i = 0; i < vc->floor_count; ++i) { |
| 216 | if (vc->floors[i].floor_type == 0) { |
| 217 | av_freep(&vc->floors[i].data.t0.map[0]); |
| 218 | av_freep(&vc->floors[i].data.t0.map[1]); |
| 219 | av_freep(&vc->floors[i].data.t0.book_list); |
| 220 | av_freep(&vc->floors[i].data.t0.lsp); |
| 221 | } else { |
| 222 | av_freep(&vc->floors[i].data.t1.list); |
| 223 | } |
| 224 | } |
| 225 | av_freep(&vc->floors); |
| 226 | |
| 227 | if (vc->mappings) |
| 228 | for (i = 0; i < vc->mapping_count; ++i) { |
| 229 | av_freep(&vc->mappings[i].magnitude); |
| 230 | av_freep(&vc->mappings[i].angle); |
| 231 | av_freep(&vc->mappings[i].mux); |
| 232 | } |
| 233 | av_freep(&vc->mappings); |
| 234 | } |
| 235 | |
| 236 | // Parse setup header ------------------------------------------------- |
| 237 | |
| 238 | // Process codebooks part |
| 239 | |
| 240 | static int vorbis_parse_setup_hdr_codebooks(vorbis_context *vc) |
| 241 | { |
| 242 | unsigned cb; |
| 243 | uint8_t *tmp_vlc_bits = NULL; |
| 244 | uint32_t *tmp_vlc_codes = NULL; |
| 245 | GetBitContext *gb = &vc->gb; |
| 246 | uint16_t *codebook_multiplicands = NULL; |
| 247 | int ret = 0; |
| 248 | |
| 249 | vc->codebook_count = get_bits(gb, 8) + 1; |
| 250 | |
| 251 | av_dlog(NULL, " Codebooks: %d \n", vc->codebook_count); |
| 252 | |
| 253 | vc->codebooks = av_mallocz(vc->codebook_count * sizeof(*vc->codebooks)); |
| 254 | tmp_vlc_bits = av_mallocz(V_MAX_VLCS * sizeof(*tmp_vlc_bits)); |
| 255 | tmp_vlc_codes = av_mallocz(V_MAX_VLCS * sizeof(*tmp_vlc_codes)); |
| 256 | codebook_multiplicands = av_malloc(V_MAX_VLCS * sizeof(*codebook_multiplicands)); |
| 257 | if (!vc->codebooks || |
| 258 | !tmp_vlc_bits || !tmp_vlc_codes || !codebook_multiplicands) { |
| 259 | ret = AVERROR(ENOMEM); |
| 260 | goto error; |
| 261 | } |
| 262 | |
| 263 | for (cb = 0; cb < vc->codebook_count; ++cb) { |
| 264 | vorbis_codebook *codebook_setup = &vc->codebooks[cb]; |
| 265 | unsigned ordered, t, entries, used_entries = 0; |
| 266 | |
| 267 | av_dlog(NULL, " %u. Codebook\n", cb); |
| 268 | |
| 269 | if (get_bits(gb, 24) != 0x564342) { |
| 270 | av_log(vc->avctx, AV_LOG_ERROR, |
| 271 | " %u. Codebook setup data corrupt.\n", cb); |
| 272 | ret = AVERROR_INVALIDDATA; |
| 273 | goto error; |
| 274 | } |
| 275 | |
| 276 | codebook_setup->dimensions=get_bits(gb, 16); |
| 277 | if (codebook_setup->dimensions > 16 || codebook_setup->dimensions == 0) { |
| 278 | av_log(vc->avctx, AV_LOG_ERROR, |
| 279 | " %u. Codebook's dimension is invalid (%d).\n", |
| 280 | cb, codebook_setup->dimensions); |
| 281 | ret = AVERROR_INVALIDDATA; |
| 282 | goto error; |
| 283 | } |
| 284 | entries = get_bits(gb, 24); |
| 285 | if (entries > V_MAX_VLCS) { |
| 286 | av_log(vc->avctx, AV_LOG_ERROR, |
| 287 | " %u. Codebook has too many entries (%u).\n", |
| 288 | cb, entries); |
| 289 | ret = AVERROR_INVALIDDATA; |
| 290 | goto error; |
| 291 | } |
| 292 | |
| 293 | ordered = get_bits1(gb); |
| 294 | |
| 295 | av_dlog(NULL, " codebook_dimensions %d, codebook_entries %u\n", |
| 296 | codebook_setup->dimensions, entries); |
| 297 | |
| 298 | if (!ordered) { |
| 299 | unsigned ce, flag; |
| 300 | unsigned sparse = get_bits1(gb); |
| 301 | |
| 302 | av_dlog(NULL, " not ordered \n"); |
| 303 | |
| 304 | if (sparse) { |
| 305 | av_dlog(NULL, " sparse \n"); |
| 306 | |
| 307 | used_entries = 0; |
| 308 | for (ce = 0; ce < entries; ++ce) { |
| 309 | flag = get_bits1(gb); |
| 310 | if (flag) { |
| 311 | tmp_vlc_bits[ce] = get_bits(gb, 5) + 1; |
| 312 | ++used_entries; |
| 313 | } else |
| 314 | tmp_vlc_bits[ce] = 0; |
| 315 | } |
| 316 | } else { |
| 317 | av_dlog(NULL, " not sparse \n"); |
| 318 | |
| 319 | used_entries = entries; |
| 320 | for (ce = 0; ce < entries; ++ce) |
| 321 | tmp_vlc_bits[ce] = get_bits(gb, 5) + 1; |
| 322 | } |
| 323 | } else { |
| 324 | unsigned current_entry = 0; |
| 325 | unsigned current_length = get_bits(gb, 5) + 1; |
| 326 | |
| 327 | av_dlog(NULL, " ordered, current length: %u\n", current_length); //FIXME |
| 328 | |
| 329 | used_entries = entries; |
| 330 | for (; current_entry < used_entries && current_length <= 32; ++current_length) { |
| 331 | unsigned i, number; |
| 332 | |
| 333 | av_dlog(NULL, " number bits: %u ", ilog(entries - current_entry)); |
| 334 | |
| 335 | number = get_bits(gb, ilog(entries - current_entry)); |
| 336 | |
| 337 | av_dlog(NULL, " number: %u\n", number); |
| 338 | |
| 339 | for (i = current_entry; i < number+current_entry; ++i) |
| 340 | if (i < used_entries) |
| 341 | tmp_vlc_bits[i] = current_length; |
| 342 | |
| 343 | current_entry+=number; |
| 344 | } |
| 345 | if (current_entry>used_entries) { |
| 346 | av_log(vc->avctx, AV_LOG_ERROR, " More codelengths than codes in codebook. \n"); |
| 347 | ret = AVERROR_INVALIDDATA; |
| 348 | goto error; |
| 349 | } |
| 350 | } |
| 351 | |
| 352 | codebook_setup->lookup_type = get_bits(gb, 4); |
| 353 | |
| 354 | av_dlog(NULL, " lookup type: %d : %s \n", codebook_setup->lookup_type, |
| 355 | codebook_setup->lookup_type ? "vq" : "no lookup"); |
| 356 | |
| 357 | // If the codebook is used for (inverse) VQ, calculate codevectors. |
| 358 | |
| 359 | if (codebook_setup->lookup_type == 1) { |
| 360 | unsigned i, j, k; |
| 361 | unsigned codebook_lookup_values = ff_vorbis_nth_root(entries, codebook_setup->dimensions); |
| 362 | |
| 363 | float codebook_minimum_value = vorbisfloat2float(get_bits_long(gb, 32)); |
| 364 | float codebook_delta_value = vorbisfloat2float(get_bits_long(gb, 32)); |
| 365 | unsigned codebook_value_bits = get_bits(gb, 4) + 1; |
| 366 | unsigned codebook_sequence_p = get_bits1(gb); |
| 367 | |
| 368 | av_dlog(NULL, " We expect %d numbers for building the codevectors. \n", |
| 369 | codebook_lookup_values); |
| 370 | av_dlog(NULL, " delta %f minmum %f \n", |
| 371 | codebook_delta_value, codebook_minimum_value); |
| 372 | |
| 373 | for (i = 0; i < codebook_lookup_values; ++i) { |
| 374 | codebook_multiplicands[i] = get_bits(gb, codebook_value_bits); |
| 375 | |
| 376 | av_dlog(NULL, " multiplicands*delta+minmum : %e \n", |
| 377 | (float)codebook_multiplicands[i] * codebook_delta_value + codebook_minimum_value); |
| 378 | av_dlog(NULL, " multiplicand %u\n", codebook_multiplicands[i]); |
| 379 | } |
| 380 | |
| 381 | // Weed out unused vlcs and build codevector vector |
| 382 | codebook_setup->codevectors = used_entries ? av_mallocz_array(used_entries, |
| 383 | codebook_setup->dimensions * |
| 384 | sizeof(*codebook_setup->codevectors)) |
| 385 | : NULL; |
| 386 | for (j = 0, i = 0; i < entries; ++i) { |
| 387 | unsigned dim = codebook_setup->dimensions; |
| 388 | |
| 389 | if (tmp_vlc_bits[i]) { |
| 390 | float last = 0.0; |
| 391 | unsigned lookup_offset = i; |
| 392 | |
| 393 | av_dlog(vc->avctx, "Lookup offset %u ,", i); |
| 394 | |
| 395 | for (k = 0; k < dim; ++k) { |
| 396 | unsigned multiplicand_offset = lookup_offset % codebook_lookup_values; |
| 397 | codebook_setup->codevectors[j * dim + k] = codebook_multiplicands[multiplicand_offset] * codebook_delta_value + codebook_minimum_value + last; |
| 398 | if (codebook_sequence_p) |
| 399 | last = codebook_setup->codevectors[j * dim + k]; |
| 400 | lookup_offset/=codebook_lookup_values; |
| 401 | } |
| 402 | tmp_vlc_bits[j] = tmp_vlc_bits[i]; |
| 403 | |
| 404 | av_dlog(vc->avctx, "real lookup offset %u, vector: ", j); |
| 405 | for (k = 0; k < dim; ++k) |
| 406 | av_dlog(vc->avctx, " %f ", |
| 407 | codebook_setup->codevectors[j * dim + k]); |
| 408 | av_dlog(vc->avctx, "\n"); |
| 409 | |
| 410 | ++j; |
| 411 | } |
| 412 | } |
| 413 | if (j != used_entries) { |
| 414 | av_log(vc->avctx, AV_LOG_ERROR, "Bug in codevector vector building code. \n"); |
| 415 | ret = AVERROR_INVALIDDATA; |
| 416 | goto error; |
| 417 | } |
| 418 | entries = used_entries; |
| 419 | } else if (codebook_setup->lookup_type >= 2) { |
| 420 | av_log(vc->avctx, AV_LOG_ERROR, "Codebook lookup type not supported. \n"); |
| 421 | ret = AVERROR_INVALIDDATA; |
| 422 | goto error; |
| 423 | } |
| 424 | |
| 425 | // Initialize VLC table |
| 426 | if (ff_vorbis_len2vlc(tmp_vlc_bits, tmp_vlc_codes, entries)) { |
| 427 | av_log(vc->avctx, AV_LOG_ERROR, " Invalid code lengths while generating vlcs. \n"); |
| 428 | ret = AVERROR_INVALIDDATA; |
| 429 | goto error; |
| 430 | } |
| 431 | codebook_setup->maxdepth = 0; |
| 432 | for (t = 0; t < entries; ++t) |
| 433 | if (tmp_vlc_bits[t] >= codebook_setup->maxdepth) |
| 434 | codebook_setup->maxdepth = tmp_vlc_bits[t]; |
| 435 | |
| 436 | if (codebook_setup->maxdepth > 3 * V_NB_BITS) |
| 437 | codebook_setup->nb_bits = V_NB_BITS2; |
| 438 | else |
| 439 | codebook_setup->nb_bits = V_NB_BITS; |
| 440 | |
| 441 | codebook_setup->maxdepth = (codebook_setup->maxdepth+codebook_setup->nb_bits - 1) / codebook_setup->nb_bits; |
| 442 | |
| 443 | if ((ret = init_vlc(&codebook_setup->vlc, codebook_setup->nb_bits, |
| 444 | entries, tmp_vlc_bits, sizeof(*tmp_vlc_bits), |
| 445 | sizeof(*tmp_vlc_bits), tmp_vlc_codes, |
| 446 | sizeof(*tmp_vlc_codes), sizeof(*tmp_vlc_codes), |
| 447 | INIT_VLC_LE))) { |
| 448 | av_log(vc->avctx, AV_LOG_ERROR, " Error generating vlc tables. \n"); |
| 449 | goto error; |
| 450 | } |
| 451 | } |
| 452 | |
| 453 | av_free(tmp_vlc_bits); |
| 454 | av_free(tmp_vlc_codes); |
| 455 | av_free(codebook_multiplicands); |
| 456 | return 0; |
| 457 | |
| 458 | // Error: |
| 459 | error: |
| 460 | av_free(tmp_vlc_bits); |
| 461 | av_free(tmp_vlc_codes); |
| 462 | av_free(codebook_multiplicands); |
| 463 | return ret; |
| 464 | } |
| 465 | |
| 466 | // Process time domain transforms part (unused in Vorbis I) |
| 467 | |
| 468 | static int vorbis_parse_setup_hdr_tdtransforms(vorbis_context *vc) |
| 469 | { |
| 470 | GetBitContext *gb = &vc->gb; |
| 471 | unsigned i, vorbis_time_count = get_bits(gb, 6) + 1; |
| 472 | |
| 473 | for (i = 0; i < vorbis_time_count; ++i) { |
| 474 | unsigned vorbis_tdtransform = get_bits(gb, 16); |
| 475 | |
| 476 | av_dlog(NULL, " Vorbis time domain transform %u: %u\n", |
| 477 | vorbis_time_count, vorbis_tdtransform); |
| 478 | |
| 479 | if (vorbis_tdtransform) { |
| 480 | av_log(vc->avctx, AV_LOG_ERROR, "Vorbis time domain transform data nonzero. \n"); |
| 481 | return AVERROR_INVALIDDATA; |
| 482 | } |
| 483 | } |
| 484 | return 0; |
| 485 | } |
| 486 | |
| 487 | // Process floors part |
| 488 | |
| 489 | static int vorbis_floor0_decode(vorbis_context *vc, |
| 490 | vorbis_floor_data *vfu, float *vec); |
| 491 | static int create_map(vorbis_context *vc, unsigned floor_number); |
| 492 | static int vorbis_floor1_decode(vorbis_context *vc, |
| 493 | vorbis_floor_data *vfu, float *vec); |
| 494 | static int vorbis_parse_setup_hdr_floors(vorbis_context *vc) |
| 495 | { |
| 496 | GetBitContext *gb = &vc->gb; |
| 497 | int i, j, k, ret; |
| 498 | |
| 499 | vc->floor_count = get_bits(gb, 6) + 1; |
| 500 | |
| 501 | vc->floors = av_mallocz(vc->floor_count * sizeof(*vc->floors)); |
| 502 | if (!vc->floors) |
| 503 | return AVERROR(ENOMEM); |
| 504 | |
| 505 | for (i = 0; i < vc->floor_count; ++i) { |
| 506 | vorbis_floor *floor_setup = &vc->floors[i]; |
| 507 | |
| 508 | floor_setup->floor_type = get_bits(gb, 16); |
| 509 | |
| 510 | av_dlog(NULL, " %d. floor type %d \n", i, floor_setup->floor_type); |
| 511 | |
| 512 | if (floor_setup->floor_type == 1) { |
| 513 | int maximum_class = -1; |
| 514 | unsigned rangebits, rangemax, floor1_values = 2; |
| 515 | |
| 516 | floor_setup->decode = vorbis_floor1_decode; |
| 517 | |
| 518 | floor_setup->data.t1.partitions = get_bits(gb, 5); |
| 519 | |
| 520 | av_dlog(NULL, " %d.floor: %d partitions \n", |
| 521 | i, floor_setup->data.t1.partitions); |
| 522 | |
| 523 | for (j = 0; j < floor_setup->data.t1.partitions; ++j) { |
| 524 | floor_setup->data.t1.partition_class[j] = get_bits(gb, 4); |
| 525 | if (floor_setup->data.t1.partition_class[j] > maximum_class) |
| 526 | maximum_class = floor_setup->data.t1.partition_class[j]; |
| 527 | |
| 528 | av_dlog(NULL, " %d. floor %d partition class %d \n", |
| 529 | i, j, floor_setup->data.t1.partition_class[j]); |
| 530 | |
| 531 | } |
| 532 | |
| 533 | av_dlog(NULL, " maximum class %d \n", maximum_class); |
| 534 | |
| 535 | for (j = 0; j <= maximum_class; ++j) { |
| 536 | floor_setup->data.t1.class_dimensions[j] = get_bits(gb, 3) + 1; |
| 537 | floor_setup->data.t1.class_subclasses[j] = get_bits(gb, 2); |
| 538 | |
| 539 | av_dlog(NULL, " %d floor %d class dim: %d subclasses %d \n", i, j, |
| 540 | floor_setup->data.t1.class_dimensions[j], |
| 541 | floor_setup->data.t1.class_subclasses[j]); |
| 542 | |
| 543 | if (floor_setup->data.t1.class_subclasses[j]) { |
| 544 | GET_VALIDATED_INDEX(floor_setup->data.t1.class_masterbook[j], 8, vc->codebook_count) |
| 545 | |
| 546 | av_dlog(NULL, " masterbook: %d \n", floor_setup->data.t1.class_masterbook[j]); |
| 547 | } |
| 548 | |
| 549 | for (k = 0; k < (1 << floor_setup->data.t1.class_subclasses[j]); ++k) { |
| 550 | int16_t bits = get_bits(gb, 8) - 1; |
| 551 | if (bits != -1) |
| 552 | VALIDATE_INDEX(bits, vc->codebook_count) |
| 553 | floor_setup->data.t1.subclass_books[j][k] = bits; |
| 554 | |
| 555 | av_dlog(NULL, " book %d. : %d \n", k, floor_setup->data.t1.subclass_books[j][k]); |
| 556 | } |
| 557 | } |
| 558 | |
| 559 | floor_setup->data.t1.multiplier = get_bits(gb, 2) + 1; |
| 560 | floor_setup->data.t1.x_list_dim = 2; |
| 561 | |
| 562 | for (j = 0; j < floor_setup->data.t1.partitions; ++j) |
| 563 | floor_setup->data.t1.x_list_dim+=floor_setup->data.t1.class_dimensions[floor_setup->data.t1.partition_class[j]]; |
| 564 | |
| 565 | floor_setup->data.t1.list = av_mallocz_array(floor_setup->data.t1.x_list_dim, |
| 566 | sizeof(*floor_setup->data.t1.list)); |
| 567 | if (!floor_setup->data.t1.list) |
| 568 | return AVERROR(ENOMEM); |
| 569 | |
| 570 | rangebits = get_bits(gb, 4); |
| 571 | rangemax = (1 << rangebits); |
| 572 | if (rangemax > vc->blocksize[1] / 2) { |
| 573 | av_log(vc->avctx, AV_LOG_ERROR, |
| 574 | "Floor value is too large for blocksize: %u (%"PRIu32")\n", |
| 575 | rangemax, vc->blocksize[1] / 2); |
| 576 | return AVERROR_INVALIDDATA; |
| 577 | } |
| 578 | floor_setup->data.t1.list[0].x = 0; |
| 579 | floor_setup->data.t1.list[1].x = rangemax; |
| 580 | |
| 581 | for (j = 0; j < floor_setup->data.t1.partitions; ++j) { |
| 582 | for (k = 0; k < floor_setup->data.t1.class_dimensions[floor_setup->data.t1.partition_class[j]]; ++k, ++floor1_values) { |
| 583 | floor_setup->data.t1.list[floor1_values].x = get_bits(gb, rangebits); |
| 584 | |
| 585 | av_dlog(NULL, " %u. floor1 Y coord. %d\n", floor1_values, |
| 586 | floor_setup->data.t1.list[floor1_values].x); |
| 587 | } |
| 588 | } |
| 589 | |
| 590 | // Precalculate order of x coordinates - needed for decode |
| 591 | if (ff_vorbis_ready_floor1_list(vc->avctx, |
| 592 | floor_setup->data.t1.list, |
| 593 | floor_setup->data.t1.x_list_dim)) { |
| 594 | return AVERROR_INVALIDDATA; |
| 595 | } |
| 596 | } else if (floor_setup->floor_type == 0) { |
| 597 | unsigned max_codebook_dim = 0; |
| 598 | |
| 599 | floor_setup->decode = vorbis_floor0_decode; |
| 600 | |
| 601 | floor_setup->data.t0.order = get_bits(gb, 8); |
| 602 | if (!floor_setup->data.t0.order) { |
| 603 | av_log(vc->avctx, AV_LOG_ERROR, "Floor 0 order is 0.\n"); |
| 604 | return AVERROR_INVALIDDATA; |
| 605 | } |
| 606 | floor_setup->data.t0.rate = get_bits(gb, 16); |
| 607 | if (!floor_setup->data.t0.rate) { |
| 608 | av_log(vc->avctx, AV_LOG_ERROR, "Floor 0 rate is 0.\n"); |
| 609 | return AVERROR_INVALIDDATA; |
| 610 | } |
| 611 | floor_setup->data.t0.bark_map_size = get_bits(gb, 16); |
| 612 | if (!floor_setup->data.t0.bark_map_size) { |
| 613 | av_log(vc->avctx, AV_LOG_ERROR, |
| 614 | "Floor 0 bark map size is 0.\n"); |
| 615 | return AVERROR_INVALIDDATA; |
| 616 | } |
| 617 | floor_setup->data.t0.amplitude_bits = get_bits(gb, 6); |
| 618 | floor_setup->data.t0.amplitude_offset = get_bits(gb, 8); |
| 619 | floor_setup->data.t0.num_books = get_bits(gb, 4) + 1; |
| 620 | |
| 621 | /* allocate mem for booklist */ |
| 622 | floor_setup->data.t0.book_list = |
| 623 | av_malloc(floor_setup->data.t0.num_books); |
| 624 | if (!floor_setup->data.t0.book_list) |
| 625 | return AVERROR(ENOMEM); |
| 626 | /* read book indexes */ |
| 627 | { |
| 628 | int idx; |
| 629 | unsigned book_idx; |
| 630 | for (idx = 0; idx < floor_setup->data.t0.num_books; ++idx) { |
| 631 | GET_VALIDATED_INDEX(book_idx, 8, vc->codebook_count) |
| 632 | floor_setup->data.t0.book_list[idx] = book_idx; |
| 633 | if (vc->codebooks[book_idx].dimensions > max_codebook_dim) |
| 634 | max_codebook_dim = vc->codebooks[book_idx].dimensions; |
| 635 | } |
| 636 | } |
| 637 | |
| 638 | if ((ret = create_map(vc, i)) < 0) |
| 639 | return ret; |
| 640 | |
| 641 | /* codebook dim is for padding if codebook dim doesn't * |
| 642 | * divide order+1 then we need to read more data */ |
| 643 | floor_setup->data.t0.lsp = |
| 644 | av_malloc_array((floor_setup->data.t0.order + 1 + max_codebook_dim), |
| 645 | sizeof(*floor_setup->data.t0.lsp)); |
| 646 | if (!floor_setup->data.t0.lsp) |
| 647 | return AVERROR(ENOMEM); |
| 648 | |
| 649 | /* debug output parsed headers */ |
| 650 | av_dlog(NULL, "floor0 order: %u\n", floor_setup->data.t0.order); |
| 651 | av_dlog(NULL, "floor0 rate: %u\n", floor_setup->data.t0.rate); |
| 652 | av_dlog(NULL, "floor0 bark map size: %u\n", |
| 653 | floor_setup->data.t0.bark_map_size); |
| 654 | av_dlog(NULL, "floor0 amplitude bits: %u\n", |
| 655 | floor_setup->data.t0.amplitude_bits); |
| 656 | av_dlog(NULL, "floor0 amplitude offset: %u\n", |
| 657 | floor_setup->data.t0.amplitude_offset); |
| 658 | av_dlog(NULL, "floor0 number of books: %u\n", |
| 659 | floor_setup->data.t0.num_books); |
| 660 | av_dlog(NULL, "floor0 book list pointer: %p\n", |
| 661 | floor_setup->data.t0.book_list); |
| 662 | { |
| 663 | int idx; |
| 664 | for (idx = 0; idx < floor_setup->data.t0.num_books; ++idx) { |
| 665 | av_dlog(NULL, " Book %d: %u\n", idx + 1, |
| 666 | floor_setup->data.t0.book_list[idx]); |
| 667 | } |
| 668 | } |
| 669 | } else { |
| 670 | av_log(vc->avctx, AV_LOG_ERROR, "Invalid floor type!\n"); |
| 671 | return AVERROR_INVALIDDATA; |
| 672 | } |
| 673 | } |
| 674 | return 0; |
| 675 | } |
| 676 | |
| 677 | // Process residues part |
| 678 | |
| 679 | static int vorbis_parse_setup_hdr_residues(vorbis_context *vc) |
| 680 | { |
| 681 | GetBitContext *gb = &vc->gb; |
| 682 | unsigned i, j, k; |
| 683 | |
| 684 | vc->residue_count = get_bits(gb, 6)+1; |
| 685 | vc->residues = av_mallocz(vc->residue_count * sizeof(*vc->residues)); |
| 686 | if (!vc->residues) |
| 687 | return AVERROR(ENOMEM); |
| 688 | |
| 689 | av_dlog(NULL, " There are %d residues. \n", vc->residue_count); |
| 690 | |
| 691 | for (i = 0; i < vc->residue_count; ++i) { |
| 692 | vorbis_residue *res_setup = &vc->residues[i]; |
| 693 | uint8_t cascade[64]; |
| 694 | unsigned high_bits, low_bits; |
| 695 | |
| 696 | res_setup->type = get_bits(gb, 16); |
| 697 | |
| 698 | av_dlog(NULL, " %u. residue type %d\n", i, res_setup->type); |
| 699 | |
| 700 | res_setup->begin = get_bits(gb, 24); |
| 701 | res_setup->end = get_bits(gb, 24); |
| 702 | res_setup->partition_size = get_bits(gb, 24) + 1; |
| 703 | /* Validations to prevent a buffer overflow later. */ |
| 704 | if (res_setup->begin>res_setup->end || |
| 705 | (res_setup->end-res_setup->begin) / res_setup->partition_size > FFMIN(V_MAX_PARTITIONS, 65535)) { |
| 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_array(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_array(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_array(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_array(vc->blocksize[1] / 2, vc->audio_channels * sizeof(*vc->channel_residues)); |
| 988 | vc->saved = av_mallocz_array(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 = -1; |
| 993 | |
| 994 | ff_mdct_init(&vc->mdct[0], bl0, 1, -1.0); |
| 995 | ff_mdct_init(&vc->mdct[1], bl1, 1, -1.0); |
| 996 | vc->fdsp = avpriv_float_dsp_alloc(vc->avctx->flags & CODEC_FLAG_BITEXACT); |
| 997 | if (!vc->fdsp) |
| 998 | return AVERROR(ENOMEM); |
| 999 | |
| 1000 | 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 ", |
| 1001 | vc->version, vc->audio_channels, vc->audio_samplerate, vc->bitrate_maximum, vc->bitrate_nominal, vc->bitrate_minimum, vc->blocksize[0], vc->blocksize[1]); |
| 1002 | |
| 1003 | /* |
| 1004 | BLK = vc->blocksize[0]; |
| 1005 | for (i = 0; i < BLK / 2; ++i) { |
| 1006 | 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))); |
| 1007 | } |
| 1008 | */ |
| 1009 | |
| 1010 | return 0; |
| 1011 | } |
| 1012 | |
| 1013 | // Process the extradata using the functions above (identification header, setup header) |
| 1014 | |
| 1015 | static av_cold int vorbis_decode_init(AVCodecContext *avctx) |
| 1016 | { |
| 1017 | vorbis_context *vc = avctx->priv_data; |
| 1018 | uint8_t *headers = avctx->extradata; |
| 1019 | int headers_len = avctx->extradata_size; |
| 1020 | uint8_t *header_start[3]; |
| 1021 | int header_len[3]; |
| 1022 | GetBitContext *gb = &vc->gb; |
| 1023 | int hdr_type, ret; |
| 1024 | |
| 1025 | vc->avctx = avctx; |
| 1026 | ff_vorbisdsp_init(&vc->dsp); |
| 1027 | ff_fmt_convert_init(&vc->fmt_conv, avctx); |
| 1028 | |
| 1029 | avctx->sample_fmt = AV_SAMPLE_FMT_FLTP; |
| 1030 | |
| 1031 | if (!headers_len) { |
| 1032 | av_log(avctx, AV_LOG_ERROR, "Extradata missing.\n"); |
| 1033 | return AVERROR_INVALIDDATA; |
| 1034 | } |
| 1035 | |
| 1036 | if ((ret = avpriv_split_xiph_headers(headers, headers_len, 30, header_start, header_len)) < 0) { |
| 1037 | av_log(avctx, AV_LOG_ERROR, "Extradata corrupt.\n"); |
| 1038 | return ret; |
| 1039 | } |
| 1040 | |
| 1041 | init_get_bits(gb, header_start[0], header_len[0]*8); |
| 1042 | hdr_type = get_bits(gb, 8); |
| 1043 | if (hdr_type != 1) { |
| 1044 | av_log(avctx, AV_LOG_ERROR, "First header is not the id header.\n"); |
| 1045 | return AVERROR_INVALIDDATA; |
| 1046 | } |
| 1047 | if ((ret = vorbis_parse_id_hdr(vc))) { |
| 1048 | av_log(avctx, AV_LOG_ERROR, "Id header corrupt.\n"); |
| 1049 | vorbis_free(vc); |
| 1050 | return ret; |
| 1051 | } |
| 1052 | |
| 1053 | init_get_bits(gb, header_start[2], header_len[2]*8); |
| 1054 | hdr_type = get_bits(gb, 8); |
| 1055 | if (hdr_type != 5) { |
| 1056 | av_log(avctx, AV_LOG_ERROR, "Third header is not the setup header.\n"); |
| 1057 | vorbis_free(vc); |
| 1058 | return AVERROR_INVALIDDATA; |
| 1059 | } |
| 1060 | if ((ret = vorbis_parse_setup_hdr(vc))) { |
| 1061 | av_log(avctx, AV_LOG_ERROR, "Setup header corrupt.\n"); |
| 1062 | vorbis_free(vc); |
| 1063 | return ret; |
| 1064 | } |
| 1065 | |
| 1066 | if (vc->audio_channels > 8) |
| 1067 | avctx->channel_layout = 0; |
| 1068 | else |
| 1069 | avctx->channel_layout = ff_vorbis_channel_layouts[vc->audio_channels - 1]; |
| 1070 | |
| 1071 | avctx->channels = vc->audio_channels; |
| 1072 | avctx->sample_rate = vc->audio_samplerate; |
| 1073 | |
| 1074 | return 0; |
| 1075 | } |
| 1076 | |
| 1077 | // Decode audiopackets ------------------------------------------------- |
| 1078 | |
| 1079 | // Read and decode floor |
| 1080 | |
| 1081 | static int vorbis_floor0_decode(vorbis_context *vc, |
| 1082 | vorbis_floor_data *vfu, float *vec) |
| 1083 | { |
| 1084 | vorbis_floor0 *vf = &vfu->t0; |
| 1085 | float *lsp = vf->lsp; |
| 1086 | unsigned amplitude, book_idx; |
| 1087 | unsigned blockflag = vc->modes[vc->mode_number].blockflag; |
| 1088 | |
| 1089 | if (!vf->amplitude_bits) |
| 1090 | return 1; |
| 1091 | |
| 1092 | amplitude = get_bits(&vc->gb, vf->amplitude_bits); |
| 1093 | if (amplitude > 0) { |
| 1094 | float last = 0; |
| 1095 | unsigned idx, lsp_len = 0; |
| 1096 | vorbis_codebook codebook; |
| 1097 | |
| 1098 | book_idx = get_bits(&vc->gb, ilog(vf->num_books)); |
| 1099 | if (book_idx >= vf->num_books) { |
| 1100 | av_log(vc->avctx, AV_LOG_ERROR, "floor0 dec: booknumber too high!\n"); |
| 1101 | book_idx = 0; |
| 1102 | } |
| 1103 | av_dlog(NULL, "floor0 dec: booknumber: %u\n", book_idx); |
| 1104 | codebook = vc->codebooks[vf->book_list[book_idx]]; |
| 1105 | /* Invalid codebook! */ |
| 1106 | if (!codebook.codevectors) |
| 1107 | return AVERROR_INVALIDDATA; |
| 1108 | |
| 1109 | while (lsp_len<vf->order) { |
| 1110 | int vec_off; |
| 1111 | |
| 1112 | av_dlog(NULL, "floor0 dec: book dimension: %d\n", codebook.dimensions); |
| 1113 | av_dlog(NULL, "floor0 dec: maximum depth: %d\n", codebook.maxdepth); |
| 1114 | /* read temp vector */ |
| 1115 | vec_off = get_vlc2(&vc->gb, codebook.vlc.table, |
| 1116 | codebook.nb_bits, codebook.maxdepth) |
| 1117 | * codebook.dimensions; |
| 1118 | av_dlog(NULL, "floor0 dec: vector offset: %d\n", vec_off); |
| 1119 | /* copy each vector component and add last to it */ |
| 1120 | for (idx = 0; idx < codebook.dimensions; ++idx) |
| 1121 | lsp[lsp_len+idx] = codebook.codevectors[vec_off+idx] + last; |
| 1122 | last = lsp[lsp_len+idx-1]; /* set last to last vector component */ |
| 1123 | |
| 1124 | lsp_len += codebook.dimensions; |
| 1125 | } |
| 1126 | /* DEBUG: output lsp coeffs */ |
| 1127 | { |
| 1128 | int idx; |
| 1129 | for (idx = 0; idx < lsp_len; ++idx) |
| 1130 | av_dlog(NULL, "floor0 dec: coeff at %d is %f\n", idx, lsp[idx]); |
| 1131 | } |
| 1132 | |
| 1133 | /* synthesize floor output vector */ |
| 1134 | { |
| 1135 | int i; |
| 1136 | int order = vf->order; |
| 1137 | float wstep = M_PI / vf->bark_map_size; |
| 1138 | |
| 1139 | for (i = 0; i < order; i++) |
| 1140 | lsp[i] = 2.0f * cos(lsp[i]); |
| 1141 | |
| 1142 | av_dlog(NULL, "floor0 synth: map_size = %"PRIu32"; m = %d; wstep = %f\n", |
| 1143 | vf->map_size[blockflag], order, wstep); |
| 1144 | |
| 1145 | i = 0; |
| 1146 | while (i < vf->map_size[blockflag]) { |
| 1147 | int j, iter_cond = vf->map[blockflag][i]; |
| 1148 | float p = 0.5f; |
| 1149 | float q = 0.5f; |
| 1150 | float two_cos_w = 2.0f * cos(wstep * iter_cond); // needed all times |
| 1151 | |
| 1152 | /* similar part for the q and p products */ |
| 1153 | for (j = 0; j + 1 < order; j += 2) { |
| 1154 | q *= lsp[j] - two_cos_w; |
| 1155 | p *= lsp[j + 1] - two_cos_w; |
| 1156 | } |
| 1157 | if (j == order) { // even order |
| 1158 | p *= p * (2.0f - two_cos_w); |
| 1159 | q *= q * (2.0f + two_cos_w); |
| 1160 | } else { // odd order |
| 1161 | q *= two_cos_w-lsp[j]; // one more time for q |
| 1162 | |
| 1163 | /* final step and square */ |
| 1164 | p *= p * (4.f - two_cos_w * two_cos_w); |
| 1165 | q *= q; |
| 1166 | } |
| 1167 | |
| 1168 | /* calculate linear floor value */ |
| 1169 | q = exp((((amplitude*vf->amplitude_offset) / |
| 1170 | (((1 << vf->amplitude_bits) - 1) * sqrt(p + q))) |
| 1171 | - vf->amplitude_offset) * .11512925f); |
| 1172 | |
| 1173 | /* fill vector */ |
| 1174 | do { |
| 1175 | vec[i] = q; ++i; |
| 1176 | } while (vf->map[blockflag][i] == iter_cond); |
| 1177 | } |
| 1178 | } |
| 1179 | } else { |
| 1180 | /* this channel is unused */ |
| 1181 | return 1; |
| 1182 | } |
| 1183 | |
| 1184 | av_dlog(NULL, " Floor0 decoded\n"); |
| 1185 | |
| 1186 | return 0; |
| 1187 | } |
| 1188 | |
| 1189 | static int vorbis_floor1_decode(vorbis_context *vc, |
| 1190 | vorbis_floor_data *vfu, float *vec) |
| 1191 | { |
| 1192 | vorbis_floor1 *vf = &vfu->t1; |
| 1193 | GetBitContext *gb = &vc->gb; |
| 1194 | uint16_t range_v[4] = { 256, 128, 86, 64 }; |
| 1195 | unsigned range = range_v[vf->multiplier - 1]; |
| 1196 | uint16_t floor1_Y[258]; |
| 1197 | uint16_t floor1_Y_final[258]; |
| 1198 | int floor1_flag[258]; |
| 1199 | unsigned partition_class, cdim, cbits, csub, cval, offset, i, j; |
| 1200 | int book, adx, ady, dy, off, predicted, err; |
| 1201 | |
| 1202 | |
| 1203 | if (!get_bits1(gb)) // silence |
| 1204 | return 1; |
| 1205 | |
| 1206 | // Read values (or differences) for the floor's points |
| 1207 | |
| 1208 | floor1_Y[0] = get_bits(gb, ilog(range - 1)); |
| 1209 | floor1_Y[1] = get_bits(gb, ilog(range - 1)); |
| 1210 | |
| 1211 | av_dlog(NULL, "floor 0 Y %d floor 1 Y %d \n", floor1_Y[0], floor1_Y[1]); |
| 1212 | |
| 1213 | offset = 2; |
| 1214 | for (i = 0; i < vf->partitions; ++i) { |
| 1215 | partition_class = vf->partition_class[i]; |
| 1216 | cdim = vf->class_dimensions[partition_class]; |
| 1217 | cbits = vf->class_subclasses[partition_class]; |
| 1218 | csub = (1 << cbits) - 1; |
| 1219 | cval = 0; |
| 1220 | |
| 1221 | av_dlog(NULL, "Cbits %u\n", cbits); |
| 1222 | |
| 1223 | if (cbits) // this reads all subclasses for this partition's class |
| 1224 | cval = get_vlc2(gb, vc->codebooks[vf->class_masterbook[partition_class]].vlc.table, |
| 1225 | vc->codebooks[vf->class_masterbook[partition_class]].nb_bits, 3); |
| 1226 | |
| 1227 | for (j = 0; j < cdim; ++j) { |
| 1228 | book = vf->subclass_books[partition_class][cval & csub]; |
| 1229 | |
| 1230 | av_dlog(NULL, "book %d Cbits %u cval %u bits:%d\n", |
| 1231 | book, cbits, cval, get_bits_count(gb)); |
| 1232 | |
| 1233 | cval = cval >> cbits; |
| 1234 | if (book > -1) { |
| 1235 | int v = get_vlc2(gb, vc->codebooks[book].vlc.table, |
| 1236 | vc->codebooks[book].nb_bits, 3); |
| 1237 | if (v < 0) |
| 1238 | return AVERROR_INVALIDDATA; |
| 1239 | floor1_Y[offset+j] = v; |
| 1240 | } else { |
| 1241 | floor1_Y[offset+j] = 0; |
| 1242 | } |
| 1243 | |
| 1244 | av_dlog(NULL, " floor(%d) = %d \n", |
| 1245 | vf->list[offset+j].x, floor1_Y[offset+j]); |
| 1246 | } |
| 1247 | offset+=cdim; |
| 1248 | } |
| 1249 | |
| 1250 | // Amplitude calculation from the differences |
| 1251 | |
| 1252 | floor1_flag[0] = 1; |
| 1253 | floor1_flag[1] = 1; |
| 1254 | floor1_Y_final[0] = floor1_Y[0]; |
| 1255 | floor1_Y_final[1] = floor1_Y[1]; |
| 1256 | |
| 1257 | for (i = 2; i < vf->x_list_dim; ++i) { |
| 1258 | unsigned val, highroom, lowroom, room, high_neigh_offs, low_neigh_offs; |
| 1259 | |
| 1260 | low_neigh_offs = vf->list[i].low; |
| 1261 | high_neigh_offs = vf->list[i].high; |
| 1262 | dy = floor1_Y_final[high_neigh_offs] - floor1_Y_final[low_neigh_offs]; // render_point begin |
| 1263 | adx = vf->list[high_neigh_offs].x - vf->list[low_neigh_offs].x; |
| 1264 | ady = FFABS(dy); |
| 1265 | err = ady * (vf->list[i].x - vf->list[low_neigh_offs].x); |
| 1266 | off = err / adx; |
| 1267 | if (dy < 0) { |
| 1268 | predicted = floor1_Y_final[low_neigh_offs] - off; |
| 1269 | } else { |
| 1270 | predicted = floor1_Y_final[low_neigh_offs] + off; |
| 1271 | } // render_point end |
| 1272 | |
| 1273 | val = floor1_Y[i]; |
| 1274 | highroom = range-predicted; |
| 1275 | lowroom = predicted; |
| 1276 | if (highroom < lowroom) { |
| 1277 | room = highroom * 2; |
| 1278 | } else { |
| 1279 | room = lowroom * 2; // SPEC misspelling |
| 1280 | } |
| 1281 | if (val) { |
| 1282 | floor1_flag[low_neigh_offs] = 1; |
| 1283 | floor1_flag[high_neigh_offs] = 1; |
| 1284 | floor1_flag[i] = 1; |
| 1285 | if (val >= room) { |
| 1286 | if (highroom > lowroom) { |
| 1287 | floor1_Y_final[i] = av_clip_uint16(val - lowroom + predicted); |
| 1288 | } else { |
| 1289 | floor1_Y_final[i] = av_clip_uint16(predicted - val + highroom - 1); |
| 1290 | } |
| 1291 | } else { |
| 1292 | if (val & 1) { |
| 1293 | floor1_Y_final[i] = av_clip_uint16(predicted - (val + 1) / 2); |
| 1294 | } else { |
| 1295 | floor1_Y_final[i] = av_clip_uint16(predicted + val / 2); |
| 1296 | } |
| 1297 | } |
| 1298 | } else { |
| 1299 | floor1_flag[i] = 0; |
| 1300 | floor1_Y_final[i] = av_clip_uint16(predicted); |
| 1301 | } |
| 1302 | |
| 1303 | av_dlog(NULL, " Decoded floor(%d) = %u / val %u\n", |
| 1304 | vf->list[i].x, floor1_Y_final[i], val); |
| 1305 | } |
| 1306 | |
| 1307 | // Curve synth - connect the calculated dots and convert from dB scale FIXME optimize ? |
| 1308 | |
| 1309 | ff_vorbis_floor1_render_list(vf->list, vf->x_list_dim, floor1_Y_final, floor1_flag, vf->multiplier, vec, vf->list[1].x); |
| 1310 | |
| 1311 | av_dlog(NULL, " Floor decoded\n"); |
| 1312 | |
| 1313 | return 0; |
| 1314 | } |
| 1315 | |
| 1316 | static av_always_inline int setup_classifs(vorbis_context *vc, |
| 1317 | vorbis_residue *vr, |
| 1318 | uint8_t *do_not_decode, |
| 1319 | unsigned ch_used, |
| 1320 | int partition_count, |
| 1321 | int ptns_to_read |
| 1322 | ) |
| 1323 | { |
| 1324 | int p, j, i; |
| 1325 | unsigned c_p_c = vc->codebooks[vr->classbook].dimensions; |
| 1326 | unsigned inverse_class = ff_inverse[vr->classifications]; |
| 1327 | unsigned temp, temp2; |
| 1328 | for (p = 0, j = 0; j < ch_used; ++j) { |
| 1329 | if (!do_not_decode[j]) { |
| 1330 | temp = get_vlc2(&vc->gb, vc->codebooks[vr->classbook].vlc.table, |
| 1331 | vc->codebooks[vr->classbook].nb_bits, 3); |
| 1332 | |
| 1333 | av_dlog(NULL, "Classword: %u\n", temp); |
| 1334 | |
| 1335 | if ((int)temp < 0) |
| 1336 | return temp; |
| 1337 | |
| 1338 | av_assert0(vr->classifications > 1); //needed for inverse[] |
| 1339 | |
| 1340 | if (temp <= 65536) { |
| 1341 | for (i = partition_count + c_p_c - 1; i >= partition_count; i--) { |
| 1342 | temp2 = (((uint64_t)temp) * inverse_class) >> 32; |
| 1343 | |
| 1344 | if (i < ptns_to_read) |
| 1345 | vr->classifs[p + i] = temp - temp2 * vr->classifications; |
| 1346 | temp = temp2; |
| 1347 | } |
| 1348 | } else { |
| 1349 | for (i = partition_count + c_p_c - 1; i >= partition_count; i--) { |
| 1350 | temp2 = temp / vr->classifications; |
| 1351 | |
| 1352 | if (i < ptns_to_read) |
| 1353 | vr->classifs[p + i] = temp - temp2 * vr->classifications; |
| 1354 | temp = temp2; |
| 1355 | } |
| 1356 | } |
| 1357 | } |
| 1358 | p += ptns_to_read; |
| 1359 | } |
| 1360 | return 0; |
| 1361 | } |
| 1362 | // Read and decode residue |
| 1363 | |
| 1364 | static av_always_inline int vorbis_residue_decode_internal(vorbis_context *vc, |
| 1365 | vorbis_residue *vr, |
| 1366 | unsigned ch, |
| 1367 | uint8_t *do_not_decode, |
| 1368 | float *vec, |
| 1369 | unsigned vlen, |
| 1370 | unsigned ch_left, |
| 1371 | int vr_type) |
| 1372 | { |
| 1373 | GetBitContext *gb = &vc->gb; |
| 1374 | unsigned c_p_c = vc->codebooks[vr->classbook].dimensions; |
| 1375 | uint8_t *classifs = vr->classifs; |
| 1376 | unsigned pass, ch_used, i, j, k, l; |
| 1377 | unsigned max_output = (ch - 1) * vlen; |
| 1378 | int ptns_to_read = vr->ptns_to_read; |
| 1379 | int libvorbis_bug = 0; |
| 1380 | |
| 1381 | if (vr_type == 2) { |
| 1382 | for (j = 1; j < ch; ++j) |
| 1383 | do_not_decode[0] &= do_not_decode[j]; // FIXME - clobbering input |
| 1384 | if (do_not_decode[0]) |
| 1385 | return 0; |
| 1386 | ch_used = 1; |
| 1387 | max_output += vr->end / ch; |
| 1388 | } else { |
| 1389 | ch_used = ch; |
| 1390 | max_output += vr->end; |
| 1391 | } |
| 1392 | |
| 1393 | if (max_output > ch_left * vlen) { |
| 1394 | if (max_output <= ch_left * vlen + vr->partition_size*ch_used/ch) { |
| 1395 | ptns_to_read--; |
| 1396 | libvorbis_bug = 1; |
| 1397 | } else { |
| 1398 | av_log(vc->avctx, AV_LOG_ERROR, "Insufficient output buffer\n"); |
| 1399 | return AVERROR_INVALIDDATA; |
| 1400 | } |
| 1401 | } |
| 1402 | |
| 1403 | av_dlog(NULL, " residue type 0/1/2 decode begin, ch: %d cpc %d \n", ch, c_p_c); |
| 1404 | |
| 1405 | for (pass = 0; pass <= vr->maxpass; ++pass) { // FIXME OPTIMIZE? |
| 1406 | int voffset, partition_count, j_times_ptns_to_read; |
| 1407 | |
| 1408 | voffset = vr->begin; |
| 1409 | for (partition_count = 0; partition_count < ptns_to_read;) { // SPEC error |
| 1410 | if (!pass) { |
| 1411 | int ret; |
| 1412 | if ((ret = setup_classifs(vc, vr, do_not_decode, ch_used, partition_count, ptns_to_read)) < 0) |
| 1413 | return ret; |
| 1414 | } |
| 1415 | for (i = 0; (i < c_p_c) && (partition_count < ptns_to_read); ++i) { |
| 1416 | for (j_times_ptns_to_read = 0, j = 0; j < ch_used; ++j) { |
| 1417 | unsigned voffs; |
| 1418 | |
| 1419 | if (!do_not_decode[j]) { |
| 1420 | unsigned vqclass = classifs[j_times_ptns_to_read + partition_count]; |
| 1421 | int vqbook = vr->books[vqclass][pass]; |
| 1422 | |
| 1423 | if (vqbook >= 0 && vc->codebooks[vqbook].codevectors) { |
| 1424 | unsigned coffs; |
| 1425 | unsigned dim = vc->codebooks[vqbook].dimensions; |
| 1426 | unsigned step = FASTDIV(vr->partition_size << 1, dim << 1); |
| 1427 | vorbis_codebook codebook = vc->codebooks[vqbook]; |
| 1428 | |
| 1429 | if (vr_type == 0) { |
| 1430 | |
| 1431 | voffs = voffset+j*vlen; |
| 1432 | for (k = 0; k < step; ++k) { |
| 1433 | coffs = get_vlc2(gb, codebook.vlc.table, codebook.nb_bits, 3) * dim; |
| 1434 | for (l = 0; l < dim; ++l) |
| 1435 | vec[voffs + k + l * step] += codebook.codevectors[coffs + l]; |
| 1436 | } |
| 1437 | } else if (vr_type == 1) { |
| 1438 | voffs = voffset + j * vlen; |
| 1439 | for (k = 0; k < step; ++k) { |
| 1440 | coffs = get_vlc2(gb, codebook.vlc.table, codebook.nb_bits, 3) * dim; |
| 1441 | for (l = 0; l < dim; ++l, ++voffs) { |
| 1442 | vec[voffs]+=codebook.codevectors[coffs+l]; |
| 1443 | |
| 1444 | av_dlog(NULL, " pass %d offs: %d curr: %f change: %f cv offs.: %d \n", |
| 1445 | pass, voffs, vec[voffs], codebook.codevectors[coffs+l], coffs); |
| 1446 | } |
| 1447 | } |
| 1448 | } else if (vr_type == 2 && ch == 2 && (voffset & 1) == 0 && (dim & 1) == 0) { // most frequent case optimized |
| 1449 | voffs = voffset >> 1; |
| 1450 | |
| 1451 | if (dim == 2) { |
| 1452 | for (k = 0; k < step; ++k) { |
| 1453 | coffs = get_vlc2(gb, codebook.vlc.table, codebook.nb_bits, 3) * 2; |
| 1454 | vec[voffs + k ] += codebook.codevectors[coffs ]; |
| 1455 | vec[voffs + k + vlen] += codebook.codevectors[coffs + 1]; |
| 1456 | } |
| 1457 | } else if (dim == 4) { |
| 1458 | for (k = 0; k < step; ++k, voffs += 2) { |
| 1459 | coffs = get_vlc2(gb, codebook.vlc.table, codebook.nb_bits, 3) * 4; |
| 1460 | vec[voffs ] += codebook.codevectors[coffs ]; |
| 1461 | vec[voffs + 1 ] += codebook.codevectors[coffs + 2]; |
| 1462 | vec[voffs + vlen ] += codebook.codevectors[coffs + 1]; |
| 1463 | vec[voffs + vlen + 1] += codebook.codevectors[coffs + 3]; |
| 1464 | } |
| 1465 | } else |
| 1466 | for (k = 0; k < step; ++k) { |
| 1467 | coffs = get_vlc2(gb, codebook.vlc.table, codebook.nb_bits, 3) * dim; |
| 1468 | for (l = 0; l < dim; l += 2, voffs++) { |
| 1469 | vec[voffs ] += codebook.codevectors[coffs + l ]; |
| 1470 | vec[voffs + vlen] += codebook.codevectors[coffs + l + 1]; |
| 1471 | |
| 1472 | av_dlog(NULL, " pass %d offs: %d curr: %f change: %f cv offs.: %d+%d \n", |
| 1473 | pass, voffset / ch + (voffs % ch) * vlen, |
| 1474 | vec[voffset / ch + (voffs % ch) * vlen], |
| 1475 | codebook.codevectors[coffs + l], coffs, l); |
| 1476 | } |
| 1477 | } |
| 1478 | |
| 1479 | } else if (vr_type == 2) { |
| 1480 | unsigned voffs_div = FASTDIV(voffset << 1, ch <<1); |
| 1481 | unsigned voffs_mod = voffset - voffs_div * ch; |
| 1482 | |
| 1483 | for (k = 0; k < step; ++k) { |
| 1484 | coffs = get_vlc2(gb, codebook.vlc.table, codebook.nb_bits, 3) * dim; |
| 1485 | for (l = 0; l < dim; ++l) { |
| 1486 | vec[voffs_div + voffs_mod * vlen] += |
| 1487 | codebook.codevectors[coffs + l]; |
| 1488 | |
| 1489 | av_dlog(NULL, " pass %d offs: %d curr: %f change: %f cv offs.: %d+%d \n", |
| 1490 | pass, voffs_div + voffs_mod * vlen, |
| 1491 | vec[voffs_div + voffs_mod * vlen], |
| 1492 | codebook.codevectors[coffs + l], coffs, l); |
| 1493 | |
| 1494 | if (++voffs_mod == ch) { |
| 1495 | voffs_div++; |
| 1496 | voffs_mod = 0; |
| 1497 | } |
| 1498 | } |
| 1499 | } |
| 1500 | } |
| 1501 | } |
| 1502 | } |
| 1503 | j_times_ptns_to_read += ptns_to_read; |
| 1504 | } |
| 1505 | ++partition_count; |
| 1506 | voffset += vr->partition_size; |
| 1507 | } |
| 1508 | } |
| 1509 | if (libvorbis_bug && !pass) { |
| 1510 | for (j = 0; j < ch_used; ++j) { |
| 1511 | if (!do_not_decode[j]) { |
| 1512 | get_vlc2(&vc->gb, vc->codebooks[vr->classbook].vlc.table, |
| 1513 | vc->codebooks[vr->classbook].nb_bits, 3); |
| 1514 | } |
| 1515 | } |
| 1516 | } |
| 1517 | } |
| 1518 | return 0; |
| 1519 | } |
| 1520 | |
| 1521 | static inline int vorbis_residue_decode(vorbis_context *vc, vorbis_residue *vr, |
| 1522 | unsigned ch, |
| 1523 | uint8_t *do_not_decode, |
| 1524 | float *vec, unsigned vlen, |
| 1525 | unsigned ch_left) |
| 1526 | { |
| 1527 | if (vr->type == 2) |
| 1528 | return vorbis_residue_decode_internal(vc, vr, ch, do_not_decode, vec, vlen, ch_left, 2); |
| 1529 | else if (vr->type == 1) |
| 1530 | return vorbis_residue_decode_internal(vc, vr, ch, do_not_decode, vec, vlen, ch_left, 1); |
| 1531 | else if (vr->type == 0) |
| 1532 | return vorbis_residue_decode_internal(vc, vr, ch, do_not_decode, vec, vlen, ch_left, 0); |
| 1533 | else { |
| 1534 | av_log(vc->avctx, AV_LOG_ERROR, " Invalid residue type while residue decode?! \n"); |
| 1535 | return AVERROR_INVALIDDATA; |
| 1536 | } |
| 1537 | } |
| 1538 | |
| 1539 | void ff_vorbis_inverse_coupling(float *mag, float *ang, intptr_t blocksize) |
| 1540 | { |
| 1541 | int i; |
| 1542 | for (i = 0; i < blocksize; i++) { |
| 1543 | if (mag[i] > 0.0) { |
| 1544 | if (ang[i] > 0.0) { |
| 1545 | ang[i] = mag[i] - ang[i]; |
| 1546 | } else { |
| 1547 | float temp = ang[i]; |
| 1548 | ang[i] = mag[i]; |
| 1549 | mag[i] += temp; |
| 1550 | } |
| 1551 | } else { |
| 1552 | if (ang[i] > 0.0) { |
| 1553 | ang[i] += mag[i]; |
| 1554 | } else { |
| 1555 | float temp = ang[i]; |
| 1556 | ang[i] = mag[i]; |
| 1557 | mag[i] -= temp; |
| 1558 | } |
| 1559 | } |
| 1560 | } |
| 1561 | } |
| 1562 | |
| 1563 | // Decode the audio packet using the functions above |
| 1564 | |
| 1565 | static int vorbis_parse_audio_packet(vorbis_context *vc, float **floor_ptr) |
| 1566 | { |
| 1567 | GetBitContext *gb = &vc->gb; |
| 1568 | FFTContext *mdct; |
| 1569 | int previous_window = vc->previous_window; |
| 1570 | unsigned mode_number, blockflag, blocksize; |
| 1571 | int i, j; |
| 1572 | uint8_t no_residue[255]; |
| 1573 | uint8_t do_not_decode[255]; |
| 1574 | vorbis_mapping *mapping; |
| 1575 | float *ch_res_ptr = vc->channel_residues; |
| 1576 | uint8_t res_chan[255]; |
| 1577 | unsigned res_num = 0; |
| 1578 | int retlen = 0; |
| 1579 | unsigned ch_left = vc->audio_channels; |
| 1580 | unsigned vlen; |
| 1581 | |
| 1582 | if (get_bits1(gb)) { |
| 1583 | av_log(vc->avctx, AV_LOG_ERROR, "Not a Vorbis I audio packet.\n"); |
| 1584 | return AVERROR_INVALIDDATA; // packet type not audio |
| 1585 | } |
| 1586 | |
| 1587 | if (vc->mode_count == 1) { |
| 1588 | mode_number = 0; |
| 1589 | } else { |
| 1590 | GET_VALIDATED_INDEX(mode_number, ilog(vc->mode_count-1), vc->mode_count) |
| 1591 | } |
| 1592 | vc->mode_number = mode_number; |
| 1593 | mapping = &vc->mappings[vc->modes[mode_number].mapping]; |
| 1594 | |
| 1595 | av_dlog(NULL, " Mode number: %u , mapping: %d , blocktype %d\n", mode_number, |
| 1596 | vc->modes[mode_number].mapping, vc->modes[mode_number].blockflag); |
| 1597 | |
| 1598 | blockflag = vc->modes[mode_number].blockflag; |
| 1599 | blocksize = vc->blocksize[blockflag]; |
| 1600 | vlen = blocksize / 2; |
| 1601 | if (blockflag) { |
| 1602 | int code = get_bits(gb, 2); |
| 1603 | if (previous_window < 0) |
| 1604 | previous_window = code>>1; |
| 1605 | } else if (previous_window < 0) |
| 1606 | previous_window = 0; |
| 1607 | |
| 1608 | memset(ch_res_ptr, 0, sizeof(float) * vc->audio_channels * vlen); //FIXME can this be removed ? |
| 1609 | for (i = 0; i < vc->audio_channels; ++i) |
| 1610 | memset(floor_ptr[i], 0, vlen * sizeof(floor_ptr[0][0])); //FIXME can this be removed ? |
| 1611 | |
| 1612 | // Decode floor |
| 1613 | |
| 1614 | for (i = 0; i < vc->audio_channels; ++i) { |
| 1615 | vorbis_floor *floor; |
| 1616 | int ret; |
| 1617 | if (mapping->submaps > 1) { |
| 1618 | floor = &vc->floors[mapping->submap_floor[mapping->mux[i]]]; |
| 1619 | } else { |
| 1620 | floor = &vc->floors[mapping->submap_floor[0]]; |
| 1621 | } |
| 1622 | |
| 1623 | ret = floor->decode(vc, &floor->data, floor_ptr[i]); |
| 1624 | |
| 1625 | if (ret < 0) { |
| 1626 | av_log(vc->avctx, AV_LOG_ERROR, "Invalid codebook in vorbis_floor_decode.\n"); |
| 1627 | return AVERROR_INVALIDDATA; |
| 1628 | } |
| 1629 | no_residue[i] = ret; |
| 1630 | } |
| 1631 | |
| 1632 | // Nonzero vector propagate |
| 1633 | |
| 1634 | for (i = mapping->coupling_steps - 1; i >= 0; --i) { |
| 1635 | if (!(no_residue[mapping->magnitude[i]] & no_residue[mapping->angle[i]])) { |
| 1636 | no_residue[mapping->magnitude[i]] = 0; |
| 1637 | no_residue[mapping->angle[i]] = 0; |
| 1638 | } |
| 1639 | } |
| 1640 | |
| 1641 | // Decode residue |
| 1642 | |
| 1643 | for (i = 0; i < mapping->submaps; ++i) { |
| 1644 | vorbis_residue *residue; |
| 1645 | unsigned ch = 0; |
| 1646 | int ret; |
| 1647 | |
| 1648 | for (j = 0; j < vc->audio_channels; ++j) { |
| 1649 | if ((mapping->submaps == 1) || (i == mapping->mux[j])) { |
| 1650 | res_chan[j] = res_num; |
| 1651 | if (no_residue[j]) { |
| 1652 | do_not_decode[ch] = 1; |
| 1653 | } else { |
| 1654 | do_not_decode[ch] = 0; |
| 1655 | } |
| 1656 | ++ch; |
| 1657 | ++res_num; |
| 1658 | } |
| 1659 | } |
| 1660 | residue = &vc->residues[mapping->submap_residue[i]]; |
| 1661 | if (ch_left < ch) { |
| 1662 | av_log(vc->avctx, AV_LOG_ERROR, "Too many channels in vorbis_floor_decode.\n"); |
| 1663 | return AVERROR_INVALIDDATA; |
| 1664 | } |
| 1665 | if (ch) { |
| 1666 | ret = vorbis_residue_decode(vc, residue, ch, do_not_decode, ch_res_ptr, vlen, ch_left); |
| 1667 | if (ret < 0) |
| 1668 | return ret; |
| 1669 | } |
| 1670 | |
| 1671 | ch_res_ptr += ch * vlen; |
| 1672 | ch_left -= ch; |
| 1673 | } |
| 1674 | |
| 1675 | if (ch_left > 0) |
| 1676 | return AVERROR_INVALIDDATA; |
| 1677 | |
| 1678 | // Inverse coupling |
| 1679 | |
| 1680 | for (i = mapping->coupling_steps - 1; i >= 0; --i) { //warning: i has to be signed |
| 1681 | float *mag, *ang; |
| 1682 | |
| 1683 | mag = vc->channel_residues+res_chan[mapping->magnitude[i]] * blocksize / 2; |
| 1684 | ang = vc->channel_residues+res_chan[mapping->angle[i]] * blocksize / 2; |
| 1685 | vc->dsp.vorbis_inverse_coupling(mag, ang, blocksize / 2); |
| 1686 | } |
| 1687 | |
| 1688 | // Dotproduct, MDCT |
| 1689 | |
| 1690 | mdct = &vc->mdct[blockflag]; |
| 1691 | |
| 1692 | for (j = vc->audio_channels-1;j >= 0; j--) { |
| 1693 | ch_res_ptr = vc->channel_residues + res_chan[j] * blocksize / 2; |
| 1694 | vc->fdsp->vector_fmul(floor_ptr[j], floor_ptr[j], ch_res_ptr, blocksize / 2); |
| 1695 | mdct->imdct_half(mdct, ch_res_ptr, floor_ptr[j]); |
| 1696 | } |
| 1697 | |
| 1698 | // Overlap/add, save data for next overlapping |
| 1699 | |
| 1700 | retlen = (blocksize + vc->blocksize[previous_window]) / 4; |
| 1701 | for (j = 0; j < vc->audio_channels; j++) { |
| 1702 | unsigned bs0 = vc->blocksize[0]; |
| 1703 | unsigned bs1 = vc->blocksize[1]; |
| 1704 | float *residue = vc->channel_residues + res_chan[j] * blocksize / 2; |
| 1705 | float *saved = vc->saved + j * bs1 / 4; |
| 1706 | float *ret = floor_ptr[j]; |
| 1707 | float *buf = residue; |
| 1708 | const float *win = vc->win[blockflag & previous_window]; |
| 1709 | |
| 1710 | if (blockflag == previous_window) { |
| 1711 | vc->fdsp->vector_fmul_window(ret, saved, buf, win, blocksize / 4); |
| 1712 | } else if (blockflag > previous_window) { |
| 1713 | vc->fdsp->vector_fmul_window(ret, saved, buf, win, bs0 / 4); |
| 1714 | memcpy(ret+bs0/2, buf+bs0/4, ((bs1-bs0)/4) * sizeof(float)); |
| 1715 | } else { |
| 1716 | memcpy(ret, saved, ((bs1 - bs0) / 4) * sizeof(float)); |
| 1717 | vc->fdsp->vector_fmul_window(ret + (bs1 - bs0) / 4, saved + (bs1 - bs0) / 4, buf, win, bs0 / 4); |
| 1718 | } |
| 1719 | memcpy(saved, buf + blocksize / 4, blocksize / 4 * sizeof(float)); |
| 1720 | } |
| 1721 | |
| 1722 | vc->previous_window = blockflag; |
| 1723 | return retlen; |
| 1724 | } |
| 1725 | |
| 1726 | // Return the decoded audio packet through the standard api |
| 1727 | |
| 1728 | static int vorbis_decode_frame(AVCodecContext *avctx, void *data, |
| 1729 | int *got_frame_ptr, AVPacket *avpkt) |
| 1730 | { |
| 1731 | const uint8_t *buf = avpkt->data; |
| 1732 | int buf_size = avpkt->size; |
| 1733 | vorbis_context *vc = avctx->priv_data; |
| 1734 | AVFrame *frame = data; |
| 1735 | GetBitContext *gb = &vc->gb; |
| 1736 | float *channel_ptrs[255]; |
| 1737 | int i, len, ret; |
| 1738 | |
| 1739 | av_dlog(NULL, "packet length %d \n", buf_size); |
| 1740 | |
| 1741 | if (*buf == 1 && buf_size > 7) { |
| 1742 | init_get_bits(gb, buf+1, buf_size*8 - 8); |
| 1743 | vorbis_free(vc); |
| 1744 | if ((ret = vorbis_parse_id_hdr(vc))) { |
| 1745 | av_log(avctx, AV_LOG_ERROR, "Id header corrupt.\n"); |
| 1746 | vorbis_free(vc); |
| 1747 | return ret; |
| 1748 | } |
| 1749 | |
| 1750 | if (vc->audio_channels > 8) |
| 1751 | avctx->channel_layout = 0; |
| 1752 | else |
| 1753 | avctx->channel_layout = ff_vorbis_channel_layouts[vc->audio_channels - 1]; |
| 1754 | |
| 1755 | avctx->channels = vc->audio_channels; |
| 1756 | avctx->sample_rate = vc->audio_samplerate; |
| 1757 | return buf_size; |
| 1758 | } |
| 1759 | |
| 1760 | if (*buf == 3 && buf_size > 7) { |
| 1761 | av_log(avctx, AV_LOG_DEBUG, "Ignoring comment header\n"); |
| 1762 | return buf_size; |
| 1763 | } |
| 1764 | |
| 1765 | if (*buf == 5 && buf_size > 7 && vc->channel_residues && !vc->modes) { |
| 1766 | init_get_bits(gb, buf+1, buf_size*8 - 8); |
| 1767 | if ((ret = vorbis_parse_setup_hdr(vc))) { |
| 1768 | av_log(avctx, AV_LOG_ERROR, "Setup header corrupt.\n"); |
| 1769 | vorbis_free(vc); |
| 1770 | return ret; |
| 1771 | } |
| 1772 | return buf_size; |
| 1773 | } |
| 1774 | |
| 1775 | if (!vc->channel_residues || !vc->modes) { |
| 1776 | av_log(avctx, AV_LOG_ERROR, "Data packet before valid headers\n"); |
| 1777 | return AVERROR_INVALIDDATA; |
| 1778 | } |
| 1779 | |
| 1780 | /* get output buffer */ |
| 1781 | frame->nb_samples = vc->blocksize[1] / 2; |
| 1782 | if ((ret = ff_get_buffer(avctx, frame, 0)) < 0) |
| 1783 | return ret; |
| 1784 | |
| 1785 | if (vc->audio_channels > 8) { |
| 1786 | for (i = 0; i < vc->audio_channels; i++) |
| 1787 | channel_ptrs[i] = (float *)frame->extended_data[i]; |
| 1788 | } else { |
| 1789 | for (i = 0; i < vc->audio_channels; i++) { |
| 1790 | int ch = ff_vorbis_channel_layout_offsets[vc->audio_channels - 1][i]; |
| 1791 | channel_ptrs[ch] = (float *)frame->extended_data[i]; |
| 1792 | } |
| 1793 | } |
| 1794 | |
| 1795 | init_get_bits(gb, buf, buf_size*8); |
| 1796 | |
| 1797 | if ((len = vorbis_parse_audio_packet(vc, channel_ptrs)) <= 0) |
| 1798 | return len; |
| 1799 | |
| 1800 | if (!vc->first_frame) { |
| 1801 | vc->first_frame = 1; |
| 1802 | *got_frame_ptr = 0; |
| 1803 | av_frame_unref(frame); |
| 1804 | return buf_size; |
| 1805 | } |
| 1806 | |
| 1807 | av_dlog(NULL, "parsed %d bytes %d bits, returned %d samples (*ch*bits) \n", |
| 1808 | get_bits_count(gb) / 8, get_bits_count(gb) % 8, len); |
| 1809 | |
| 1810 | frame->nb_samples = len; |
| 1811 | *got_frame_ptr = 1; |
| 1812 | |
| 1813 | return buf_size; |
| 1814 | } |
| 1815 | |
| 1816 | // Close decoder |
| 1817 | |
| 1818 | static av_cold int vorbis_decode_close(AVCodecContext *avctx) |
| 1819 | { |
| 1820 | vorbis_context *vc = avctx->priv_data; |
| 1821 | |
| 1822 | vorbis_free(vc); |
| 1823 | |
| 1824 | return 0; |
| 1825 | } |
| 1826 | |
| 1827 | static av_cold void vorbis_decode_flush(AVCodecContext *avctx) |
| 1828 | { |
| 1829 | vorbis_context *vc = avctx->priv_data; |
| 1830 | |
| 1831 | if (vc->saved) { |
| 1832 | memset(vc->saved, 0, (vc->blocksize[1] / 4) * vc->audio_channels * |
| 1833 | sizeof(*vc->saved)); |
| 1834 | } |
| 1835 | vc->previous_window = -1; |
| 1836 | vc->first_frame = 0; |
| 1837 | } |
| 1838 | |
| 1839 | AVCodec ff_vorbis_decoder = { |
| 1840 | .name = "vorbis", |
| 1841 | .long_name = NULL_IF_CONFIG_SMALL("Vorbis"), |
| 1842 | .type = AVMEDIA_TYPE_AUDIO, |
| 1843 | .id = AV_CODEC_ID_VORBIS, |
| 1844 | .priv_data_size = sizeof(vorbis_context), |
| 1845 | .init = vorbis_decode_init, |
| 1846 | .close = vorbis_decode_close, |
| 1847 | .decode = vorbis_decode_frame, |
| 1848 | .flush = vorbis_decode_flush, |
| 1849 | .capabilities = CODEC_CAP_DR1, |
| 1850 | .channel_layouts = ff_vorbis_channel_layouts, |
| 1851 | .sample_fmts = (const enum AVSampleFormat[]) { AV_SAMPLE_FMT_FLTP, |
| 1852 | AV_SAMPLE_FMT_NONE }, |
| 1853 | }; |