115 static const char overread_err[] =
"Input buffer exhausted before END element found\n";
120 for (i = 0; i < tags; i++) {
121 int syn_ele =
layout[i][0];
123 sum += (1 + (syn_ele ==
TYPE_CPE)) *
143 int type,
int id,
int *channels)
148 if (!ac->
che[type][
id]) {
161 if (ac->
che[type][
id])
171 int type,
id, ch, ret;
174 for (type = 0; type < 4; type++) {
192 for (ch = 0; ch < avctx->
channels; ch++) {
208 uint8_t (*layout_map)[3],
int offset, uint64_t left,
209 uint64_t right,
int pos)
211 if (layout_map[offset][0] ==
TYPE_CPE) {
213 .av_position = left | right,
215 .elem_id = layout_map[offset][1],
223 .elem_id = layout_map[offset][1],
227 .av_position = right,
229 .elem_id = layout_map[offset + 1][1],
239 int num_pos_channels = 0;
243 for (i = *current; i < tags; i++) {
244 if (layout_map[i][2] != pos)
254 num_pos_channels += 2;
265 return num_pos_channels;
270 int i, n, total_non_cc_elements;
272 int num_front_channels, num_side_channels, num_back_channels;
281 if (num_front_channels < 0)
285 if (num_side_channels < 0)
289 if (num_back_channels < 0)
293 if (num_front_channels & 1) {
297 .elem_id = layout_map[i][1],
301 num_front_channels--;
303 if (num_front_channels >= 4) {
308 num_front_channels -= 2;
310 if (num_front_channels >= 2) {
315 num_front_channels -= 2;
317 while (num_front_channels >= 2) {
322 num_front_channels -= 2;
325 if (num_side_channels >= 2) {
330 num_side_channels -= 2;
332 while (num_side_channels >= 2) {
337 num_side_channels -= 2;
340 while (num_back_channels >= 4) {
345 num_back_channels -= 2;
347 if (num_back_channels >= 2) {
352 num_back_channels -= 2;
354 if (num_back_channels) {
358 .elem_id = layout_map[i][1],
369 .elem_id = layout_map[i][1],
378 .elem_id = layout_map[i][1],
385 total_non_cc_elements = n = i;
388 for (i = 1; i < n; i++)
389 if (e2c_vec[i - 1].av_position > e2c_vec[i].av_position) {
397 for (i = 0; i < total_non_cc_elements; i++) {
398 layout_map[i][0] = e2c_vec[i].
syn_ele;
399 layout_map[i][1] = e2c_vec[i].
elem_id;
401 if (e2c_vec[i].av_position != UINT64_MAX) {
414 ac->
oc[0] = ac->
oc[1];
425 ac->
oc[1] = ac->
oc[0];
438 uint8_t layout_map[MAX_ELEM_ID * 4][3],
int tags,
439 enum OCStatus oc_type,
int get_new_frame)
442 int i, channels = 0, ret;
446 memcpy(ac->
oc[1].
layout_map, layout_map, tags *
sizeof(layout_map[0]));
454 for (i = 0; i < tags; i++) {
455 int type = layout_map[i][0];
456 int id = layout_map[i][1];
457 int position = layout_map[i][2];
464 if (ac->
oc[1].
m4ac.
ps == 1 && channels == 2) {
496 if (channel_config < 1 || channel_config > 7) {
498 "invalid default channel configuration (%d)\n",
504 *tags *
sizeof(*layout_map));
518 uint8_t layout_map[MAX_ELEM_ID*4][3];
523 &layout_map_tags, 2) < 0)
535 uint8_t layout_map[MAX_ELEM_ID * 4][3];
540 &layout_map_tags, 1) < 0)
626 layout_map[0][0] = syn_ele;
628 layout_map[0][2] = type;
642 int num_front, num_side, num_back, num_lfe, num_assoc_data, num_cc;
652 "Sample rate index in program config element does not "
653 "match the sample rate index configured by the container.\n");
709 int extension_flag, ret;
710 uint8_t layout_map[MAX_ELEM_ID*4][3];
726 if (channel_config == 0) {
728 tags =
decode_pce(avctx, m4ac, layout_map, gb);
733 &tags, channel_config)))
739 }
else if (m4ac->
sbr == 1 && m4ac->
ps == -1)
745 if (extension_flag) {
796 sync_extension)) < 0)
800 "invalid sampling rate index %d\n",
817 "Audio object type %s%d is not supported.\n",
818 m4ac->
sbr == 1 ?
"SBR+" :
"",
824 "AOT %d chan config %d sampling index %d (%d) SBR %d PS %d\n",
841 union {
unsigned u;
int s; } v = { previous_val * 1664525u + 1013904223 };
864 if (92017 <= rate)
return 0;
865 else if (75132 <= rate)
return 1;
866 else if (55426 <= rate)
return 2;
867 else if (46009 <= rate)
return 3;
868 else if (37566 <= rate)
return 4;
869 else if (27713 <= rate)
return 5;
870 else if (23004 <= rate)
return 6;
871 else if (18783 <= rate)
return 7;
872 else if (13856 <= rate)
return 8;
873 else if (11502 <= rate)
return 9;
874 else if (9391 <= rate)
return 10;
885 #define AAC_INIT_VLC_STATIC(num, size) \
886 INIT_VLC_STATIC(&vlc_spectral[num], 8, ff_aac_spectral_sizes[num], \
887 ff_aac_spectral_bits[num], sizeof(ff_aac_spectral_bits[num][0]), \
888 sizeof(ff_aac_spectral_bits[num][0]), \
889 ff_aac_spectral_codes[num], sizeof(ff_aac_spectral_codes[num][0]), \
890 sizeof(ff_aac_spectral_codes[num][0]), \
911 uint8_t layout_map[MAX_ELEM_ID*4][3];
1017 "Invalid Predictor Reset Group.\n");
1060 for (i = 0; i < 7; i++) {
1088 "Prediction is not allowed in AAC-LC.\n");
1099 "Number of scalefactor bands in group (%d) "
1100 "exceeds limit (%d).\n",
1124 while (k < ics->max_sfb) {
1127 int sect_band_type =
get_bits(gb, 4);
1128 if (sect_band_type == 12) {
1133 sect_len_incr =
get_bits(gb, bits);
1134 sect_end += sect_len_incr;
1139 if (sect_end > ics->
max_sfb) {
1141 "Number of bands (%d) exceeds limit (%d).\n",
1145 }
while (sect_len_incr == (1 << bits) - 1);
1146 for (; k < sect_end; k++) {
1147 band_type [idx] = sect_band_type;
1148 band_type_run_end[idx++] = sect_end;
1166 unsigned int global_gain,
1169 int band_type_run_end[120])
1172 int offset[3] = { global_gain, global_gain - 90, 0 };
1176 for (i = 0; i < ics->
max_sfb;) {
1177 int run_end = band_type_run_end[idx];
1178 if (band_type[idx] ==
ZERO_BT) {
1179 for (; i < run_end; i++, idx++)
1183 for (; i < run_end; i++, idx++) {
1184 offset[2] +=
get_vlc2(gb, vlc_scalefactors.
table, 7, 3) - 60;
1185 clipped_offset = av_clip(offset[2], -155, 100);
1186 if (offset[2] != clipped_offset) {
1188 "position clipped (%d -> %d).\nIf you heard an "
1189 "audible artifact, there may be a bug in the "
1190 "decoder. ", offset[2], clipped_offset);
1194 }
else if (band_type[idx] ==
NOISE_BT) {
1195 for (; i < run_end; i++, idx++) {
1196 if (noise_flag-- > 0)
1197 offset[1] +=
get_bits(gb, 9) - 256;
1199 offset[1] +=
get_vlc2(gb, vlc_scalefactors.
table, 7, 3) - 60;
1200 clipped_offset = av_clip(offset[1], -100, 155);
1201 if (offset[1] != clipped_offset) {
1203 "(%d -> %d).\nIf you heard an audible "
1204 "artifact, there may be a bug in the decoder. ",
1205 offset[1], clipped_offset);
1210 for (; i < run_end; i++, idx++) {
1211 offset[0] +=
get_vlc2(gb, vlc_scalefactors.
table, 7, 3) - 60;
1212 if (offset[0] > 255
U) {
1214 "Scalefactor (%d) out of range.\n", offset[0]);
1229 const uint16_t *swb_offset,
int num_swb)
1234 if (pulse_swb >= num_swb)
1236 pulse->
pos[0] = swb_offset[pulse_swb];
1238 if (pulse->
pos[0] > 1023)
1241 for (i = 1; i < pulse->
num_pulse; i++) {
1243 if (pulse->
pos[i] > 1023)
1258 int w, filt, i, coef_len, coef_res, coef_compress;
1265 for (filt = 0; filt < tns->
n_filt[w]; filt++) {
1269 if ((tns->
order[w][filt] =
get_bits(gb, 5 - 2 * is8)) > tns_max_order) {
1271 "TNS filter order %d is greater than maximum %d.\n",
1272 tns->
order[w][filt], tns_max_order);
1273 tns->
order[w][filt] = 0;
1276 if (tns->
order[w][filt]) {
1279 coef_len = coef_res + 3 - coef_compress;
1280 tmp2_idx = 2 * coef_compress + coef_res;
1282 for (i = 0; i < tns->
order[w][filt]; i++)
1302 if (ms_present == 1) {
1307 }
else if (ms_present == 2) {
1313 static inline float *
VMUL2(
float *dst,
const float *v,
unsigned idx,
1317 *dst++ = v[idx & 15] * s;
1318 *dst++ = v[idx>>4 & 15] * s;
1324 static inline float *
VMUL4(
float *dst,
const float *v,
unsigned idx,
1328 *dst++ = v[idx & 3] * s;
1329 *dst++ = v[idx>>2 & 3] * s;
1330 *dst++ = v[idx>>4 & 3] * s;
1331 *dst++ = v[idx>>6 & 3] * s;
1337 static inline float *
VMUL2S(
float *dst,
const float *v,
unsigned idx,
1338 unsigned sign,
const float *
scale)
1343 s0.
i ^= sign >> 1 << 31;
1346 *dst++ = v[idx & 15] * s0.
f;
1347 *dst++ = v[idx>>4 & 15] * s1.
f;
1354 static inline float *
VMUL4S(
float *dst,
const float *v,
unsigned idx,
1355 unsigned sign,
const float *
scale)
1357 unsigned nz = idx >> 12;
1361 t.
i = s.
i ^ (sign & 1
U<<31);
1362 *dst++ = v[idx & 3] * t.
f;
1364 sign <<= nz & 1; nz >>= 1;
1365 t.
i = s.
i ^ (sign & 1
U<<31);
1366 *dst++ = v[idx>>2 & 3] * t.
f;
1368 sign <<= nz & 1; nz >>= 1;
1369 t.
i = s.
i ^ (sign & 1
U<<31);
1370 *dst++ = v[idx>>4 & 3] * t.
f;
1373 t.
i = s.
i ^ (sign & 1
U<<31);
1374 *dst++ = v[idx>>6 & 3] * t.
f;
1394 int pulse_present,
const Pulse *pulse,
1398 int i, k,
g, idx = 0;
1401 float *coef_base = coef;
1404 memset(coef + g * 128 + offsets[ics->
max_sfb], 0,
1405 sizeof(
float) * (c - offsets[ics->
max_sfb]));
1410 for (i = 0; i < ics->
max_sfb; i++, idx++) {
1411 const unsigned cbt_m1 = band_type[idx] - 1;
1412 float *cfo = coef + offsets[
i];
1413 int off_len = offsets[i + 1] - offsets[
i];
1417 for (group = 0; group < g_len; group++, cfo+=128) {
1418 memset(cfo, 0, off_len *
sizeof(
float));
1420 }
else if (cbt_m1 ==
NOISE_BT - 1) {
1421 for (group = 0; group < g_len; group++, cfo+=128) {
1425 for (k = 0; k < off_len; k++) {
1431 scale = sf[idx] / sqrtf(band_energy);
1440 switch (cbt_m1 >> 1) {
1442 for (group = 0; group < g_len; group++, cfo+=128) {
1452 cb_idx = cb_vector_idx[code];
1453 cf =
VMUL4(cf, vq, cb_idx, sf + idx);
1459 for (group = 0; group < g_len; group++, cfo+=128) {
1471 cb_idx = cb_vector_idx[code];
1472 nnz = cb_idx >> 8 & 15;
1475 cf =
VMUL4S(cf, vq, cb_idx, bits, sf + idx);
1481 for (group = 0; group < g_len; group++, cfo+=128) {
1491 cb_idx = cb_vector_idx[code];
1492 cf =
VMUL2(cf, vq, cb_idx, sf + idx);
1499 for (group = 0; group < g_len; group++, cfo+=128) {
1511 cb_idx = cb_vector_idx[code];
1512 nnz = cb_idx >> 8 & 15;
1513 sign = nnz ?
SHOW_UBITS(
re, gb, nnz) << (cb_idx >> 12) : 0;
1515 cf =
VMUL2S(cf, vq, cb_idx, sign, sf + idx);
1521 for (group = 0; group < g_len; group++, cfo+=128) {
1523 uint32_t *icf = (uint32_t *) cf;
1542 cb_idx = cb_vector_idx[code];
1548 for (j = 0; j < 2; j++) {
1570 unsigned v = ((
const uint32_t*)vq)[cb_idx & 15];
1571 *icf++ = (bits & 1
U<<31) | v;
1588 if (pulse_present) {
1590 for (i = 0; i < pulse->
num_pulse; i++) {
1591 float co = coef_base[ pulse->
pos[
i] ];
1592 while (offsets[idx + 1] <= pulse->
pos[i])
1594 if (band_type[idx] !=
NOISE_BT && sf[idx]) {
1595 float ico = -pulse->
amp[
i];
1598 ico = co / sqrtf(sqrtf(fabsf(co))) + (co > 0 ? -ico : ico);
1600 coef_base[ pulse->
pos[
i] ] =
cbrtf(fabsf(ico)) * ico * sf[idx];
1611 tmp.
i = (tmp.
i + 0x00008000
U) & 0xFFFF0000U;
1619 tmp.
i = (tmp.
i + 0x00007FFF
U + (tmp.
i & 0x00010000
U >> 16)) & 0xFFFF0000
U;
1627 pun.
i &= 0xFFFF0000
U;
1634 const float a = 0.953125;
1635 const float alpha = 0.90625;
1639 float r0 = ps->
r0, r1 = ps->
r1;
1640 float cor0 = ps->
cor0, cor1 = ps->
cor1;
1641 float var0 = ps->
var0, var1 = ps->
var1;
1643 k1 = var0 > 1 ? cor0 *
flt16_even(a / var0) : 0;
1644 k2 = var1 > 1 ? cor1 *
flt16_even(a / var1) : 0;
1679 k < sce->ics.swb_offset[sfb + 1];
1707 float *out = sce->
coeffs;
1708 int global_gain, pulse_present = 0;
1718 if (!common_window && !scale_flag) {
1735 "Pulse tool not allowed in eight short sequence.\n");
1740 "Pulse data corrupt or invalid.\n");
1770 int g,
i, group, idx = 0;
1773 for (i = 0; i < ics->
max_sfb; i++, idx++) {
1777 for (group = 0; group < ics->
group_len[
g]; group++) {
1779 ch1 + group * 128 + offsets[i],
1780 offsets[i+1] - offsets[i]);
1803 int g, group,
i, idx = 0;
1807 for (i = 0; i < ics->
max_sfb;) {
1811 for (; i < bt_run_end; i++, idx++) {
1812 c = -1 + 2 * (sce1->
band_type[idx] - 14);
1814 c *= 1 - 2 * cpe->
ms_mask[idx];
1815 scale = c * sce1->
sf[idx];
1816 for (group = 0; group < ics->
group_len[
g]; group++)
1818 coef0 + group * 128 + offsets[i],
1820 offsets[i + 1] - offsets[i]);
1824 idx += bt_run_end -
i;
1840 int i, ret, common_window, ms_present = 0;
1843 if (common_window) {
1854 if (ms_present == 3) {
1857 }
else if (ms_present)
1860 if ((ret =
decode_ics(ac, &cpe->
ch[0], gb, common_window, 0)))
1862 if ((ret =
decode_ics(ac, &cpe->
ch[1], gb, common_window, 0)))
1865 if (common_window) {
1879 1.09050773266525765921,
1880 1.18920711500272106672,
1915 scale = cce_scale[
get_bits(gb, 2)];
1920 for (c = 0; c < num_gain; c++) {
1924 float gain_cache = 1.;
1927 gain = cge ?
get_vlc2(gb, vlc_scalefactors.
table, 7, 3) - 60: 0;
1928 gain_cache =
powf(scale, -gain);
1931 coup->
gain[c][0] = gain_cache;
1934 for (sfb = 0; sfb < sce->
ics.
max_sfb; sfb++, idx++) {
1945 gain_cache =
powf(scale, -t) * s;
1948 coup->
gain[c][idx] = gain_cache;
1966 int num_excl_chan = 0;
1969 for (i = 0; i < 7; i++)
1973 return num_excl_chan / 7;
1985 int drc_num_bands = 1;
2006 for (i = 0; i < drc_num_bands; i++) {
2019 for (i = 0; i < drc_num_bands; i++) {
2089 int bottom, top, order, start, end,
size, inc;
2095 for (filt = 0; filt < tns->
n_filt[w]; filt++) {
2098 order = tns->
order[w][filt];
2107 if ((size = end - start) <= 0)
2119 for (m = 0; m <
size; m++, start += inc)
2120 for (i = 1; i <=
FFMIN(m, order); i++)
2121 coef[start] -= coef[start - i * inc] * lpc[i - 1];
2124 for (m = 0; m <
size; m++, start += inc) {
2125 tmp[0] = coef[start];
2126 for (i = 1; i <=
FFMIN(m, order); i++)
2127 coef[start] += tmp[i] * lpc[i - 1];
2128 for (i = order; i > 0; i--)
2129 tmp[i] = tmp[i - 1];
2151 memset(in, 0, 448 *
sizeof(
float));
2158 memset(in + 1024 + 576, 0, 448 *
sizeof(
float));
2173 float *predTime = sce->
ret;
2175 int16_t num_samples = 2048;
2177 if (ltp->
lag < 1024)
2178 num_samples = ltp->
lag + 1024;
2179 for (i = 0; i < num_samples; i++)
2181 memset(&predTime[i], 0, (2048 - i) *
sizeof(
float));
2190 for (i = offsets[sfb]; i < offsets[sfb + 1]; i++)
2191 sce->
coeffs[i] += predFreq[i];
2201 float *saved = sce->
saved;
2202 float *saved_ltp = sce->
coeffs;
2208 memcpy(saved_ltp, saved, 512 *
sizeof(
float));
2209 memset(saved_ltp + 576, 0, 448 *
sizeof(
float));
2211 for (i = 0; i < 64; i++)
2212 saved_ltp[i + 512] = ac->
buf_mdct[1023 - i] * swindow[63 - i];
2214 memcpy(saved_ltp, ac->
buf_mdct + 512, 448 *
sizeof(
float));
2215 memset(saved_ltp + 576, 0, 448 *
sizeof(
float));
2217 for (i = 0; i < 64; i++)
2218 saved_ltp[i + 512] = ac->
buf_mdct[1023 - i] * swindow[63 - i];
2221 for (i = 0; i < 512; i++)
2222 saved_ltp[i + 512] = ac->
buf_mdct[1023 - i] * lwindow[511 - i];
2237 float *out = sce->
ret;
2238 float *saved = sce->
saved;
2243 float *temp = ac->
temp;
2248 for (i = 0; i < 1024; i += 128)
2263 memcpy( out, saved, 448 *
sizeof(
float));
2271 memcpy( out + 448 + 4*128, temp, 64 *
sizeof(
float));
2274 memcpy( out + 576, buf + 64, 448 *
sizeof(
float));
2280 memcpy( saved, temp + 64, 64 *
sizeof(
float));
2284 memcpy( saved + 448, buf + 7*128 + 64, 64 *
sizeof(
float));
2286 memcpy( saved, buf + 512, 448 *
sizeof(
float));
2287 memcpy( saved + 448, buf + 7*128 + 64, 64 *
sizeof(
float));
2289 memcpy( saved, buf + 512, 512 *
sizeof(
float));
2304 float *dest = target->
coeffs;
2305 const float *src = cce->
ch[0].
coeffs;
2306 int g,
i, group, k, idx = 0;
2309 "Dependent coupling is not supported together with LTP\n");
2313 for (i = 0; i < ics->
max_sfb; i++, idx++) {
2316 for (group = 0; group < ics->
group_len[
g]; group++) {
2317 for (k = offsets[i]; k < offsets[i + 1]; k++) {
2319 dest[group * 128 + k] += gain * src[group * 128 + k];
2340 const float *src = cce->
ch[0].
ret;
2341 float *dest = target->
ret;
2344 for (i = 0; i <
len; i++)
2345 dest[i] += gain * src[i];
2368 if (coup->
type[c] == type && coup->
id_select[c] == elem_id) {
2370 apply_coupling_method(ac, &cc->
ch[0], cce, index);
2375 apply_coupling_method(ac, &cc->
ch[1], cce, index++);
2389 for (type = 3; type >= 0; type--) {
2433 uint8_t layout_map[MAX_ELEM_ID*4][3];
2434 int layout_map_tags;
2476 int samples = 0, multiplier, audio_found = 0, pce_found = 0;
2502 if (!(che=
get_che(ac, elem_type, elem_id))) {
2504 elem_type, elem_id);
2511 switch (elem_type) {
2537 uint8_t layout_map[MAX_ELEM_ID*4][3];
2547 "Not evaluating a further program_config_element as this construct is dubious at best.\n");
2575 elem_type_prev = elem_type;
2590 samples <<= multiplier;
2598 if (ac->
oc[1].
status && audio_found) {
2611 int *got_frame_ptr,
AVPacket *avpkt)
2615 int buf_size = avpkt->
size;
2620 int new_extradata_size;
2623 &new_extradata_size);
2625 if (new_extradata) {
2632 memcpy(avctx->
extradata, new_extradata, new_extradata_size);
2649 for (buf_offset = buf_consumed; buf_offset < buf_size; buf_offset++)
2650 if (buf[buf_offset])
2653 return buf_size > buf_offset ? buf_consumed : buf_size;
2662 for (type = 0; type < 4; type++) {
2663 if (ac->
che[type][i])
2676 #define LOAS_SYNC_WORD 0x2b7
2702 int sync_extension = 0;
2703 int bits_consumed, esize;
2711 if (config_start_bit % 8) {
2713 "Non-byte-aligned audio-specific config", 1);
2719 gb->
buffer + (config_start_bit / 8),
2720 asclen, sync_extension);
2722 if (bits_consumed < 0)
2731 esize = (bits_consumed+7) / 8;
2746 return bits_consumed;
2752 int ret, audio_mux_version =
get_bits(gb, 1);
2755 if (audio_mux_version)
2760 if (audio_mux_version)
2768 "Multiple programs", 1);
2777 "Multiple layers", 1);
2782 if (!audio_mux_version) {
2813 if (audio_mux_version) {
2836 int mux_slot_length = 0;
2839 mux_slot_length += tmp;
2840 }
while (tmp == 255);
2841 return mux_slot_length;
2857 if (!use_same_mux) {
2862 "no decoder config found\n");
2870 }
else if (mux_slot_length_bytes * 8 + 256 <
get_bits_left(gb)) {
2872 "frame length mismatch %d << %d\n",
2882 int *got_frame_ptr,
AVPacket *avpkt)
2897 if (muxlength > avpkt->
size)
2921 "ADTS header detected, probably as result of configuration "