40 #define ENVELOPE_ADJUSTMENT_OFFSET 2
41 #define NOISE_FLOOR_OFFSET 6.0f
75 { 60, 60, 24, 24, 31, 31, 12, 12, 31, 12 };
78 #define SBR_INIT_VLC_STATIC(num, size) \
79 INIT_VLC_STATIC(&vlc_sbr[num], 9, sbr_tmp[num].table_size / sbr_tmp[num].elem_size, \
80 sbr_tmp[num].sbr_bits , 1, 1, \
81 sbr_tmp[num].sbr_codes, sbr_tmp[num].elem_size, sbr_tmp[num].elem_size, \
84 #define SBR_VLC_ROW(name) \
85 { name ## _codes, name ## _bits, sizeof(name ## _codes), sizeof(name ## _codes[0]) }
91 const void *sbr_codes, *sbr_bits;
92 const unsigned int table_size, elem_size;
118 for (n = 1; n < 320; n++)
123 for (n = 0; n < 320; n++)
132 sbr->
kx[0] = sbr->
kx[1] = 32;
153 return *(
const int16_t *)a - *(
const int16_t *)
b;
156 static inline int in_table_int16(
const int16_t *table,
int last_el, int16_t needle)
159 for (i = 0; i <= last_el; i++)
160 if (table[i] == needle)
170 static const float bands_warped[3] = { 1.32715174233856803909f,
171 1.18509277094158210129f,
172 1.11987160404675912501f };
173 const float lim_bands_per_octave_warped = bands_warped[sbr->
bs_limiter_bands - 1];
174 int16_t patch_borders[7];
177 patch_borders[0] = sbr->
kx[1];
184 memcpy(sbr->
f_tablelim + sbr->
n[0] + 1, patch_borders + 1,
185 (sbr->
num_patches - 1) *
sizeof(patch_borders[0]));
192 while (out < sbr->f_tablelim + sbr->
n_lim) {
193 if (*in >= *out * lim_bands_per_octave_warped) {
195 }
else if (*in == *out ||
216 uint8_t bs_header_extra_1;
217 uint8_t bs_header_extra_2;
235 if (bs_header_extra_1) {
249 if (bs_header_extra_2) {
269 int i,
min = array[0];
270 for (i = 1; i < nel; i++)
271 min =
FFMIN(array[i], min);
275 static void make_bands(int16_t* bands,
int start,
int stop,
int num_bands)
277 int k, previous, present;
280 base = powf((
float)stop / start, 1.0
f / num_bands);
284 for (k = 0; k < num_bands-1; k++) {
287 bands[k] = present - previous;
290 bands[num_bands-1] = stop - previous;
300 if (bs_xover_band >= n_master) {
302 "Invalid bitstream, crossover band index beyond array bounds: %d\n",
313 unsigned int temp, max_qmf_subbands;
314 unsigned int start_min, stop_min;
316 const int8_t *sbr_offset_ptr;
342 case 44100:
case 48000:
case 64000:
345 case 88200:
case 96000:
case 128000:
case 176400:
case 192000:
350 "Unsupported sample rate for SBR: %d\n", sbr->
sample_rate);
357 sbr->
k[2] = stop_min;
361 sbr->
k[2] += stop_dk[k];
363 sbr->
k[2] = 2*sbr->
k[0];
365 sbr->
k[2] = 3*sbr->
k[0];
371 sbr->
k[2] =
FFMIN(64, sbr->
k[2]);
375 max_qmf_subbands = 48;
377 max_qmf_subbands = 35;
379 max_qmf_subbands = 32;
381 if (sbr->
k[2] - sbr->
k[0] > max_qmf_subbands) {
383 "Invalid bitstream, too many QMF subbands: %d\n", sbr->
k[2] - sbr->
k[0]);
391 sbr->
n_master = ((sbr->
k[2] - sbr->
k[0] + (dk&2)) >> dk) << 1;
395 for (k = 1; k <= sbr->
n_master; k++)
398 k2diff = sbr->
k[2] - sbr->
k[0] - sbr->
n_master * dk;
407 for (k = 1; k <= sbr->
n_master; k++)
412 int two_regions, num_bands_0;
413 int vdk0_max, vdk1_min;
416 if (49 * sbr->
k[2] > 110 * sbr->
k[0]) {
418 sbr->
k[1] = 2 * sbr->
k[0];
421 sbr->
k[1] = sbr->
k[2];
424 num_bands_0 =
lrintf(half_bands *
log2f(sbr->
k[1] / (
float)sbr->
k[0])) * 2;
426 if (num_bands_0 <= 0) {
436 vdk0_max = vk0[num_bands_0];
439 for (k = 1; k <= num_bands_0; k++) {
449 float invwarp = spectrum->
bs_alter_scale ? 0.76923076923076923077f
451 int num_bands_1 =
lrintf(half_bands * invwarp *
452 log2f(sbr->
k[2] / (
float)sbr->
k[1])) * 2;
458 if (vdk1_min < vdk0_max) {
461 change =
FFMIN(vdk0_max - vk1[1], (vk1[num_bands_1] - vk1[1]) >> 1);
463 vk1[num_bands_1] -= change;
469 for (k = 1; k <= num_bands_1; k++) {
477 sbr->
n_master = num_bands_0 + num_bands_1;
481 (num_bands_0 + 1) *
sizeof(sbr->
f_master[0]));
482 memcpy(&sbr->
f_master[num_bands_0 + 1], vk1 + 1,
483 num_bands_1 *
sizeof(sbr->
f_master[0]));
501 int usb = sbr->
kx[1];
506 if (goal_sb < sbr->kx[1] + sbr->
m[1]) {
507 for (k = 0; sbr->
f_master[k] < goal_sb; k++) ;
513 for (i = k; i == k || sb > (sbr->
k[0] - 1 + msb - odd); i--) {
515 odd = (sb + sbr->
k[0]) & 1;
539 }
while (sb != sbr->
kx[1] + sbr->
m[1]);
553 sbr->
n[0] = (sbr->
n[1] + 1) >> 1;
556 (sbr->
n[1] + 1) *
sizeof(sbr->
f_master[0]));
561 if (sbr->
kx[1] + sbr->
m[1] > 64) {
563 "Stop frequency border too high: %d\n", sbr->
kx[1] + sbr->
m[1]);
566 if (sbr->
kx[1] > 32) {
572 temp = sbr->
n[1] & 1;
573 for (k = 1; k <= sbr->
n[0]; k++)
577 log2f(sbr->
k[2] / (
float)sbr->
kx[1])));
585 for (k = 1; k <= sbr->
n_q; k++) {
586 temp += (sbr->
n[0] - temp) / (sbr->
n_q + 1 - k);
605 for (i = 0; i < elements; i++) {
619 unsigned bs_pointer = 0;
621 int abs_bord_trail = 16;
622 int num_rel_lead, num_rel_trail;
623 unsigned bs_num_env_old = ch_data->
bs_num_env;
638 "Invalid bitstream, too many SBR envelopes in FIXFIX type SBR frame: %d\n",
643 ch_data->
t_env[0] = 0;
646 abs_bord_trail = (abs_bord_trail + (ch_data->
bs_num_env >> 1)) /
648 for (i = 0; i < num_rel_lead; i++)
649 ch_data->
t_env[i + 1] = ch_data->
t_env[i] + abs_bord_trail;
659 ch_data->
t_env[0] = 0;
662 for (i = 0; i < num_rel_trail; i++)
677 for (i = 0; i < num_rel_lead; i++)
689 ch_data->
bs_num_env = num_rel_lead + num_rel_trail + 1;
693 "Invalid bitstream, too many SBR envelopes in VARVAR type SBR frame: %d\n",
700 for (i = 0; i < num_rel_lead; i++)
702 for (i = 0; i < num_rel_trail; i++)
714 "Invalid bitstream, bs_pointer points to a middle noise border outside the time borders table: %d\n",
720 if (ch_data->
t_env[i-1] > ch_data->
t_env[i]) {
728 ch_data->
t_q[0] = ch_data->
t_env[0];
739 else if (bs_pointer == 1)
742 idx = bs_pointer - 1;
744 ch_data->
t_q[1] = ch_data->
t_env[idx];
747 ch_data->
e_a[0] = -(ch_data->
e_a[1] != bs_num_env_old);
748 ch_data->
e_a[1] = -1;
752 ch_data->
e_a[1] = bs_pointer - 1;
766 memcpy(dst->
t_q, src->
t_q,
sizeof(dst->
t_q));
771 dst->
e_a[1] = src->
e_a[1];
789 for (i = 0; i < sbr->
n_q; i++)
798 VLC_TYPE (*t_huff)[2], (*f_huff)[2];
801 const int odd = sbr->
n[1] & 1;
837 for (j = 0; j < sbr->
n[ch_data->
bs_freq_res[i + 1]]; j++)
840 for (j = 0; j < sbr->
n[ch_data->
bs_freq_res[i + 1]]; j++) {
845 for (j = 0; j < sbr->
n[ch_data->
bs_freq_res[i + 1]]; j++) {
846 k = j ? 2*j - odd : 0;
852 for (j = 1; j < sbr->
n[ch_data->
bs_freq_res[i + 1]]; j++)
866 VLC_TYPE (*t_huff)[2], (*f_huff)[2];
884 for (j = 0; j < sbr->
n_q; j++)
888 for (j = 1; j < sbr->
n_q; j++)
900 int bs_extension_id,
int *num_bits_left)
902 switch (bs_extension_id) {
905 av_log(ac->
avctx,
AV_LOG_ERROR,
"Parametric Stereo signaled to be not-present but was found in the bitstream.\n");
1005 "Invalid bitstream - cannot apply SBR to element type %d\n", id_aac);
1010 int num_bits_left =
get_bits(gb, 4);
1011 if (num_bits_left == 15)
1014 num_bits_left <<= 3;
1015 while (num_bits_left > 7) {
1019 if (num_bits_left < 0) {
1022 if (num_bits_left > 0)
1037 "SBR reset failed. Switching SBR to pure upsampling mode.\n");
1053 unsigned int num_sbr_bits = 0, num_align_bits;
1054 unsigned bytes_read;
1071 sbr->
kx[0] = sbr->
kx[1];
1072 sbr->
m[0] = sbr->
m[1];
1084 num_align_bits = ((cnt << 3) - 4 - num_sbr_bits) & 7;
1085 bytes_read = ((num_sbr_bits + num_align_bits + 4) >> 3);
1087 if (bytes_read > cnt) {
1089 "Expected to read %d SBR bytes actually read %d.\n", cnt, bytes_read);
1107 float fac = temp1 / (1.0f + temp2);
1113 for (k = 0; k < sbr->
n_q; k++) {
1116 float fac = temp1 / (1.0f + temp2);
1122 for (ch = 0; ch < (id_aac ==
TYPE_CPE) + 1; ch++) {
1129 for (k = 0; k < sbr->
n_q; k++)
1143 float z[320],
float W[2][32][32][2])
1146 memcpy(W[0], W[1],
sizeof(W[0]));
1147 memcpy(x , x+1024, (320-32)*
sizeof(x[0]));
1148 memcpy(x+288, in, 1024*
sizeof(x[0]));
1149 for (i = 0; i < 32; i++) {
1152 for (k = 0; k < 64; k++) {
1153 float f = z[k] + z[k + 64] + z[k + 128] + z[k + 192] + z[k + 256];
1158 for (k = 1; k < 32; k++) {
1159 z[64+2*k-1] = z[ k];
1160 z[64+2*k ] = -z[64-k];
1165 for (k = 0; k < 32; k++) {
1166 W[1][i][k][0] = -z[63-k];
1167 W[1][i][k][1] = z[k];
1178 float *out,
float X[2][38][64],
1179 float mdct_buf[2][64],
1180 float *
v0,
int *v_off,
const unsigned int div)
1184 const int step = 128 >> div;
1186 for (i = 0; i < 32; i++) {
1187 if (*v_off < step) {
1188 int saved_samples = (1280 - 128) >> div;
1196 for (n = 0; n < 32; n++) {
1197 X[0][i][ n] = -X[0][i][n];
1198 X[0][i][32+n] = X[1][i][31-n];
1200 mdct->
imdct_half(mdct, mdct_buf[0], X[0][i]);
1201 for (n = 0; n < 32; n++) {
1202 v[ n] = mdct_buf[0][63 - 2*n];
1203 v[63 - n] = -mdct_buf[0][62 - 2*n];
1206 for (n = 1; n < 64; n+=2) {
1207 X[1][i][n] = -X[1][i][n];
1209 mdct->
imdct_half(mdct, mdct_buf[0], X[0][i]);
1210 mdct->
imdct_half(mdct, mdct_buf[1], X[1][i]);
1211 for (n = 0; n < 64; n++) {
1212 v[ n] = -mdct_buf[0][63 - n] + mdct_buf[1][ n ];
1213 v[127 - n] = mdct_buf[0][63 - n] + mdct_buf[1][ n ];
1217 dsp->
vector_fmul_add(out, v + ( 192 >> div), sbr_qmf_window + ( 64 >> div), out , 64 >> div);
1218 dsp->
vector_fmul_add(out, v + ( 256 >> div), sbr_qmf_window + (128 >> div), out , 64 >> div);
1219 dsp->
vector_fmul_add(out, v + ( 448 >> div), sbr_qmf_window + (192 >> div), out , 64 >> div);
1220 dsp->
vector_fmul_add(out, v + ( 512 >> div), sbr_qmf_window + (256 >> div), out , 64 >> div);
1221 dsp->
vector_fmul_add(out, v + ( 704 >> div), sbr_qmf_window + (320 >> div), out , 64 >> div);
1222 dsp->
vector_fmul_add(out, v + ( 768 >> div), sbr_qmf_window + (384 >> div), out , 64 >> div);
1223 dsp->
vector_fmul_add(out, v + ( 960 >> div), sbr_qmf_window + (448 >> div), out , 64 >> div);
1224 dsp->
vector_fmul_add(out, v + (1024 >> div), sbr_qmf_window + (512 >> div), out , 64 >> div);
1225 dsp->
vector_fmul_add(out, v + (1216 >> div), sbr_qmf_window + (576 >> div), out , 64 >> div);
1233 float real_sum = 0.0f;
1234 float imag_sum = 0.0f;
1236 for (i = 1; i < 38; i++) {
1237 real_sum += x[i][0] * x[i+lag][0] + x[i][1] * x[i+lag][1];
1238 imag_sum += x[i][0] * x[i+lag][1] - x[i][1] * x[i+lag][0];
1240 phi[2-lag][1][0] = real_sum + x[ 0][0] * x[lag][0] + x[ 0][1] * x[lag][1];
1241 phi[2-lag][1][1] = imag_sum + x[ 0][0] * x[lag][1] - x[ 0][1] * x[lag][0];
1243 phi[0][0][0] = real_sum + x[38][0] * x[39][0] + x[38][1] * x[39][1];
1244 phi[0][0][1] = imag_sum + x[38][0] * x[39][1] - x[38][1] * x[39][0];
1247 for (i = 1; i < 38; i++) {
1248 real_sum += x[i][0] * x[i][0] + x[i][1] * x[i][1];
1250 phi[2][1][0] = real_sum + x[ 0][0] * x[ 0][0] + x[ 0][1] * x[ 0][1];
1251 phi[1][0][0] = real_sum + x[38][0] * x[38][0] + x[38][1] * x[38][1];
1260 const float X_low[32][40][2],
int k0)
1263 for (k = 0; k < k0; k++) {
1264 float phi[3][2][2], dk;
1270 dk = phi[2][1][0] * phi[1][0][0] -
1271 (phi[1][1][0] * phi[1][1][0] + phi[1][1][1] * phi[1][1][1]) / 1.000001
f;
1277 float temp_real, temp_im;
1278 temp_real = phi[0][0][0] * phi[1][1][0] -
1279 phi[0][0][1] * phi[1][1][1] -
1280 phi[0][1][0] * phi[1][0][0];
1281 temp_im = phi[0][0][0] * phi[1][1][1] +
1282 phi[0][0][1] * phi[1][1][0] -
1283 phi[0][1][1] * phi[1][0][0];
1285 alpha1[k][0] = temp_real / dk;
1286 alpha1[k][1] = temp_im / dk;
1289 if (!phi[1][0][0]) {
1293 float temp_real, temp_im;
1294 temp_real = phi[0][0][0] + alpha1[k][0] * phi[1][1][0] +
1295 alpha1[k][1] * phi[1][1][1];
1296 temp_im = phi[0][0][1] + alpha1[k][1] * phi[1][1][0] -
1297 alpha1[k][0] * phi[1][1][1];
1299 alpha0[k][0] = -temp_real / phi[1][0][0];
1300 alpha0[k][1] = -temp_im / phi[1][0][0];
1303 if (alpha1[k][0] * alpha1[k][0] + alpha1[k][1] * alpha1[k][1] >= 16.0
f ||
1304 alpha0[k][0] * alpha0[k][0] + alpha0[k][1] * alpha0[k][1] >= 16.0
f) {
1318 static const float bw_tab[] = { 0.0f, 0.75f, 0.9f, 0.98f };
1320 for (i = 0; i < sbr->
n_q; i++) {
1326 if (new_bw < ch_data->bw_array[i]) {
1327 new_bw = 0.75f * new_bw + 0.25f * ch_data->
bw_array[i];
1329 new_bw = 0.90625f * new_bw + 0.09375f * ch_data->
bw_array[i];
1330 ch_data->
bw_array[i] = new_bw < 0.015625f ? 0.0f : new_bw;
1336 float X_low[32][40][2],
const float W[2][32][32][2])
1339 const int t_HFGen = 8;
1341 memset(X_low, 0, 32*
sizeof(*X_low));
1342 for (k = 0; k < sbr->
kx[1]; k++) {
1343 for (i = t_HFGen; i < i_f + t_HFGen; i++) {
1344 X_low[k][i][0] = W[1][i - t_HFGen][k][0];
1345 X_low[k][i][1] = W[1][i - t_HFGen][k][1];
1348 for (k = 0; k < sbr->
kx[0]; k++) {
1349 for (i = 0; i < t_HFGen; i++) {
1350 X_low[k][i][0] = W[0][i + i_f - t_HFGen][k][0];
1351 X_low[k][i][1] = W[0][i + i_f - t_HFGen][k][1];
1359 float X_high[64][40][2],
const float X_low[32][40][2],
1360 const float (*alpha0)[2],
const float (*alpha1)[2],
1361 const float bw_array[5],
const uint8_t *t_env,
1377 "ERROR : no subband found for frequency %d\n", k);
1381 alpha[0] = alpha1[p][0] * bw_array[
g] * bw_array[
g];
1382 alpha[1] = alpha1[p][1] * bw_array[
g] * bw_array[
g];
1383 alpha[2] = alpha0[p][0] * bw_array[
g];
1384 alpha[3] = alpha0[p][1] * bw_array[
g];
1386 for (i = 2 * t_env[0]; i < 2 * t_env[bs_num_env]; i++) {
1389 X_low[p][idx - 2][0] * alpha[0] -
1390 X_low[p][idx - 2][1] * alpha[1] +
1391 X_low[p][idx - 1][0] * alpha[2] -
1392 X_low[p][idx - 1][1] * alpha[3] +
1395 X_low[p][idx - 2][1] * alpha[0] +
1396 X_low[p][idx - 2][0] * alpha[1] +
1397 X_low[p][idx - 1][1] * alpha[2] +
1398 X_low[p][idx - 1][0] * alpha[3] +
1403 if (k < sbr->m[1] + sbr->
kx[1])
1404 memset(X_high + k, 0, (sbr->
m[1] + sbr->
kx[1] - k) *
sizeof(*X_high));
1411 const float X_low[32][40][2],
const float Y[2][38][64][2],
1417 memset(X, 0, 2*
sizeof(*X));
1418 for (k = 0; k < sbr->
kx[0]; k++) {
1419 for (i = 0; i < i_Temp; i++) {
1424 for (; k < sbr->
kx[0] + sbr->
m[0]; k++) {
1425 for (i = 0; i < i_Temp; i++) {
1426 X[0][i][k] = Y[0][i + i_f][k][0];
1427 X[1][i][k] = Y[0][i + i_f][k][1];
1431 for (k = 0; k < sbr->
kx[1]; k++) {
1432 for (i = i_Temp; i < 38; i++) {
1437 for (; k < sbr->
kx[1] + sbr->
m[1]; k++) {
1438 for (i = i_Temp; i < i_f; i++) {
1439 X[0][i][k] = Y[1][i][k][0];
1440 X[1][i][k] = Y[1][i][k][1];
1456 const unsigned int ilim = sbr->
n[ch_data->
bs_freq_res[e + 1]];
1460 for (i = 0; i < ilim; i++)
1461 for (m = table[i]; m < table[i + 1]; m++)
1466 for (i = 0; i < sbr->
n_q; i++)
1467 for (m = sbr->
f_tablenoise[i]; m < sbr->f_tablenoise[i + 1]; m++)
1470 for (i = 0; i < sbr->
n[1]; i++) {
1472 const unsigned int m_midpoint =
1476 (e >= e_a[1] || (ch_data->
s_indexmapped[0][m_midpoint - sbr->
kx[1]] == 1));
1480 for (i = 0; i < ilim; i++) {
1481 int additional_sinusoid_present = 0;
1482 for (m = table[i]; m < table[i + 1]; m++) {
1484 additional_sinusoid_present = 1;
1488 memset(&sbr->
s_mapped[e][table[i] - sbr->
kx[1]], additional_sinusoid_present,
1489 (table[i + 1] - table[i]) *
sizeof(sbr->
s_mapped[e][0]));
1504 const float recip_env_size = 0.5f / (ch_data->
t_env[e + 1] - ch_data->
t_env[e]);
1508 for (m = 0; m < sbr->
m[1]; m++) {
1511 for (i = ilb; i < iub; i++) {
1512 sum += X_high[m + sbr->
kx[1]][i][0] * X_high[m + sbr->
kx[1]][i][0] +
1513 X_high[m + sbr->
kx[1]][i][1] * X_high[m + sbr->
kx[1]][i][1];
1515 e_curr[e][m] = sum * recip_env_size;
1522 const int env_size = 2 * (ch_data->
t_env[e + 1] - ch_data->
t_env[e]);
1527 for (p = 0; p < sbr->
n[ch_data->
bs_freq_res[e + 1]]; p++) {
1529 const int den = env_size * (table[p + 1] - table[p]);
1531 for (k = table[p]; k < table[p + 1]; k++) {
1532 for (i = ilb; i < iub; i++) {
1533 sum += X_high[k][i][0] * X_high[k][i][0] +
1534 X_high[k][i][1] * X_high[k][i][1];
1538 for (k = table[p]; k < table[p + 1]; k++) {
1539 e_curr[e][k - sbr->
kx[1]] = sum;
1551 SBRData *ch_data,
const int e_a[2])
1555 static const float limgain[4] = { 0.70795, 1.0, 1.41254, 10000000000 };
1558 int delta = !((e == e_a[1]) || (e == e_a[0]));
1559 for (k = 0; k < sbr->
n_lim; k++) {
1560 float gain_boost, gain_max;
1561 float sum[2] = { 0.0f, 0.0f };
1562 for (m = sbr->
f_tablelim[k] - sbr->
kx[1]; m < sbr->f_tablelim[k + 1] - sbr->
kx[1]; m++) {
1564 sbr->
q_m[e][m] = sqrtf(temp * sbr->
q_mapped[e][m]);
1568 ((1.0f + sbr->
e_curr[e][m]) *
1569 (1.0f + sbr->
q_mapped[e][m] * delta)));
1572 ((1.0f + sbr->
e_curr[e][m]) *
1576 for (m = sbr->
f_tablelim[k] - sbr->
kx[1]; m < sbr->f_tablelim[k + 1] - sbr->
kx[1]; m++) {
1578 sum[1] += sbr->
e_curr[e][m];
1580 gain_max = limgain[sbr->
bs_limiter_gains] * sqrtf((FLT_EPSILON + sum[0]) / (FLT_EPSILON + sum[1]));
1581 gain_max =
FFMIN(100000.
f, gain_max);
1582 for (m = sbr->
f_tablelim[k] - sbr->
kx[1]; m < sbr->f_tablelim[k + 1] - sbr->
kx[1]; m++) {
1583 float q_m_max = sbr->
q_m[e][m] * gain_max / sbr->
gain[e][m];
1587 sum[0] = sum[1] = 0.0f;
1588 for (m = sbr->
f_tablelim[k] - sbr->
kx[1]; m < sbr->f_tablelim[k + 1] - sbr->
kx[1]; m++) {
1591 + sbr->
s_m[e][m] * sbr->
s_m[e][m]
1592 + (delta && !sbr->
s_m[e][m]) * sbr->
q_m[e][m] * sbr->
q_m[e][m];
1594 gain_boost = sqrtf((FLT_EPSILON + sum[0]) / (FLT_EPSILON + sum[1]));
1595 gain_boost =
FFMIN(1.584893192
f, gain_boost);
1596 for (m = sbr->
f_tablelim[k] - sbr->
kx[1]; m < sbr->f_tablelim[k + 1] - sbr->
kx[1]; m++) {
1597 sbr->
gain[e][m] *= gain_boost;
1598 sbr->
q_m[e][m] *= gain_boost;
1599 sbr->
s_m[e][m] *= gain_boost;
1612 const int kx = sbr->
kx[1];
1613 const int m_max = sbr->
m[1];
1614 static const float h_smooth[5] = {
1621 static const int8_t phi[2][4] = {
1625 float (*g_temp)[48] = ch_data->
g_temp, (*q_temp)[48] = ch_data->
q_temp;
1628 memcpy(Y[0], Y[1],
sizeof(Y[0]));
1631 for (i = 0; i < h_SL; i++) {
1632 memcpy(g_temp[i + 2*ch_data->
t_env[0]], sbr->
gain[0], m_max *
sizeof(sbr->
gain[0][0]));
1633 memcpy(q_temp[i + 2*ch_data->
t_env[0]], sbr->
q_m[0], m_max *
sizeof(sbr->
q_m[0][0]));
1641 for (i = 2 * ch_data->
t_env[e]; i < 2 * ch_data->t_env[e + 1]; i++) {
1642 memcpy(g_temp[h_SL + i], sbr->
gain[e], m_max *
sizeof(sbr->
gain[0][0]));
1643 memcpy(q_temp[h_SL + i], sbr->
q_m[e], m_max *
sizeof(sbr->
q_m[0][0]));
1648 for (i = 2 * ch_data->
t_env[e]; i < 2 * ch_data->t_env[e + 1]; i++) {
1649 int phi_sign = (1 - 2*(kx & 1));
1651 if (h_SL && e != e_a[0] && e != e_a[1]) {
1652 for (m = 0; m < m_max; m++) {
1653 const int idx1 = i + h_SL;
1654 float g_filt = 0.0f;
1655 for (j = 0; j <= h_SL; j++)
1656 g_filt += g_temp[idx1 - j][m] * h_smooth[j];
1657 Y[1][i][m + kx][0] =
1659 Y[1][i][m + kx][1] =
1663 for (m = 0; m < m_max; m++) {
1664 const float g_filt = g_temp[i + h_SL][m];
1665 Y[1][i][m + kx][0] =
1667 Y[1][i][m + kx][1] =
1672 if (e != e_a[0] && e != e_a[1]) {
1673 for (m = 0; m < m_max; m++) {
1674 indexnoise = (indexnoise + 1) & 0x1ff;
1675 if (sbr->
s_m[e][m]) {
1676 Y[1][i][m + kx][0] +=
1677 sbr->
s_m[e][m] * phi[0][indexsine];
1678 Y[1][i][m + kx][1] +=
1679 sbr->
s_m[e][m] * (phi[1][indexsine] * phi_sign);
1683 const int idx1 = i + h_SL;
1685 for (j = 0; j <= h_SL; j++)
1686 q_filt += q_temp[idx1 - j][m] * h_smooth[j];
1688 q_filt = q_temp[i][m];
1690 Y[1][i][m + kx][0] +=
1692 Y[1][i][m + kx][1] +=
1693 q_filt * sbr_noise_table[indexnoise][1];
1695 phi_sign = -phi_sign;
1698 indexnoise = (indexnoise + m_max) & 0x1ff;
1699 for (m = 0; m < m_max; m++) {
1700 Y[1][i][m + kx][0] +=
1701 sbr->
s_m[e][m] * phi[0][indexsine];
1702 Y[1][i][m + kx][1] +=
1703 sbr->
s_m[e][m] * (phi[1][indexsine] * phi_sign);
1704 phi_sign = -phi_sign;
1707 indexsine = (indexsine + 1) & 3;
1719 int nch = (id_aac ==
TYPE_CPE) ? 2 : 1;
1724 for (ch = 0; ch < nch; ch++) {
1753 memcpy(sbr->
X[1], sbr->
X[0],
sizeof(sbr->
X[0]));