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++) {
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;
314 for (i = 0; i <
len; i++, in++) {
316 out[0][i][0] = temp[6][0];
317 out[0][i][1] = temp[6][1];
318 out[1][i][0] = temp[7][0];
319 out[1][i][1] = temp[7][1];
320 out[2][i][0] = temp[0][0];
321 out[2][i][1] = temp[0][1];
322 out[3][i][0] = temp[1][0];
323 out[3][i][1] = temp[1][1];
324 out[4][i][0] = temp[2][0] + temp[5][0];
325 out[4][i][1] = temp[2][1] + temp[5][1];
326 out[5][i][0] = temp[3][0] + temp[4][0];
327 out[5][i][1] = temp[3][1] + temp[4][1];
335 for (i = 0; i <
len; i++, in++) {
341 float in[5][44][2],
float L[2][38][64],
345 for (i = 0; i < 5; i++) {
346 for (j = 0; j < 38; j++) {
347 in[i][j+6][0] = L[0][j][i];
348 in[i][j+6][1] = L[1][j][i];
365 for (i = 0; i < 5; i++) {
366 memcpy(in[i], in[i]+32, 6 *
sizeof(in[i][0]));
371 float in[91][32][2],
int is34,
int len)
375 for (n = 0; n <
len; n++) {
376 memset(out[0][n], 0, 5*
sizeof(out[0][n][0]));
377 memset(out[1][n], 0, 5*
sizeof(out[1][n][0]));
378 for (i = 0; i < 12; i++) {
379 out[0][n][0] += in[ i][n][0];
380 out[1][n][0] += in[ i][n][1];
382 for (i = 0; i < 8; i++) {
383 out[0][n][1] += in[12+i][n][0];
384 out[1][n][1] += in[12+i][n][1];
386 for (i = 0; i < 4; i++) {
387 out[0][n][2] += in[20+i][n][0];
388 out[1][n][2] += in[20+i][n][1];
389 out[0][n][3] += in[24+i][n][0];
390 out[1][n][3] += in[24+i][n][1];
391 out[0][n][4] += in[28+i][n][0];
392 out[1][n][4] += in[28+i][n][1];
397 for (n = 0; n <
len; n++) {
398 out[0][n][0] = in[0][n][0] + in[1][n][0] + in[2][n][0] +
399 in[3][n][0] + in[4][n][0] + in[5][n][0];
400 out[1][n][0] = in[0][n][1] + in[1][n][1] + in[2][n][1] +
401 in[3][n][1] + in[4][n][1] + in[5][n][1];
402 out[0][n][1] = in[6][n][0] + in[7][n][0];
403 out[1][n][1] = in[6][n][1] + in[7][n][1];
404 out[0][n][2] = in[8][n][0] + in[9][n][0];
405 out[1][n][2] = in[8][n][1] + in[9][n][1];
412 #define DECAY_SLOPE 0.05f
434 for (; b >= 0; b--) {
435 par_mapped[2*b+1] = par_mapped[2*
b] = par[
b];
441 par_mapped[ 0] = (2*par[ 0] + par[ 1]) / 3;
442 par_mapped[ 1] = ( par[ 1] + 2*par[ 2]) / 3;
443 par_mapped[ 2] = (2*par[ 3] + par[ 4]) / 3;
444 par_mapped[ 3] = ( par[ 4] + 2*par[ 5]) / 3;
445 par_mapped[ 4] = ( par[ 6] + par[ 7]) / 2;
446 par_mapped[ 5] = ( par[ 8] + par[ 9]) / 2;
447 par_mapped[ 6] = par[10];
448 par_mapped[ 7] = par[11];
449 par_mapped[ 8] = ( par[12] + par[13]) / 2;
450 par_mapped[ 9] = ( par[14] + par[15]) / 2;
451 par_mapped[10] = par[16];
453 par_mapped[11] = par[17];
454 par_mapped[12] = par[18];
455 par_mapped[13] = par[19];
456 par_mapped[14] = ( par[20] + par[21]) / 2;
457 par_mapped[15] = ( par[22] + par[23]) / 2;
458 par_mapped[16] = ( par[24] + par[25]) / 2;
459 par_mapped[17] = ( par[26] + par[27]) / 2;
460 par_mapped[18] = ( par[28] + par[29] + par[30] + par[31]) / 4;
461 par_mapped[19] = ( par[32] + par[33]) / 2;
467 par[ 0] = (2*par[ 0] + par[ 1]) * 0.33333333f;
468 par[ 1] = ( par[ 1] + 2*par[ 2]) * 0.33333333f;
469 par[ 2] = (2*par[ 3] + par[ 4]) * 0.33333333f;
470 par[ 3] = ( par[ 4] + 2*par[ 5]) * 0.33333333f;
471 par[ 4] = ( par[ 6] + par[ 7]) * 0.5f;
472 par[ 5] = ( par[ 8] + par[ 9]) * 0.5f;
475 par[ 8] = ( par[12] + par[13]) * 0.5f;
476 par[ 9] = ( par[14] + par[15]) * 0.5f;
481 par[14] = ( par[20] + par[21]) * 0.5f;
482 par[15] = ( par[22] + par[23]) * 0.5f;
483 par[16] = ( par[24] + par[25]) * 0.5f;
484 par[17] = ( par[26] + par[27]) * 0.5f;
485 par[18] = ( par[28] + par[29] + par[30] + par[31]) * 0.25f;
486 par[19] = ( par[32] + par[33]) * 0.5f;
492 par_mapped[33] = par[9];
493 par_mapped[32] = par[9];
494 par_mapped[31] = par[9];
495 par_mapped[30] = par[9];
496 par_mapped[29] = par[9];
497 par_mapped[28] = par[9];
498 par_mapped[27] = par[8];
499 par_mapped[26] = par[8];
500 par_mapped[25] = par[8];
501 par_mapped[24] = par[8];
502 par_mapped[23] = par[7];
503 par_mapped[22] = par[7];
504 par_mapped[21] = par[7];
505 par_mapped[20] = par[7];
506 par_mapped[19] = par[6];
507 par_mapped[18] = par[6];
508 par_mapped[17] = par[5];
509 par_mapped[16] = par[5];
513 par_mapped[15] = par[4];
514 par_mapped[14] = par[4];
515 par_mapped[13] = par[4];
516 par_mapped[12] = par[4];
517 par_mapped[11] = par[3];
518 par_mapped[10] = par[3];
519 par_mapped[ 9] = par[2];
520 par_mapped[ 8] = par[2];
521 par_mapped[ 7] = par[2];
522 par_mapped[ 6] = par[2];
523 par_mapped[ 5] = par[1];
524 par_mapped[ 4] = par[1];
525 par_mapped[ 3] = par[1];
526 par_mapped[ 2] = par[0];
527 par_mapped[ 1] = par[0];
528 par_mapped[ 0] = par[0];
534 par_mapped[33] = par[19];
535 par_mapped[32] = par[19];
536 par_mapped[31] = par[18];
537 par_mapped[30] = par[18];
538 par_mapped[29] = par[18];
539 par_mapped[28] = par[18];
540 par_mapped[27] = par[17];
541 par_mapped[26] = par[17];
542 par_mapped[25] = par[16];
543 par_mapped[24] = par[16];
544 par_mapped[23] = par[15];
545 par_mapped[22] = par[15];
546 par_mapped[21] = par[14];
547 par_mapped[20] = par[14];
548 par_mapped[19] = par[13];
549 par_mapped[18] = par[12];
550 par_mapped[17] = par[11];
552 par_mapped[16] = par[10];
553 par_mapped[15] = par[ 9];
554 par_mapped[14] = par[ 9];
555 par_mapped[13] = par[ 8];
556 par_mapped[12] = par[ 8];
557 par_mapped[11] = par[ 7];
558 par_mapped[10] = par[ 6];
559 par_mapped[ 9] = par[ 5];
560 par_mapped[ 8] = par[ 5];
561 par_mapped[ 7] = par[ 4];
562 par_mapped[ 6] = par[ 4];
563 par_mapped[ 5] = par[ 3];
564 par_mapped[ 4] = (par[ 2] + par[ 3]) / 2;
565 par_mapped[ 3] = par[ 2];
566 par_mapped[ 2] = par[ 1];
567 par_mapped[ 1] = (par[ 0] + par[ 1]) / 2;
568 par_mapped[ 0] = par[ 0];
602 par[ 4] = (par[ 2] + par[ 3]) * 0.5f;
605 par[ 1] = (par[ 0] + par[ 1]) * 0.5f;
619 const float peak_decay_factor = 0.76592833836465f;
620 const float transient_impact = 1.5f;
621 const float a_smooth = 0.25f;
625 memset(power, 0, 34 *
sizeof(*power));
635 for (k = 0; k < NR_BANDS[is34]; k++) {
641 for (i = 0; i < NR_PAR_BANDS[is34]; i++) {
642 for (n = n0; n < nL; n++) {
643 float decayed_peak = peak_decay_factor * peak_decay_nrg[i];
645 peak_decay_nrg[i] =
FFMAX(decayed_peak, power[i][n]);
646 power_smooth[i] += a_smooth * (power[i][n] - power_smooth[i]);
647 peak_decay_diff_smooth[i] += a_smooth * (peak_decay_nrg[i] - power[i][n] - peak_decay_diff_smooth[i]);
648 denom = transient_impact * peak_decay_diff_smooth[i];
649 transient_gain[i][n] = (denom > power_smooth[i]) ?
650 power_smooth[i] / denom : 1.0f;
662 for (k = 0; k < NR_ALLPASS_BANDS[is34]; k++) {
664 float g_decay_slope = 1.f -
DECAY_SLOPE * (k - DECAY_CUTOFF[is34]);
665 g_decay_slope = av_clipf(g_decay_slope, 0.f, 1.f);
666 memcpy(delay[k], delay[k]+nL,
PS_MAX_DELAY*
sizeof(delay[k][0]));
669 memcpy(ap_delay[k][m], ap_delay[k][m]+
numQMFSlots, 5*
sizeof(ap_delay[k][m][0]));
673 transient_gain[b], g_decay_slope, nL - n0);
675 for (; k < SHORT_DELAY_BAND[is34]; k++) {
677 memcpy(delay[k], delay[k]+nL,
PS_MAX_DELAY*
sizeof(delay[k][0]));
681 transient_gain[i], nL - n0);
683 for (; k < NR_BANDS[is34]; k++) {
685 memcpy(delay[k], delay[k]+nL,
PS_MAX_DELAY*
sizeof(delay[k][0]));
689 transient_gain[i], nL - n0);
694 int8_t (*par)[PS_MAX_NR_IIDICC],
695 int num_par,
int num_env,
int full)
699 if (num_par == 20 || num_par == 11) {
700 for (e = 0; e < num_env; e++) {
703 }
else if (num_par == 10 || num_par == 5) {
704 for (e = 0; e < num_env; e++) {
712 static void remap20(int8_t (**p_par_mapped)[PS_MAX_NR_IIDICC],
713 int8_t (*par)[PS_MAX_NR_IIDICC],
714 int num_par,
int num_env,
int full)
718 if (num_par == 34 || num_par == 17) {
719 for (e = 0; e < num_env; e++) {
722 }
else if (num_par == 10 || num_par == 5) {
723 for (e = 0; e < num_env; e++) {
754 memcpy(H11[0][0], H11[0][ps->
num_env_old], PS_MAX_NR_IIDICC*
sizeof(H11[0][0][0]));
755 memcpy(H11[1][0], H11[1][ps->
num_env_old], PS_MAX_NR_IIDICC*
sizeof(H11[1][0][0]));
756 memcpy(H12[0][0], H12[0][ps->
num_env_old], PS_MAX_NR_IIDICC*
sizeof(H12[0][0][0]));
757 memcpy(H12[1][0], H12[1][ps->
num_env_old], PS_MAX_NR_IIDICC*
sizeof(H12[1][0][0]));
758 memcpy(H21[0][0], H21[0][ps->
num_env_old], PS_MAX_NR_IIDICC*
sizeof(H21[0][0][0]));
759 memcpy(H21[1][0], H21[1][ps->
num_env_old], PS_MAX_NR_IIDICC*
sizeof(H21[1][0][0]));
760 memcpy(H22[0][0], H22[0][ps->
num_env_old], PS_MAX_NR_IIDICC*
sizeof(H22[0][0][0]));
761 memcpy(H22[1][0], H22[1][ps->
num_env_old], PS_MAX_NR_IIDICC*
sizeof(H22[1][0][0]));
803 for (e = 0; e < ps->
num_env; e++) {
804 for (b = 0; b < NR_PAR_BANDS[is34]; b++) {
805 float h11, h12, h21, h22;
806 h11 = H_LUT[iid_mapped[e][
b] + 7 + 23 * ps->
iid_quant][icc_mapped[e][
b]][0];
807 h12 = H_LUT[iid_mapped[e][
b] + 7 + 23 * ps->
iid_quant][icc_mapped[e][
b]][1];
808 h21 = H_LUT[iid_mapped[e][
b] + 7 + 23 * ps->
iid_quant][icc_mapped[e][
b]][2];
809 h22 = H_LUT[iid_mapped[e][
b] + 7 + 23 * ps->
iid_quant][icc_mapped[e][
b]][3];
813 float h11i, h12i, h21i, h22i;
814 float ipd_adj_re, ipd_adj_im;
815 int opd_idx = opd_hist[
b] * 8 + opd_mapped[e][
b];
816 int ipd_idx = ipd_hist[
b] * 8 + ipd_mapped[e][
b];
821 opd_hist[
b] = opd_idx & 0x3F;
822 ipd_hist[
b] = ipd_idx & 0x3F;
824 ipd_adj_re = opd_re*ipd_re + opd_im*ipd_im;
825 ipd_adj_im = opd_im*ipd_re - opd_re*ipd_im;
828 h12i = h12 * ipd_adj_im;
829 h12 = h12 * ipd_adj_re;
832 h22i = h22 * ipd_adj_im;
833 h22 = h22 * ipd_adj_re;
834 H11[1][e+1][
b] = h11i;
835 H12[1][e+1][
b] = h12i;
836 H21[1][e+1][
b] = h21i;
837 H22[1][e+1][
b] = h22i;
839 H11[0][e+1][
b] = h11;
840 H12[0][e+1][
b] = h12;
841 H21[0][e+1][
b] = h21;
842 H22[0][e+1][
b] = h22;
844 for (k = 0; k < NR_BANDS[is34]; k++) {
849 float width = 1.f / (stop - start);
851 h[0][0] = H11[0][e][
b];
852 h[0][1] = H12[0][e][
b];
853 h[0][2] = H21[0][e][
b];
854 h[0][3] = H22[0][e][
b];
857 if ((is34 && k <= 13 && k >= 9) || (!is34 && k <= 1)) {
858 h[1][0] = -H11[1][e][
b];
859 h[1][1] = -H12[1][e][
b];
860 h[1][2] = -H21[1][e][
b];
861 h[1][3] = -H22[1][e][
b];
863 h[1][0] = H11[1][e][
b];
864 h[1][1] = H12[1][e][
b];
865 h[1][2] = H21[1][e][
b];
866 h[1][3] = H22[1][e][
b];
870 h_step[0][0] = (H11[0][e+1][
b] - h[0][0]) * width;
871 h_step[0][1] = (H12[0][e+1][
b] - h[0][1]) * width;
872 h_step[0][2] = (H21[0][e+1][
b] - h[0][2]) * width;
873 h_step[0][3] = (H22[0][e+1][
b] - h[0][3]) * width;
875 h_step[1][0] = (H11[1][e+1][
b] - h[1][0]) * width;
876 h_step[1][1] = (H12[1][e+1][
b] - h[1][1]) * width;
877 h_step[1][2] = (H21[1][e+1][
b] - h[1][2]) * width;
878 h_step[1][3] = (H22[1][e+1][
b] - h[1][3]) * width;
881 l[k] + start + 1,
r[k] + start + 1,
882 h, h_step, stop - start);
894 top += NR_BANDS[is34] - 64;
895 memset(ps->
delay+top, 0, (NR_BANDS[is34] - top)*
sizeof(ps->
delay[0]));
896 if (top < NR_ALLPASS_BANDS[is34])
897 memset(ps->
ap_delay + top, 0, (NR_ALLPASS_BANDS[is34] - top)*
sizeof(ps->
ap_delay[0]));
908 #define PS_INIT_VLC_STATIC(num, size) \
909 INIT_VLC_STATIC(&vlc_ps[num], 9, ps_tmp[num].table_size / ps_tmp[num].elem_size, \
910 ps_tmp[num].ps_bits, 1, 1, \
911 ps_tmp[num].ps_codes, ps_tmp[num].elem_size, ps_tmp[num].elem_size, \
914 #define PS_VLC_ROW(name) \
915 { name ## _codes, name ## _bits, sizeof(name ## _codes), sizeof(name ## _codes[0]) }
919 static const struct {
920 const void *ps_codes, *ps_bits;
921 const unsigned int table_size, elem_size;
static void map_idx_20_to_34(int8_t *par_mapped, const int8_t *par, int full)
int ff_ps_apply(AVCodecContext *avctx, PSContext *ps, float L[2][38][64], float R[2][38][64], int top)
static float Q_fract_allpass[2][50][3][2]
static unsigned int get_bits(GetBitContext *s, int n)
Read 1-25 bits.
av_cold void ff_psdsp_init(PSDSPContext *s)
static void skip_bits_long(GetBitContext *s, int n)
void(* hybrid_analysis_ileave)(float(*out)[32][2], float L[2][38][64], int i, int len)
float H12[2][PS_MAX_NUM_ENV+1][PS_MAX_NR_IIDICC]
static void hybrid_analysis(PSDSPContext *dsp, float out[91][32][2], float in[5][44][2], float L[2][38][64], int is34, int len)
static const int NR_BANDS[]
Number of frequency bands that can be addressed by the sub subband index, k.
float delay[PS_MAX_SSB][PS_QMF_TIME_SLOTS+PS_MAX_DELAY][2]
static const int8_t nr_iidopd_par_tab[]
static const int8_t huff_offset[]
static void map_idx_10_to_20(int8_t *par_mapped, const int8_t *par, int full)
Table 8.46.
static void map_idx_34_to_20(int8_t *par_mapped, const int8_t *par, int full)
void(* mul_pair_single)(float(*dst)[2], float(*src0)[2], float *src1, int n)
static float pd_im_smooth[8 *8 *8]
static const int8_t k_to_i_34[]
Table 8.49.
float peak_decay_diff_smooth[34]
static float f34_1_8[8][8][2]
#define PS_INIT_VLC_STATIC(num, size)
#define READ_PAR_DATA(PAR, OFFSET, MASK, ERR_CONDITION)
void(* stereo_interpolate[2])(float(*l)[2], float(*r)[2], float h[2][4], float h_step[2][4], int len)
static void hybrid6_cx(PSDSPContext *dsp, float(*in)[2], float(*out)[32][2], const float(*filter)[8][2], int len)
Split one subband into 6 subsubbands with a complex filter.
static void decorrelation(PSContext *ps, float(*out)[32][2], const float(*s)[32][2], int is34)
static float f20_0_8[8][8][2]
static void hybrid4_8_12_cx(PSDSPContext *dsp, float(*in)[2], float(*out)[32][2], const float(*filter)[8][2], int N, int len)
static void remap34(int8_t(**p_par_mapped)[PS_MAX_NR_IIDICC], int8_t(*par)[PS_MAX_NR_IIDICC], int num_par, int num_env, int full)
static int get_bits_count(const GetBitContext *s)
static float phi_fract[2][50][2]
bitstream reader API header.
const uint8_t ff_log2_tab[256]
#define LOCAL_ALIGNED_16(t, v,...)
static void map_idx_10_to_34(int8_t *par_mapped, const int8_t *par, int full)
#define DECAY_SLOPE
All-pass filter decay slope.
int8_t icc_par[PS_MAX_NUM_ENV][PS_MAX_NR_IIDICC]
Inter-Channel Coherence Parameters.
static float f34_2_4[4][8][2]
void(* hybrid_synthesis_deint)(float out[2][38][64], float(*in)[32][2], int i, int len)
static int read_iid_data(AVCodecContext *avctx, GetBitContext *gb, PSContext *ps, int8_t(*iid)[34], int table_idx, int e, int dt)
\ * Read Inter-channel Intensity Difference/Inter-Channel Coherence/ \ * Inter-channel Phase Differen...
static const int SHORT_DELAY_BAND[]
First stereo band using the short one sample delay.
#define PS_BASELINE
Operate in Baseline PS mode.
av_cold void ff_ps_ctx_init(PSContext *ps)
void(* hybrid_analysis)(float(*out)[2], float(*in)[2], const float(*filter)[8][2], int stride, int n)
void av_log(void *avcl, int level, const char *fmt,...)
static int read_icc_data(AVCodecContext *avctx, GetBitContext *gb, PSContext *ps, int8_t(*icc)[34], int table_idx, int e, int dt)
\ * Read Inter-channel Intensity Difference/Inter-Channel Coherence/ \ * Inter-channel Phase Differen...
float H11[2][PS_MAX_NUM_ENV+1][PS_MAX_NR_IIDICC]
static void filter(MpegAudioContext *s, int ch, const short *samples, int incr)
static const float g1_Q2[]
static float HA[46][8][4]
static int read_ipdopd_data(AVCodecContext *avctx, GetBitContext *gb, PSContext *ps, int8_t(*ipdopd)[34], int table_idx, int e, int dt)
\ * Read Inter-channel Intensity Difference/Inter-Channel Coherence/ \ * Inter-channel Phase Differen...
int8_t ipd_hist[PS_MAX_NR_IIDICC]
static float pd_re_smooth[8 *8 *8]
static const int huff_iid[]
void(* decorrelate)(float(*out)[2], float(*delay)[2], float(*ap_delay)[PS_QMF_TIME_SLOTS+PS_MAX_AP_DELAY][2], const float phi_fract[2], float(*Q_fract)[2], const float *transient_gain, float g_decay_slope, int len)
int8_t ipd_par[PS_MAX_NUM_ENV][PS_MAX_NR_IIDICC]
Inter-channel Phase Difference Parameters.
static const int8_t nr_iidicc_par_tab[]
float ap_delay[PS_MAX_AP_BANDS][PS_AP_LINKS][PS_QMF_TIME_SLOTS+PS_MAX_AP_DELAY][2]
av_cold void ff_ps_init(void)
static const int NR_ALLPASS_BANDS[]
Number of all-pass filer bands.
static const int8_t num_env_tab[2][4]
main external API structure.
static unsigned int get_bits1(GetBitContext *s)
int8_t iid_par[PS_MAX_NUM_ENV][PS_MAX_NR_IIDICC]
Inter-channel Intensity Difference Parameters.
static void skip_bits1(GetBitContext *s)
static const int NR_PAR_BANDS[]
Number of frequency bands that can be addressed by the parameter index, b(k)
static void map_val_34_to_20(float par[PS_MAX_NR_IIDICC])
static void skip_bits(GetBitContext *s, int n)
static void hybrid_synthesis(PSDSPContext *dsp, float out[2][38][64], float in[91][32][2], int is34, int len)
static void ps_tableinit(void)
static float HB[46][8][4]
int8_t opd_hist[PS_MAX_NR_IIDICC]
int border_position[PS_MAX_NUM_ENV+1]
static void hybrid2_re(float(*in)[2], float(*out)[32][2], const float filter[8], int len, int reverse)
Split one subband into 2 subsubbands with a symmetric real filter.
static float f34_0_12[12][8][2]
common internal and external API header
static const int8_t k_to_i_20[]
Table 8.48.
float H22[2][PS_MAX_NUM_ENV+1][PS_MAX_NR_IIDICC]
static void ipdopd_reset(int8_t *opd_hist, int8_t *ipd_hist)
#define PS_QMF_TIME_SLOTS
static void remap20(int8_t(**p_par_mapped)[PS_MAX_NR_IIDICC], int8_t(*par)[PS_MAX_NR_IIDICC], int num_par, int num_env, int full)
void(* add_squares)(float *dst, const float(*src)[2], int n)
static const int DECAY_CUTOFF[]
Start frequency band for the all-pass filter decay slope.
static int ps_read_extension_data(GetBitContext *gb, PSContext *ps, int ps_extension_id)
static void stereo_processing(PSContext *ps, float(*l)[32][2], float(*r)[32][2], int is34)
int ff_ps_read_data(AVCodecContext *avctx, GetBitContext *gb_host, PSContext *ps, int bits_left)
int8_t opd_par[PS_MAX_NUM_ENV][PS_MAX_NR_IIDICC]
Overall Phase Difference Parameters.
static void map_val_20_to_34(float par[PS_MAX_NR_IIDICC])
float H21[2][PS_MAX_NUM_ENV+1][PS_MAX_NR_IIDICC]