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)
195 if (che_pos[type][
id]) {
196 if (!ac->
che[type][
id]) {
209 if (ac->
che[type][
id])
227 int channel_config,
enum OCStatus oc_type)
230 int i, type, channels = 0, ret;
232 if (new_che_pos != che_pos)
233 memcpy(che_pos, new_che_pos, 4 * MAX_ELEM_ID *
sizeof(new_che_pos[0][0]));
235 if (channel_config) {
244 memset(ac->
tag_che_map, 0, 4 * MAX_ELEM_ID *
sizeof(ac->
che[0][0]));
258 for (type = 0; type < 4; type++) {
305 int num_front, num_side, num_back, num_lfe, num_assoc_data, num_cc, sampling_index;
312 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");
362 if (channel_config < 1 || channel_config > 7) {
379 if (channel_config != 2)
381 if (channel_config > 1)
383 if (channel_config == 4)
385 if (channel_config > 4)
386 new_che_pos[
TYPE_CPE][(channel_config == 7) + 1]
388 if (channel_config > 5)
390 if (channel_config == 7)
410 int extension_flag, ret;
425 memset(new_che_pos, 0, 4 *
MAX_ELEM_ID *
sizeof(new_che_pos[0][0]));
426 if (channel_config == 0) {
428 if ((ret =
decode_pce(avctx, m4ac, new_che_pos, gb)))
437 }
else if (m4ac->
sbr == 1 && m4ac->
ps == -1)
443 if (extension_flag) {
479 const uint8_t *
data,
int bit_size,
514 av_dlog(avctx,
"AOT %d chan config %d sampling index %d (%d) SBR %d PS %d\n",
530 return previous_val * 1664525 + 1013904223;
552 if (92017 <= rate)
return 0;
553 else if (75132 <= rate)
return 1;
554 else if (55426 <= rate)
return 2;
555 else if (46009 <= rate)
return 3;
556 else if (37566 <= rate)
return 4;
557 else if (27713 <= rate)
return 5;
558 else if (23004 <= rate)
return 6;
559 else if (18783 <= rate)
return 7;
560 else if (13856 <= rate)
return 8;
561 else if (11502 <= rate)
return 9;
562 else if (9391 <= rate)
return 10;
573 #define AAC_INIT_VLC_STATIC(num, size) \
574 INIT_VLC_STATIC(&vlc_spectral[num], 8, ff_aac_spectral_sizes[num], \
575 ff_aac_spectral_bits[num], sizeof( ff_aac_spectral_bits[num][0]), sizeof( ff_aac_spectral_bits[num][0]), \
576 ff_aac_spectral_codes[num], sizeof(ff_aac_spectral_codes[num][0]), sizeof(ff_aac_spectral_codes[num][0]), \
582 float output_scale_factor;
621 output_scale_factor = 1.0 / 32768.0;
624 output_scale_factor = 1.0;
740 for (i = 0; i < 7; i++) {
778 "Number of scalefactor bands in group (%d) exceeds limit (%d).\n",
802 while (k < ics->max_sfb) {
803 uint8_t sect_end = k;
805 int sect_band_type =
get_bits(gb, 4);
806 if (sect_band_type == 12) {
812 sect_end += sect_len_incr;
819 "Number of bands (%d) exceeds limit (%d).\n",
823 }
while (sect_len_incr == (1 << bits) - 1);
824 for (; k < sect_end; k++) {
825 band_type [idx] = sect_band_type;
826 band_type_run_end[idx++] = sect_end;
844 unsigned int global_gain,
847 int band_type_run_end[120])
850 int offset[3] = { global_gain, global_gain - 90, 0 };
853 static const char *sf_str[3] = {
"Global gain",
"Noise gain",
"Intensity stereo position" };
855 for (i = 0; i < ics->
max_sfb;) {
856 int run_end = band_type_run_end[idx];
857 if (band_type[idx] ==
ZERO_BT) {
858 for (; i < run_end; i++, idx++)
861 for (; i < run_end; i++, idx++) {
862 offset[2] +=
get_vlc2(gb, vlc_scalefactors.
table, 7, 3) - 60;
863 clipped_offset = av_clip(offset[2], -155, 100);
864 if (offset[2] != clipped_offset) {
866 "position clipped (%d -> %d).\nIf you heard an "
867 "audible artifact, there may be a bug in the "
868 "decoder. ", offset[2], clipped_offset);
872 }
else if (band_type[idx] ==
NOISE_BT) {
873 for (; i < run_end; i++, idx++) {
874 if (noise_flag-- > 0)
877 offset[1] +=
get_vlc2(gb, vlc_scalefactors.
table, 7, 3) - 60;
878 clipped_offset = av_clip(offset[1], -100, 155);
879 if (offset[1] != clipped_offset) {
881 "(%d -> %d).\nIf you heard an audible "
882 "artifact, there may be a bug in the decoder. ",
883 offset[1], clipped_offset);
888 for (; i < run_end; i++, idx++) {
889 offset[0] +=
get_vlc2(gb, vlc_scalefactors.
table, 7, 3) - 60;
890 if (offset[0] > 255
U) {
892 "%s (%d) out of range.\n", sf_str[0], offset[0]);
907 const uint16_t *swb_offset,
int num_swb)
912 if (pulse_swb >= num_swb)
914 pulse->
pos[0] = swb_offset[pulse_swb];
916 if (pulse->
pos[0] > 1023)
921 if (pulse->
pos[i] > 1023)
936 int w, filt, i, coef_len, coef_res, coef_compress;
943 for (filt = 0; filt < tns->
n_filt[w]; filt++) {
947 if ((tns->
order[w][filt] =
get_bits(gb, 5 - 2 * is8)) > tns_max_order) {
949 tns->
order[w][filt], tns_max_order);
950 tns->
order[w][filt] = 0;
953 if (tns->
order[w][filt]) {
956 coef_len = coef_res + 3 - coef_compress;
957 tmp2_idx = 2 * coef_compress + coef_res;
959 for (i = 0; i < tns->
order[w][filt]; i++)
979 if (ms_present == 1) {
982 }
else if (ms_present == 2) {
988 static inline float *
VMUL2(
float *dst,
const float *
v,
unsigned idx,
992 *dst++ = v[idx & 15] * s;
993 *dst++ = v[idx>>4 & 15] * s;
999 static inline float *
VMUL4(
float *dst,
const float *
v,
unsigned idx,
1003 *dst++ = v[idx & 3] * s;
1004 *dst++ = v[idx>>2 & 3] * s;
1005 *dst++ = v[idx>>4 & 3] * s;
1006 *dst++ = v[idx>>6 & 3] * s;
1012 static inline float *
VMUL2S(
float *dst,
const float *
v,
unsigned idx,
1013 unsigned sign,
const float *
scale)
1018 s0.
i ^= sign >> 1 << 31;
1021 *dst++ = v[idx & 15] * s0.
f;
1022 *dst++ = v[idx>>4 & 15] * s1.
f;
1029 static inline float *
VMUL4S(
float *dst,
const float *
v,
unsigned idx,
1030 unsigned sign,
const float *
scale)
1032 unsigned nz = idx >> 12;
1036 t.
i = s.
i ^ (sign & 1
U<<31);
1037 *dst++ = v[idx & 3] * t.
f;
1039 sign <<= nz & 1; nz >>= 1;
1040 t.
i = s.
i ^ (sign & 1
U<<31);
1041 *dst++ = v[idx>>2 & 3] * t.
f;
1043 sign <<= nz & 1; nz >>= 1;
1044 t.
i = s.
i ^ (sign & 1
U<<31);
1045 *dst++ = v[idx>>4 & 3] * t.
f;
1047 sign <<= nz & 1; nz >>= 1;
1048 t.
i = s.
i ^ (sign & 1
U<<31);
1049 *dst++ = v[idx>>6 & 3] * t.
f;
1069 int pulse_present,
const Pulse *pulse,
1073 int i, k,
g, idx = 0;
1076 float *coef_base = coef;
1079 memset(coef + g * 128 + offsets[ics->
max_sfb], 0,
sizeof(
float) * (c - offsets[ics->
max_sfb]));
1084 for (i = 0; i < ics->
max_sfb; i++, idx++) {
1085 const unsigned cbt_m1 = band_type[idx] - 1;
1086 float *cfo = coef + offsets[
i];
1087 int off_len = offsets[i + 1] - offsets[
i];
1091 for (group = 0; group < g_len; group++, cfo+=128) {
1092 memset(cfo, 0, off_len *
sizeof(
float));
1094 }
else if (cbt_m1 ==
NOISE_BT - 1) {
1095 for (group = 0; group < g_len; group++, cfo+=128) {
1099 for (k = 0; k < off_len; k++) {
1105 scale = sf[idx] / sqrtf(band_energy);
1114 switch (cbt_m1 >> 1) {
1116 for (group = 0; group < g_len; group++, cfo+=128) {
1126 cb_idx = cb_vector_idx[code];
1127 cf =
VMUL4(cf, vq, cb_idx, sf + idx);
1133 for (group = 0; group < g_len; group++, cfo+=128) {
1145 cb_idx = cb_vector_idx[code];
1146 nnz = cb_idx >> 8 & 15;
1149 cf =
VMUL4S(cf, vq, cb_idx, bits, sf + idx);
1155 for (group = 0; group < g_len; group++, cfo+=128) {
1165 cb_idx = cb_vector_idx[code];
1166 cf =
VMUL2(cf, vq, cb_idx, sf + idx);
1173 for (group = 0; group < g_len; group++, cfo+=128) {
1185 cb_idx = cb_vector_idx[code];
1186 nnz = cb_idx >> 8 & 15;
1187 sign = nnz ?
SHOW_UBITS(
re, gb, nnz) << (cb_idx >> 12) : 0;
1189 cf =
VMUL2S(cf, vq, cb_idx, sign, sf + idx);
1195 for (group = 0; group < g_len; group++, cfo+=128) {
1197 uint32_t *icf = (uint32_t *) cf;
1216 cb_idx = cb_vector_idx[code];
1222 for (j = 0; j < 2; j++) {
1230 b = 31 - av_log2(~b);
1244 unsigned v = ((
const uint32_t*)vq)[cb_idx & 15];
1245 *icf++ = (bits & 1
U<<31) | v;
1262 if (pulse_present) {
1264 for (i = 0; i < pulse->
num_pulse; i++) {
1265 float co = coef_base[ pulse->
pos[
i] ];
1266 while (offsets[idx + 1] <= pulse->
pos[i])
1268 if (band_type[idx] !=
NOISE_BT && sf[idx]) {
1269 float ico = -pulse->
amp[
i];
1272 ico = co / sqrtf(sqrtf(fabsf(co))) + (co > 0 ? -ico : ico);
1274 coef_base[ pulse->
pos[
i] ] = cbrtf(fabsf(ico)) * ico * sf[idx];
1285 tmp.
i = (tmp.
i + 0x00008000
U) & 0xFFFF0000U;
1293 tmp.
i = (tmp.
i + 0x00007FFF
U + (tmp.
i & 0x00010000
U >> 16)) & 0xFFFF0000
U;
1301 pun.
i &= 0xFFFF0000
U;
1308 const float a = 0.953125;
1309 const float alpha = 0.90625;
1313 float r0 = ps->
r0, r1 = ps->
r1;
1314 float cor0 = ps->
cor0, cor1 = ps->
cor1;
1315 float var0 = ps->
var0, var1 = ps->
var1;
1317 k1 = var0 > 1 ? cor0 *
flt16_even(a / var0) : 0;
1318 k2 = var1 > 1 ? cor1 *
flt16_even(a / var1) : 0;
1350 for (k = sce->
ics.
swb_offset[sfb]; k < sce->ics.swb_offset[sfb + 1]; k++) {
1375 float *out = sce->
coeffs;
1376 int global_gain, pulse_present = 0;
1385 if (!common_window && !scale_flag) {
1432 int g,
i, group, idx = 0;
1435 for (i = 0; i < ics->
max_sfb; i++, idx++) {
1438 for (group = 0; group < ics->
group_len[
g]; group++) {
1440 ch1 + group * 128 + offsets[i],
1441 offsets[i+1] - offsets[i]);
1463 int g, group,
i, idx = 0;
1467 for (i = 0; i < ics->
max_sfb;) {
1470 for (; i < bt_run_end; i++, idx++) {
1471 c = -1 + 2 * (sce1->
band_type[idx] - 14);
1473 c *= 1 - 2 * cpe->
ms_mask[idx];
1474 scale = c * sce1->
sf[idx];
1475 for (group = 0; group < ics->
group_len[
g]; group++)
1477 coef0 + group * 128 + offsets[i],
1479 offsets[i + 1] - offsets[i]);
1483 idx += bt_run_end -
i;
1499 int i, ret, common_window, ms_present = 0;
1502 if (common_window) {
1512 if (ms_present == 3) {
1515 }
else if (ms_present)
1518 if ((ret =
decode_ics(ac, &cpe->
ch[0], gb, common_window, 0)))
1520 if ((ret =
decode_ics(ac, &cpe->
ch[1], gb, common_window, 0)))
1523 if (common_window) {
1537 1.09050773266525765921,
1538 1.18920711500272106672,
1578 for (c = 0; c < num_gain; c++) {
1582 float gain_cache = 1.;
1585 gain = cge ?
get_vlc2(gb, vlc_scalefactors.
table, 7, 3) - 60: 0;
1586 gain_cache = powf(scale, -gain);
1589 coup->
gain[c][0] = gain_cache;
1592 for (sfb = 0; sfb < sce->
ics.
max_sfb; sfb++, idx++) {
1603 gain_cache = powf(scale, -t) * s;
1606 coup->
gain[c][idx] = gain_cache;
1624 int num_excl_chan = 0;
1627 for (i = 0; i < 7; i++)
1631 return num_excl_chan / 7;
1645 int drc_num_bands = 1;
1666 for (i = 0; i < drc_num_bands; i++) {
1679 for (i = 0; i < drc_num_bands; i++) {
1748 int bottom, top, order, start, end,
size, inc;
1754 for (filt = 0; filt < tns->
n_filt[w]; filt++) {
1757 order = tns->
order[w][filt];
1766 if ((size = end - start) <= 0)
1778 for (m = 0; m <
size; m++, start += inc)
1779 for (i = 1; i <=
FFMIN(m, order); i++)
1780 coef[start] -= coef[start - i * inc] * lpc[i - 1];
1783 for (m = 0; m <
size; m++, start += inc) {
1784 tmp[0] = coef[start];
1785 for (i = 1; i <=
FFMIN(m, order); i++)
1786 coef[start] += tmp[i] * lpc[i - 1];
1787 for (i = order; i > 0; i--)
1788 tmp[i] = tmp[i - 1];
1810 memset(in, 0, 448 *
sizeof(
float));
1817 memset(in + 1024 + 576, 0, 448 *
sizeof(
float));
1832 float *predTime = sce->
ret;
1834 int16_t num_samples = 2048;
1836 if (ltp->
lag < 1024)
1837 num_samples = ltp->
lag + 1024;
1838 for (i = 0; i < num_samples; i++)
1840 memset(&predTime[i], 0, (2048 - i) *
sizeof(
float));
1849 for (i = offsets[sfb]; i < offsets[sfb + 1]; i++)
1850 sce->
coeffs[i] += predFreq[i];
1860 float *saved = sce->
saved;
1861 float *saved_ltp = sce->
coeffs;
1867 memcpy(saved_ltp, saved, 512 *
sizeof(
float));
1868 memset(saved_ltp + 576, 0, 448 *
sizeof(
float));
1870 for (i = 0; i < 64; i++)
1871 saved_ltp[i + 512] = ac->
buf_mdct[1023 - i] * swindow[63 - i];
1873 memcpy(saved_ltp, ac->
buf_mdct + 512, 448 *
sizeof(
float));
1874 memset(saved_ltp + 576, 0, 448 *
sizeof(
float));
1876 for (i = 0; i < 64; i++)
1877 saved_ltp[i + 512] = ac->
buf_mdct[1023 - i] * swindow[63 - i];
1880 for (i = 0; i < 512; i++)
1881 saved_ltp[i + 512] = ac->
buf_mdct[1023 - i] * lwindow[511 - i];
1896 float *out = sce->
ret;
1897 float *saved = sce->
saved;
1902 float *temp = ac->
temp;
1907 for (i = 0; i < 1024; i += 128)
1922 memcpy( out, saved, 448 *
sizeof(
float));
1930 memcpy( out + 448 + 4*128, temp, 64 *
sizeof(
float));
1933 memcpy( out + 576, buf + 64, 448 *
sizeof(
float));
1939 memcpy( saved, temp + 64, 64 *
sizeof(
float));
1943 memcpy( saved + 448, buf + 7*128 + 64, 64 *
sizeof(
float));
1945 memcpy( saved, buf + 512, 448 *
sizeof(
float));
1946 memcpy( saved + 448, buf + 7*128 + 64, 64 *
sizeof(
float));
1948 memcpy( saved, buf + 512, 512 *
sizeof(
float));
1963 float *dest = target->
coeffs;
1964 const float *src = cce->
ch[0].
coeffs;
1965 int g,
i, group, k, idx = 0;
1968 "Dependent coupling is not supported together with LTP\n");
1972 for (i = 0; i < ics->
max_sfb; i++, idx++) {
1975 for (group = 0; group < ics->
group_len[
g]; group++) {
1976 for (k = offsets[i]; k < offsets[i + 1]; k++) {
1978 dest[group * 128 + k] += gain * src[group * 128 + k];
1999 const float *src = cce->
ch[0].
ret;
2000 float *dest = target->
ret;
2003 for (i = 0; i <
len; i++)
2004 dest[i] += gain * src[i];
2027 if (coup->
type[c] == type && coup->
id_select[c] == elem_id) {
2029 apply_coupling_method(ac, &cc->
ch[0], cce, index);
2034 apply_coupling_method(ac, &cc->
ch[1], cce, index++);
2048 for (type = 3; type >= 0; type--) {
2097 memset(new_che_pos, 0, 4 *
MAX_ELEM_ID *
sizeof(new_che_pos[0][0]));
2135 int samples = 0, multiplier, audio_found = 0;
2154 if (!(che=
get_che(ac, elem_type, elem_id))) {
2156 elem_type, elem_id);
2162 switch (elem_type) {
2189 memset(new_che_pos, 0, 4 *
MAX_ELEM_ID *
sizeof(new_che_pos[0][0]));
2194 "Not evaluating a further program_config_element as this construct is dubious at best.\n");
2218 elem_type_prev = elem_type;
2232 samples <<= multiplier;
2266 int *got_frame_ptr,
AVPacket *avpkt)
2269 const uint8_t *buf = avpkt->
data;
2270 int buf_size = avpkt->
size;
2275 int new_extradata_size;
2278 &new_extradata_size);
2280 if (new_extradata) {
2287 memcpy(avctx->
extradata, new_extradata, new_extradata_size);
2300 for (buf_offset = buf_consumed; buf_offset < buf_size; buf_offset++)
2301 if (buf[buf_offset])
2304 return buf_size > buf_offset ? buf_consumed : buf_size;
2313 for (type = 0; type < 4; type++) {
2314 if (ac->
che[type][i])
2327 #define LOAS_SYNC_WORD 0x2b7
2353 int sync_extension = 0;
2354 int bits_consumed, esize;
2362 if (config_start_bit % 8) {
2364 "config not byte aligned.\n", 1);
2370 gb->
buffer + (config_start_bit / 8),
2371 asclen, sync_extension);
2373 if (bits_consumed < 0)
2382 esize = (bits_consumed+7) / 8;
2397 return bits_consumed;
2403 int ret, audio_mux_version =
get_bits(gb, 1);
2406 if (audio_mux_version)
2411 if (audio_mux_version)
2419 "multiple programs are not supported\n", 1);
2428 "multiple layers are not supported\n", 1);
2433 if (!audio_mux_version) {
2464 if (audio_mux_version) {
2487 int mux_slot_length = 0;
2490 mux_slot_length += tmp;
2491 }
while (tmp == 255);
2492 return mux_slot_length;
2507 uint8_t use_same_mux =
get_bits(gb, 1);
2508 if (!use_same_mux) {
2513 "no decoder config found\n");
2521 }
else if (mux_slot_length_bytes * 8 + 256 <
get_bits_left(gb)) {
2523 "frame length mismatch %d << %d\n",
2533 int *got_frame_ptr,
AVPacket *avpkt)
2547 if (muxlength > avpkt->
size)
2568 "ADTS header detected, probably as result of configuration "