115 static const char overread_err[] =
"Input buffer exhausted before END element found\n";
168 int i, type, sum = 0;
170 for (type = 0; type < 4; type++) {
193 int type,
int id,
int *channels)
197 if (che_pos[type][
id]) {
198 if (!ac->
che[type][
id]) {
211 if (ac->
che[type][
id])
229 int channel_config,
enum OCStatus oc_type)
232 int i, type, channels = 0, ret;
234 if (new_che_pos != che_pos)
235 memcpy(che_pos, new_che_pos, 4 * MAX_ELEM_ID *
sizeof(new_che_pos[0][0]));
237 if (channel_config) {
246 memset(ac->
tag_che_map, 0, 4 * MAX_ELEM_ID *
sizeof(ac->
che[0][0]));
260 for (type = 0; type < 4; type++) {
307 int num_front, num_side, num_back, num_lfe, num_assoc_data, num_cc, sampling_index;
314 av_log(avctx,
AV_LOG_WARNING,
"Sample rate index in program config element does not match the sample rate index configured by the container.\n");
364 if (channel_config < 1 || channel_config > 7) {
381 if (channel_config != 2)
383 if (channel_config > 1)
385 if (channel_config == 4)
387 if (channel_config > 4)
388 new_che_pos[
TYPE_CPE][(channel_config == 7) + 1]
390 if (channel_config > 5)
392 if (channel_config == 7)
412 int extension_flag, ret;
427 memset(new_che_pos, 0, 4 *
MAX_ELEM_ID *
sizeof(new_che_pos[0][0]));
428 if (channel_config == 0) {
430 if ((ret =
decode_pce(avctx, m4ac, new_che_pos, gb)))
439 }
else if (m4ac->
sbr == 1 && m4ac->
ps == -1)
445 if (extension_flag) {
481 const uint8_t *
data,
int bit_size,
518 av_dlog(avctx,
"AOT %d chan config %d sampling index %d (%d) SBR %d PS %d\n",
534 return previous_val * 1664525 + 1013904223;
556 if (92017 <= rate)
return 0;
557 else if (75132 <= rate)
return 1;
558 else if (55426 <= rate)
return 2;
559 else if (46009 <= rate)
return 3;
560 else if (37566 <= rate)
return 4;
561 else if (27713 <= rate)
return 5;
562 else if (23004 <= rate)
return 6;
563 else if (18783 <= rate)
return 7;
564 else if (13856 <= rate)
return 8;
565 else if (11502 <= rate)
return 9;
566 else if (9391 <= rate)
return 10;
577 #define AAC_INIT_VLC_STATIC(num, size) \
578 INIT_VLC_STATIC(&vlc_spectral[num], 8, ff_aac_spectral_sizes[num], \
579 ff_aac_spectral_bits[num], sizeof( ff_aac_spectral_bits[num][0]), sizeof( ff_aac_spectral_bits[num][0]), \
580 ff_aac_spectral_codes[num], sizeof(ff_aac_spectral_codes[num][0]), sizeof(ff_aac_spectral_codes[num][0]), \
587 float output_scale_factor;
626 output_scale_factor = 1.0 / 32768.0;
629 output_scale_factor = 1.0;
745 for (i = 0; i < 7; i++) {
783 "Number of scalefactor bands in group (%d) exceeds limit (%d).\n",
807 while (k < ics->max_sfb) {
808 uint8_t sect_end = k;
810 int sect_band_type =
get_bits(gb, 4);
811 if (sect_band_type == 12) {
817 sect_end += sect_len_incr;
824 "Number of bands (%d) exceeds limit (%d).\n",
828 }
while (sect_len_incr == (1 << bits) - 1);
829 for (; k < sect_end; k++) {
830 band_type [idx] = sect_band_type;
831 band_type_run_end[idx++] = sect_end;
849 unsigned int global_gain,
852 int band_type_run_end[120])
855 int offset[3] = { global_gain, global_gain - 90, 0 };
858 static const char *sf_str[3] = {
"Global gain",
"Noise gain",
"Intensity stereo position" };
860 for (i = 0; i < ics->
max_sfb;) {
861 int run_end = band_type_run_end[idx];
862 if (band_type[idx] ==
ZERO_BT) {
863 for (; i < run_end; i++, idx++)
866 for (; i < run_end; i++, idx++) {
867 offset[2] +=
get_vlc2(gb, vlc_scalefactors.
table, 7, 3) - 60;
868 clipped_offset = av_clip(offset[2], -155, 100);
869 if (offset[2] != clipped_offset) {
871 "position clipped (%d -> %d).\nIf you heard an "
872 "audible artifact, there may be a bug in the "
873 "decoder. ", offset[2], clipped_offset);
877 }
else if (band_type[idx] ==
NOISE_BT) {
878 for (; i < run_end; i++, idx++) {
879 if (noise_flag-- > 0)
882 offset[1] +=
get_vlc2(gb, vlc_scalefactors.
table, 7, 3) - 60;
883 clipped_offset = av_clip(offset[1], -100, 155);
884 if (offset[1] != clipped_offset) {
886 "(%d -> %d).\nIf you heard an audible "
887 "artifact, there may be a bug in the decoder. ",
888 offset[1], clipped_offset);
893 for (; i < run_end; i++, idx++) {
894 offset[0] +=
get_vlc2(gb, vlc_scalefactors.
table, 7, 3) - 60;
895 if (offset[0] > 255
U) {
897 "%s (%d) out of range.\n", sf_str[0], offset[0]);
912 const uint16_t *swb_offset,
int num_swb)
917 if (pulse_swb >= num_swb)
919 pulse->
pos[0] = swb_offset[pulse_swb];
921 if (pulse->
pos[0] > 1023)
926 if (pulse->
pos[i] > 1023)
941 int w, filt, i, coef_len, coef_res, coef_compress;
948 for (filt = 0; filt < tns->
n_filt[w]; filt++) {
952 if ((tns->
order[w][filt] =
get_bits(gb, 5 - 2 * is8)) > tns_max_order) {
954 tns->
order[w][filt], tns_max_order);
955 tns->
order[w][filt] = 0;
958 if (tns->
order[w][filt]) {
961 coef_len = coef_res + 3 - coef_compress;
962 tmp2_idx = 2 * coef_compress + coef_res;
964 for (i = 0; i < tns->
order[w][filt]; i++)
984 if (ms_present == 1) {
987 }
else if (ms_present == 2) {
993 static inline float *
VMUL2(
float *dst,
const float *
v,
unsigned idx,
997 *dst++ = v[idx & 15] * s;
998 *dst++ = v[idx>>4 & 15] * s;
1004 static inline float *
VMUL4(
float *dst,
const float *
v,
unsigned idx,
1008 *dst++ = v[idx & 3] * s;
1009 *dst++ = v[idx>>2 & 3] * s;
1010 *dst++ = v[idx>>4 & 3] * s;
1011 *dst++ = v[idx>>6 & 3] * s;
1017 static inline float *
VMUL2S(
float *dst,
const float *
v,
unsigned idx,
1018 unsigned sign,
const float *
scale)
1023 s0.
i ^= sign >> 1 << 31;
1026 *dst++ = v[idx & 15] * s0.
f;
1027 *dst++ = v[idx>>4 & 15] * s1.
f;
1034 static inline float *
VMUL4S(
float *dst,
const float *
v,
unsigned idx,
1035 unsigned sign,
const float *
scale)
1037 unsigned nz = idx >> 12;
1041 t.
i = s.
i ^ (sign & 1
U<<31);
1042 *dst++ = v[idx & 3] * t.
f;
1044 sign <<= nz & 1; nz >>= 1;
1045 t.
i = s.
i ^ (sign & 1
U<<31);
1046 *dst++ = v[idx>>2 & 3] * t.
f;
1048 sign <<= nz & 1; nz >>= 1;
1049 t.
i = s.
i ^ (sign & 1
U<<31);
1050 *dst++ = v[idx>>4 & 3] * t.
f;
1052 sign <<= nz & 1; nz >>= 1;
1053 t.
i = s.
i ^ (sign & 1
U<<31);
1054 *dst++ = v[idx>>6 & 3] * t.
f;
1074 int pulse_present,
const Pulse *pulse,
1078 int i, k,
g, idx = 0;
1081 float *coef_base = coef;
1084 memset(coef + g * 128 + offsets[ics->
max_sfb], 0,
sizeof(
float) * (c - offsets[ics->
max_sfb]));
1089 for (i = 0; i < ics->
max_sfb; i++, idx++) {
1090 const unsigned cbt_m1 = band_type[idx] - 1;
1091 float *cfo = coef + offsets[
i];
1092 int off_len = offsets[i + 1] - offsets[
i];
1096 for (group = 0; group < g_len; group++, cfo+=128) {
1097 memset(cfo, 0, off_len *
sizeof(
float));
1099 }
else if (cbt_m1 ==
NOISE_BT - 1) {
1100 for (group = 0; group < g_len; group++, cfo+=128) {
1104 for (k = 0; k < off_len; k++) {
1110 scale = sf[idx] / sqrtf(band_energy);
1119 switch (cbt_m1 >> 1) {
1121 for (group = 0; group < g_len; group++, cfo+=128) {
1131 cb_idx = cb_vector_idx[code];
1132 cf =
VMUL4(cf, vq, cb_idx, sf + idx);
1138 for (group = 0; group < g_len; group++, cfo+=128) {
1150 cb_idx = cb_vector_idx[code];
1151 nnz = cb_idx >> 8 & 15;
1154 cf =
VMUL4S(cf, vq, cb_idx, bits, sf + idx);
1160 for (group = 0; group < g_len; group++, cfo+=128) {
1170 cb_idx = cb_vector_idx[code];
1171 cf =
VMUL2(cf, vq, cb_idx, sf + idx);
1178 for (group = 0; group < g_len; group++, cfo+=128) {
1190 cb_idx = cb_vector_idx[code];
1191 nnz = cb_idx >> 8 & 15;
1192 sign = nnz ?
SHOW_UBITS(
re, gb, nnz) << (cb_idx >> 12) : 0;
1194 cf =
VMUL2S(cf, vq, cb_idx, sign, sf + idx);
1200 for (group = 0; group < g_len; group++, cfo+=128) {
1202 uint32_t *icf = (uint32_t *) cf;
1221 cb_idx = cb_vector_idx[code];
1227 for (j = 0; j < 2; j++) {
1235 b = 31 - av_log2(~b);
1249 unsigned v = ((
const uint32_t*)vq)[cb_idx & 15];
1250 *icf++ = (bits & 1
U<<31) | v;
1267 if (pulse_present) {
1269 for (i = 0; i < pulse->
num_pulse; i++) {
1270 float co = coef_base[ pulse->
pos[
i] ];
1271 while (offsets[idx + 1] <= pulse->
pos[i])
1273 if (band_type[idx] !=
NOISE_BT && sf[idx]) {
1274 float ico = -pulse->
amp[
i];
1277 ico = co / sqrtf(sqrtf(fabsf(co))) + (co > 0 ? -ico : ico);
1279 coef_base[ pulse->
pos[
i] ] = cbrtf(fabsf(ico)) * ico * sf[idx];
1290 tmp.
i = (tmp.
i + 0x00008000
U) & 0xFFFF0000U;
1298 tmp.
i = (tmp.
i + 0x00007FFF
U + (tmp.
i & 0x00010000
U >> 16)) & 0xFFFF0000
U;
1306 pun.
i &= 0xFFFF0000
U;
1313 const float a = 0.953125;
1314 const float alpha = 0.90625;
1318 float r0 = ps->
r0, r1 = ps->
r1;
1319 float cor0 = ps->
cor0, cor1 = ps->
cor1;
1320 float var0 = ps->
var0, var1 = ps->
var1;
1322 k1 = var0 > 1 ? cor0 *
flt16_even(a / var0) : 0;
1323 k2 = var1 > 1 ? cor1 *
flt16_even(a / var1) : 0;
1355 for (k = sce->
ics.
swb_offset[sfb]; k < sce->ics.swb_offset[sfb + 1]; k++) {
1380 float *out = sce->
coeffs;
1381 int global_gain, pulse_present = 0;
1391 if (!common_window && !scale_flag) {
1440 int g,
i, group, idx = 0;
1443 for (i = 0; i < ics->
max_sfb; i++, idx++) {
1446 for (group = 0; group < ics->
group_len[
g]; group++) {
1448 ch1 + group * 128 + offsets[i],
1449 offsets[i+1] - offsets[i]);
1471 int g, group,
i, idx = 0;
1475 for (i = 0; i < ics->
max_sfb;) {
1478 for (; i < bt_run_end; i++, idx++) {
1479 c = -1 + 2 * (sce1->
band_type[idx] - 14);
1481 c *= 1 - 2 * cpe->
ms_mask[idx];
1482 scale = c * sce1->
sf[idx];
1483 for (group = 0; group < ics->
group_len[
g]; group++)
1485 coef0 + group * 128 + offsets[i],
1487 offsets[i + 1] - offsets[i]);
1491 idx += bt_run_end -
i;
1507 int i, ret, common_window, ms_present = 0;
1510 if (common_window) {
1520 if (ms_present == 3) {
1523 }
else if (ms_present)
1526 if ((ret =
decode_ics(ac, &cpe->
ch[0], gb, common_window, 0)))
1528 if ((ret =
decode_ics(ac, &cpe->
ch[1], gb, common_window, 0)))
1531 if (common_window) {
1545 1.09050773266525765921,
1546 1.18920711500272106672,
1586 for (c = 0; c < num_gain; c++) {
1590 float gain_cache = 1.;
1593 gain = cge ?
get_vlc2(gb, vlc_scalefactors.
table, 7, 3) - 60: 0;
1594 gain_cache = powf(scale, -gain);
1597 coup->
gain[c][0] = gain_cache;
1600 for (sfb = 0; sfb < sce->
ics.
max_sfb; sfb++, idx++) {
1611 gain_cache = powf(scale, -t) * s;
1614 coup->
gain[c][idx] = gain_cache;
1632 int num_excl_chan = 0;
1635 for (i = 0; i < 7; i++)
1639 return num_excl_chan / 7;
1653 int drc_num_bands = 1;
1674 for (i = 0; i < drc_num_bands; i++) {
1687 for (i = 0; i < drc_num_bands; i++) {
1756 int bottom, top, order, start, end,
size, inc;
1762 for (filt = 0; filt < tns->
n_filt[w]; filt++) {
1765 order = tns->
order[w][filt];
1774 if ((size = end - start) <= 0)
1786 for (m = 0; m <
size; m++, start += inc)
1787 for (i = 1; i <=
FFMIN(m, order); i++)
1788 coef[start] -= coef[start - i * inc] * lpc[i - 1];
1791 for (m = 0; m <
size; m++, start += inc) {
1792 tmp[0] = coef[start];
1793 for (i = 1; i <=
FFMIN(m, order); i++)
1794 coef[start] += tmp[i] * lpc[i - 1];
1795 for (i = order; i > 0; i--)
1796 tmp[i] = tmp[i - 1];
1818 memset(in, 0, 448 *
sizeof(
float));
1825 memset(in + 1024 + 576, 0, 448 *
sizeof(
float));
1840 float *predTime = sce->
ret;
1842 int16_t num_samples = 2048;
1844 if (ltp->
lag < 1024)
1845 num_samples = ltp->
lag + 1024;
1846 for (i = 0; i < num_samples; i++)
1848 memset(&predTime[i], 0, (2048 - i) *
sizeof(
float));
1857 for (i = offsets[sfb]; i < offsets[sfb + 1]; i++)
1858 sce->
coeffs[i] += predFreq[i];
1868 float *saved = sce->
saved;
1869 float *saved_ltp = sce->
coeffs;
1875 memcpy(saved_ltp, saved, 512 *
sizeof(
float));
1876 memset(saved_ltp + 576, 0, 448 *
sizeof(
float));
1878 for (i = 0; i < 64; i++)
1879 saved_ltp[i + 512] = ac->
buf_mdct[1023 - i] * swindow[63 - i];
1881 memcpy(saved_ltp, ac->
buf_mdct + 512, 448 *
sizeof(
float));
1882 memset(saved_ltp + 576, 0, 448 *
sizeof(
float));
1884 for (i = 0; i < 64; i++)
1885 saved_ltp[i + 512] = ac->
buf_mdct[1023 - i] * swindow[63 - i];
1888 for (i = 0; i < 512; i++)
1889 saved_ltp[i + 512] = ac->
buf_mdct[1023 - i] * lwindow[511 - i];
1904 float *out = sce->
ret;
1905 float *saved = sce->
saved;
1910 float *temp = ac->
temp;
1915 for (i = 0; i < 1024; i += 128)
1930 memcpy( out, saved, 448 *
sizeof(
float));
1938 memcpy( out + 448 + 4*128, temp, 64 *
sizeof(
float));
1941 memcpy( out + 576, buf + 64, 448 *
sizeof(
float));
1947 memcpy( saved, temp + 64, 64 *
sizeof(
float));
1951 memcpy( saved + 448, buf + 7*128 + 64, 64 *
sizeof(
float));
1953 memcpy( saved, buf + 512, 448 *
sizeof(
float));
1954 memcpy( saved + 448, buf + 7*128 + 64, 64 *
sizeof(
float));
1956 memcpy( saved, buf + 512, 512 *
sizeof(
float));
1971 float *dest = target->
coeffs;
1972 const float *src = cce->
ch[0].
coeffs;
1973 int g,
i, group, k, idx = 0;
1976 "Dependent coupling is not supported together with LTP\n");
1980 for (i = 0; i < ics->
max_sfb; i++, idx++) {
1983 for (group = 0; group < ics->
group_len[
g]; group++) {
1984 for (k = offsets[i]; k < offsets[i + 1]; k++) {
1986 dest[group * 128 + k] += gain * src[group * 128 + k];
2007 const float *src = cce->
ch[0].
ret;
2008 float *dest = target->
ret;
2011 for (i = 0; i <
len; i++)
2012 dest[i] += gain * src[i];
2035 if (coup->
type[c] == type && coup->
id_select[c] == elem_id) {
2037 apply_coupling_method(ac, &cc->
ch[0], cce, index);
2042 apply_coupling_method(ac, &cc->
ch[1], cce, index++);
2056 for (type = 3; type >= 0; type--) {
2105 memset(new_che_pos, 0, 4 *
MAX_ELEM_ID *
sizeof(new_che_pos[0][0]));
2143 int samples = 0, multiplier, audio_found = 0;
2162 if (!(che=
get_che(ac, elem_type, elem_id))) {
2164 elem_type, elem_id);
2170 switch (elem_type) {
2197 memset(new_che_pos, 0, 4 *
MAX_ELEM_ID *
sizeof(new_che_pos[0][0]));
2202 "Not evaluating a further program_config_element as this construct is dubious at best.\n");
2226 elem_type_prev = elem_type;
2240 samples <<= multiplier;
2274 int *got_frame_ptr,
AVPacket *avpkt)
2277 const uint8_t *buf = avpkt->
data;
2278 int buf_size = avpkt->
size;
2283 int new_extradata_size;
2286 &new_extradata_size);
2288 if (new_extradata) {
2295 memcpy(avctx->
extradata, new_extradata, new_extradata_size);
2309 for (buf_offset = buf_consumed; buf_offset < buf_size; buf_offset++)
2310 if (buf[buf_offset])
2313 return buf_size > buf_offset ? buf_consumed : buf_size;
2322 for (type = 0; type < 4; type++) {
2323 if (ac->
che[type][i])
2336 #define LOAS_SYNC_WORD 0x2b7
2362 int sync_extension = 0;
2363 int bits_consumed, esize;
2371 if (config_start_bit % 8) {
2373 "config not byte aligned.\n", 1);
2379 gb->
buffer + (config_start_bit / 8),
2380 asclen, sync_extension);
2382 if (bits_consumed < 0)
2391 esize = (bits_consumed+7) / 8;
2406 return bits_consumed;
2412 int ret, audio_mux_version =
get_bits(gb, 1);
2415 if (audio_mux_version)
2420 if (audio_mux_version)
2428 "multiple programs are not supported\n", 1);
2437 "multiple layers are not supported\n", 1);
2442 if (!audio_mux_version) {
2473 if (audio_mux_version) {
2496 int mux_slot_length = 0;
2499 mux_slot_length += tmp;
2500 }
while (tmp == 255);
2501 return mux_slot_length;
2516 uint8_t use_same_mux =
get_bits(gb, 1);
2517 if (!use_same_mux) {
2522 "no decoder config found\n");
2530 }
else if (mux_slot_length_bytes * 8 + 256 <
get_bits_left(gb)) {
2532 "frame length mismatch %d << %d\n",
2542 int *got_frame_ptr,
AVPacket *avpkt)
2557 if (muxlength > avpkt->
size)
2578 "ADTS header detected, probably as result of configuration "