130 float c_re,
float c_im,
131 float mag_totall,
float mag_totalr,
132 float fl_phase,
float fr_phase,
133 float bl_phase,
float br_phase,
134 float sl_phase,
float sr_phase,
139 float c_re,
float c_im,
140 float lfe_re,
float lfe_im,
141 float mag_totall,
float mag_totalr,
142 float fl_phase,
float fr_phase,
143 float bl_phase,
float br_phase,
144 float sl_phase,
float sr_phase,
189 for (
int ch = 0; ch <
s->nb_in_channels &&
s->level_in >= 0.f; ch++)
190 s->input_levels[ch] =
s->level_in;
196 s->input_levels[ch] =
s->f_i[n];
204 for (
int ch = 0; ch <
s->nb_out_channels &&
s->level_out >= 0.f; ch++)
205 s->output_levels[ch] =
s->level_out;
211 s->output_levels[ch] =
s->f_o[n];
224 s->nb_in_channels =
inlink->ch_layout.nb_channels;
226 for (
int ch = 0; ch <
inlink->ch_layout.nb_channels; ch++) {
230 0,
s->win_size, &
scale, 0);
236 if (!
s->input_levels)
253 s->lowcut = 1.f *
s->lowcutf / (
inlink->sample_rate * 0.5) * (
s->win_size / 2);
254 s->highcut = 1.f *
s->highcutf / (
inlink->sample_rate * 0.5) * (
s->win_size / 2);
274 1,
s->win_size, &iscale, 0);
280 if (!
s->output_levels)
292 if (!
s->overlap_buffer || !
s->output || !
s->output_out || !
s->output_mag ||
293 !
s->output_ph || !
s->factors || !
s->sfactors)
296 s->rdft_size =
s->win_size / 2 + 1;
300 s->l_phase =
av_calloc(
s->rdft_size,
sizeof(*
s->l_phase));
301 s->r_phase =
av_calloc(
s->rdft_size,
sizeof(*
s->r_phase));
303 s->c_phase =
av_calloc(
s->rdft_size,
sizeof(*
s->c_phase));
304 s->mag_total =
av_calloc(
s->rdft_size,
sizeof(*
s->mag_total));
305 s->lfe_mag =
av_calloc(
s->rdft_size,
sizeof(*
s->lfe_mag));
306 s->lfe_phase =
av_calloc(
s->rdft_size,
sizeof(*
s->lfe_phase));
307 if (!
s->x_pos || !
s->y_pos || !
s->l_phase || !
s->r_phase || !
s->lfe_phase ||
308 !
s->c_phase || !
s->mag_total || !
s->lfe_mag || !
s->c_mag)
324 #define MIN_MAG_SUM 0.00000001f
328 float reference,
r,
a;
333 reference = angle *
M_PIf / 180.f;
360 r = focus > 0.f ? 1.f -
powf(1.
f -
r, 1.
f + focus * 20.
f) :
powf(
r, 1.
f - focus * 20.
f);
374 static inline void get_lfe(
int output_lfe,
int n,
float lowcut,
float highcut,
375 float *lfe_mag,
float c_mag,
float *mag_total,
int lfe_mode)
377 if (output_lfe && n < highcut) {
378 *lfe_mag = n < lowcut ? 1.f : .5f*(1.f+
cosf(
M_PIf*(lowcut-n)/(lowcut-highcut)));
381 *mag_total -= *lfe_mag;
388 dst[2 * n ] = mag * cosf(ph); \
389 dst[2 * n + 1] = mag * sinf(ph);
394 float *
factor = (
float *)
s->factors->extended_data[ch];
395 const float f_x =
s->f_x[
sc_map[chan >= 0 ? chan : 0]];
396 const float f_y =
s->f_y[
sc_map[chan >= 0 ? chan : 0]];
397 const int rdft_size =
s->rdft_size;
398 const float *x =
s->x_pos;
399 const float *y =
s->y_pos;
403 for (
int n = 0; n < rdft_size; n++)
407 for (
int n = 0; n < rdft_size; n++)
411 for (
int n = 0; n < rdft_size; n++)
415 for (
int n = 0; n < rdft_size; n++)
419 for (
int n = 0; n < rdft_size; n++)
423 for (
int n = 0; n < rdft_size; n++)
427 for (
int n = 0; n < rdft_size; n++)
431 for (
int n = 0; n < rdft_size; n++)
435 for (
int n = 0; n < rdft_size; n++)
439 for (
int n = 0; n < rdft_size; n++)
448 float *sfactor = (
float *)
s->sfactors->extended_data[ch];
449 float *
factor = (
float *)
s->factors->extended_data[ch];
450 float *omag = (
float *)
s->output_mag->extended_data[ch];
451 float *oph = (
float *)
s->output_ph->extended_data[ch];
452 float *dst = (
float *)
s->output->extended_data[ch];
453 const int rdft_size =
s->rdft_size;
454 const float smooth =
s->smooth;
457 for (
int n = 0; n < rdft_size; n++)
463 for (
int n = 0; n < rdft_size; n++)
466 for (
int n = 0; n < rdft_size; n++) {
467 const float mag = omag[n];
468 const float ph = oph[n];
477 float *omag = (
float *)
s->output_mag->extended_data[ch];
478 float *oph = (
float *)
s->output_ph->extended_data[ch];
479 const float *mag_total =
s->mag_total;
480 const int rdft_size =
s->rdft_size;
481 const float *c_phase =
s->c_phase;
482 const float *l_phase =
s->l_phase;
483 const float *r_phase =
s->r_phase;
484 const float *lfe_mag =
s->lfe_mag;
485 const float *c_mag =
s->c_mag;
489 memcpy(omag, c_mag, rdft_size *
sizeof(*omag));
492 memcpy(omag, lfe_mag, rdft_size *
sizeof(*omag));
501 memcpy(omag, mag_total, rdft_size *
sizeof(*omag));
511 memcpy(oph, c_phase, rdft_size *
sizeof(*oph));
516 memcpy(oph, l_phase, rdft_size *
sizeof(*oph));
521 memcpy(oph, r_phase, rdft_size *
sizeof(*oph));
544 float *omag = (
float *)
s->output_mag->extended_data[ch];
545 float *oph = (
float *)
s->output_ph->extended_data[ch];
546 const float *mag_total =
s->mag_total;
547 const float *lfe_phase =
s->lfe_phase;
548 const int rdft_size =
s->rdft_size;
549 const float *c_phase =
s->c_phase;
550 const float *l_phase =
s->l_phase;
551 const float *r_phase =
s->r_phase;
552 const float *lfe_mag =
s->lfe_mag;
553 const float *c_mag =
s->c_mag;
566 memcpy(omag, c_mag, rdft_size *
sizeof(*omag));
569 memcpy(omag, lfe_mag, rdft_size *
sizeof(*omag));
578 memcpy(omag, mag_total, rdft_size *
sizeof(*omag));
586 memcpy(oph, lfe_phase, rdft_size *
sizeof(*oph));
590 memcpy(oph, c_phase, rdft_size *
sizeof(*oph));
595 memcpy(oph, l_phase, rdft_size *
sizeof(*oph));
600 memcpy(oph, r_phase, rdft_size *
sizeof(*oph));
629 float c_re,
float c_im,
630 float mag_totall,
float mag_totalr,
631 float fl_phase,
float fr_phase,
632 float bl_phase,
float br_phase,
633 float sl_phase,
float sr_phase,
638 float fl_mag, fr_mag, ls_mag, rs_mag, lb_mag, rb_mag;
639 float *dstc, *dstl, *dstr, *dstls, *dstrs, *dstlb, *dstrb, *dstlfe;
640 float lfe_mag, c_phase, mag_total = (mag_totall + mag_totalr) * 0.5
f;
643 dstl = (
float *)
s->output->extended_data[0];
644 dstr = (
float *)
s->output->extended_data[1];
645 dstc = (
float *)
s->output->extended_data[2];
646 dstlfe = (
float *)
s->output->extended_data[3];
647 dstlb = (
float *)
s->output->extended_data[4];
648 dstrb = (
float *)
s->output->extended_data[5];
649 dstls = (
float *)
s->output->extended_data[6];
650 dstrs = (
float *)
s->output->extended_data[7];
652 c_phase =
atan2f(c_im, c_re);
654 get_lfe(
s->output_lfe, n,
s->lowcut,
s->highcut, &lfe_mag, hypotf(c_re, c_im), &mag_total,
s->lfe_mode);
663 dstl[2 * n ] = fl_mag *
cosf(fl_phase);
664 dstl[2 * n + 1] = fl_mag *
sinf(fl_phase);
666 dstr[2 * n ] = fr_mag *
cosf(fr_phase);
667 dstr[2 * n + 1] = fr_mag *
sinf(fr_phase);
670 dstc[2 * n + 1] = c_im;
672 dstlfe[2 * n ] = lfe_mag *
cosf(c_phase);
673 dstlfe[2 * n + 1] = lfe_mag *
sinf(c_phase);
675 dstlb[2 * n ] = lb_mag *
cosf(bl_phase);
676 dstlb[2 * n + 1] = lb_mag *
sinf(bl_phase);
678 dstrb[2 * n ] = rb_mag *
cosf(br_phase);
679 dstrb[2 * n + 1] = rb_mag *
sinf(br_phase);
681 dstls[2 * n ] = ls_mag *
cosf(sl_phase);
682 dstls[2 * n + 1] = ls_mag *
sinf(sl_phase);
684 dstrs[2 * n ] = rs_mag *
cosf(sr_phase);
685 dstrs[2 * n + 1] = rs_mag *
sinf(sr_phase);
689 float c_re,
float c_im,
690 float lfe_re,
float lfe_im,
691 float mag_totall,
float mag_totalr,
692 float fl_phase,
float fr_phase,
693 float bl_phase,
float br_phase,
694 float sl_phase,
float sr_phase,
699 float fl_mag, fr_mag, ls_mag, rs_mag, lb_mag, rb_mag;
700 float *dstc, *dstl, *dstr, *dstls, *dstrs, *dstlb, *dstrb, *dstlfe;
703 dstl = (
float *)
s->output->extended_data[0];
704 dstr = (
float *)
s->output->extended_data[1];
705 dstc = (
float *)
s->output->extended_data[2];
706 dstlfe = (
float *)
s->output->extended_data[3];
707 dstlb = (
float *)
s->output->extended_data[4];
708 dstrb = (
float *)
s->output->extended_data[5];
709 dstls = (
float *)
s->output->extended_data[6];
710 dstrs = (
float *)
s->output->extended_data[7];
719 dstl[2 * n ] = fl_mag *
cosf(fl_phase);
720 dstl[2 * n + 1] = fl_mag *
sinf(fl_phase);
722 dstr[2 * n ] = fr_mag *
cosf(fr_phase);
723 dstr[2 * n + 1] = fr_mag *
sinf(fr_phase);
726 dstc[2 * n + 1] = c_im;
728 dstlfe[2 * n ] = lfe_re;
729 dstlfe[2 * n + 1] = lfe_im;
731 dstlb[2 * n ] = lb_mag *
cosf(bl_phase);
732 dstlb[2 * n + 1] = lb_mag *
sinf(bl_phase);
734 dstrb[2 * n ] = rb_mag *
cosf(br_phase);
735 dstrb[2 * n + 1] = rb_mag *
sinf(br_phase);
737 dstls[2 * n ] = ls_mag *
cosf(sl_phase);
738 dstls[2 * n + 1] = ls_mag *
sinf(sl_phase);
740 dstrs[2 * n ] = rs_mag *
cosf(sr_phase);
741 dstrs[2 * n + 1] = rs_mag *
sinf(sr_phase);
747 const float *srcl = (
const float *)
s->input->extended_data[0];
748 const float *srcr = (
const float *)
s->input->extended_data[1];
749 const int output_lfe =
s->output_lfe &&
s->create_lfe;
750 const int rdft_size =
s->rdft_size;
751 const int lfe_mode =
s->lfe_mode;
752 const float highcut =
s->highcut;
753 const float lowcut =
s->lowcut;
754 const float angle =
s->angle;
755 const float focus =
s->focus;
756 float *magtotal =
s->mag_total;
757 float *lfemag =
s->lfe_mag;
758 float *lphase =
s->l_phase;
759 float *rphase =
s->r_phase;
760 float *cphase =
s->c_phase;
761 float *cmag =
s->c_mag;
762 float *xpos =
s->x_pos;
763 float *ypos =
s->y_pos;
765 for (
int n = 0; n < rdft_size; n++) {
766 float l_re = srcl[2 * n], r_re = srcr[2 * n];
767 float l_im = srcl[2 * n + 1], r_im = srcr[2 * n + 1];
768 float c_phase =
atan2f(l_im + r_im, l_re + r_re);
769 float l_mag = hypotf(l_re, l_im);
770 float r_mag = hypotf(r_re, r_im);
771 float mag_total = hypotf(l_mag, r_mag);
772 float l_phase =
atan2f(l_im, l_re);
773 float r_phase =
atan2f(r_im, r_re);
774 float phase_dif =
fabsf(l_phase - r_phase);
775 float mag_sum = l_mag + r_mag;
776 float c_mag = mag_sum * 0.5f;
780 mag_dif = (l_mag - r_mag) / mag_sum;
781 if (phase_dif >
M_PIf)
782 phase_dif = 2.f *
M_PIf - phase_dif;
787 get_lfe(output_lfe, n, lowcut, highcut, &lfemag[n], c_mag, &mag_total, lfe_mode);
795 magtotal[n] = mag_total;
802 const float *srcl = (
const float *)
s->input->extended_data[0];
803 const float *srcr = (
const float *)
s->input->extended_data[1];
804 const float *srclfe = (
const float *)
s->input->extended_data[2];
805 const int rdft_size =
s->rdft_size;
806 const float angle =
s->angle;
807 const float focus =
s->focus;
808 float *magtotal =
s->mag_total;
809 float *lfephase =
s->lfe_phase;
810 float *lfemag =
s->lfe_mag;
811 float *lphase =
s->l_phase;
812 float *rphase =
s->r_phase;
813 float *cphase =
s->c_phase;
814 float *cmag =
s->c_mag;
815 float *xpos =
s->x_pos;
816 float *ypos =
s->y_pos;
818 for (
int n = 0; n < rdft_size; n++) {
819 float l_re = srcl[2 * n], r_re = srcr[2 * n];
820 float l_im = srcl[2 * n + 1], r_im = srcr[2 * n + 1];
821 float lfe_re = srclfe[2 * n], lfe_im = srclfe[2 * n + 1];
822 float c_phase =
atan2f(l_im + r_im, l_re + r_re);
823 float l_mag = hypotf(l_re, l_im);
824 float r_mag = hypotf(r_re, r_im);
825 float lfe_mag = hypotf(lfe_re, lfe_im);
826 float lfe_phase =
atan2f(lfe_im, lfe_re);
827 float mag_total = hypotf(l_mag, r_mag);
828 float l_phase =
atan2f(l_im, l_re);
829 float r_phase =
atan2f(r_im, r_re);
830 float phase_dif =
fabsf(l_phase - r_phase);
831 float mag_sum = l_mag + r_mag;
832 float c_mag = mag_sum * 0.5f;
836 mag_dif = (l_mag - r_mag) / mag_sum;
837 if (phase_dif >
M_PIf)
838 phase_dif = 2.f *
M_PIf - phase_dif;
851 lfephase[n] = lfe_phase;
852 magtotal[n] = mag_total;
859 const float *srcl = (
const float *)
s->input->extended_data[0];
860 const float *srcr = (
const float *)
s->input->extended_data[1];
861 const float *srcc = (
const float *)
s->input->extended_data[2];
862 const int output_lfe =
s->output_lfe &&
s->create_lfe;
863 const int rdft_size =
s->rdft_size;
864 const int lfe_mode =
s->lfe_mode;
865 const float highcut =
s->highcut;
866 const float lowcut =
s->lowcut;
867 const float angle =
s->angle;
868 const float focus =
s->focus;
869 float *magtotal =
s->mag_total;
870 float *lfemag =
s->lfe_mag;
871 float *lphase =
s->l_phase;
872 float *rphase =
s->r_phase;
873 float *cphase =
s->c_phase;
874 float *cmag =
s->c_mag;
875 float *xpos =
s->x_pos;
876 float *ypos =
s->y_pos;
878 for (
int n = 0; n < rdft_size; n++) {
879 float l_re = srcl[2 * n], r_re = srcr[2 * n];
880 float l_im = srcl[2 * n + 1], r_im = srcr[2 * n + 1];
881 float c_re = srcc[2 * n], c_im = srcc[2 * n + 1];
882 float c_phase =
atan2f(c_im, c_re);
883 float c_mag = hypotf(c_re, c_im);
884 float l_mag = hypotf(l_re, l_im);
885 float r_mag = hypotf(r_re, r_im);
886 float mag_total = hypotf(l_mag, r_mag);
887 float l_phase =
atan2f(l_im, l_re);
888 float r_phase =
atan2f(r_im, r_re);
889 float phase_dif =
fabsf(l_phase - r_phase);
890 float mag_sum = l_mag + r_mag;
894 mag_dif = (l_mag - r_mag) / mag_sum;
895 if (phase_dif >
M_PIf)
896 phase_dif = 2.f *
M_PIf - phase_dif;
901 get_lfe(output_lfe, n, lowcut, highcut, &lfemag[n], c_mag, &mag_total, lfe_mode);
909 magtotal[n] = mag_total;
916 const int rdft_size =
s->rdft_size;
917 float *srcl, *srcr, *srcc, *srcsl, *srcsr;
920 srcl = (
float *)
s->input->extended_data[0];
921 srcr = (
float *)
s->input->extended_data[1];
922 srcc = (
float *)
s->input->extended_data[2];
923 srcsl = (
float *)
s->input->extended_data[3];
924 srcsr = (
float *)
s->input->extended_data[4];
926 for (n = 0; n < rdft_size; n++) {
927 float fl_re = srcl[2 * n], fr_re = srcr[2 * n];
928 float fl_im = srcl[2 * n + 1], fr_im = srcr[2 * n + 1];
929 float c_re = srcc[2 * n], c_im = srcc[2 * n + 1];
930 float sl_re = srcsl[2 * n], sl_im = srcsl[2 * n + 1];
931 float sr_re = srcsr[2 * n], sr_im = srcsr[2 * n + 1];
932 float fl_mag = hypotf(fl_re, fl_im);
933 float fr_mag = hypotf(fr_re, fr_im);
934 float fl_phase =
atan2f(fl_im, fl_re);
935 float fr_phase =
atan2f(fr_im, fr_re);
936 float sl_mag = hypotf(sl_re, sl_im);
937 float sr_mag = hypotf(sr_re, sr_im);
938 float sl_phase =
atan2f(sl_im, sl_re);
939 float sr_phase =
atan2f(sr_im, sr_re);
940 float phase_difl =
fabsf(fl_phase - sl_phase);
941 float phase_difr =
fabsf(fr_phase - sr_phase);
942 float magl_sum = fl_mag + sl_mag;
943 float magr_sum = fr_mag + sr_mag;
946 float mag_totall = hypotf(fl_mag, sl_mag);
947 float mag_totalr = hypotf(fr_mag, sr_mag);
948 float bl_phase =
atan2f(fl_im + sl_im, fl_re + sl_re);
949 float br_phase =
atan2f(fr_im + sr_im, fr_re + sr_re);
953 if (phase_difl >
M_PIf)
954 phase_difl = 2.f *
M_PIf - phase_difl;
956 if (phase_difr >
M_PIf)
957 phase_difr = 2.f *
M_PIf - phase_difr;
962 s->upmix_5_0(
ctx, c_re, c_im,
963 mag_totall, mag_totalr,
974 const int rdft_size =
s->rdft_size;
975 float *srcl, *srcr, *srcc, *srclfe, *srcsl, *srcsr;
978 srcl = (
float *)
s->input->extended_data[0];
979 srcr = (
float *)
s->input->extended_data[1];
980 srcc = (
float *)
s->input->extended_data[2];
981 srclfe = (
float *)
s->input->extended_data[3];
982 srcsl = (
float *)
s->input->extended_data[4];
983 srcsr = (
float *)
s->input->extended_data[5];
985 for (n = 0; n < rdft_size; n++) {
986 float fl_re = srcl[2 * n], fr_re = srcr[2 * n];
987 float fl_im = srcl[2 * n + 1], fr_im = srcr[2 * n + 1];
988 float c_re = srcc[2 * n], c_im = srcc[2 * n + 1];
989 float lfe_re = srclfe[2 * n], lfe_im = srclfe[2 * n + 1];
990 float sl_re = srcsl[2 * n], sl_im = srcsl[2 * n + 1];
991 float sr_re = srcsr[2 * n], sr_im = srcsr[2 * n + 1];
992 float fl_mag = hypotf(fl_re, fl_im);
993 float fr_mag = hypotf(fr_re, fr_im);
994 float fl_phase =
atan2f(fl_im, fl_re);
995 float fr_phase =
atan2f(fr_im, fr_re);
996 float sl_mag = hypotf(sl_re, sl_im);
997 float sr_mag = hypotf(sr_re, sr_im);
998 float sl_phase =
atan2f(sl_im, sl_re);
999 float sr_phase =
atan2f(sr_im, sr_re);
1000 float phase_difl =
fabsf(fl_phase - sl_phase);
1001 float phase_difr =
fabsf(fr_phase - sr_phase);
1002 float magl_sum = fl_mag + sl_mag;
1003 float magr_sum = fr_mag + sr_mag;
1004 float mag_difl = magl_sum <
MIN_MAG_SUM ?
FFDIFFSIGN(fl_mag, sl_mag) : (fl_mag - sl_mag) / magl_sum;
1005 float mag_difr = magr_sum <
MIN_MAG_SUM ?
FFDIFFSIGN(fr_mag, sr_mag) : (fr_mag - sr_mag) / magr_sum;
1006 float mag_totall = hypotf(fl_mag, sl_mag);
1007 float mag_totalr = hypotf(fr_mag, sr_mag);
1008 float bl_phase =
atan2f(fl_im + sl_im, fl_re + sl_re);
1009 float br_phase =
atan2f(fr_im + sr_im, fr_re + sr_re);
1013 if (phase_difl >
M_PIf)
1014 phase_difl = 2.f *
M_PIf - phase_difl;
1016 if (phase_difr >
M_PIf)
1017 phase_difr = 2.f *
M_PIf - phase_difr;
1022 s->upmix_5_1(
ctx, c_re, c_im, lfe_re, lfe_im,
1023 mag_totall, mag_totalr,
1034 const int rdft_size =
s->rdft_size;
1035 float *srcl, *srcr, *srcc, *srclfe, *srcbl, *srcbr;
1038 srcl = (
float *)
s->input->extended_data[0];
1039 srcr = (
float *)
s->input->extended_data[1];
1040 srcc = (
float *)
s->input->extended_data[2];
1041 srclfe = (
float *)
s->input->extended_data[3];
1042 srcbl = (
float *)
s->input->extended_data[4];
1043 srcbr = (
float *)
s->input->extended_data[5];
1045 for (n = 0; n < rdft_size; n++) {
1046 float fl_re = srcl[2 * n], fr_re = srcr[2 * n];
1047 float fl_im = srcl[2 * n + 1], fr_im = srcr[2 * n + 1];
1048 float c_re = srcc[2 * n], c_im = srcc[2 * n + 1];
1049 float lfe_re = srclfe[2 * n], lfe_im = srclfe[2 * n + 1];
1050 float bl_re = srcbl[2 * n], bl_im = srcbl[2 * n + 1];
1051 float br_re = srcbr[2 * n], br_im = srcbr[2 * n + 1];
1052 float fl_mag = hypotf(fl_re, fl_im);
1053 float fr_mag = hypotf(fr_re, fr_im);
1054 float fl_phase =
atan2f(fl_im, fl_re);
1055 float fr_phase =
atan2f(fr_im, fr_re);
1056 float bl_mag = hypotf(bl_re, bl_im);
1057 float br_mag = hypotf(br_re, br_im);
1058 float bl_phase =
atan2f(bl_im, bl_re);
1059 float br_phase =
atan2f(br_im, br_re);
1060 float phase_difl =
fabsf(fl_phase - bl_phase);
1061 float phase_difr =
fabsf(fr_phase - br_phase);
1062 float magl_sum = fl_mag + bl_mag;
1063 float magr_sum = fr_mag + br_mag;
1064 float mag_difl = magl_sum <
MIN_MAG_SUM ?
FFDIFFSIGN(fl_mag, bl_mag) : (fl_mag - bl_mag) / magl_sum;
1065 float mag_difr = magr_sum <
MIN_MAG_SUM ?
FFDIFFSIGN(fr_mag, br_mag) : (fr_mag - br_mag) / magr_sum;
1066 float mag_totall = hypotf(fl_mag, bl_mag);
1067 float mag_totalr = hypotf(fr_mag, br_mag);
1068 float sl_phase =
atan2f(fl_im + bl_im, fl_re + bl_re);
1069 float sr_phase =
atan2f(fr_im + br_im, fr_re + br_re);
1073 if (phase_difl >
M_PIf)
1074 phase_difl = 2.f *
M_PIf - phase_difl;
1076 if (phase_difr >
M_PIf)
1077 phase_difr = 2.f *
M_PIf - phase_difr;
1082 s->upmix_5_1(
ctx, c_re, c_im, lfe_re, lfe_im,
1083 mag_totall, mag_totalr,
1095 if (
s->all_x >= 0.f)
1096 for (
int n = 0; n <
SC_NB; n++)
1097 s->f_x[n] =
s->all_x;
1099 if (
s->all_y >= 0.f)
1100 for (
int n = 0; n <
SC_NB; n++)
1101 s->f_y[n] =
s->all_y;
1108 int64_t in_channel_layout, out_channel_layout;
1109 char in_name[128], out_name[128];
1112 if (
s->lowcutf >=
s->highcutf) {
1114 s->lowcutf,
s->highcutf);
1119 s->in_ch_layout.u.mask : 0;
1121 s->out_ch_layout.u.mask : 0;
1126 switch (in_channel_layout) {
1141 switch (out_channel_layout) {
1151 switch (out_channel_layout) {
1161 switch (out_channel_layout) {
1178 s->window_func_lut =
av_calloc(
s->win_size,
sizeof(*
s->window_func_lut));
1179 if (!
s->window_func_lut)
1183 if (
s->overlap == 1)
1184 s->overlap = overlap;
1186 for (
int i = 0;
i <
s->win_size;
i++)
1187 s->window_func_lut[
i] =
sqrtf(
s->window_func_lut[
i] /
s->win_size);
1188 s->hop_size =
FFMAX(1,
s->win_size * (1. -
s->overlap));
1191 float max = 0.f, *temp_lut =
av_calloc(
s->win_size,
sizeof(*temp_lut));
1195 for (
int j = 0; j <
s->win_size; j +=
s->hop_size) {
1196 for (
int i = 0;
i <
s->win_size;
i++)
1197 temp_lut[(
i + j) %
s->win_size] +=
s->window_func_lut[
i];
1200 for (
int i = 0;
i <
s->win_size;
i++)
1204 s->win_gain = 1.f / (
max *
sqrtf(
s->win_size));
1215 float *
src = (
float *)
s->input_in->extended_data[ch];
1216 float *
win = (
float *)
s->window->extended_data[ch];
1217 const float *window_func_lut =
s->window_func_lut;
1218 const int offset =
s->win_size -
s->hop_size;
1219 const float level_in =
s->input_levels[ch];
1220 const int win_size =
s->win_size;
1226 for (
int n = 0; n < win_size; n++)
1227 win[n] =
src[n] * window_func_lut[n] * level_in;
1229 s->tx_fn(
s->rdft[ch], (
float *)
s->input->extended_data[ch],
win,
sizeof(
float));
1240 for (
int ch = start; ch < end; ch++)
1249 const float level_out =
s->output_levels[ch] *
s->win_gain;
1250 const float *window_func_lut =
s->window_func_lut;
1251 const int win_size =
s->win_size;
1254 dst = (
float *)
s->output_out->extended_data[ch];
1255 ptr = (
float *)
s->overlap_buffer->extended_data[ch];
1256 s->itx_fn(
s->irdft[ch], dst, (
float *)
s->output->extended_data[ch],
sizeof(
AVComplexFloat));
1258 memmove(
s->overlap_buffer->extended_data[ch],
1259 s->overlap_buffer->extended_data[ch] +
s->hop_size *
sizeof(
float),
1260 s->win_size *
sizeof(
float));
1261 memset(
s->overlap_buffer->extended_data[ch] +
s->win_size *
sizeof(
float),
1262 0,
s->hop_size *
sizeof(
float));
1264 for (
int n = 0; n < win_size; n++)
1265 ptr[n] += dst[n] * window_func_lut[n] * level_out;
1267 ptr = (
float *)
s->overlap_buffer->extended_data[ch];
1268 dst = (
float *)
out->extended_data[ch];
1269 memcpy(dst, ptr,
s->hop_size *
sizeof(
float));
1278 const int start = (
out->ch_layout.nb_channels * jobnr) / nb_jobs;
1279 const int end = (
out->ch_layout.nb_channels * (jobnr+1)) / nb_jobs;
1281 for (
int ch = start; ch < end; ch++) {
1368 for (
int ch = 0; ch <
s->nb_in_channels; ch++)
1370 for (
int ch = 0; ch <
s->nb_out_channels; ch++)
1390 char *res,
int res_len,
int flags)
1399 s->hop_size =
FFMAX(1,
s->win_size * (1. -
s->overlap));
1408 #define OFFSET(x) offsetof(AudioSurroundContext, x)
1409 #define FLAGS AV_OPT_FLAG_AUDIO_PARAM|AV_OPT_FLAG_FILTERING_PARAM
1410 #define TFLAGS AV_OPT_FLAG_AUDIO_PARAM|AV_OPT_FLAG_FILTERING_PARAM|AV_OPT_FLAG_RUNTIME_PARAM
1422 {
"sub",
"subtract LFE channel with others", 0,
AV_OPT_TYPE_CONST, {.i64=1}, 0, 1,
TFLAGS, .unit =
"lfe_mode" },
1490 .priv_class = &surround_class,