32 #define BITSTREAM_READER_LE 
   47 #define V_MAX_VLCS (1 << 16) 
   48 #define V_MAX_PARTITIONS (1 << 20) 
   62 struct vorbis_context_s;
 
   88             int16_t       subclass_books[16][8];
 
  103     int16_t       books[64][8];
 
  126 typedef struct vorbis_context_s {
 
  142     uint32_t      blocksize[2];
 
  163     (13.1f * atan(0.00074f * (x)) + 2.24f * atan(1.85e-8f * (x) * (x)) + 1e-4f * (x)) 
  165 static const char idx_err_str[] = 
"Index value %d out of range (0 - %d) for %s at %s:%i\n";
 
  166 #define VALIDATE_INDEX(idx, limit) \ 
  168         av_log(vc->avccontext, AV_LOG_ERROR,\ 
  170                (int)(idx), (int)(limit - 1), #idx, __FILE__, __LINE__);\ 
  171         return AVERROR_INVALIDDATA;\ 
  173 #define GET_VALIDATED_INDEX(idx, bits, limit) \ 
  175         idx = get_bits(gb, bits);\ 
  176         VALIDATE_INDEX(idx, limit)\ 
  181     double mant = val & 0x1fffff;
 
  182     long exp    = (val & 0x7fe00000
L) >> 21;
 
  183     if (val & 0x80000000)
 
  185     return ldexp(mant, exp - 20 - 768);
 
  244     uint32_t *tmp_vlc_codes;
 
  246     uint16_t *codebook_multiplicands;
 
  260         unsigned ordered, 
t, entries, used_entries = 0;
 
  266                    " %u. Codebook setup data corrupt.\n", cb);
 
  274                    " %u. Codebook's dimension is invalid (%d).\n",
 
  282                    " %u. Codebook has too many entries (%u).\n",
 
  290         av_dlog(
NULL, 
" codebook_dimensions %d, codebook_entries %u\n",
 
  303                 for (ce = 0; ce < entries; ++ce) {
 
  306                         tmp_vlc_bits[ce] = 
get_bits(gb, 5) + 1;
 
  309                         tmp_vlc_bits[ce] = 0;
 
  314                 used_entries = entries;
 
  315                 for (ce = 0; ce < entries; ++ce)
 
  316                     tmp_vlc_bits[ce] = 
get_bits(gb, 5) + 1;
 
  319             unsigned current_entry  = 0;
 
  320             unsigned current_length = 
get_bits(gb, 5) + 1;
 
  322             av_dlog(
NULL, 
" ordered, current length: %u\n", current_length);  
 
  324             used_entries = entries;
 
  325             for (; current_entry < used_entries && current_length <= 32; ++current_length) {
 
  334                 for (i = current_entry; i < number+current_entry; ++i)
 
  335                     if (i < used_entries)
 
  336                         tmp_vlc_bits[i] = current_length;
 
  338                 current_entry+=number;
 
  340             if (current_entry>used_entries) {
 
  360             unsigned codebook_value_bits = 
get_bits(gb, 4) + 1;
 
  361             unsigned codebook_sequence_p = 
get_bits1(gb);
 
  363             av_dlog(
NULL, 
" We expect %d numbers for building the codevectors. \n",
 
  364                     codebook_lookup_values);
 
  366                     codebook_delta_value, codebook_minimum_value);
 
  368             for (i = 0; i < codebook_lookup_values; ++i) {
 
  369                 codebook_multiplicands[i] = 
get_bits(gb, codebook_value_bits);
 
  371                 av_dlog(
NULL, 
" multiplicands*delta+minmum : %e \n",
 
  372                         (
float)codebook_multiplicands[i] * codebook_delta_value + codebook_minimum_value);
 
  373                 av_dlog(
NULL, 
" multiplicand %u\n", codebook_multiplicands[i]);
 
  381             for (j = 0, i = 0; i < entries; ++i) {
 
  384                 if (tmp_vlc_bits[i]) {
 
  386                     unsigned lookup_offset = i;
 
  390                     for (k = 0; k < 
dim; ++k) {
 
  391                         unsigned multiplicand_offset = lookup_offset % codebook_lookup_values;
 
  392                         codebook_setup->
codevectors[j * dim + k] = codebook_multiplicands[multiplicand_offset] * codebook_delta_value + codebook_minimum_value + last;
 
  393                         if (codebook_sequence_p)
 
  395                         lookup_offset/=codebook_lookup_values;
 
  397                     tmp_vlc_bits[j] = tmp_vlc_bits[i];
 
  400                     for (k = 0; k < 
dim; ++k)
 
  408             if (j != used_entries) {
 
  413             entries = used_entries;
 
  427         for (t = 0; t < entries; ++
t)
 
  428             if (tmp_vlc_bits[t] >= codebook_setup->
maxdepth)
 
  429                 codebook_setup->
maxdepth = tmp_vlc_bits[
t];
 
  439                             entries, tmp_vlc_bits, 
sizeof(*tmp_vlc_bits),
 
  440                             sizeof(*tmp_vlc_bits), tmp_vlc_codes,
 
  441                             sizeof(*tmp_vlc_codes), 
sizeof(*tmp_vlc_codes),
 
  450     av_free(codebook_multiplicands);
 
  457     av_free(codebook_multiplicands);
 
  466     unsigned i, vorbis_time_count = 
get_bits(gb, 6) + 1;
 
  468     for (i = 0; i < vorbis_time_count; ++i) {
 
  469         unsigned vorbis_tdtransform = 
get_bits(gb, 16);
 
  471         av_dlog(
NULL, 
" Vorbis time domain transform %u: %u\n",
 
  472                 vorbis_time_count, vorbis_tdtransform);
 
  474         if (vorbis_tdtransform) {
 
  506             int maximum_class = -1;
 
  507             unsigned rangebits, rangemax, floor1_values = 2;
 
  521                 av_dlog(
NULL, 
" %d. floor %d partition class %d \n",
 
  526             av_dlog(
NULL, 
" maximum class %d \n", maximum_class);
 
  528             for (j = 0; j <= maximum_class; ++j) {
 
  532                 av_dlog(
NULL, 
" %d floor %d class dim: %d subclasses %d \n", i, j,
 
  563             rangemax = (1 << rangebits);
 
  566                        "Floor value is too large for blocksize: %u (%"PRIu32
")\n",
 
  577                     av_dlog(
NULL, 
" %u. floor1 Y coord. %d\n", floor1_values,
 
  589             unsigned max_codebook_dim = 0;
 
  596                        "Floor 0 order is 0.\n");
 
  602                        "Floor 0 rate is 0.\n");
 
  608                        "Floor 0 bark map size is 0.\n");
 
  685         unsigned high_bits, low_bits;
 
  695         if (res_setup->
begin>res_setup->
end ||
 
  699                    "partition out of bounds: type, begin, end, size, blocksize: %"PRIu16
", %"PRIu32
", %"PRIu32
", %u, %"PRIu32
"\n",
 
  716         av_dlog(
NULL, 
"    begin %d end %d part.size %d classif.s %d classbook %d \n",
 
  725             cascade[j] = (high_bits << 3) + low_bits;
 
  732             for (k = 0; k < 8; ++k) {
 
  733                 if (cascade[j]&(1 << k)) {
 
  736                     av_dlog(
NULL, 
"     %u class cascade depth %u book: %d\n",
 
  737                             j, k, res_setup->
books[j][k]);
 
  742                     res_setup->
books[j][k] = -1;
 
  780                                                        sizeof(*mapping_setup->
angle));
 
  797         if (mapping_setup->
submaps>1) {
 
  799                                             sizeof(*mapping_setup->
mux));
 
  804         for (j = 0; j < mapping_setup->
submaps; ++j) {
 
  809             av_dlog(
NULL, 
"   %u mapping %u submap : floor %d, residue %d\n", i, j,
 
  827     for (blockflag = 0; blockflag < 2; ++blockflag) {
 
  829         floors[floor_number].
data.
t0.
map[blockflag] =
 
  832         map =  floors[floor_number].
data.
t0.
map[blockflag];
 
  833         vf  = &floors[floor_number].
data.
t0;
 
  835         for (idx = 0; idx < n; ++idx) {
 
  836             map[idx] = floor(
BARK((vf->rate * idx) / (2.0f * n)) *
 
  837                              (vf->bark_map_size / 
BARK(vf->rate / 2.0f)));
 
  838             if (vf->bark_map_size-1 < map[idx])
 
  839                 map[idx] = vf->bark_map_size - 1;
 
  842         vf->map_size[blockflag] = n;
 
  845     for (idx = 0; idx <= n; ++idx) {
 
  846         av_dlog(
NULL, 
"floor0 map: map at pos %d is %d\n", idx, map[idx]);
 
  868         av_dlog(
NULL, 
" %u mode: blockflag %d, windowtype %d, transformtype %d, mapping %d\n",
 
  951     if (bl0 > 13 || bl0 < 6 || bl1 > 13 || bl1 < 6 || bl1 < bl0) {
 
  972     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 ",
 
 1016     if (hdr_type != 1) {
 
 1028     if (hdr_type != 5) {
 
 1061     float *lsp = vf->lsp;
 
 1062     unsigned amplitude, book_idx;
 
 1065     if (!vf->amplitude_bits)
 
 1068     amplitude = 
get_bits(&vc->
gb, vf->amplitude_bits);
 
 1069     if (amplitude > 0) {
 
 1071         unsigned idx, lsp_len = 0;
 
 1075         if (book_idx >= vf->num_books) {
 
 1077                     "floor0 dec: booknumber too high!\n");
 
 1080         av_dlog(
NULL, 
"floor0 dec: booknumber: %u\n", book_idx);
 
 1081         codebook = vc->
codebooks[vf->book_list[book_idx]];
 
 1086         while (lsp_len<vf->order) {
 
 1095             av_dlog(
NULL, 
"floor0 dec: vector offset: %d\n", vec_off);
 
 1097             for (idx = 0; idx < codebook.
dimensions; ++idx)
 
 1098                 lsp[lsp_len+idx] = codebook.
codevectors[vec_off+idx] + last;
 
 1099             last = lsp[lsp_len+idx-1]; 
 
 1106             for (idx = 0; idx < lsp_len; ++idx)
 
 1107                 av_dlog(
NULL, 
"floor0 dec: coeff at %d is %f\n", idx, lsp[idx]);
 
 1113             int order = vf->order;
 
 1114             float wstep = 
M_PI / vf->bark_map_size;
 
 1116             for (i = 0; i < order; i++)
 
 1117                 lsp[i] = 2.0f * cos(lsp[i]);
 
 1119             av_dlog(
NULL, 
"floor0 synth: map_size = %"PRIu32
"; m = %d; wstep = %f\n",
 
 1120                     vf->map_size[blockflag], order, wstep);
 
 1123             while (i < vf->map_size[blockflag]) {
 
 1124                 int j, iter_cond = vf->map[blockflag][i];
 
 1127                 float two_cos_w = 2.0f * cos(wstep * iter_cond); 
 
 1130                 for (j = 0; j + 1 < order; j += 2) {
 
 1131                     q *= lsp[j]     - two_cos_w;
 
 1132                     p *= lsp[j + 1] - two_cos_w;
 
 1135                     p *= p * (2.0f - two_cos_w);
 
 1136                     q *= q * (2.0f + two_cos_w);
 
 1138                     q *= two_cos_w-lsp[j]; 
 
 1141                     p *= p * (4.f - two_cos_w * two_cos_w);
 
 1146                 q = exp((((amplitude*vf->amplitude_offset) /
 
 1147                           (((1 << vf->amplitude_bits) - 1) * sqrt(p + q)))
 
 1148                          - vf->amplitude_offset) * .11512925f);
 
 1153                 } 
while (vf->map[blockflag][i] == iter_cond);
 
 1171     uint16_t range_v[4] = { 256, 128, 86, 64 };
 
 1172     unsigned range = range_v[vf->multiplier - 1];
 
 1173     uint16_t floor1_Y[258];
 
 1174     uint16_t floor1_Y_final[258];
 
 1175     int floor1_flag[258];
 
 1176     unsigned partition_class, cdim, cbits, csub, cval, 
offset, i, j;
 
 1177     int book, adx, ady, dy, 
off, predicted, err;
 
 1188     av_dlog(
NULL, 
"floor 0 Y %d floor 1 Y %d \n", floor1_Y[0], floor1_Y[1]);
 
 1191     for (i = 0; i < vf->partitions; ++i) {
 
 1192         partition_class = vf->partition_class[i];
 
 1193         cdim   = vf->class_dimensions[partition_class];
 
 1194         cbits  = vf->class_subclasses[partition_class];
 
 1195         csub = (1 << cbits) - 1;
 
 1204         for (j = 0; j < cdim; ++j) {
 
 1205             book = vf->subclass_books[partition_class][cval & csub];
 
 1207             av_dlog(
NULL, 
"book %d Cbits %u cval %u  bits:%d\n",
 
 1210             cval = cval >> cbits;
 
 1215                 floor1_Y[offset+j] = 0;
 
 1219                     vf->list[offset+j].x, floor1_Y[offset+j]);
 
 1228     floor1_Y_final[0] = floor1_Y[0];
 
 1229     floor1_Y_final[1] = floor1_Y[1];
 
 1231     for (i = 2; i < vf->x_list_dim; ++i) {
 
 1232         unsigned val, highroom, lowroom, room, high_neigh_offs, low_neigh_offs;
 
 1234         low_neigh_offs  = vf->list[i].low;
 
 1235         high_neigh_offs = vf->list[i].high;
 
 1236         dy  = floor1_Y_final[high_neigh_offs] - floor1_Y_final[low_neigh_offs];  
 
 1237         adx = vf->list[high_neigh_offs].x - vf->list[low_neigh_offs].x;
 
 1239         err = ady * (vf->list[i].x - vf->list[low_neigh_offs].x);
 
 1242             predicted = floor1_Y_final[low_neigh_offs] - 
off;
 
 1244             predicted = floor1_Y_final[low_neigh_offs] + 
off;
 
 1248         highroom = range-predicted;
 
 1249         lowroom  = predicted;
 
 1250         if (highroom < lowroom) {
 
 1251             room = highroom * 2;
 
 1256             floor1_flag[low_neigh_offs]  = 1;
 
 1257             floor1_flag[high_neigh_offs] = 1;
 
 1260                 if (highroom > lowroom) {
 
 1261                     floor1_Y_final[i] = av_clip_uint16(val - lowroom + predicted);
 
 1263                     floor1_Y_final[i] = av_clip_uint16(predicted - val + highroom - 1);
 
 1267                     floor1_Y_final[i] = av_clip_uint16(predicted - (val + 1) / 2);
 
 1269                     floor1_Y_final[i] = av_clip_uint16(predicted + val / 2);
 
 1274             floor1_Y_final[i] = av_clip_uint16(predicted);
 
 1277         av_dlog(
NULL, 
" Decoded floor(%d) = %u / val %u\n",
 
 1278                 vf->list[i].x, floor1_Y_final[i], val);
 
 1305     unsigned pass, ch_used, i, j, k, l;
 
 1306     unsigned max_output = (ch - 1) * vlen;
 
 1309         for (j = 1; j < ch; ++j)
 
 1310             do_not_decode[0] &= do_not_decode[j];  
 
 1311         if (do_not_decode[0])
 
 1314         max_output += vr->
end / ch;
 
 1317         max_output += vr->
end;
 
 1320     if (max_output > ch_left * vlen) {
 
 1325     av_dlog(
NULL, 
" residue type 0/1/2 decode begin, ch: %d  cpc %d  \n", ch, c_p_c);
 
 1328         uint16_t voffset, partition_count, j_times_ptns_to_read;
 
 1330         voffset = vr->
begin;
 
 1331         for (partition_count = 0; partition_count < ptns_to_read;) {  
 
 1334                 for (j_times_ptns_to_read = 0, j = 0; j < ch_used; ++j) {
 
 1335                     if (!do_not_decode[j]) {
 
 1342                         for (i = 0; i < c_p_c; ++i) {
 
 1345                             temp2 = (((uint64_t)temp) * inverse_class) >> 32;
 
 1346                             if (partition_count + c_p_c - 1 - i < ptns_to_read)
 
 1347                                 classifs[j_times_ptns_to_read + partition_count + c_p_c - 1 - i] = temp - temp2 * vr->
classifications;
 
 1351                     j_times_ptns_to_read += ptns_to_read;
 
 1354             for (i = 0; (i < c_p_c) && (partition_count < ptns_to_read); ++i) {
 
 1355                 for (j_times_ptns_to_read = 0, j = 0; j < ch_used; ++j) {
 
 1358                     if (!do_not_decode[j]) {
 
 1359                         unsigned vqclass = classifs[j_times_ptns_to_read + partition_count];
 
 1370                                 voffs = voffset+j*vlen;
 
 1371                                 for (k = 0; k < step; ++k) {
 
 1373                                     for (l = 0; l < 
dim; ++l)
 
 1374                                         vec[voffs + k + l * step] += codebook.
codevectors[coffs + l];
 
 1376                             } 
else if (vr_type == 1) {
 
 1377                                 voffs = voffset + j * vlen;
 
 1378                                 for (k = 0; k < step; ++k) {
 
 1380                                     for (l = 0; l < 
dim; ++l, ++voffs) {
 
 1383                                         av_dlog(
NULL, 
" pass %d offs: %d curr: %f change: %f cv offs.: %d  \n",
 
 1384                                                 pass, voffs, vec[voffs], codebook.
codevectors[coffs+l], coffs);
 
 1387                             } 
else if (vr_type == 2 && ch == 2 && (voffset & 1) == 0 && (dim & 1) == 0) { 
 
 1388                                 voffs = voffset >> 1;
 
 1391                                     for (k = 0; k < step; ++k) {
 
 1394                                         vec[voffs + k + vlen] += codebook.
codevectors[coffs + 1];
 
 1396                                 } 
else if (dim == 4) {
 
 1397                                     for (k = 0; k < step; ++k, voffs += 2) {
 
 1400                                         vec[voffs + 1       ] += codebook.
codevectors[coffs + 2];
 
 1401                                         vec[voffs + vlen    ] += codebook.
codevectors[coffs + 1];
 
 1402                                         vec[voffs + vlen + 1] += codebook.
codevectors[coffs + 3];
 
 1405                                 for (k = 0; k < step; ++k) {
 
 1407                                     for (l = 0; l < 
dim; l += 2, voffs++) {
 
 1409                                         vec[voffs + vlen] += codebook.
codevectors[coffs + l + 1];
 
 1411                                         av_dlog(
NULL, 
" pass %d offs: %d curr: %f change: %f cv offs.: %d+%d  \n",
 
 1412                                                 pass, voffset / ch + (voffs % ch) * vlen,
 
 1413                                                 vec[voffset / ch + (voffs % ch) * vlen],
 
 1418                             } 
else if (vr_type == 2) {
 
 1419                                 unsigned voffs_div = 
FASTDIV(voffset << 1, ch <<1);
 
 1420                                 unsigned voffs_mod = voffset - voffs_div * ch;
 
 1422                                 for (k = 0; k < step; ++k) {
 
 1424                                     for (l = 0; l < 
dim; ++l) {
 
 1425                                         vec[voffs_div + voffs_mod * vlen] +=
 
 1428                                         av_dlog(
NULL, 
" pass %d offs: %d curr: %f change: %f cv offs.: %d+%d  \n",
 
 1429                                                 pass, voffs_div + voffs_mod * vlen,
 
 1430                                                 vec[voffs_div + voffs_mod * vlen],
 
 1433                                         if (++voffs_mod == ch) {
 
 1442                     j_times_ptns_to_read += ptns_to_read;
 
 1455                                         float *vec, 
unsigned vlen,
 
 1460     else if (vr->
type == 1)
 
 1462     else if (vr->
type == 0)
 
 1473     for (i = 0;  i < blocksize;  i++) {
 
 1476                 ang[i] = mag[i] - ang[i];
 
 1478                 float temp = ang[i];
 
 1486                 float temp = ang[i];
 
 1501     unsigned mode_number, blockflag, blocksize;
 
 1508     unsigned res_num = 0;
 
 1526     av_dlog(
NULL, 
" Mode number: %u , mapping: %d , blocktype %d\n", mode_number,
 
 1531     vlen = blocksize / 2;
 
 1537     memset(ch_res_ptr,   0, 
sizeof(
float) * vc->
audio_channels * vlen); 
 
 1539         memset(floor_ptr[i], 0, vlen * 
sizeof(floor_ptr[0][0])); 
 
 1552         ret = floor->
decode(vc, &floor->
data, floor_ptr[i]);
 
 1558         no_residue[i] = ret;
 
 1564         if (!(no_residue[mapping->
magnitude[i]] & no_residue[mapping->
angle[i]])) {
 
 1566             no_residue[mapping->
angle[i]]     = 0;
 
 1572     for (i = 0; i < mapping->
submaps; ++i) {
 
 1578             if ((mapping->
submaps == 1) || (i == mapping->
mux[j])) {
 
 1579                 res_chan[j] = res_num;
 
 1580                 if (no_residue[j]) {
 
 1581                     do_not_decode[ch] = 1;
 
 1583                     do_not_decode[ch] = 0;
 
 1600         ch_res_ptr += ch * vlen;
 
 1619     mdct = &vc->
mdct[blockflag];
 
 1623         vc->
fdsp.
vector_fmul(floor_ptr[j], floor_ptr[j], ch_res_ptr, blocksize / 2);
 
 1624         mdct->
imdct_half(mdct, ch_res_ptr, floor_ptr[j]);
 
 1629     retlen = (blocksize + vc->
blocksize[previous_window]) / 4;
 
 1634         float *saved      = vc->
saved + j * bs1 / 4;
 
 1635         float *ret        = floor_ptr[j];
 
 1636         float *buf        = residue;
 
 1637         const float *win  = vc->
win[blockflag & previous_window];
 
 1639         if (blockflag == previous_window) {
 
 1641         } 
else if (blockflag > previous_window) {
 
 1643             memcpy(ret+bs0/2, buf+bs0/4, ((bs1-bs0)/4) * 
sizeof(
float));
 
 1645             memcpy(ret, saved, ((bs1 - bs0) / 4) * 
sizeof(
float));
 
 1648         memcpy(saved, buf + blocksize / 4, blocksize / 4 * 
sizeof(
float));
 
 1658                                int *got_frame_ptr, 
AVPacket *avpkt)
 
 1661     int buf_size       = avpkt->
size;
 
 1664     float *channel_ptrs[255];
 
 1669     if (*buf == 1 && buf_size > 7) {
 
 1688     if (*buf == 3 && buf_size > 7) {
 
 1736     av_dlog(
NULL, 
"parsed %d bytes %d bits, returned %d samples (*ch*bits) \n",
 
 1763                              sizeof(*vc->
saved));