Lines Matching defs:sbr
30 #include "sbr.h"
139 static void sbr_turnoff(SpectralBandReplication *sbr) {
140 sbr->start = 0;
142 sbr->kx[1] = 32; //Typo in spec, kx' inits to 32
143 sbr->m[1] = 0;
145 sbr->data[0].e_a[1] = sbr->data[1].e_a[1] = -1;
146 memset(&sbr->spectrum_params, -1, sizeof(SpectrumParameters));
149 av_cold void ff_aac_sbr_ctx_init(AACContext *ac, SpectralBandReplication *sbr)
151 if(sbr->mdct.mdct_bits)
153 sbr->kx[0] = sbr->kx[1];
154 sbr_turnoff(sbr);
155 sbr->data[0].synthesis_filterbank_samples_offset = SBR_SYNTHESIS_BUF_SIZE - (1280 - 128);
156 sbr->data[1].synthesis_filterbank_samples_offset = SBR_SYNTHESIS_BUF_SIZE - (1280 - 128);
160 ff_mdct_init(&sbr->mdct, 7, 1, 1.0 / (64 * 32768.0));
161 ff_mdct_init(&sbr->mdct_ana, 7, 1, -2.0 * 32768.0);
162 ff_ps_ctx_init(&sbr->ps);
163 ff_sbrdsp_init(&sbr->dsp);
164 aacsbr_func_ptr_init(&sbr->c);
167 av_cold void ff_aac_sbr_ctx_close(SpectralBandReplication *sbr)
169 ff_mdct_end(&sbr->mdct);
170 ff_mdct_end(&sbr->mdct_ana);
188 static void sbr_make_f_tablelim(SpectralBandReplication *sbr)
191 if (sbr->bs_limiter_bands > 0) {
195 const float lim_bands_per_octave_warped = bands_warped[sbr->bs_limiter_bands - 1];
197 uint16_t *in = sbr->f_tablelim + 1, *out = sbr->f_tablelim;
199 patch_borders[0] = sbr->kx[1];
200 for (k = 1; k <= sbr->num_patches; k++)
201 patch_borders[k] = patch_borders[k-1] + sbr->patch_num_subbands[k-1];
203 memcpy(sbr->f_tablelim, sbr->f_tablelow,
204 (sbr->n[0] + 1) * sizeof(sbr->f_tablelow[0]));
205 if (sbr->num_patches > 1)
206 memcpy(sbr->f_tablelim + sbr->n[0] + 1, patch_borders + 1,
207 (sbr->num_patches - 1) * sizeof(patch_borders[0]));
209 qsort(sbr->f_tablelim, sbr->num_patches + sbr->n[0],
210 sizeof(sbr->f_tablelim[0]),
213 sbr->n_lim = sbr->n[0] + sbr->num_patches - 1;
214 while (out < sbr->f_tablelim + sbr->n_lim) {
218 !in_table_int16(patch_borders, sbr->num_patches, *in)) {
220 sbr->n_lim--;
221 } else if (!in_table_int16(patch_borders, sbr->num_patches, *out)) {
223 sbr->n_lim--;
229 sbr->f_tablelim[0] = sbr->f_tablelow[0];
230 sbr->f_tablelim[1] = sbr->f_tablelow[sbr->n[0]];
231 sbr->n_lim = 1;
235 static unsigned int read_sbr_header(SpectralBandReplication *sbr, GetBitContext *gb)
240 int old_bs_limiter_bands = sbr->bs_limiter_bands;
243 sbr->start = 1;
246 memcpy(&old_spectrum_params, &sbr->spectrum_params, sizeof(SpectrumParameters));
248 sbr->bs_amp_res_header = get_bits1(gb);
249 sbr->spectrum_params.bs_start_freq = get_bits(gb, 4);
250 sbr->spectrum_params.bs_stop_freq = get_bits(gb, 4);
251 sbr->spectrum_params.bs_xover_band = get_bits(gb, 3);
258 sbr->spectrum_params.bs_freq_scale = get_bits(gb, 2);
259 sbr->spectrum_params.bs_alter_scale = get_bits1(gb);
260 sbr->spectrum_params.bs_noise_bands = get_bits(gb, 2);
262 sbr->spectrum_params.bs_freq_scale = 2;
263 sbr->spectrum_params.bs_alter_scale = 1;
264 sbr->spectrum_params.bs_noise_bands = 2;
268 if (memcmp(&old_spectrum_params, &sbr->spectrum_params, sizeof(SpectrumParameters)))
269 sbr->reset = 1;
272 sbr->bs_limiter_bands = get_bits(gb, 2);
273 sbr->bs_limiter_gains = get_bits(gb, 2);
274 sbr->bs_interpol_freq = get_bits1(gb);
275 sbr->bs_smoothing_mode = get_bits1(gb);
277 sbr->bs_limiter_bands = 2;
278 sbr->bs_limiter_gains = 2;
279 sbr->bs_interpol_freq = 1;
280 sbr->bs_smoothing_mode = 1;
283 if (sbr->bs_limiter_bands != old_bs_limiter_bands && !sbr->reset)
284 sbr_make_f_tablelim(sbr);
332 static int sbr_make_f_master(AACContext *ac, SpectralBandReplication *sbr,
341 if (sbr->sample_rate < 32000) {
343 } else if (sbr->sample_rate < 64000) {
348 switch (sbr->sample_rate) {
369 "Unsupported sample rate for SBR: %d\n", sbr->sample_rate);
373 start_min = ((temp << 7) + (sbr->sample_rate >> 1)) / sbr->sample_rate;
374 stop_min = ((temp << 8) + (sbr->sample_rate >> 1)) / sbr->sample_rate;
376 sbr->k[0] = start_min + sbr_offset_ptr[spectrum->bs_start_freq];
379 sbr->k[2] = stop_min;
383 sbr->k[2] += stop_dk[k];
385 sbr->k[2] = 2*sbr->k[0];
387 sbr->k[2] = 3*sbr->k[0];
393 sbr->k[2] = FFMIN(64, sbr->k[2]);
396 if (sbr->sample_rate <= 32000) {
398 } else if (sbr->sample_rate == 44100) {
400 } else if (sbr->sample_rate >= 48000)
405 if (sbr->k[2] - sbr->k[0] > max_qmf_subbands) {
407 "Invalid bitstream, too many QMF subbands: %d\n", sbr->k[2] - sbr->k[0]);
415 sbr->n_master = ((sbr->k[2] - sbr->k[0] + (dk&2)) >> dk) << 1;
416 if (check_n_master(ac->avctx, sbr->n_master, sbr->spectrum_params.bs_xover_band))
419 for (k = 1; k <= sbr->n_master; k++)
420 sbr->f_master[k] = dk;
422 k2diff = sbr->k[2] - sbr->k[0] - sbr->n_master * dk;
424 sbr->f_master[1]--;
425 sbr->f_master[2]-= (k2diff < -1);
427 sbr->f_master[sbr->n_master]++;
430 sbr->f_master[0] = sbr->k[0];
431 for (k = 1; k <= sbr->n_master; k++)
432 sbr->f_master[k] += sbr->f_master[k - 1];
440 if (49 * sbr->k[2] > 110 * sbr->k[0]) {
442 sbr->k[1] = 2 * sbr->k[0];
445 sbr->k[1] = sbr->k[2];
448 num_bands_0 = lrintf(half_bands * log2f(sbr->k[1] / (float)sbr->k[0])) * 2;
457 make_bands(vk0+1, sbr->k[0], sbr->k[1], num_bands_0);
462 vk0[0] = sbr->k[0];
476 log2f(sbr->k[2] / (float)sbr->k[1])) * 2;
478 make_bands(vk1+1, sbr->k[1], sbr->k[2], num_bands_1);
492 vk1[0] = sbr->k[1];
501 sbr->n_master = num_bands_0 + num_bands_1;
502 if (check_n_master(ac->avctx, sbr->n_master, sbr->spectrum_params.bs_xover_band))
504 memcpy(&sbr->f_master[0], vk0,
505 (num_bands_0 + 1) * sizeof(sbr->f_master[0]));
506 memcpy(&sbr->f_master[num_bands_0 + 1], vk1 + 1,
507 num_bands_1 * sizeof(sbr->f_master[0]));
510 sbr->n_master = num_bands_0;
511 if (check_n_master(ac->avctx, sbr->n_master, sbr->spectrum_params.bs_xover_band))
513 memcpy(sbr->f_master, vk0, (num_bands_0 + 1) * sizeof(sbr->f_master[0]));
521 static int sbr_hf_calc_npatches(AACContext *ac, SpectralBandReplication *sbr)
524 int msb = sbr->k[0];
525 int usb = sbr->kx[1];
526 int goal_sb = ((1000 << 11) + (sbr->sample_rate >> 1)) / sbr->sample_rate;
528 sbr->num_patches = 0;
530 if (goal_sb < sbr->kx[1] + sbr->m[1]) {
531 for (k = 0; sbr->f_master[k] < goal_sb; k++) ;
533 k = sbr->n_master;
537 for (i = k; i == k || sb > (sbr->k[0] - 1 + msb - odd); i--) {
538 sb = sbr->f_master[i];
539 odd = (sb + sbr->k[0]) & 1;
546 if (sbr->num_patches > 5) {
547 av_log(ac->avctx, AV_LOG_ERROR, "Too many patches: %d\n", sbr->num_patches);
551 sbr->patch_num_subbands[sbr->num_patches] = FFMAX(sb - usb, 0);
552 sbr->patch_start_subband[sbr->num_patches] = sbr->k[0] - odd - sbr->patch_num_subbands[sbr->num_patches];
554 if (sbr->patch_num_subbands[sbr->num_patches] > 0) {
557 sbr->num_patches++;
559 msb = sbr->kx[1];
561 if (sbr->f_master[k] - sb < 3)
562 k = sbr->n_master;
563 } while (sb != sbr->kx[1] + sbr->m[1]);
565 if (sbr->num_patches > 1 && sbr->patch_num_subbands[sbr->num_patches-1] < 3)
566 sbr->num_patches--;
572 static int sbr_make_f_derived(AACContext *ac, SpectralBandReplication *sbr)
576 sbr->n[1] = sbr->n_master - sbr->spectrum_params.bs_xover_band;
577 sbr->n[0] = (sbr->n[1] + 1) >> 1;
579 memcpy(sbr->f_tablehigh, &sbr->f_master[sbr->spectrum_params.bs_xover_band],
580 (sbr->n[1] + 1) * sizeof(sbr->f_master[0]));
581 sbr->m[1] = sbr->f_tablehigh[sbr->n[1]] - sbr->f_tablehigh[0];
582 sbr->kx[1] = sbr->f_tablehigh[0];
585 if (sbr->kx[1] + sbr->m[1] > 64) {
587 "Stop frequency border too high: %d\n", sbr->kx[1] + sbr->m[1]);
590 if (sbr->kx[1] > 32) {
591 av_log(ac->avctx, AV_LOG_ERROR, "Start frequency border too high: %d\n", sbr->kx[1]);
595 sbr->f_tablelow[0] = sbr->f_tablehigh[0];
596 temp = sbr->n[1] & 1;
597 for (k = 1; k <= sbr->n[0]; k++)
598 sbr->f_tablelow[k] = sbr->f_tablehigh[2 * k - temp];
600 sbr->n_q = FFMAX(1, lrintf(sbr->spectrum_params.bs_noise_bands *
601 log2f(sbr->k[2] / (float)sbr->kx[1]))); // 0 <= bs_noise_bands <= 3
602 if (sbr->n_q > 5) {
603 av_log(ac->avctx, AV_LOG_ERROR, "Too many noise floor scale factors: %d\n", sbr->n_q);
607 sbr->f_tablenoise[0] = sbr->f_tablelow[0];
609 for (k = 1; k <= sbr->n_q; k++) {
610 temp += (sbr->n[0] - temp) / (sbr->n_q + 1 - k);
611 sbr->f_tablenoise[k] = sbr->f_tablelow[temp];
614 if (sbr_hf_calc_npatches(ac, sbr) < 0)
617 sbr_make_f_tablelim(sbr);
619 sbr->data[0].f_indexnoise = 0;
620 sbr->data[1].f_indexnoise = 0;
639 static int read_sbr_grid(AACContext *ac, SpectralBandReplication *sbr,
650 ch_data->bs_amp_res = sbr->bs_amp_res_header;
799 static void read_sbr_dtdf(SpectralBandReplication *sbr, GetBitContext *gb,
807 static void read_sbr_invf(SpectralBandReplication *sbr, GetBitContext *gb,
813 for (i = 0; i < sbr->n_q; i++)
817 static void read_sbr_envelope(SpectralBandReplication *sbr, GetBitContext *gb,
824 const int delta = (ch == 1 && sbr->bs_coupling == 1) + 1;
825 const int odd = sbr->n[1] & 1;
827 if (sbr->bs_coupling && ch) {
861 for (j = 0; j < sbr->n[ch_data->bs_freq_res[i + 1]]; j++)
864 for (j = 0; j < sbr->n[ch_data->bs_freq_res[i + 1]]; j++) {
869 for (j = 0; j < sbr->n[ch_data->bs_freq_res[i + 1]]; j++) {
876 for (j = 1; j < sbr->n[ch_data->bs_freq_res[i + 1]]; j++)
886 static void read_sbr_noise(SpectralBandReplication *sbr, GetBitContext *gb,
892 int delta = (ch == 1 && sbr->bs_coupling == 1) + 1;
894 if (sbr->bs_coupling && ch) {
908 for (j = 0; j < sbr->n_q; j++)
912 for (j = 1; j < sbr->n_q; j++)
922 static void read_sbr_extension(AACContext *ac, SpectralBandReplication *sbr,
934 *num_bits_left -= ff_ps_read_data(ac->avctx, gb, &sbr->ps, *num_bits_left);
954 SpectralBandReplication *sbr,
960 if (read_sbr_grid(ac, sbr, gb, &sbr->data[0]))
962 read_sbr_dtdf(sbr, gb, &sbr->data[0]);
963 read_sbr_invf(sbr, gb, &sbr->data[0]);
964 read_sbr_envelope(sbr, gb, &sbr->data[0], 0);
965 read_sbr_noise(sbr, gb, &sbr->data[0], 0);
967 if ((sbr->data[0].bs_add_harmonic_flag = get_bits1(gb)))
968 get_bits1_vector(gb, sbr->data[0].bs_add_harmonic, sbr->n[1]);
974 SpectralBandReplication *sbr,
980 if ((sbr->bs_coupling = get_bits1(gb))) {
981 if (read_sbr_grid(ac, sbr, gb, &sbr->data[0]))
983 copy_sbr_grid(&sbr->data[1], &sbr->data[0]);
984 read_sbr_dtdf(sbr, gb, &sbr->data[0]);
985 read_sbr_dtdf(sbr, gb, &sbr->data[1]);
986 read_sbr_invf(sbr, gb, &sbr->data[0]);
987 memcpy(sbr->data[1].bs_invf_mode[1], sbr->data[1].bs_invf_mode[0], sizeof(sbr->data[1].bs_invf_mode[0]));
988 memcpy(sbr->data[1].bs_invf_mode[0], sbr->data[0].bs_invf_mode[0], sizeof(sbr->data[1].bs_invf_mode[0]));
989 read_sbr_envelope(sbr, gb, &sbr->data[0], 0);
990 read_sbr_noise(sbr, gb, &sbr->data[0], 0);
991 read_sbr_envelope(sbr, gb, &sbr->data[1], 1);
992 read_sbr_noise(sbr, gb, &sbr->data[1], 1);
994 if (read_sbr_grid(ac, sbr, gb, &sbr->data[0]) ||
995 read_sbr_grid(ac, sbr, gb, &sbr->data[1]))
997 read_sbr_dtdf(sbr, gb, &sbr->data[0]);
998 read_sbr_dtdf(sbr, gb, &sbr->data[1]);
999 read_sbr_invf(sbr, gb, &sbr->data[0]);
1000 read_sbr_invf(sbr, gb, &sbr->data[1]);
1001 read_sbr_envelope(sbr, gb, &sbr->data[0], 0);
1002 read_sbr_envelope(sbr, gb, &sbr->data[1], 1);
1003 read_sbr_noise(sbr, gb, &sbr->data[0], 0);
1004 read_sbr_noise(sbr, gb, &sbr->data[1], 1);
1007 if ((sbr->data[0].bs_add_harmonic_flag = get_bits1(gb)))
1008 get_bits1_vector(gb, sbr->data[0].bs_add_harmonic, sbr->n[1]);
1009 if ((sbr->data[1].bs_add_harmonic_flag = get_bits1(gb)))
1010 get_bits1_vector(gb, sbr->data[1].bs_add_harmonic, sbr->n[1]);
1015 static unsigned int read_sbr_data(AACContext *ac, SpectralBandReplication *sbr,
1021 if (read_sbr_single_channel_element(ac, sbr, gb)) {
1022 sbr_turnoff(sbr);
1026 if (read_sbr_channel_pair_element(ac, sbr, gb)) {
1027 sbr_turnoff(sbr);
1033 sbr_turnoff(sbr);
1044 read_sbr_extension(ac, sbr, gb, get_bits(gb, 2), &num_bits_left); // bs_extension_id
1056 static void sbr_reset(AACContext *ac, SpectralBandReplication *sbr)
1059 err = sbr_make_f_master(ac, sbr, &sbr->spectrum_params);
1061 err = sbr_make_f_derived(ac, sbr);
1065 sbr_turnoff(sbr);
1077 int ff_decode_sbr_extension(AACContext *ac, SpectralBandReplication *sbr,
1085 sbr->reset = 0;
1087 if (!sbr->sample_rate)
1088 sbr->sample_rate = 2 * ac->oc[1].m4ac.sample_rate; //TODO use the nominal sample rate for arbitrary sample rate support
1098 sbr->kx[0] = sbr->kx[1];
1099 sbr->m[0] = sbr->m[1];
1100 sbr->kx_and_m_pushed = 1;
1104 num_sbr_bits += read_sbr_header(sbr, gb);
1106 if (sbr->reset)
1107 sbr_reset(ac, sbr);
1109 if (sbr->start)
1110 num_sbr_bits += read_sbr_data(ac, sbr, gb, id_aac);
1123 static void sbr_dequant(SpectralBandReplication *sbr, int id_aac)
1128 if (id_aac == TYPE_CPE && sbr->bs_coupling) {
1129 float alpha = sbr->data[0].bs_amp_res ? 1.0f : 0.5f;
1130 float pan_offset = sbr->data[0].bs_amp_res ? 12.0f : 24.0f;
1131 for (e = 1; e <= sbr->data[0].bs_num_env; e++) {
1132 for (k = 0; k < sbr->n[sbr->data[0].bs_freq_res[e]]; k++) {
1133 float temp1 = exp2f(sbr->data[0].env_facs[e][k] * alpha + 7.0f);
1134 float temp2 = exp2f((pan_offset - sbr->data[1].env_facs[e][k]) * alpha);
1141 sbr->data[0].env_facs[e][k] = fac;
1142 sbr->data[1].env_facs[e][k] = fac * temp2;
1145 for (e = 1; e <= sbr->data[0].bs_num_noise; e++) {
1146 for (k = 0; k < sbr->n_q; k++) {
1147 float temp1 = exp2f(NOISE_FLOOR_OFFSET - sbr->data[0].noise_facs[e][k] + 1);
1148 float temp2 = exp2f(12 - sbr->data[1].noise_facs[e][k]);
1155 sbr->data[0].noise_facs[e][k] = fac;
1156 sbr->data[1].noise_facs[e][k] = fac * temp2;
1161 float alpha = sbr->data[ch].bs_amp_res ? 1.0f : 0.5f;
1162 for (e = 1; e <= sbr->data[ch].bs_num_env; e++)
1163 for (k = 0; k < sbr->n[sbr->data[ch].bs_freq_res[e]]; k++){
1164 sbr->data[ch].env_facs[e][k] =
1165 exp2f(alpha * sbr->data[ch].env_facs[e][k] + 6.0f);
1166 if (sbr->data[ch].env_facs[e][k] > 1E20) {
1168 sbr->data[ch].env_facs[e][k] = 1;
1172 for (e = 1; e <= sbr->data[ch].bs_num_noise; e++)
1173 for (k = 0; k < sbr->n_q; k++)
1174 sbr->data[ch].noise_facs[e][k] =
1175 exp2f(NOISE_FLOOR_OFFSET - sbr->data[ch].noise_facs[e][k]);
1317 static void sbr_chirp(SpectralBandReplication *sbr, SBRData *ch_data)
1323 for (i = 0; i < sbr->n_q; i++) {
1338 static int sbr_lf_gen(AACContext *ac, SpectralBandReplication *sbr,
1346 for (k = 0; k < sbr->kx[1]; k++) {
1353 for (k = 0; k < sbr->kx[0]; k++) {
1363 static int sbr_hf_gen(AACContext *ac, SpectralBandReplication *sbr,
1371 int k = sbr->kx[1];
1372 for (j = 0; j < sbr->num_patches; j++) {
1373 for (x = 0; x < sbr->patch_num_subbands[j]; x++, k++) {
1374 const int p = sbr->patch_start_subband[j] + x;
1375 while (g <= sbr->n_q && k >= sbr->f_tablenoise[g])
1385 sbr->dsp.hf_gen(X_high[k] + ENVELOPE_ADJUSTMENT_OFFSET,
1391 if (k < sbr->m[1] + sbr->kx[1])
1392 memset(X_high + k, 0, (sbr->m[1] + sbr->kx[1] - k) * sizeof(*X_high));
1398 static int sbr_x_gen(SpectralBandReplication *sbr, float X[2][38][64],
1404 const int i_Temp = FFMAX(2*sbr->data[ch].t_env_num_env_old - i_f, 0);
1406 for (k = 0; k < sbr->kx[0]; k++) {
1412 for (; k < sbr->kx[0] + sbr->m[0]; k++) {
1419 for (k = 0; k < sbr->kx[1]; k++) {
1425 for (; k < sbr->kx[1] + sbr->m[1]; k++) {
1437 static int sbr_mapping(AACContext *ac, SpectralBandReplication *sbr,
1444 const unsigned int ilim = sbr->n[ch_data->bs_freq_res[e + 1]];
1445 uint16_t *table = ch_data->bs_freq_res[e + 1] ? sbr->f_tablehigh : sbr->f_tablelow;
1448 if (sbr->kx[1] != table[0]) {
1451 sbr_turnoff(sbr);
1456 sbr->e_origmapped[e][m - sbr->kx[1]] = ch_data->env_facs[e+1][i];
1460 for (i = 0; i < sbr->n_q; i++)
1461 for (m = sbr->f_tablenoise[i]; m < sbr->f_tablenoise[i + 1]; m++)
1462 sbr->q_mapped[e][m - sbr->kx[1]] = ch_data->noise_facs[k+1][i];
1464 for (i = 0; i < sbr->n[1]; i++) {
1467 (sbr->f_tablehigh[i] + sbr->f_tablehigh[i + 1]) >> 1;
1469 ch_data->s_indexmapped[e + 1][m_midpoint - sbr->kx[1]] = ch_data->bs_add_harmonic[i] *
1470 (e >= e_a[1] || (ch_data->s_indexmapped[0][m_midpoint - sbr->kx[1]] == 1));
1477 if (ch_data->s_indexmapped[e + 1][m - sbr->kx[1]]) {
1482 memset(&sbr->s_mapped[e][table[i] - sbr->kx[1]], additional_sinusoid_present,
1483 (table[i + 1] - table[i]) * sizeof(sbr->s_mapped[e][0]));
1493 SpectralBandReplication *sbr, SBRData *ch_data)
1496 int kx1 = sbr->kx[1];
1498 if (sbr->bs_interpol_freq) {
1504 for (m = 0; m < sbr->m[1]; m++) {
1505 float sum = sbr->dsp.sum_square(X_high[m+kx1] + ilb, iub - ilb);
1516 const uint16_t *table = ch_data->bs_freq_res[e + 1] ? sbr->f_tablehigh : sbr->f_tablelow;
1518 for (p = 0; p < sbr->n[ch_data->bs_freq_res[e + 1]]; p++) {
1523 sum += sbr->dsp.sum_square(X_high[k] + ilb, iub - ilb);
1538 static void sbr_gain_calc(AACContext *ac, SpectralBandReplication *sbr,
1547 for (k = 0; k < sbr->n_lim; k++) {
1550 for (m = sbr->f_tablelim[k] - sbr->kx[1]; m < sbr->f_tablelim[k + 1] - sbr->kx[1]; m++) {
1551 const float temp = sbr->e_origmapped[e][m] / (1.0f + sbr->q_mapped[e][m]);
1552 sbr->q_m[e][m] = sqrtf(temp * sbr->q_mapped[e][m]);
1553 sbr->s_m[e][m] = sqrtf(temp * ch_data->s_indexmapped[e + 1][m]);
1554 if (!sbr->s_mapped[e][m]) {
1555 sbr->gain[e][m] = sqrtf(sbr->e_origmapped[e][m] /
1556 ((1.0f + sbr->e_curr[e][m]) *
1557 (1.0f + sbr->q_mapped[e][m] * delta)));
1559 sbr->gain[e][m] = sqrtf(sbr->e_origmapped[e][m] * sbr->q_mapped[e][m] /
1560 ((1.0f + sbr->e_curr[e][m]) *
1561 (1.0f + sbr->q_mapped[e][m])));
1564 for (m = sbr->f_tablelim[k] - sbr->kx[1]; m < sbr->f_tablelim[k + 1] - sbr->kx[1]; m++) {
1565 sum[0] += sbr->e_origmapped[e][m];
1566 sum[1] += sbr->e_curr[e][m];
1568 gain_max = limgain[sbr->bs_limiter_gains] * sqrtf((FLT_EPSILON + sum[0]) / (FLT_EPSILON + sum[1]));
1570 for (m = sbr->f_tablelim[k] - sbr->kx[1]; m < sbr->f_tablelim[k + 1] - sbr->kx[1]; m++) {
1571 float q_m_max = sbr->q_m[e][m] * gain_max / sbr->gain[e][m];
1572 sbr->q_m[e][m] = FFMIN(sbr->q_m[e][m], q_m_max);
1573 sbr->gain[e][m] = FFMIN(sbr->gain[e][m], gain_max);
1576 for (m = sbr->f_tablelim[k] - sbr->kx[1]; m < sbr->f_tablelim[k + 1] - sbr->kx[1]; m++) {
1577 sum[0] += sbr->e_origmapped[e][m];
1578 sum[1] += sbr->e_curr[e][m] * sbr->gain[e][m] * sbr->gain[e][m]
1579 + sbr->s_m[e][m] * sbr->s_m[e][m]
1580 + (delta && !sbr->s_m[e][m]) * sbr->q_m[e][m] * sbr->q_m[e][m];
1584 for (m = sbr->f_tablelim[k] - sbr->kx[1]; m < sbr->f_tablelim[k + 1] - sbr->kx[1]; m++) {
1585 sbr->gain[e][m] *= gain_boost;
1586 sbr->q_m[e][m] *= gain_boost;
1587 sbr->s_m[e][m] *= gain_boost;
1596 SpectralBandReplication *sbr, SBRData *ch_data,
1600 const int h_SL = 4 * !sbr->bs_smoothing_mode;
1601 const int kx = sbr->kx[1];
1602 const int m_max = sbr->m[1];
1614 if (sbr->reset) {
1616 memcpy(g_temp[i + 2*ch_data->t_env[0]], sbr->gain[0], m_max * sizeof(sbr->gain[0][0]));
1617 memcpy(q_temp[i + 2*ch_data->t_env[0]], sbr->q_m[0], m_max * sizeof(sbr->q_m[0][0]));
1626 memcpy(g_temp[h_SL + i], sbr->gain[e], m_max * sizeof(sbr->gain[0][0]));
1627 memcpy(q_temp[h_SL + i], sbr->q_m[e], m_max * sizeof(sbr->q_m[0][0]));
1654 sbr->dsp.hf_g_filt(Y1[i] + kx, X_high + kx, g_filt, m_max,
1658 sbr->dsp.hf_apply_noise[indexsine](Y1[i] + kx, sbr->s_m[e],
1666 float *in = sbr->s_m[e];
1682 void ff_sbr_apply(AACContext *ac, SpectralBandReplication *sbr, int id_aac,
1685 int downsampled = ac->oc[1].m4ac.ext_sample_rate < sbr->sample_rate;
1690 if (!sbr->kx_and_m_pushed) {
1691 sbr->kx[0] = sbr->kx[1];
1692 sbr->m[0] = sbr->m[1];
1694 sbr->kx_and_m_pushed = 0;
1697 if (sbr->start) {
1698 sbr_dequant(sbr, id_aac);
1702 sbr_qmf_analysis(&ac->fdsp, &sbr->mdct_ana, &sbr->dsp, ch ? R : L, sbr->data[ch].analysis_filterbank_samples,
1703 (float*)sbr->qmf_filter_scratch,
1704 sbr->data[ch].W, sbr->data[ch].Ypos);
1705 sbr->c.sbr_lf_gen(ac, sbr, sbr->X_low,
1706 (const float (*)[32][32][2]) sbr->data[ch].W,
1707 sbr->data[ch].Ypos);
1708 sbr->data[ch].Ypos ^= 1;
1709 if (sbr->start) {
1710 sbr->c.sbr_hf_inverse_filter(&sbr->dsp, sbr->alpha0, sbr->alpha1,
1711 (const float (*)[40][2]) sbr->X_low, sbr->k[0]);
1712 sbr_chirp(sbr, &sbr->data[ch]);
1713 sbr_hf_gen(ac, sbr, sbr->X_high,
1714 (const float (*)[40][2]) sbr->X_low,
1715 (const float (*)[2]) sbr->alpha0,
1716 (const float (*)[2]) sbr->alpha1,
1717 sbr->data[ch].bw_array, sbr->data[ch].t_env,
1718 sbr->data[ch].bs_num_env);
1721 err = sbr_mapping(ac, sbr, &sbr->data[ch], sbr->data[ch].e_a);
1723 sbr_env_estimate(sbr->e_curr, sbr->X_high, sbr, &sbr->data[ch]);
1724 sbr_gain_calc(ac, sbr, &sbr->data[ch], sbr->data[ch].e_a);
1725 sbr->c.sbr_hf_assemble(sbr->data[ch].Y[sbr->data[ch].Ypos],
1726 (const float (*)[40][2]) sbr->X_high,
1727 sbr, &sbr->data[ch],
1728 sbr->data[ch].e_a);
1733 sbr->c.sbr_x_gen(sbr, sbr->X[ch],
1734 (const float (*)[64][2]) sbr->data[ch].Y[1-sbr->data[ch].Ypos],
1735 (const float (*)[64][2]) sbr->data[ch].Y[ sbr->data[ch].Ypos],
1736 (const float (*)[40][2]) sbr->X_low, ch);
1740 if (sbr->ps.start) {
1741 ff_ps_apply(ac->avctx, &sbr->ps, sbr->X[0], sbr->X[1], sbr->kx[1] + sbr->m[1]);
1743 memcpy(sbr->X[1], sbr->X[0], sizeof(sbr->X[0]));
1748 sbr_qmf_synthesis(&sbr->mdct, &sbr->dsp, &ac->fdsp,
1749 L, sbr->X[0], sbr->qmf_filter_scratch,
1750 sbr->data[0].synthesis_filterbank_samples,
1751 &sbr->data[0].synthesis_filterbank_samples_offset,
1754 sbr_qmf_synthesis(&sbr->mdct, &sbr->dsp, &ac->fdsp,
1755 R, sbr->X[1], sbr->qmf_filter_scratch,
1756 sbr->data[1].synthesis_filterbank_samples,
1757 &sbr->data[1].synthesis_filterbank_samples_offset,