129 float c_re,
float c_im,
130 float mag_totall,
float mag_totalr,
131 float fl_phase,
float fr_phase,
132 float bl_phase,
float br_phase,
133 float sl_phase,
float sr_phase,
138 float c_re,
float c_im,
139 float lfe_re,
float lfe_im,
140 float mag_totall,
float mag_totalr,
141 float fl_phase,
float fr_phase,
142 float bl_phase,
float br_phase,
143 float sl_phase,
float sr_phase,
188 for (
int ch = 0; ch <
s->nb_in_channels &&
s->level_in >= 0.f; ch++)
189 s->input_levels[ch] =
s->level_in;
195 s->input_levels[ch] =
s->f_i[n];
203 for (
int ch = 0; ch <
s->nb_out_channels &&
s->level_out >= 0.f; ch++)
204 s->output_levels[ch] =
s->level_out;
210 s->output_levels[ch] =
s->f_o[n];
223 s->nb_in_channels =
inlink->ch_layout.nb_channels;
225 for (
int ch = 0; ch <
inlink->ch_layout.nb_channels; ch++) {
229 0,
s->win_size, &
scale, 0);
235 if (!
s->input_levels)
252 s->lowcut = 1.f *
s->lowcutf / (
inlink->sample_rate * 0.5) * (
s->win_size / 2);
253 s->highcut = 1.f *
s->highcutf / (
inlink->sample_rate * 0.5) * (
s->win_size / 2);
273 1,
s->win_size, &iscale, 0);
279 if (!
s->output_levels)
291 if (!
s->overlap_buffer || !
s->output || !
s->output_out || !
s->output_mag ||
292 !
s->output_ph || !
s->factors || !
s->sfactors)
295 s->rdft_size =
s->win_size / 2 + 1;
299 s->l_phase =
av_calloc(
s->rdft_size,
sizeof(*
s->l_phase));
300 s->r_phase =
av_calloc(
s->rdft_size,
sizeof(*
s->r_phase));
302 s->c_phase =
av_calloc(
s->rdft_size,
sizeof(*
s->c_phase));
303 s->mag_total =
av_calloc(
s->rdft_size,
sizeof(*
s->mag_total));
304 s->lfe_mag =
av_calloc(
s->rdft_size,
sizeof(*
s->lfe_mag));
305 s->lfe_phase =
av_calloc(
s->rdft_size,
sizeof(*
s->lfe_phase));
306 if (!
s->x_pos || !
s->y_pos || !
s->l_phase || !
s->r_phase || !
s->lfe_phase ||
307 !
s->c_phase || !
s->mag_total || !
s->lfe_mag || !
s->c_mag)
323 #define MIN_MAG_SUM 0.00000001f
327 float reference,
r,
a;
332 reference = angle *
M_PIf / 180.f;
359 r = focus > 0.f ? 1.f -
powf(1.
f -
r, 1.
f + focus * 20.
f) :
powf(
r, 1.
f - focus * 20.
f);
373 static inline void get_lfe(
int output_lfe,
int n,
float lowcut,
float highcut,
374 float *lfe_mag,
float c_mag,
float *mag_total,
int lfe_mode)
376 if (output_lfe && n < highcut) {
377 *lfe_mag = n < lowcut ? 1.f : .5f*(1.f+
cosf(
M_PIf*(lowcut-n)/(lowcut-highcut)));
380 *mag_total -= *lfe_mag;
387 dst[2 * n ] = mag * cosf(ph); \
388 dst[2 * n + 1] = mag * sinf(ph);
393 float *
factor = (
float *)
s->factors->extended_data[ch];
394 const float f_x =
s->f_x[
sc_map[chan >= 0 ? chan : 0]];
395 const float f_y =
s->f_y[
sc_map[chan >= 0 ? chan : 0]];
396 const int rdft_size =
s->rdft_size;
397 const float *x =
s->x_pos;
398 const float *y =
s->y_pos;
402 for (
int n = 0; n < rdft_size; n++)
406 for (
int n = 0; n < rdft_size; n++)
410 for (
int n = 0; n < rdft_size; n++)
414 for (
int n = 0; n < rdft_size; n++)
418 for (
int n = 0; n < rdft_size; n++)
422 for (
int n = 0; n < rdft_size; n++)
426 for (
int n = 0; n < rdft_size; n++)
430 for (
int n = 0; n < rdft_size; n++)
434 for (
int n = 0; n < rdft_size; n++)
438 for (
int n = 0; n < rdft_size; n++)
447 float *sfactor = (
float *)
s->sfactors->extended_data[ch];
448 float *
factor = (
float *)
s->factors->extended_data[ch];
449 float *omag = (
float *)
s->output_mag->extended_data[ch];
450 float *oph = (
float *)
s->output_ph->extended_data[ch];
451 float *dst = (
float *)
s->output->extended_data[ch];
452 const int rdft_size =
s->rdft_size;
453 const float smooth =
s->smooth;
456 for (
int n = 0; n < rdft_size; n++)
462 for (
int n = 0; n < rdft_size; n++)
465 for (
int n = 0; n < rdft_size; n++) {
466 const float mag = omag[n];
467 const float ph = oph[n];
476 float *omag = (
float *)
s->output_mag->extended_data[ch];
477 float *oph = (
float *)
s->output_ph->extended_data[ch];
478 const float *mag_total =
s->mag_total;
479 const int rdft_size =
s->rdft_size;
480 const float *c_phase =
s->c_phase;
481 const float *l_phase =
s->l_phase;
482 const float *r_phase =
s->r_phase;
483 const float *lfe_mag =
s->lfe_mag;
484 const float *c_mag =
s->c_mag;
488 memcpy(omag, c_mag, rdft_size *
sizeof(*omag));
491 memcpy(omag, lfe_mag, rdft_size *
sizeof(*omag));
500 memcpy(omag, mag_total, rdft_size *
sizeof(*omag));
510 memcpy(oph, c_phase, rdft_size *
sizeof(*oph));
515 memcpy(oph, l_phase, rdft_size *
sizeof(*oph));
520 memcpy(oph, r_phase, rdft_size *
sizeof(*oph));
543 float *omag = (
float *)
s->output_mag->extended_data[ch];
544 float *oph = (
float *)
s->output_ph->extended_data[ch];
545 const float *mag_total =
s->mag_total;
546 const float *lfe_phase =
s->lfe_phase;
547 const int rdft_size =
s->rdft_size;
548 const float *c_phase =
s->c_phase;
549 const float *l_phase =
s->l_phase;
550 const float *r_phase =
s->r_phase;
551 const float *lfe_mag =
s->lfe_mag;
552 const float *c_mag =
s->c_mag;
565 memcpy(omag, c_mag, rdft_size *
sizeof(*omag));
568 memcpy(omag, lfe_mag, rdft_size *
sizeof(*omag));
577 memcpy(omag, mag_total, rdft_size *
sizeof(*omag));
585 memcpy(oph, lfe_phase, rdft_size *
sizeof(*oph));
589 memcpy(oph, c_phase, rdft_size *
sizeof(*oph));
594 memcpy(oph, l_phase, rdft_size *
sizeof(*oph));
599 memcpy(oph, r_phase, rdft_size *
sizeof(*oph));
628 float c_re,
float c_im,
629 float mag_totall,
float mag_totalr,
630 float fl_phase,
float fr_phase,
631 float bl_phase,
float br_phase,
632 float sl_phase,
float sr_phase,
637 float fl_mag, fr_mag, ls_mag, rs_mag, lb_mag, rb_mag;
638 float *dstc, *dstl, *dstr, *dstls, *dstrs, *dstlb, *dstrb, *dstlfe;
639 float lfe_mag, c_phase, mag_total = (mag_totall + mag_totalr) * 0.5
f;
642 dstl = (
float *)
s->output->extended_data[0];
643 dstr = (
float *)
s->output->extended_data[1];
644 dstc = (
float *)
s->output->extended_data[2];
645 dstlfe = (
float *)
s->output->extended_data[3];
646 dstlb = (
float *)
s->output->extended_data[4];
647 dstrb = (
float *)
s->output->extended_data[5];
648 dstls = (
float *)
s->output->extended_data[6];
649 dstrs = (
float *)
s->output->extended_data[7];
651 c_phase =
atan2f(c_im, c_re);
653 get_lfe(
s->output_lfe, n,
s->lowcut,
s->highcut, &lfe_mag, hypotf(c_re, c_im), &mag_total,
s->lfe_mode);
662 dstl[2 * n ] = fl_mag *
cosf(fl_phase);
663 dstl[2 * n + 1] = fl_mag *
sinf(fl_phase);
665 dstr[2 * n ] = fr_mag *
cosf(fr_phase);
666 dstr[2 * n + 1] = fr_mag *
sinf(fr_phase);
669 dstc[2 * n + 1] = c_im;
671 dstlfe[2 * n ] = lfe_mag *
cosf(c_phase);
672 dstlfe[2 * n + 1] = lfe_mag *
sinf(c_phase);
674 dstlb[2 * n ] = lb_mag *
cosf(bl_phase);
675 dstlb[2 * n + 1] = lb_mag *
sinf(bl_phase);
677 dstrb[2 * n ] = rb_mag *
cosf(br_phase);
678 dstrb[2 * n + 1] = rb_mag *
sinf(br_phase);
680 dstls[2 * n ] = ls_mag *
cosf(sl_phase);
681 dstls[2 * n + 1] = ls_mag *
sinf(sl_phase);
683 dstrs[2 * n ] = rs_mag *
cosf(sr_phase);
684 dstrs[2 * n + 1] = rs_mag *
sinf(sr_phase);
688 float c_re,
float c_im,
689 float lfe_re,
float lfe_im,
690 float mag_totall,
float mag_totalr,
691 float fl_phase,
float fr_phase,
692 float bl_phase,
float br_phase,
693 float sl_phase,
float sr_phase,
698 float fl_mag, fr_mag, ls_mag, rs_mag, lb_mag, rb_mag;
699 float *dstc, *dstl, *dstr, *dstls, *dstrs, *dstlb, *dstrb, *dstlfe;
702 dstl = (
float *)
s->output->extended_data[0];
703 dstr = (
float *)
s->output->extended_data[1];
704 dstc = (
float *)
s->output->extended_data[2];
705 dstlfe = (
float *)
s->output->extended_data[3];
706 dstlb = (
float *)
s->output->extended_data[4];
707 dstrb = (
float *)
s->output->extended_data[5];
708 dstls = (
float *)
s->output->extended_data[6];
709 dstrs = (
float *)
s->output->extended_data[7];
718 dstl[2 * n ] = fl_mag *
cosf(fl_phase);
719 dstl[2 * n + 1] = fl_mag *
sinf(fl_phase);
721 dstr[2 * n ] = fr_mag *
cosf(fr_phase);
722 dstr[2 * n + 1] = fr_mag *
sinf(fr_phase);
725 dstc[2 * n + 1] = c_im;
727 dstlfe[2 * n ] = lfe_re;
728 dstlfe[2 * n + 1] = lfe_im;
730 dstlb[2 * n ] = lb_mag *
cosf(bl_phase);
731 dstlb[2 * n + 1] = lb_mag *
sinf(bl_phase);
733 dstrb[2 * n ] = rb_mag *
cosf(br_phase);
734 dstrb[2 * n + 1] = rb_mag *
sinf(br_phase);
736 dstls[2 * n ] = ls_mag *
cosf(sl_phase);
737 dstls[2 * n + 1] = ls_mag *
sinf(sl_phase);
739 dstrs[2 * n ] = rs_mag *
cosf(sr_phase);
740 dstrs[2 * n + 1] = rs_mag *
sinf(sr_phase);
746 const float *srcl = (
const float *)
s->input->extended_data[0];
747 const float *srcr = (
const float *)
s->input->extended_data[1];
748 const int output_lfe =
s->output_lfe &&
s->create_lfe;
749 const int rdft_size =
s->rdft_size;
750 const int lfe_mode =
s->lfe_mode;
751 const float highcut =
s->highcut;
752 const float lowcut =
s->lowcut;
753 const float angle =
s->angle;
754 const float focus =
s->focus;
755 float *magtotal =
s->mag_total;
756 float *lfemag =
s->lfe_mag;
757 float *lphase =
s->l_phase;
758 float *rphase =
s->r_phase;
759 float *cphase =
s->c_phase;
760 float *cmag =
s->c_mag;
761 float *xpos =
s->x_pos;
762 float *ypos =
s->y_pos;
764 for (
int n = 0; n < rdft_size; n++) {
765 float l_re = srcl[2 * n], r_re = srcr[2 * n];
766 float l_im = srcl[2 * n + 1], r_im = srcr[2 * n + 1];
767 float c_phase =
atan2f(l_im + r_im, l_re + r_re);
768 float l_mag = hypotf(l_re, l_im);
769 float r_mag = hypotf(r_re, r_im);
770 float mag_total = hypotf(l_mag, r_mag);
771 float l_phase =
atan2f(l_im, l_re);
772 float r_phase =
atan2f(r_im, r_re);
773 float phase_dif =
fabsf(l_phase - r_phase);
774 float mag_sum = l_mag + r_mag;
775 float c_mag = mag_sum * 0.5f;
779 mag_dif = (l_mag - r_mag) / mag_sum;
780 if (phase_dif >
M_PIf)
781 phase_dif = 2.f *
M_PIf - phase_dif;
786 get_lfe(output_lfe, n, lowcut, highcut, &lfemag[n], c_mag, &mag_total, lfe_mode);
794 magtotal[n] = mag_total;
801 const float *srcl = (
const float *)
s->input->extended_data[0];
802 const float *srcr = (
const float *)
s->input->extended_data[1];
803 const float *srclfe = (
const float *)
s->input->extended_data[2];
804 const int rdft_size =
s->rdft_size;
805 const float angle =
s->angle;
806 const float focus =
s->focus;
807 float *magtotal =
s->mag_total;
808 float *lfephase =
s->lfe_phase;
809 float *lfemag =
s->lfe_mag;
810 float *lphase =
s->l_phase;
811 float *rphase =
s->r_phase;
812 float *cphase =
s->c_phase;
813 float *cmag =
s->c_mag;
814 float *xpos =
s->x_pos;
815 float *ypos =
s->y_pos;
817 for (
int n = 0; n < rdft_size; n++) {
818 float l_re = srcl[2 * n], r_re = srcr[2 * n];
819 float l_im = srcl[2 * n + 1], r_im = srcr[2 * n + 1];
820 float lfe_re = srclfe[2 * n], lfe_im = srclfe[2 * n + 1];
821 float c_phase =
atan2f(l_im + r_im, l_re + r_re);
822 float l_mag = hypotf(l_re, l_im);
823 float r_mag = hypotf(r_re, r_im);
824 float lfe_mag = hypotf(lfe_re, lfe_im);
825 float lfe_phase =
atan2f(lfe_im, lfe_re);
826 float mag_total = hypotf(l_mag, r_mag);
827 float l_phase =
atan2f(l_im, l_re);
828 float r_phase =
atan2f(r_im, r_re);
829 float phase_dif =
fabsf(l_phase - r_phase);
830 float mag_sum = l_mag + r_mag;
831 float c_mag = mag_sum * 0.5f;
835 mag_dif = (l_mag - r_mag) / mag_sum;
836 if (phase_dif >
M_PIf)
837 phase_dif = 2.f *
M_PIf - phase_dif;
850 lfephase[n] = lfe_phase;
851 magtotal[n] = mag_total;
858 const float *srcl = (
const float *)
s->input->extended_data[0];
859 const float *srcr = (
const float *)
s->input->extended_data[1];
860 const float *srcc = (
const float *)
s->input->extended_data[2];
861 const int output_lfe =
s->output_lfe &&
s->create_lfe;
862 const int rdft_size =
s->rdft_size;
863 const int lfe_mode =
s->lfe_mode;
864 const float highcut =
s->highcut;
865 const float lowcut =
s->lowcut;
866 const float angle =
s->angle;
867 const float focus =
s->focus;
868 float *magtotal =
s->mag_total;
869 float *lfemag =
s->lfe_mag;
870 float *lphase =
s->l_phase;
871 float *rphase =
s->r_phase;
872 float *cphase =
s->c_phase;
873 float *cmag =
s->c_mag;
874 float *xpos =
s->x_pos;
875 float *ypos =
s->y_pos;
877 for (
int n = 0; n < rdft_size; n++) {
878 float l_re = srcl[2 * n], r_re = srcr[2 * n];
879 float l_im = srcl[2 * n + 1], r_im = srcr[2 * n + 1];
880 float c_re = srcc[2 * n], c_im = srcc[2 * n + 1];
881 float c_phase =
atan2f(c_im, c_re);
882 float c_mag = hypotf(c_re, c_im);
883 float l_mag = hypotf(l_re, l_im);
884 float r_mag = hypotf(r_re, r_im);
885 float mag_total = hypotf(l_mag, r_mag);
886 float l_phase =
atan2f(l_im, l_re);
887 float r_phase =
atan2f(r_im, r_re);
888 float phase_dif =
fabsf(l_phase - r_phase);
889 float mag_sum = l_mag + r_mag;
893 mag_dif = (l_mag - r_mag) / mag_sum;
894 if (phase_dif >
M_PIf)
895 phase_dif = 2.f *
M_PIf - phase_dif;
900 get_lfe(output_lfe, n, lowcut, highcut, &lfemag[n], c_mag, &mag_total, lfe_mode);
908 magtotal[n] = mag_total;
915 const int rdft_size =
s->rdft_size;
916 float *srcl, *srcr, *srcc, *srcsl, *srcsr;
919 srcl = (
float *)
s->input->extended_data[0];
920 srcr = (
float *)
s->input->extended_data[1];
921 srcc = (
float *)
s->input->extended_data[2];
922 srcsl = (
float *)
s->input->extended_data[3];
923 srcsr = (
float *)
s->input->extended_data[4];
925 for (n = 0; n < rdft_size; n++) {
926 float fl_re = srcl[2 * n], fr_re = srcr[2 * n];
927 float fl_im = srcl[2 * n + 1], fr_im = srcr[2 * n + 1];
928 float c_re = srcc[2 * n], c_im = srcc[2 * n + 1];
929 float sl_re = srcsl[2 * n], sl_im = srcsl[2 * n + 1];
930 float sr_re = srcsr[2 * n], sr_im = srcsr[2 * n + 1];
931 float fl_mag = hypotf(fl_re, fl_im);
932 float fr_mag = hypotf(fr_re, fr_im);
933 float fl_phase =
atan2f(fl_im, fl_re);
934 float fr_phase =
atan2f(fr_im, fr_re);
935 float sl_mag = hypotf(sl_re, sl_im);
936 float sr_mag = hypotf(sr_re, sr_im);
937 float sl_phase =
atan2f(sl_im, sl_re);
938 float sr_phase =
atan2f(sr_im, sr_re);
939 float phase_difl =
fabsf(fl_phase - sl_phase);
940 float phase_difr =
fabsf(fr_phase - sr_phase);
941 float magl_sum = fl_mag + sl_mag;
942 float magr_sum = fr_mag + sr_mag;
945 float mag_totall = hypotf(fl_mag, sl_mag);
946 float mag_totalr = hypotf(fr_mag, sr_mag);
947 float bl_phase =
atan2f(fl_im + sl_im, fl_re + sl_re);
948 float br_phase =
atan2f(fr_im + sr_im, fr_re + sr_re);
952 if (phase_difl >
M_PIf)
953 phase_difl = 2.f *
M_PIf - phase_difl;
955 if (phase_difr >
M_PIf)
956 phase_difr = 2.f *
M_PIf - phase_difr;
961 s->upmix_5_0(
ctx, c_re, c_im,
962 mag_totall, mag_totalr,
973 const int rdft_size =
s->rdft_size;
974 float *srcl, *srcr, *srcc, *srclfe, *srcsl, *srcsr;
977 srcl = (
float *)
s->input->extended_data[0];
978 srcr = (
float *)
s->input->extended_data[1];
979 srcc = (
float *)
s->input->extended_data[2];
980 srclfe = (
float *)
s->input->extended_data[3];
981 srcsl = (
float *)
s->input->extended_data[4];
982 srcsr = (
float *)
s->input->extended_data[5];
984 for (n = 0; n < rdft_size; n++) {
985 float fl_re = srcl[2 * n], fr_re = srcr[2 * n];
986 float fl_im = srcl[2 * n + 1], fr_im = srcr[2 * n + 1];
987 float c_re = srcc[2 * n], c_im = srcc[2 * n + 1];
988 float lfe_re = srclfe[2 * n], lfe_im = srclfe[2 * n + 1];
989 float sl_re = srcsl[2 * n], sl_im = srcsl[2 * n + 1];
990 float sr_re = srcsr[2 * n], sr_im = srcsr[2 * n + 1];
991 float fl_mag = hypotf(fl_re, fl_im);
992 float fr_mag = hypotf(fr_re, fr_im);
993 float fl_phase =
atan2f(fl_im, fl_re);
994 float fr_phase =
atan2f(fr_im, fr_re);
995 float sl_mag = hypotf(sl_re, sl_im);
996 float sr_mag = hypotf(sr_re, sr_im);
997 float sl_phase =
atan2f(sl_im, sl_re);
998 float sr_phase =
atan2f(sr_im, sr_re);
999 float phase_difl =
fabsf(fl_phase - sl_phase);
1000 float phase_difr =
fabsf(fr_phase - sr_phase);
1001 float magl_sum = fl_mag + sl_mag;
1002 float magr_sum = fr_mag + sr_mag;
1003 float mag_difl = magl_sum <
MIN_MAG_SUM ?
FFDIFFSIGN(fl_mag, sl_mag) : (fl_mag - sl_mag) / magl_sum;
1004 float mag_difr = magr_sum <
MIN_MAG_SUM ?
FFDIFFSIGN(fr_mag, sr_mag) : (fr_mag - sr_mag) / magr_sum;
1005 float mag_totall = hypotf(fl_mag, sl_mag);
1006 float mag_totalr = hypotf(fr_mag, sr_mag);
1007 float bl_phase =
atan2f(fl_im + sl_im, fl_re + sl_re);
1008 float br_phase =
atan2f(fr_im + sr_im, fr_re + sr_re);
1012 if (phase_difl >
M_PIf)
1013 phase_difl = 2.f *
M_PIf - phase_difl;
1015 if (phase_difr >
M_PIf)
1016 phase_difr = 2.f *
M_PIf - phase_difr;
1021 s->upmix_5_1(
ctx, c_re, c_im, lfe_re, lfe_im,
1022 mag_totall, mag_totalr,
1033 const int rdft_size =
s->rdft_size;
1034 float *srcl, *srcr, *srcc, *srclfe, *srcbl, *srcbr;
1037 srcl = (
float *)
s->input->extended_data[0];
1038 srcr = (
float *)
s->input->extended_data[1];
1039 srcc = (
float *)
s->input->extended_data[2];
1040 srclfe = (
float *)
s->input->extended_data[3];
1041 srcbl = (
float *)
s->input->extended_data[4];
1042 srcbr = (
float *)
s->input->extended_data[5];
1044 for (n = 0; n < rdft_size; n++) {
1045 float fl_re = srcl[2 * n], fr_re = srcr[2 * n];
1046 float fl_im = srcl[2 * n + 1], fr_im = srcr[2 * n + 1];
1047 float c_re = srcc[2 * n], c_im = srcc[2 * n + 1];
1048 float lfe_re = srclfe[2 * n], lfe_im = srclfe[2 * n + 1];
1049 float bl_re = srcbl[2 * n], bl_im = srcbl[2 * n + 1];
1050 float br_re = srcbr[2 * n], br_im = srcbr[2 * n + 1];
1051 float fl_mag = hypotf(fl_re, fl_im);
1052 float fr_mag = hypotf(fr_re, fr_im);
1053 float fl_phase =
atan2f(fl_im, fl_re);
1054 float fr_phase =
atan2f(fr_im, fr_re);
1055 float bl_mag = hypotf(bl_re, bl_im);
1056 float br_mag = hypotf(br_re, br_im);
1057 float bl_phase =
atan2f(bl_im, bl_re);
1058 float br_phase =
atan2f(br_im, br_re);
1059 float phase_difl =
fabsf(fl_phase - bl_phase);
1060 float phase_difr =
fabsf(fr_phase - br_phase);
1061 float magl_sum = fl_mag + bl_mag;
1062 float magr_sum = fr_mag + br_mag;
1063 float mag_difl = magl_sum <
MIN_MAG_SUM ?
FFDIFFSIGN(fl_mag, bl_mag) : (fl_mag - bl_mag) / magl_sum;
1064 float mag_difr = magr_sum <
MIN_MAG_SUM ?
FFDIFFSIGN(fr_mag, br_mag) : (fr_mag - br_mag) / magr_sum;
1065 float mag_totall = hypotf(fl_mag, bl_mag);
1066 float mag_totalr = hypotf(fr_mag, br_mag);
1067 float sl_phase =
atan2f(fl_im + bl_im, fl_re + bl_re);
1068 float sr_phase =
atan2f(fr_im + br_im, fr_re + br_re);
1072 if (phase_difl >
M_PIf)
1073 phase_difl = 2.f *
M_PIf - phase_difl;
1075 if (phase_difr >
M_PIf)
1076 phase_difr = 2.f *
M_PIf - phase_difr;
1081 s->upmix_5_1(
ctx, c_re, c_im, lfe_re, lfe_im,
1082 mag_totall, mag_totalr,
1094 if (
s->all_x >= 0.f)
1095 for (
int n = 0; n <
SC_NB; n++)
1096 s->f_x[n] =
s->all_x;
1098 if (
s->all_y >= 0.f)
1099 for (
int n = 0; n <
SC_NB; n++)
1100 s->f_y[n] =
s->all_y;
1107 int64_t in_channel_layout, out_channel_layout;
1108 char in_name[128], out_name[128];
1111 if (
s->lowcutf >=
s->highcutf) {
1113 s->lowcutf,
s->highcutf);
1118 s->in_ch_layout.u.mask : 0;
1120 s->out_ch_layout.u.mask : 0;
1125 switch (in_channel_layout) {
1140 switch (out_channel_layout) {
1150 switch (out_channel_layout) {
1160 switch (out_channel_layout) {
1177 s->window_func_lut =
av_calloc(
s->win_size,
sizeof(*
s->window_func_lut));
1178 if (!
s->window_func_lut)
1182 if (
s->overlap == 1)
1183 s->overlap = overlap;
1185 for (
int i = 0;
i <
s->win_size;
i++)
1186 s->window_func_lut[
i] =
sqrtf(
s->window_func_lut[
i] /
s->win_size);
1187 s->hop_size =
FFMAX(1,
s->win_size * (1. -
s->overlap));
1190 float max = 0.f, *temp_lut =
av_calloc(
s->win_size,
sizeof(*temp_lut));
1194 for (
int j = 0; j <
s->win_size; j +=
s->hop_size) {
1195 for (
int i = 0;
i <
s->win_size;
i++)
1196 temp_lut[(
i + j) %
s->win_size] +=
s->window_func_lut[
i];
1199 for (
int i = 0;
i <
s->win_size;
i++)
1203 s->win_gain = 1.f / (
max *
sqrtf(
s->win_size));
1214 float *
src = (
float *)
s->input_in->extended_data[ch];
1215 float *
win = (
float *)
s->window->extended_data[ch];
1216 const float *window_func_lut =
s->window_func_lut;
1217 const int offset =
s->win_size -
s->hop_size;
1218 const float level_in =
s->input_levels[ch];
1219 const int win_size =
s->win_size;
1225 for (
int n = 0; n < win_size; n++)
1226 win[n] =
src[n] * window_func_lut[n] * level_in;
1228 s->tx_fn(
s->rdft[ch], (
float *)
s->input->extended_data[ch],
win,
sizeof(
float));
1239 for (
int ch = start; ch < end; ch++)
1248 const float level_out =
s->output_levels[ch] *
s->win_gain;
1249 const float *window_func_lut =
s->window_func_lut;
1250 const int win_size =
s->win_size;
1253 dst = (
float *)
s->output_out->extended_data[ch];
1254 ptr = (
float *)
s->overlap_buffer->extended_data[ch];
1255 s->itx_fn(
s->irdft[ch], dst, (
float *)
s->output->extended_data[ch],
sizeof(
AVComplexFloat));
1257 memmove(
s->overlap_buffer->extended_data[ch],
1258 s->overlap_buffer->extended_data[ch] +
s->hop_size *
sizeof(
float),
1259 s->win_size *
sizeof(
float));
1260 memset(
s->overlap_buffer->extended_data[ch] +
s->win_size *
sizeof(
float),
1261 0,
s->hop_size *
sizeof(
float));
1263 for (
int n = 0; n < win_size; n++)
1264 ptr[n] += dst[n] * window_func_lut[n] * level_out;
1266 ptr = (
float *)
s->overlap_buffer->extended_data[ch];
1267 dst = (
float *)
out->extended_data[ch];
1268 memcpy(dst, ptr,
s->hop_size *
sizeof(
float));
1277 const int start = (
out->ch_layout.nb_channels * jobnr) / nb_jobs;
1278 const int end = (
out->ch_layout.nb_channels * (jobnr+1)) / nb_jobs;
1280 for (
int ch = start; ch < end; ch++) {
1367 for (
int ch = 0; ch <
s->nb_in_channels; ch++)
1369 for (
int ch = 0; ch <
s->nb_out_channels; ch++)
1389 char *res,
int res_len,
int flags)
1398 s->hop_size =
FFMAX(1,
s->win_size * (1. -
s->overlap));
1407 #define OFFSET(x) offsetof(AudioSurroundContext, x)
1408 #define FLAGS AV_OPT_FLAG_AUDIO_PARAM|AV_OPT_FLAG_FILTERING_PARAM
1409 #define TFLAGS AV_OPT_FLAG_AUDIO_PARAM|AV_OPT_FLAG_FILTERING_PARAM|AV_OPT_FLAG_RUNTIME_PARAM
1421 {
"sub",
"subtract LFE channel with others", 0,
AV_OPT_TYPE_CONST, {.i64=1}, 0, 1,
TFLAGS, .unit =
"lfe_mode" },
1489 .priv_class = &surround_class,