36 #define numQMFSlots 32 //numTimeSlots * RATE
44 10, 20, 34, 10, 20, 34,
73 #define READ_PAR_DATA(PAR, OFFSET, MASK, ERR_CONDITION) \
86 static int read_ ## PAR ## _data(AVCodecContext *avctx, GetBitContext *gb, PSContext *ps, \
87 int8_t (*PAR)[PS_MAX_NR_IIDICC], int table_idx, int e, int dt) \
89 int b, num = ps->nr_ ## PAR ## _par; \
90 VLC_TYPE (*vlc_table)[2] = vlc_ps[table_idx].table; \
92 int e_prev = e ? e - 1 : ps->num_env_old - 1; \
93 e_prev = FFMAX(e_prev, 0); \
94 for (b = 0; b < num; b++) { \
95 int val = PAR[e_prev][b] + get_vlc2(gb, vlc_table, 9, 3) - OFFSET; \
96 if (MASK) val &= MASK; \
103 for (b = 0; b < num; b++) { \
104 val += get_vlc2(gb, vlc_table, 9, 3) - OFFSET; \
105 if (MASK) val &= MASK; \
113 av_log(avctx, AV_LOG_ERROR, "illegal "#PAR"\n"); \
130 if (ps->enable_ipdopd) {
131 for (e = 0; e < ps->num_env; e++) {
142 static void ipdopd_reset(int8_t *opd_hist, int8_t *ipd_hist)
192 for (e = 1; e <= ps->
num_env; e++)
195 for (e = 1; e <= ps->
num_env; e++)
199 for (e = 0; e < ps->
num_env; e++) {
208 for (e = 0; e < ps->
num_env; e++) {
223 int ps_extension_id =
get_bits(gb, 2);
239 if (source >= 0 && source != ps->
num_env) {
271 if (bits_consumed <= bits_left) {
273 return bits_consumed;
275 av_log(avctx,
AV_LOG_ERROR,
"Expected to read %d PS bits actually read %d.\n", bits_left, bits_consumed);
291 for (i = 0; i <
len; i++,
in++) {
292 float re_in = filter[6] *
in[6][0];
294 float im_in = filter[6] * in[6][1];
296 for (j = 0; j < 6; j += 2) {
297 re_op += filter[j+1] * (in[j+1][0] + in[12-j-1][0]);
298 im_op += filter[j+1] * (in[j+1][1] + in[12-j-1][1]);
300 out[ reverse][i][0] = re_in + re_op;
301 out[ reverse][i][1] = im_in + im_op;
302 out[!reverse][i][0] = re_in - re_op;
303 out[!reverse][i][1] = im_in - im_op;
315 for (i = 0; i <
len; i++,
in++) {
317 out[0][i][0] = temp[6][0];
318 out[0][i][1] = temp[6][1];
319 out[1][i][0] = temp[7][0];
320 out[1][i][1] = temp[7][1];
321 out[2][i][0] = temp[0][0];
322 out[2][i][1] = temp[0][1];
323 out[3][i][0] = temp[1][0];
324 out[3][i][1] = temp[1][1];
325 out[4][i][0] = temp[2][0] + temp[5][0];
326 out[4][i][1] = temp[2][1] + temp[5][1];
327 out[5][i][0] = temp[3][0] + temp[4][0];
328 out[5][i][1] = temp[3][1] + temp[4][1];
333 float (*
in)[2],
float (*
out)[32][2],
338 for (i = 0; i <
len; i++,
in++) {
344 float in[5][44][2],
float L[2][38][64],
348 for (i = 0; i < 5; i++) {
349 for (j = 0; j < 38; j++) {
350 in[i][j+6][0] = L[0][j][i];
351 in[i][j+6][1] = L[1][j][i];
368 for (i = 0; i < 5; i++) {
369 memcpy(in[i], in[i]+32, 6 *
sizeof(in[i][0]));
374 float in[91][32][2],
int is34,
int len)
378 for (n = 0; n <
len; n++) {
379 memset(out[0][n], 0, 5*
sizeof(out[0][n][0]));
380 memset(out[1][n], 0, 5*
sizeof(out[1][n][0]));
381 for (i = 0; i < 12; i++) {
382 out[0][n][0] += in[ i][n][0];
383 out[1][n][0] += in[ i][n][1];
385 for (i = 0; i < 8; i++) {
386 out[0][n][1] += in[12+i][n][0];
387 out[1][n][1] += in[12+i][n][1];
389 for (i = 0; i < 4; i++) {
390 out[0][n][2] += in[20+i][n][0];
391 out[1][n][2] += in[20+i][n][1];
392 out[0][n][3] += in[24+i][n][0];
393 out[1][n][3] += in[24+i][n][1];
394 out[0][n][4] += in[28+i][n][0];
395 out[1][n][4] += in[28+i][n][1];
400 for (n = 0; n <
len; n++) {
401 out[0][n][0] = in[0][n][0] + in[1][n][0] + in[2][n][0] +
402 in[3][n][0] + in[4][n][0] + in[5][n][0];
403 out[1][n][0] = in[0][n][1] + in[1][n][1] + in[2][n][1] +
404 in[3][n][1] + in[4][n][1] + in[5][n][1];
405 out[0][n][1] = in[6][n][0] + in[7][n][0];
406 out[1][n][1] = in[6][n][1] + in[7][n][1];
407 out[0][n][2] = in[8][n][0] + in[9][n][0];
408 out[1][n][2] = in[8][n][1] + in[9][n][1];
415 #define DECAY_SLOPE 0.05f
419 static const int NR_BANDS[] = { 71, 91 };
428 static void map_idx_10_to_20(int8_t *par_mapped,
const int8_t *par,
int full)
437 for (; b >= 0; b--) {
438 par_mapped[2*b+1] = par_mapped[2*
b] = par[
b];
442 static void map_idx_34_to_20(int8_t *par_mapped,
const int8_t *par,
int full)
444 par_mapped[ 0] = (2*par[ 0] + par[ 1]) / 3;
445 par_mapped[ 1] = ( par[ 1] + 2*par[ 2]) / 3;
446 par_mapped[ 2] = (2*par[ 3] + par[ 4]) / 3;
447 par_mapped[ 3] = ( par[ 4] + 2*par[ 5]) / 3;
448 par_mapped[ 4] = ( par[ 6] + par[ 7]) / 2;
449 par_mapped[ 5] = ( par[ 8] + par[ 9]) / 2;
450 par_mapped[ 6] = par[10];
451 par_mapped[ 7] = par[11];
452 par_mapped[ 8] = ( par[12] + par[13]) / 2;
453 par_mapped[ 9] = ( par[14] + par[15]) / 2;
454 par_mapped[10] = par[16];
456 par_mapped[11] = par[17];
457 par_mapped[12] = par[18];
458 par_mapped[13] = par[19];
459 par_mapped[14] = ( par[20] + par[21]) / 2;
460 par_mapped[15] = ( par[22] + par[23]) / 2;
461 par_mapped[16] = ( par[24] + par[25]) / 2;
462 par_mapped[17] = ( par[26] + par[27]) / 2;
463 par_mapped[18] = ( par[28] + par[29] + par[30] + par[31]) / 4;
464 par_mapped[19] = ( par[32] + par[33]) / 2;
470 par[ 0] = (2*par[ 0] + par[ 1]) * 0.33333333f;
471 par[ 1] = ( par[ 1] + 2*par[ 2]) * 0.33333333f;
472 par[ 2] = (2*par[ 3] + par[ 4]) * 0.33333333f;
473 par[ 3] = ( par[ 4] + 2*par[ 5]) * 0.33333333f;
474 par[ 4] = ( par[ 6] + par[ 7]) * 0.5f;
475 par[ 5] = ( par[ 8] + par[ 9]) * 0.5f;
478 par[ 8] = ( par[12] + par[13]) * 0.5f;
479 par[ 9] = ( par[14] + par[15]) * 0.5f;
484 par[14] = ( par[20] + par[21]) * 0.5f;
485 par[15] = ( par[22] + par[23]) * 0.5f;
486 par[16] = ( par[24] + par[25]) * 0.5f;
487 par[17] = ( par[26] + par[27]) * 0.5f;
488 par[18] = ( par[28] + par[29] + par[30] + par[31]) * 0.25f;
489 par[19] = ( par[32] + par[33]) * 0.5f;
492 static void map_idx_10_to_34(int8_t *par_mapped,
const int8_t *par,
int full)
495 par_mapped[33] = par[9];
496 par_mapped[32] = par[9];
497 par_mapped[31] = par[9];
498 par_mapped[30] = par[9];
499 par_mapped[29] = par[9];
500 par_mapped[28] = par[9];
501 par_mapped[27] = par[8];
502 par_mapped[26] = par[8];
503 par_mapped[25] = par[8];
504 par_mapped[24] = par[8];
505 par_mapped[23] = par[7];
506 par_mapped[22] = par[7];
507 par_mapped[21] = par[7];
508 par_mapped[20] = par[7];
509 par_mapped[19] = par[6];
510 par_mapped[18] = par[6];
511 par_mapped[17] = par[5];
512 par_mapped[16] = par[5];
516 par_mapped[15] = par[4];
517 par_mapped[14] = par[4];
518 par_mapped[13] = par[4];
519 par_mapped[12] = par[4];
520 par_mapped[11] = par[3];
521 par_mapped[10] = par[3];
522 par_mapped[ 9] = par[2];
523 par_mapped[ 8] = par[2];
524 par_mapped[ 7] = par[2];
525 par_mapped[ 6] = par[2];
526 par_mapped[ 5] = par[1];
527 par_mapped[ 4] = par[1];
528 par_mapped[ 3] = par[1];
529 par_mapped[ 2] = par[0];
530 par_mapped[ 1] = par[0];
531 par_mapped[ 0] = par[0];
534 static void map_idx_20_to_34(int8_t *par_mapped,
const int8_t *par,
int full)
537 par_mapped[33] = par[19];
538 par_mapped[32] = par[19];
539 par_mapped[31] = par[18];
540 par_mapped[30] = par[18];
541 par_mapped[29] = par[18];
542 par_mapped[28] = par[18];
543 par_mapped[27] = par[17];
544 par_mapped[26] = par[17];
545 par_mapped[25] = par[16];
546 par_mapped[24] = par[16];
547 par_mapped[23] = par[15];
548 par_mapped[22] = par[15];
549 par_mapped[21] = par[14];
550 par_mapped[20] = par[14];
551 par_mapped[19] = par[13];
552 par_mapped[18] = par[12];
553 par_mapped[17] = par[11];
555 par_mapped[16] = par[10];
556 par_mapped[15] = par[ 9];
557 par_mapped[14] = par[ 9];
558 par_mapped[13] = par[ 8];
559 par_mapped[12] = par[ 8];
560 par_mapped[11] = par[ 7];
561 par_mapped[10] = par[ 6];
562 par_mapped[ 9] = par[ 5];
563 par_mapped[ 8] = par[ 5];
564 par_mapped[ 7] = par[ 4];
565 par_mapped[ 6] = par[ 4];
566 par_mapped[ 5] = par[ 3];
567 par_mapped[ 4] = (par[ 2] + par[ 3]) / 2;
568 par_mapped[ 3] = par[ 2];
569 par_mapped[ 2] = par[ 1];
570 par_mapped[ 1] = (par[ 0] + par[ 1]) / 2;
571 par_mapped[ 0] = par[ 0];
605 par[ 4] = (par[ 2] + par[ 3]) * 0.5f;
608 par[ 1] = (par[ 0] + par[ 1]) * 0.5f;
622 const float peak_decay_factor = 0.76592833836465f;
623 const float transient_impact = 1.5f;
624 const float a_smooth = 0.25f;
628 memset(power, 0, 34 *
sizeof(*power));
638 for (k = 0; k < NR_BANDS[is34]; k++) {
644 for (i = 0; i < NR_PAR_BANDS[is34]; i++) {
645 for (n = n0; n < nL; n++) {
646 float decayed_peak = peak_decay_factor * peak_decay_nrg[i];
648 peak_decay_nrg[i] =
FFMAX(decayed_peak, power[i][n]);
649 power_smooth[i] += a_smooth * (power[i][n] - power_smooth[i]);
650 peak_decay_diff_smooth[i] += a_smooth * (peak_decay_nrg[i] - power[i][n] - peak_decay_diff_smooth[i]);
651 denom = transient_impact * peak_decay_diff_smooth[i];
652 transient_gain[i][n] = (denom > power_smooth[i]) ?
653 power_smooth[i] / denom : 1.0f;
665 for (k = 0; k < NR_ALLPASS_BANDS[is34]; k++) {
667 float g_decay_slope = 1.f -
DECAY_SLOPE * (k - DECAY_CUTOFF[is34]);
668 g_decay_slope = av_clipf(g_decay_slope, 0.f, 1.f);
669 memcpy(delay[k], delay[k]+nL,
PS_MAX_DELAY*
sizeof(delay[k][0]));
672 memcpy(ap_delay[k][m], ap_delay[k][m]+
numQMFSlots, 5*
sizeof(ap_delay[k][m][0]));
677 transient_gain[b], g_decay_slope, nL - n0);
679 for (; k < SHORT_DELAY_BAND[is34]; k++) {
681 memcpy(delay[k], delay[k]+nL,
PS_MAX_DELAY*
sizeof(delay[k][0]));
685 transient_gain[i], nL - n0);
687 for (; k < NR_BANDS[is34]; k++) {
689 memcpy(delay[k], delay[k]+nL,
PS_MAX_DELAY*
sizeof(delay[k][0]));
693 transient_gain[i], nL - n0);
698 int8_t (*par)[PS_MAX_NR_IIDICC],
699 int num_par,
int num_env,
int full)
703 if (num_par == 20 || num_par == 11) {
704 for (e = 0; e < num_env; e++) {
707 }
else if (num_par == 10 || num_par == 5) {
708 for (e = 0; e < num_env; e++) {
716 static void remap20(int8_t (**p_par_mapped)[PS_MAX_NR_IIDICC],
717 int8_t (*par)[PS_MAX_NR_IIDICC],
718 int num_par,
int num_env,
int full)
722 if (num_par == 34 || num_par == 17) {
723 for (e = 0; e < num_env; e++) {
726 }
else if (num_par == 10 || num_par == 5) {
727 for (e = 0; e < num_env; e++) {
758 memcpy(H11[0][0], H11[0][ps->
num_env_old], PS_MAX_NR_IIDICC*
sizeof(H11[0][0][0]));
759 memcpy(H11[1][0], H11[1][ps->
num_env_old], PS_MAX_NR_IIDICC*
sizeof(H11[1][0][0]));
760 memcpy(H12[0][0], H12[0][ps->
num_env_old], PS_MAX_NR_IIDICC*
sizeof(H12[0][0][0]));
761 memcpy(H12[1][0], H12[1][ps->
num_env_old], PS_MAX_NR_IIDICC*
sizeof(H12[1][0][0]));
762 memcpy(H21[0][0], H21[0][ps->
num_env_old], PS_MAX_NR_IIDICC*
sizeof(H21[0][0][0]));
763 memcpy(H21[1][0], H21[1][ps->
num_env_old], PS_MAX_NR_IIDICC*
sizeof(H21[1][0][0]));
764 memcpy(H22[0][0], H22[0][ps->
num_env_old], PS_MAX_NR_IIDICC*
sizeof(H22[0][0][0]));
765 memcpy(H22[1][0], H22[1][ps->
num_env_old], PS_MAX_NR_IIDICC*
sizeof(H22[1][0][0]));
807 for (e = 0; e < ps->
num_env; e++) {
808 for (b = 0; b < NR_PAR_BANDS[is34]; b++) {
809 float h11, h12, h21, h22;
810 h11 = H_LUT[iid_mapped[e][
b] + 7 + 23 * ps->
iid_quant][icc_mapped[e][
b]][0];
811 h12 = H_LUT[iid_mapped[e][
b] + 7 + 23 * ps->
iid_quant][icc_mapped[e][
b]][1];
812 h21 = H_LUT[iid_mapped[e][
b] + 7 + 23 * ps->
iid_quant][icc_mapped[e][
b]][2];
813 h22 = H_LUT[iid_mapped[e][
b] + 7 + 23 * ps->
iid_quant][icc_mapped[e][
b]][3];
817 float h11i, h12i, h21i, h22i;
818 float ipd_adj_re, ipd_adj_im;
819 int opd_idx = opd_hist[
b] * 8 + opd_mapped[e][
b];
820 int ipd_idx = ipd_hist[
b] * 8 + ipd_mapped[e][
b];
825 opd_hist[
b] = opd_idx & 0x3F;
826 ipd_hist[
b] = ipd_idx & 0x3F;
828 ipd_adj_re = opd_re*ipd_re + opd_im*ipd_im;
829 ipd_adj_im = opd_im*ipd_re - opd_re*ipd_im;
832 h12i = h12 * ipd_adj_im;
833 h12 = h12 * ipd_adj_re;
836 h22i = h22 * ipd_adj_im;
837 h22 = h22 * ipd_adj_re;
838 H11[1][e+1][
b] = h11i;
839 H12[1][e+1][
b] = h12i;
840 H21[1][e+1][
b] = h21i;
841 H22[1][e+1][
b] = h22i;
843 H11[0][e+1][
b] = h11;
844 H12[0][e+1][
b] = h12;
845 H21[0][e+1][
b] = h21;
846 H22[0][e+1][
b] = h22;
848 for (k = 0; k < NR_BANDS[is34]; k++) {
853 float width = 1.f / (stop - start);
855 h[0][0] = H11[0][e][
b];
856 h[0][1] = H12[0][e][
b];
857 h[0][2] = H21[0][e][
b];
858 h[0][3] = H22[0][e][
b];
861 if ((is34 && k <= 13 && k >= 9) || (!is34 && k <= 1)) {
862 h[1][0] = -H11[1][e][
b];
863 h[1][1] = -H12[1][e][
b];
864 h[1][2] = -H21[1][e][
b];
865 h[1][3] = -H22[1][e][
b];
867 h[1][0] = H11[1][e][
b];
868 h[1][1] = H12[1][e][
b];
869 h[1][2] = H21[1][e][
b];
870 h[1][3] = H22[1][e][
b];
874 h_step[0][0] = (H11[0][e+1][
b] - h[0][0]) * width;
875 h_step[0][1] = (H12[0][e+1][
b] - h[0][1]) * width;
876 h_step[0][2] = (H21[0][e+1][
b] - h[0][2]) * width;
877 h_step[0][3] = (H22[0][e+1][
b] - h[0][3]) * width;
879 h_step[1][0] = (H11[1][e+1][
b] - h[1][0]) * width;
880 h_step[1][1] = (H12[1][e+1][
b] - h[1][1]) * width;
881 h_step[1][2] = (H21[1][e+1][
b] - h[1][2]) * width;
882 h_step[1][3] = (H22[1][e+1][
b] - h[1][3]) * width;
885 l[k] + start + 1,
r[k] + start + 1,
886 h, h_step, stop - start);
898 top += NR_BANDS[is34] - 64;
899 memset(ps->
delay+top, 0, (NR_BANDS[is34] - top)*
sizeof(ps->
delay[0]));
900 if (top < NR_ALLPASS_BANDS[is34])
904 decorrelation(ps, Rbuf, (
const float (*)[32][2]) Lbuf, is34);
912 #define PS_INIT_VLC_STATIC(num, size) \
913 INIT_VLC_STATIC(&vlc_ps[num], 9, ps_tmp[num].table_size / ps_tmp[num].elem_size, \
914 ps_tmp[num].ps_bits, 1, 1, \
915 ps_tmp[num].ps_codes, ps_tmp[num].elem_size, ps_tmp[num].elem_size, \
918 #define PS_VLC_ROW(name) \
919 { name ## _codes, name ## _bits, sizeof(name ## _codes), sizeof(name ## _codes[0]) }
923 static const struct {
924 const void *ps_codes, *ps_bits;
925 const unsigned int table_size, elem_size;