Lines Matching refs:a
376 static void cftfsub(picoos_int32 n, PICOFFTSG_FFTTYPE *a);
377 static void cftbsub(picoos_int32 n, PICOFFTSG_FFTTYPE *a);
378 static void rftfsub(picoos_int32 n, PICOFFTSG_FFTTYPE *a);
379 static void rftbsub(picoos_int32 n, PICOFFTSG_FFTTYPE *a);
381 static void cftfsub(picoos_int32 n, PICOFFTSG_FFTTYPE *a);
382 static void cftbsub(picoos_int32 n, PICOFFTSG_FFTTYPE *a);
383 static void rftfsub(picoos_int32 n, PICOFFTSG_FFTTYPE *a);
384 static void rftbsub(picoos_int32 n, PICOFFTSG_FFTTYPE *a);
385 static void dctsub(picoos_int32 n, PICOFFTSG_FFTTYPE *a);
386 static void dctsub4(picoos_int32 n, PICOFFTSG_FFTTYPE *a);
388 static void ddct(picoos_int32 n, picoos_int32 isgn, PICOFFTSG_FFTTYPE *a);
389 static void bitrv1(picoos_int32 n, PICOFFTSG_FFTTYPE *a);
391 static void bitrv2(picoos_int32 n, PICOFFTSG_FFTTYPE *a);
392 static void bitrv216(PICOFFTSG_FFTTYPE *a);
393 static void bitrv208(PICOFFTSG_FFTTYPE *a);
394 static void cftmdl1(picoos_int32 n, PICOFFTSG_FFTTYPE *a);
395 static void cftrec4(picoos_int32 n, PICOFFTSG_FFTTYPE *a);
396 static void cftleaf(picoos_int32 n, picoos_int32 isplt, PICOFFTSG_FFTTYPE *a);
397 static void cftfx41(picoos_int32 n, PICOFFTSG_FFTTYPE *a);
398 static void cftf161(PICOFFTSG_FFTTYPE *a);
399 static void cftf081(PICOFFTSG_FFTTYPE *a);
400 static void cftf040(PICOFFTSG_FFTTYPE *a);
401 static void cftx020(PICOFFTSG_FFTTYPE *a);
403 void bitrv2conj(picoos_int32 n, PICOFFTSG_FFTTYPE *a);
404 void bitrv216neg(PICOFFTSG_FFTTYPE *a);
405 void bitrv208neg(PICOFFTSG_FFTTYPE *a);
406 void cftb1st(picoos_int32 n, PICOFFTSG_FFTTYPE *a);
407 void cftrec4(picoos_int32 n, PICOFFTSG_FFTTYPE *a);
408 void cftleaf(picoos_int32 n, picoos_int32 isplt, PICOFFTSG_FFTTYPE *a);
409 void cftfx41(picoos_int32 n, PICOFFTSG_FFTTYPE *a);
410 void cftf161(PICOFFTSG_FFTTYPE *a);
411 void cftf081(PICOFFTSG_FFTTYPE *a);
412 void cftb040(PICOFFTSG_FFTTYPE *a);
413 void cftx020(PICOFFTSG_FFTTYPE *a);
415 static picoos_int32 cfttree(picoos_int32 n, picoos_int32 j, picoos_int32 k, PICOFFTSG_FFTTYPE *a);
416 static void cftleaf(picoos_int32 n, picoos_int32 isplt, PICOFFTSG_FFTTYPE *a);
417 static void cftmdl1(picoos_int32 n, PICOFFTSG_FFTTYPE *a);
419 static void cftmdl1(picoos_int32 n, PICOFFTSG_FFTTYPE *a);
420 static void cftmdl2(picoos_int32 n, PICOFFTSG_FFTTYPE *a);
422 static void cftmdl1(picoos_int32 n, PICOFFTSG_FFTTYPE *a);
423 static void cftmdl2(picoos_int32 n, PICOFFTSG_FFTTYPE *a);
424 static void cftf161(PICOFFTSG_FFTTYPE *a);
425 static void cftf162(PICOFFTSG_FFTTYPE *a);
426 static void cftf081(PICOFFTSG_FFTTYPE *a);
427 static void cftf082(PICOFFTSG_FFTTYPE *a);
429 static void cftf161(PICOFFTSG_FFTTYPE *a);
430 static void cftf162(PICOFFTSG_FFTTYPE *a);
431 static void cftf081(PICOFFTSG_FFTTYPE *a);
432 static void cftf082(PICOFFTSG_FFTTYPE *a);
437 void rdft(picoos_int32 n, picoos_int32 isgn, PICOFFTSG_FFTTYPE *a) in rdft() argument
443 cftfsub(n, a); in rdft()
444 rftfsub(n, a); in rdft()
446 cftfsub(n, a); in rdft()
448 xi = a[0] - a[1]; in rdft()
449 a[0] += a[1]; in rdft()
450 a[1] = xi; in rdft()
452 a[1] = (a[0] - a[1]) / 2; in rdft()
453 a[0] -= a[1]; in rdft()
455 rftbsub(n, a); in rdft()
456 cftbsub(n, a); in rdft()
458 cftbsub(n, a); in rdft()
468 PICOFFTSG_FFTTYPE a,b, E; in norm_result() local
472 a = (norm_window[nI]>>18) * ((tmpX[nI]>0) ? tmpX[nI]>>11 : -((-tmpX[nI])>>11)); in norm_result()
473 tmpX[nI] = a; in norm_result()
474 b = (a>=0?a:-a) >> 18; in norm_result()
486 void ddct(picoos_int32 n, picoos_int32 isgn, PICOFFTSG_FFTTYPE *a) in ddct() argument
492 xr = a[n - 1]; in ddct()
494 a[j + 1] = a[j] - a[j - 1]; in ddct()
495 a[j] += a[j - 1]; in ddct()
497 a[1] = a[0] - xr; in ddct()
498 a[0] += xr; in ddct()
500 rftbsub(n, a); in ddct()
501 cftbsub(n, a); in ddct()
503 cftbsub(n, a); in ddct()
507 dctsub(n, a); in ddct()
509 dctsub4(n, a); in ddct()
513 cftfsub(n, a); in ddct()
514 rftfsub(n, a); in ddct()
516 cftfsub(n, a); in ddct()
520 xr = a[0] - a[1]; in ddct()
521 a[0] += a[1]; in ddct()
523 a[j - 1] = a[j] - a[j + 1]; in ddct()
524 a[j] += a[j + 1]; in ddct()
526 a[n - 1] = xr; in ddct()
530 void dfct_nmf(picoos_int32 n, picoos_int32 *a) in dfct_nmf() argument
539 xr = a[j] + a[k]; in dfct_nmf()
540 a[j] -= a[k]; in dfct_nmf()
541 a[k] = xr; in dfct_nmf()
543 an = a[n]; in dfct_nmf()
545 ddct(m, 1, a); in dfct_nmf()
547 bitrv1(m, a); in dfct_nmf()
550 xi = a[m]; in dfct_nmf()
551 a[m] = a[0]; in dfct_nmf()
552 a[0] = an - xi; in dfct_nmf()
555 aj = a + 1; ak = a + k; amj = aj + m; amk = ak + m; in dfct_nmf()
573 xi = a[1]; in dfct_nmf()
574 a[1] = a[0]; in dfct_nmf()
575 a[0] = an + xi; in dfct_nmf()
576 a[n] = an - xi; in dfct_nmf()
578 bitrv1(n, a); in dfct_nmf()
609 static void cftfsub(picoos_int32 n, PICOFFTSG_FFTTYPE *a) in cftfsub() argument
614 cftmdl1(n, a); in cftfsub()
616 cftrec4(n, a); in cftfsub()
618 cftleaf(n, 1, a); in cftfsub()
620 cftfx41(n, a); in cftfsub()
622 bitrv2(n, a); in cftfsub()
624 cftf161(a); in cftfsub()
625 bitrv216(a); in cftfsub()
627 cftf081(a); in cftfsub()
628 bitrv208(a); in cftfsub()
631 cftf040(a); in cftfsub()
633 cftx020(a); in cftfsub()
638 void cftbsub(picoos_int32 n, PICOFFTSG_FFTTYPE *a) in cftbsub() argument
642 cftb1st(n, a); in cftbsub()
644 cftrec4(n, a); in cftbsub()
646 cftleaf(n, 1, a); in cftbsub()
648 cftfx41(n, a); in cftbsub()
650 bitrv2conj(n, a); in cftbsub()
652 cftf161(a); in cftbsub()
653 bitrv216neg(a); in cftbsub()
655 cftf081(a); in cftbsub()
656 bitrv208neg(a); in cftbsub()
659 cftb040(a); in cftbsub()
661 cftx020(a); in cftbsub()
668 void bitrv2(picoos_int32 n, PICOFFTSG_FFTTYPE *a) in bitrv2() argument
684 xr = a[j]; in bitrv2()
685 xi = a[j + 1]; in bitrv2()
686 yr = a[k]; in bitrv2()
687 yi = a[k + 1]; in bitrv2()
688 a[j] = yr; in bitrv2()
689 a[j + 1] = yi; in bitrv2()
690 a[k] = xr; in bitrv2()
691 a[k + 1] = xi; in bitrv2()
694 xr = a[j1]; in bitrv2()
695 xi = a[j1 + 1]; in bitrv2()
696 yr = a[k1]; in bitrv2()
697 yi = a[k1 + 1]; in bitrv2()
698 a[j1] = yr; in bitrv2()
699 a[j1 + 1] = yi; in bitrv2()
700 a[k1] = xr; in bitrv2()
701 a[k1 + 1] = xi; in bitrv2()
704 xr = a[j1]; in bitrv2()
705 xi = a[j1 + 1]; in bitrv2()
706 yr = a[k1]; in bitrv2()
707 yi = a[k1 + 1]; in bitrv2()
708 a[j1] = yr; in bitrv2()
709 a[j1 + 1] = yi; in bitrv2()
710 a[k1] = xr; in bitrv2()
711 a[k1 + 1] = xi; in bitrv2()
714 xr = a[j1]; in bitrv2()
715 xi = a[j1 + 1]; in bitrv2()
716 yr = a[k1]; in bitrv2()
717 yi = a[k1 + 1]; in bitrv2()
718 a[j1] = yr; in bitrv2()
719 a[j1 + 1] = yi; in bitrv2()
720 a[k1] = xr; in bitrv2()
721 a[k1 + 1] = xi; in bitrv2()
724 xr = a[j1]; in bitrv2()
725 xi = a[j1 + 1]; in bitrv2()
726 yr = a[k1]; in bitrv2()
727 yi = a[k1 + 1]; in bitrv2()
728 a[j1] = yr; in bitrv2()
729 a[j1 + 1] = yi; in bitrv2()
730 a[k1] = xr; in bitrv2()
731 a[k1 + 1] = xi; in bitrv2()
734 xr = a[j1]; in bitrv2()
735 xi = a[j1 + 1]; in bitrv2()
736 yr = a[k1]; in bitrv2()
737 yi = a[k1 + 1]; in bitrv2()
738 a[j1] = yr; in bitrv2()
739 a[j1 + 1] = yi; in bitrv2()
740 a[k1] = xr; in bitrv2()
741 a[k1 + 1] = xi; in bitrv2()
744 xr = a[j1]; in bitrv2()
745 xi = a[j1 + 1]; in bitrv2()
746 yr = a[k1]; in bitrv2()
747 yi = a[k1 + 1]; in bitrv2()
748 a[j1] = yr; in bitrv2()
749 a[j1 + 1] = yi; in bitrv2()
750 a[k1] = xr; in bitrv2()
751 a[k1 + 1] = xi; in bitrv2()
754 xr = a[j1]; in bitrv2()
755 xi = a[j1 + 1]; in bitrv2()
756 yr = a[k1]; in bitrv2()
757 yi = a[k1 + 1]; in bitrv2()
758 a[j1] = yr; in bitrv2()
759 a[j1 + 1] = yi; in bitrv2()
760 a[k1] = xr; in bitrv2()
761 a[k1 + 1] = xi; in bitrv2()
764 xr = a[j1]; in bitrv2()
765 xi = a[j1 + 1]; in bitrv2()
766 yr = a[k1]; in bitrv2()
767 yi = a[k1 + 1]; in bitrv2()
768 a[j1] = yr; in bitrv2()
769 a[j1 + 1] = yi; in bitrv2()
770 a[k1] = xr; in bitrv2()
771 a[k1 + 1] = xi; in bitrv2()
774 xr = a[j1]; in bitrv2()
775 xi = a[j1 + 1]; in bitrv2()
776 yr = a[k1]; in bitrv2()
777 yi = a[k1 + 1]; in bitrv2()
778 a[j1] = yr; in bitrv2()
779 a[j1 + 1] = yi; in bitrv2()
780 a[k1] = xr; in bitrv2()
781 a[k1 + 1] = xi; in bitrv2()
784 xr = a[j1]; in bitrv2()
785 xi = a[j1 + 1]; in bitrv2()
786 yr = a[k1]; in bitrv2()
787 yi = a[k1 + 1]; in bitrv2()
788 a[j1] = yr; in bitrv2()
789 a[j1 + 1] = yi; in bitrv2()
790 a[k1] = xr; in bitrv2()
791 a[k1 + 1] = xi; in bitrv2()
794 xr = a[j1]; in bitrv2()
795 xi = a[j1 + 1]; in bitrv2()
796 yr = a[k1]; in bitrv2()
797 yi = a[k1 + 1]; in bitrv2()
798 a[j1] = yr; in bitrv2()
799 a[j1 + 1] = yi; in bitrv2()
800 a[k1] = xr; in bitrv2()
801 a[k1 + 1] = xi; in bitrv2()
804 xr = a[j1]; in bitrv2()
805 xi = a[j1 + 1]; in bitrv2()
806 yr = a[k1]; in bitrv2()
807 yi = a[k1 + 1]; in bitrv2()
808 a[j1] = yr; in bitrv2()
809 a[j1 + 1] = yi; in bitrv2()
810 a[k1] = xr; in bitrv2()
811 a[k1 + 1] = xi; in bitrv2()
814 xr = a[j1]; in bitrv2()
815 xi = a[j1 + 1]; in bitrv2()
816 yr = a[k1]; in bitrv2()
817 yi = a[k1 + 1]; in bitrv2()
818 a[j1] = yr; in bitrv2()
819 a[j1 + 1] = yi; in bitrv2()
820 a[k1] = xr; in bitrv2()
821 a[k1 + 1] = xi; in bitrv2()
824 xr = a[j1]; in bitrv2()
825 xi = a[j1 + 1]; in bitrv2()
826 yr = a[k1]; in bitrv2()
827 yi = a[k1 + 1]; in bitrv2()
828 a[j1] = yr; in bitrv2()
829 a[j1 + 1] = yi; in bitrv2()
830 a[k1] = xr; in bitrv2()
831 a[k1 + 1] = xi; in bitrv2()
834 xr = a[j1]; in bitrv2()
835 xi = a[j1 + 1]; in bitrv2()
836 yr = a[k1]; in bitrv2()
837 yi = a[k1 + 1]; in bitrv2()
838 a[j1] = yr; in bitrv2()
839 a[j1 + 1] = yi; in bitrv2()
840 a[k1] = xr; in bitrv2()
841 a[k1 + 1] = xi; in bitrv2()
849 xr = a[j1]; in bitrv2()
850 xi = a[j1 + 1]; in bitrv2()
851 yr = a[k1]; in bitrv2()
852 yi = a[k1 + 1]; in bitrv2()
853 a[j1] = yr; in bitrv2()
854 a[j1 + 1] = yi; in bitrv2()
855 a[k1] = xr; in bitrv2()
856 a[k1 + 1] = xi; in bitrv2()
859 xr = a[j1]; in bitrv2()
860 xi = a[j1 + 1]; in bitrv2()
861 yr = a[k1]; in bitrv2()
862 yi = a[k1 + 1]; in bitrv2()
863 a[j1] = yr; in bitrv2()
864 a[j1 + 1] = yi; in bitrv2()
865 a[k1] = xr; in bitrv2()
866 a[k1 + 1] = xi; in bitrv2()
869 xr = a[j1]; in bitrv2()
870 xi = a[j1 + 1]; in bitrv2()
871 yr = a[k1]; in bitrv2()
872 yi = a[k1 + 1]; in bitrv2()
873 a[j1] = yr; in bitrv2()
874 a[j1 + 1] = yi; in bitrv2()
875 a[k1] = xr; in bitrv2()
876 a[k1 + 1] = xi; in bitrv2()
879 xr = a[j1]; in bitrv2()
880 xi = a[j1 + 1]; in bitrv2()
881 yr = a[k1]; in bitrv2()
882 yi = a[k1 + 1]; in bitrv2()
883 a[j1] = yr; in bitrv2()
884 a[j1 + 1] = yi; in bitrv2()
885 a[k1] = xr; in bitrv2()
886 a[k1 + 1] = xi; in bitrv2()
889 xr = a[j1]; in bitrv2()
890 xi = a[j1 + 1]; in bitrv2()
891 yr = a[k1]; in bitrv2()
892 yi = a[k1 + 1]; in bitrv2()
893 a[j1] = yr; in bitrv2()
894 a[j1 + 1] = yi; in bitrv2()
895 a[k1] = xr; in bitrv2()
896 a[k1 + 1] = xi; in bitrv2()
899 xr = a[j1]; in bitrv2()
900 xi = a[j1 + 1]; in bitrv2()
901 yr = a[k1]; in bitrv2()
902 yi = a[k1 + 1]; in bitrv2()
903 a[j1] = yr; in bitrv2()
904 a[j1 + 1] = yi; in bitrv2()
905 a[k1] = xr; in bitrv2()
906 a[k1 + 1] = xi; in bitrv2()
916 xr = a[j]; in bitrv2()
917 xi = a[j + 1]; in bitrv2()
918 yr = a[k]; in bitrv2()
919 yi = a[k + 1]; in bitrv2()
920 a[j] = yr; in bitrv2()
921 a[j + 1] = yi; in bitrv2()
922 a[k] = xr; in bitrv2()
923 a[k + 1] = xi; in bitrv2()
926 xr = a[j1]; in bitrv2()
927 xi = a[j1 + 1]; in bitrv2()
928 yr = a[k1]; in bitrv2()
929 yi = a[k1 + 1]; in bitrv2()
930 a[j1] = yr; in bitrv2()
931 a[j1 + 1] = yi; in bitrv2()
932 a[k1] = xr; in bitrv2()
933 a[k1 + 1] = xi; in bitrv2()
936 xr = a[j1]; in bitrv2()
937 xi = a[j1 + 1]; in bitrv2()
938 yr = a[k1]; in bitrv2()
939 yi = a[k1 + 1]; in bitrv2()
940 a[j1] = yr; in bitrv2()
941 a[j1 + 1] = yi; in bitrv2()
942 a[k1] = xr; in bitrv2()
943 a[k1 + 1] = xi; in bitrv2()
946 xr = a[j1]; in bitrv2()
947 xi = a[j1 + 1]; in bitrv2()
948 yr = a[k1]; in bitrv2()
949 yi = a[k1 + 1]; in bitrv2()
950 a[j1] = yr; in bitrv2()
951 a[j1 + 1] = yi; in bitrv2()
952 a[k1] = xr; in bitrv2()
953 a[k1 + 1] = xi; in bitrv2()
956 xr = a[j1]; in bitrv2()
957 xi = a[j1 + 1]; in bitrv2()
958 yr = a[k1]; in bitrv2()
959 yi = a[k1 + 1]; in bitrv2()
960 a[j1] = yr; in bitrv2()
961 a[j1 + 1] = yi; in bitrv2()
962 a[k1] = xr; in bitrv2()
963 a[k1 + 1] = xi; in bitrv2()
966 xr = a[j1]; in bitrv2()
967 xi = a[j1 + 1]; in bitrv2()
968 yr = a[k1]; in bitrv2()
969 yi = a[k1 + 1]; in bitrv2()
970 a[j1] = yr; in bitrv2()
971 a[j1 + 1] = yi; in bitrv2()
972 a[k1] = xr; in bitrv2()
973 a[k1 + 1] = xi; in bitrv2()
976 xr = a[j1]; in bitrv2()
977 xi = a[j1 + 1]; in bitrv2()
978 yr = a[k1]; in bitrv2()
979 yi = a[k1 + 1]; in bitrv2()
980 a[j1] = yr; in bitrv2()
981 a[j1 + 1] = yi; in bitrv2()
982 a[k1] = xr; in bitrv2()
983 a[k1 + 1] = xi; in bitrv2()
986 xr = a[j1]; in bitrv2()
987 xi = a[j1 + 1]; in bitrv2()
988 yr = a[k1]; in bitrv2()
989 yi = a[k1 + 1]; in bitrv2()
990 a[j1] = yr; in bitrv2()
991 a[j1 + 1] = yi; in bitrv2()
992 a[k1] = xr; in bitrv2()
993 a[k1 + 1] = xi; in bitrv2()
1001 xr = a[j1]; in bitrv2()
1002 xi = a[j1 + 1]; in bitrv2()
1003 yr = a[k1]; in bitrv2()
1004 yi = a[k1 + 1]; in bitrv2()
1005 a[j1] = yr; in bitrv2()
1006 a[j1 + 1] = yi; in bitrv2()
1007 a[k1] = xr; in bitrv2()
1008 a[k1 + 1] = xi; in bitrv2()
1011 xr = a[j1]; in bitrv2()
1012 xi = a[j1 + 1]; in bitrv2()
1013 yr = a[k1]; in bitrv2()
1014 yi = a[k1 + 1]; in bitrv2()
1015 a[j1] = yr; in bitrv2()
1016 a[j1 + 1] = yi; in bitrv2()
1017 a[k1] = xr; in bitrv2()
1018 a[k1 + 1] = xi; in bitrv2()
1027 void bitrv2conj(picoos_int32 n, PICOFFTSG_FFTTYPE *a) in bitrv2conj() argument
1044 xr = a[j]; in bitrv2conj()
1045 xi = -a[j + 1]; in bitrv2conj()
1046 yr = a[k]; in bitrv2conj()
1047 yi = -a[k + 1]; in bitrv2conj()
1048 a[j] = yr; in bitrv2conj()
1049 a[j + 1] = yi; in bitrv2conj()
1050 a[k] = xr; in bitrv2conj()
1051 a[k + 1] = xi; in bitrv2conj()
1054 xr = a[j1]; in bitrv2conj()
1055 xi = -a[j1 + 1]; in bitrv2conj()
1056 yr = a[k1]; in bitrv2conj()
1057 yi = -a[k1 + 1]; in bitrv2conj()
1058 a[j1] = yr; in bitrv2conj()
1059 a[j1 + 1] = yi; in bitrv2conj()
1060 a[k1] = xr; in bitrv2conj()
1061 a[k1 + 1] = xi; in bitrv2conj()
1064 xr = a[j1]; in bitrv2conj()
1065 xi = -a[j1 + 1]; in bitrv2conj()
1066 yr = a[k1]; in bitrv2conj()
1067 yi = -a[k1 + 1]; in bitrv2conj()
1068 a[j1] = yr; in bitrv2conj()
1069 a[j1 + 1] = yi; in bitrv2conj()
1070 a[k1] = xr; in bitrv2conj()
1071 a[k1 + 1] = xi; in bitrv2conj()
1074 xr = a[j1]; in bitrv2conj()
1075 xi = -a[j1 + 1]; in bitrv2conj()
1076 yr = a[k1]; in bitrv2conj()
1077 yi = -a[k1 + 1]; in bitrv2conj()
1078 a[j1] = yr; in bitrv2conj()
1079 a[j1 + 1] = yi; in bitrv2conj()
1080 a[k1] = xr; in bitrv2conj()
1081 a[k1 + 1] = xi; in bitrv2conj()
1084 xr = a[j1]; in bitrv2conj()
1085 xi = -a[j1 + 1]; in bitrv2conj()
1086 yr = a[k1]; in bitrv2conj()
1087 yi = -a[k1 + 1]; in bitrv2conj()
1088 a[j1] = yr; in bitrv2conj()
1089 a[j1 + 1] = yi; in bitrv2conj()
1090 a[k1] = xr; in bitrv2conj()
1091 a[k1 + 1] = xi; in bitrv2conj()
1094 xr = a[j1]; in bitrv2conj()
1095 xi = -a[j1 + 1]; in bitrv2conj()
1096 yr = a[k1]; in bitrv2conj()
1097 yi = -a[k1 + 1]; in bitrv2conj()
1098 a[j1] = yr; in bitrv2conj()
1099 a[j1 + 1] = yi; in bitrv2conj()
1100 a[k1] = xr; in bitrv2conj()
1101 a[k1 + 1] = xi; in bitrv2conj()
1104 xr = a[j1]; in bitrv2conj()
1105 xi = -a[j1 + 1]; in bitrv2conj()
1106 yr = a[k1]; in bitrv2conj()
1107 yi = -a[k1 + 1]; in bitrv2conj()
1108 a[j1] = yr; in bitrv2conj()
1109 a[j1 + 1] = yi; in bitrv2conj()
1110 a[k1] = xr; in bitrv2conj()
1111 a[k1 + 1] = xi; in bitrv2conj()
1114 xr = a[j1]; in bitrv2conj()
1115 xi = -a[j1 + 1]; in bitrv2conj()
1116 yr = a[k1]; in bitrv2conj()
1117 yi = -a[k1 + 1]; in bitrv2conj()
1118 a[j1] = yr; in bitrv2conj()
1119 a[j1 + 1] = yi; in bitrv2conj()
1120 a[k1] = xr; in bitrv2conj()
1121 a[k1 + 1] = xi; in bitrv2conj()
1124 xr = a[j1]; in bitrv2conj()
1125 xi = -a[j1 + 1]; in bitrv2conj()
1126 yr = a[k1]; in bitrv2conj()
1127 yi = -a[k1 + 1]; in bitrv2conj()
1128 a[j1] = yr; in bitrv2conj()
1129 a[j1 + 1] = yi; in bitrv2conj()
1130 a[k1] = xr; in bitrv2conj()
1131 a[k1 + 1] = xi; in bitrv2conj()
1134 xr = a[j1]; in bitrv2conj()
1135 xi = -a[j1 + 1]; in bitrv2conj()
1136 yr = a[k1]; in bitrv2conj()
1137 yi = -a[k1 + 1]; in bitrv2conj()
1138 a[j1] = yr; in bitrv2conj()
1139 a[j1 + 1] = yi; in bitrv2conj()
1140 a[k1] = xr; in bitrv2conj()
1141 a[k1 + 1] = xi; in bitrv2conj()
1144 xr = a[j1]; in bitrv2conj()
1145 xi = -a[j1 + 1]; in bitrv2conj()
1146 yr = a[k1]; in bitrv2conj()
1147 yi = -a[k1 + 1]; in bitrv2conj()
1148 a[j1] = yr; in bitrv2conj()
1149 a[j1 + 1] = yi; in bitrv2conj()
1150 a[k1] = xr; in bitrv2conj()
1151 a[k1 + 1] = xi; in bitrv2conj()
1154 xr = a[j1]; in bitrv2conj()
1155 xi = -a[j1 + 1]; in bitrv2conj()
1156 yr = a[k1]; in bitrv2conj()
1157 yi = -a[k1 + 1]; in bitrv2conj()
1158 a[j1] = yr; in bitrv2conj()
1159 a[j1 + 1] = yi; in bitrv2conj()
1160 a[k1] = xr; in bitrv2conj()
1161 a[k1 + 1] = xi; in bitrv2conj()
1164 xr = a[j1]; in bitrv2conj()
1165 xi = -a[j1 + 1]; in bitrv2conj()
1166 yr = a[k1]; in bitrv2conj()
1167 yi = -a[k1 + 1]; in bitrv2conj()
1168 a[j1] = yr; in bitrv2conj()
1169 a[j1 + 1] = yi; in bitrv2conj()
1170 a[k1] = xr; in bitrv2conj()
1171 a[k1 + 1] = xi; in bitrv2conj()
1174 xr = a[j1]; in bitrv2conj()
1175 xi = -a[j1 + 1]; in bitrv2conj()
1176 yr = a[k1]; in bitrv2conj()
1177 yi = -a[k1 + 1]; in bitrv2conj()
1178 a[j1] = yr; in bitrv2conj()
1179 a[j1 + 1] = yi; in bitrv2conj()
1180 a[k1] = xr; in bitrv2conj()
1181 a[k1 + 1] = xi; in bitrv2conj()
1184 xr = a[j1]; in bitrv2conj()
1185 xi = -a[j1 + 1]; in bitrv2conj()
1186 yr = a[k1]; in bitrv2conj()
1187 yi = -a[k1 + 1]; in bitrv2conj()
1188 a[j1] = yr; in bitrv2conj()
1189 a[j1 + 1] = yi; in bitrv2conj()
1190 a[k1] = xr; in bitrv2conj()
1191 a[k1 + 1] = xi; in bitrv2conj()
1194 xr = a[j1]; in bitrv2conj()
1195 xi = -a[j1 + 1]; in bitrv2conj()
1196 yr = a[k1]; in bitrv2conj()
1197 yi = -a[k1 + 1]; in bitrv2conj()
1198 a[j1] = yr; in bitrv2conj()
1199 a[j1 + 1] = yi; in bitrv2conj()
1200 a[k1] = xr; in bitrv2conj()
1201 a[k1 + 1] = xi; in bitrv2conj()
1209 a[j1 - 1] = -a[j1 - 1]; in bitrv2conj()
1210 xr = a[j1]; in bitrv2conj()
1211 xi = -a[j1 + 1]; in bitrv2conj()
1212 yr = a[k1]; in bitrv2conj()
1213 yi = -a[k1 + 1]; in bitrv2conj()
1214 a[j1] = yr; in bitrv2conj()
1215 a[j1 + 1] = yi; in bitrv2conj()
1216 a[k1] = xr; in bitrv2conj()
1217 a[k1 + 1] = xi; in bitrv2conj()
1218 a[k1 + 3] = -a[k1 + 3]; in bitrv2conj()
1221 xr = a[j1]; in bitrv2conj()
1222 xi = -a[j1 + 1]; in bitrv2conj()
1223 yr = a[k1]; in bitrv2conj()
1224 yi = -a[k1 + 1]; in bitrv2conj()
1225 a[j1] = yr; in bitrv2conj()
1226 a[j1 + 1] = yi; in bitrv2conj()
1227 a[k1] = xr; in bitrv2conj()
1228 a[k1 + 1] = xi; in bitrv2conj()
1231 xr = a[j1]; in bitrv2conj()
1232 xi = -a[j1 + 1]; in bitrv2conj()
1233 yr = a[k1]; in bitrv2conj()
1234 yi = -a[k1 + 1]; in bitrv2conj()
1235 a[j1] = yr; in bitrv2conj()
1236 a[j1 + 1] = yi; in bitrv2conj()
1237 a[k1] = xr; in bitrv2conj()
1238 a[k1 + 1] = xi; in bitrv2conj()
1241 xr = a[j1]; in bitrv2conj()
1242 xi = -a[j1 + 1]; in bitrv2conj()
1243 yr = a[k1]; in bitrv2conj()
1244 yi = -a[k1 + 1]; in bitrv2conj()
1245 a[j1] = yr; in bitrv2conj()
1246 a[j1 + 1] = yi; in bitrv2conj()
1247 a[k1] = xr; in bitrv2conj()
1248 a[k1 + 1] = xi; in bitrv2conj()
1251 xr = a[j1]; in bitrv2conj()
1252 xi = -a[j1 + 1]; in bitrv2conj()
1253 yr = a[k1]; in bitrv2conj()
1254 yi = -a[k1 + 1]; in bitrv2conj()
1255 a[j1] = yr; in bitrv2conj()
1256 a[j1 + 1] = yi; in bitrv2conj()
1257 a[k1] = xr; in bitrv2conj()
1258 a[k1 + 1] = xi; in bitrv2conj()
1261 a[j1 - 1] = -a[j1 - 1]; in bitrv2conj()
1262 xr = a[j1]; in bitrv2conj()
1263 xi = -a[j1 + 1]; in bitrv2conj()
1264 yr = a[k1]; in bitrv2conj()
1265 yi = -a[k1 + 1]; in bitrv2conj()
1266 a[j1] = yr; in bitrv2conj()
1267 a[j1 + 1] = yi; in bitrv2conj()
1268 a[k1] = xr; in bitrv2conj()
1269 a[k1 + 1] = xi; in bitrv2conj()
1270 a[k1 + 3] = -a[k1 + 3]; in bitrv2conj()
1280 xr = a[j]; in bitrv2conj()
1281 xi = -a[j + 1]; in bitrv2conj()
1282 yr = a[k]; in bitrv2conj()
1283 yi = -a[k + 1]; in bitrv2conj()
1284 a[j] = yr; in bitrv2conj()
1285 a[j + 1] = yi; in bitrv2conj()
1286 a[k] = xr; in bitrv2conj()
1287 a[k + 1] = xi; in bitrv2conj()
1290 xr = a[j1]; in bitrv2conj()
1291 xi = -a[j1 + 1]; in bitrv2conj()
1292 yr = a[k1]; in bitrv2conj()
1293 yi = -a[k1 + 1]; in bitrv2conj()
1294 a[j1] = yr; in bitrv2conj()
1295 a[j1 + 1] = yi; in bitrv2conj()
1296 a[k1] = xr; in bitrv2conj()
1297 a[k1 + 1] = xi; in bitrv2conj()
1300 xr = a[j1]; in bitrv2conj()
1301 xi = -a[j1 + 1]; in bitrv2conj()
1302 yr = a[k1]; in bitrv2conj()
1303 yi = -a[k1 + 1]; in bitrv2conj()
1304 a[j1] = yr; in bitrv2conj()
1305 a[j1 + 1] = yi; in bitrv2conj()
1306 a[k1] = xr; in bitrv2conj()
1307 a[k1 + 1] = xi; in bitrv2conj()
1310 xr = a[j1]; in bitrv2conj()
1311 xi = -a[j1 + 1]; in bitrv2conj()
1312 yr = a[k1]; in bitrv2conj()
1313 yi = -a[k1 + 1]; in bitrv2conj()
1314 a[j1] = yr; in bitrv2conj()
1315 a[j1 + 1] = yi; in bitrv2conj()
1316 a[k1] = xr; in bitrv2conj()
1317 a[k1 + 1] = xi; in bitrv2conj()
1320 xr = a[j1]; in bitrv2conj()
1321 xi = -a[j1 + 1]; in bitrv2conj()
1322 yr = a[k1]; in bitrv2conj()
1323 yi = -a[k1 + 1]; in bitrv2conj()
1324 a[j1] = yr; in bitrv2conj()
1325 a[j1 + 1] = yi; in bitrv2conj()
1326 a[k1] = xr; in bitrv2conj()
1327 a[k1 + 1] = xi; in bitrv2conj()
1330 xr = a[j1]; in bitrv2conj()
1331 xi = -a[j1 + 1]; in bitrv2conj()
1332 yr = a[k1]; in bitrv2conj()
1333 yi = -a[k1 + 1]; in bitrv2conj()
1334 a[j1] = yr; in bitrv2conj()
1335 a[j1 + 1] = yi; in bitrv2conj()
1336 a[k1] = xr; in bitrv2conj()
1337 a[k1 + 1] = xi; in bitrv2conj()
1340 xr = a[j1]; in bitrv2conj()
1341 xi = -a[j1 + 1]; in bitrv2conj()
1342 yr = a[k1]; in bitrv2conj()
1343 yi = -a[k1 + 1]; in bitrv2conj()
1344 a[j1] = yr; in bitrv2conj()
1345 a[j1 + 1] = yi; in bitrv2conj()
1346 a[k1] = xr; in bitrv2conj()
1347 a[k1 + 1] = xi; in bitrv2conj()
1350 xr = a[j1]; in bitrv2conj()
1351 xi = -a[j1 + 1]; in bitrv2conj()
1352 yr = a[k1]; in bitrv2conj()
1353 yi = -a[k1 + 1]; in bitrv2conj()
1354 a[j1] = yr; in bitrv2conj()
1355 a[j1 + 1] = yi; in bitrv2conj()
1356 a[k1] = xr; in bitrv2conj()
1357 a[k1 + 1] = xi; in bitrv2conj()
1365 a[j1 - 1] = -a[j1 - 1]; in bitrv2conj()
1366 xr = a[j1]; in bitrv2conj()
1367 xi = -a[j1 + 1]; in bitrv2conj()
1368 yr = a[k1]; in bitrv2conj()
1369 yi = -a[k1 + 1]; in bitrv2conj()
1370 a[j1] = yr; in bitrv2conj()
1371 a[j1 + 1] = yi; in bitrv2conj()
1372 a[k1] = xr; in bitrv2conj()
1373 a[k1 + 1] = xi; in bitrv2conj()
1374 a[k1 + 3] = -a[k1 + 3]; in bitrv2conj()
1377 a[j1 - 1] = -a[j1 - 1]; in bitrv2conj()
1378 xr = a[j1]; in bitrv2conj()
1379 xi = -a[j1 + 1]; in bitrv2conj()
1380 yr = a[k1]; in bitrv2conj()
1381 yi = -a[k1 + 1]; in bitrv2conj()
1382 a[j1] = yr; in bitrv2conj()
1383 a[j1 + 1] = yi; in bitrv2conj()
1384 a[k1] = xr; in bitrv2conj()
1385 a[k1 + 1] = xi; in bitrv2conj()
1386 a[k1 + 3] = -a[k1 + 3]; in bitrv2conj()
1395 void bitrv216(PICOFFTSG_FFTTYPE *a) in bitrv216() argument
1401 x1r = a[2]; in bitrv216()
1402 x1i = a[3]; in bitrv216()
1403 x2r = a[4]; in bitrv216()
1404 x2i = a[5]; in bitrv216()
1405 x3r = a[6]; in bitrv216()
1406 x3i = a[7]; in bitrv216()
1407 x4r = a[8]; in bitrv216()
1408 x4i = a[9]; in bitrv216()
1409 x5r = a[10]; in bitrv216()
1410 x5i = a[11]; in bitrv216()
1411 x7r = a[14]; in bitrv216()
1412 x7i = a[15]; in bitrv216()
1413 x8r = a[16]; in bitrv216()
1414 x8i = a[17]; in bitrv216()
1415 x10r = a[20]; in bitrv216()
1416 x10i = a[21]; in bitrv216()
1417 x11r = a[22]; in bitrv216()
1418 x11i = a[23]; in bitrv216()
1419 x12r = a[24]; in bitrv216()
1420 x12i = a[25]; in bitrv216()
1421 x13r = a[26]; in bitrv216()
1422 x13i = a[27]; in bitrv216()
1423 x14r = a[28]; in bitrv216()
1424 x14i = a[29]; in bitrv216()
1425 a[2] = x8r; in bitrv216()
1426 a[3] = x8i; in bitrv216()
1427 a[4] = x4r; in bitrv216()
1428 a[5] = x4i; in bitrv216()
1429 a[6] = x12r; in bitrv216()
1430 a[7] = x12i; in bitrv216()
1431 a[8] = x2r; in bitrv216()
1432 a[9] = x2i; in bitrv216()
1433 a[10] = x10r; in bitrv216()
1434 a[11] = x10i; in bitrv216()
1435 a[14] = x14r; in bitrv216()
1436 a[15] = x14i; in bitrv216()
1437 a[16] = x1r; in bitrv216()
1438 a[17] = x1i; in bitrv216()
1439 a[20] = x5r; in bitrv216()
1440 a[21] = x5i; in bitrv216()
1441 a[22] = x13r; in bitrv216()
1442 a[23] = x13i; in bitrv216()
1443 a[24] = x3r; in bitrv216()
1444 a[25] = x3i; in bitrv216()
1445 a[26] = x11r; in bitrv216()
1446 a[27] = x11i; in bitrv216()
1447 a[28] = x7r; in bitrv216()
1448 a[29] = x7i; in bitrv216()
1452 void bitrv216neg(PICOFFTSG_FFTTYPE *a) in bitrv216neg() argument
1459 x1r = a[2]; in bitrv216neg()
1460 x1i = a[3]; in bitrv216neg()
1461 x2r = a[4]; in bitrv216neg()
1462 x2i = a[5]; in bitrv216neg()
1463 x3r = a[6]; in bitrv216neg()
1464 x3i = a[7]; in bitrv216neg()
1465 x4r = a[8]; in bitrv216neg()
1466 x4i = a[9]; in bitrv216neg()
1467 x5r = a[10]; in bitrv216neg()
1468 x5i = a[11]; in bitrv216neg()
1469 x6r = a[12]; in bitrv216neg()
1470 x6i = a[13]; in bitrv216neg()
1471 x7r = a[14]; in bitrv216neg()
1472 x7i = a[15]; in bitrv216neg()
1473 x8r = a[16]; in bitrv216neg()
1474 x8i = a[17]; in bitrv216neg()
1475 x9r = a[18]; in bitrv216neg()
1476 x9i = a[19]; in bitrv216neg()
1477 x10r = a[20]; in bitrv216neg()
1478 x10i = a[21]; in bitrv216neg()
1479 x11r = a[22]; in bitrv216neg()
1480 x11i = a[23]; in bitrv216neg()
1481 x12r = a[24]; in bitrv216neg()
1482 x12i = a[25]; in bitrv216neg()
1483 x13r = a[26]; in bitrv216neg()
1484 x13i = a[27]; in bitrv216neg()
1485 x14r = a[28]; in bitrv216neg()
1486 x14i = a[29]; in bitrv216neg()
1487 x15r = a[30]; in bitrv216neg()
1488 x15i = a[31]; in bitrv216neg()
1489 a[2] = x15r; in bitrv216neg()
1490 a[3] = x15i; in bitrv216neg()
1491 a[4] = x7r; in bitrv216neg()
1492 a[5] = x7i; in bitrv216neg()
1493 a[6] = x11r; in bitrv216neg()
1494 a[7] = x11i; in bitrv216neg()
1495 a[8] = x3r; in bitrv216neg()
1496 a[9] = x3i; in bitrv216neg()
1497 a[10] = x13r; in bitrv216neg()
1498 a[11] = x13i; in bitrv216neg()
1499 a[12] = x5r; in bitrv216neg()
1500 a[13] = x5i; in bitrv216neg()
1501 a[14] = x9r; in bitrv216neg()
1502 a[15] = x9i; in bitrv216neg()
1503 a[16] = x1r; in bitrv216neg()
1504 a[17] = x1i; in bitrv216neg()
1505 a[18] = x14r; in bitrv216neg()
1506 a[19] = x14i; in bitrv216neg()
1507 a[20] = x6r; in bitrv216neg()
1508 a[21] = x6i; in bitrv216neg()
1509 a[22] = x10r; in bitrv216neg()
1510 a[23] = x10i; in bitrv216neg()
1511 a[24] = x2r; in bitrv216neg()
1512 a[25] = x2i; in bitrv216neg()
1513 a[26] = x12r; in bitrv216neg()
1514 a[27] = x12i; in bitrv216neg()
1515 a[28] = x4r; in bitrv216neg()
1516 a[29] = x4i; in bitrv216neg()
1517 a[30] = x8r; in bitrv216neg()
1518 a[31] = x8i; in bitrv216neg()
1522 void bitrv208(PICOFFTSG_FFTTYPE *a) in bitrv208() argument
1526 x1r = a[2]; in bitrv208()
1527 x1i = a[3]; in bitrv208()
1528 x3r = a[6]; in bitrv208()
1529 x3i = a[7]; in bitrv208()
1530 x4r = a[8]; in bitrv208()
1531 x4i = a[9]; in bitrv208()
1532 x6r = a[12]; in bitrv208()
1533 x6i = a[13]; in bitrv208()
1534 a[2] = x4r; in bitrv208()
1535 a[3] = x4i; in bitrv208()
1536 a[6] = x6r; in bitrv208()
1537 a[7] = x6i; in bitrv208()
1538 a[8] = x1r; in bitrv208()
1539 a[9] = x1i; in bitrv208()
1540 a[12] = x3r; in bitrv208()
1541 a[13] = x3i; in bitrv208()
1545 void bitrv208neg(PICOFFTSG_FFTTYPE *a) in bitrv208neg() argument
1550 x1r = a[2]; in bitrv208neg()
1551 x1i = a[3]; in bitrv208neg()
1552 x2r = a[4]; in bitrv208neg()
1553 x2i = a[5]; in bitrv208neg()
1554 x3r = a[6]; in bitrv208neg()
1555 x3i = a[7]; in bitrv208neg()
1556 x4r = a[8]; in bitrv208neg()
1557 x4i = a[9]; in bitrv208neg()
1558 x5r = a[10]; in bitrv208neg()
1559 x5i = a[11]; in bitrv208neg()
1560 x6r = a[12]; in bitrv208neg()
1561 x6i = a[13]; in bitrv208neg()
1562 x7r = a[14]; in bitrv208neg()
1563 x7i = a[15]; in bitrv208neg()
1564 a[2] = x7r; in bitrv208neg()
1565 a[3] = x7i; in bitrv208neg()
1566 a[4] = x3r; in bitrv208neg()
1567 a[5] = x3i; in bitrv208neg()
1568 a[6] = x5r; in bitrv208neg()
1569 a[7] = x5i; in bitrv208neg()
1570 a[8] = x1r; in bitrv208neg()
1571 a[9] = x1i; in bitrv208neg()
1572 a[10] = x6r; in bitrv208neg()
1573 a[11] = x6i; in bitrv208neg()
1574 a[12] = x2r; in bitrv208neg()
1575 a[13] = x2i; in bitrv208neg()
1576 a[14] = x4r; in bitrv208neg()
1577 a[15] = x4i; in bitrv208neg()
1581 void bitrv1(picoos_int32 n, PICOFFTSG_FFTTYPE *a) in bitrv1() argument
1586 x = a[1]; in bitrv1()
1587 a[1] = a[nh]; in bitrv1()
1588 a[nh] = x; in bitrv1()
1596 x = a[j1]; in bitrv1()
1597 a[j1] = a[k1]; in bitrv1()
1598 a[k1] = x; in bitrv1()
1606 x = a[j]; in bitrv1()
1607 a[j] = a[k]; in bitrv1()
1608 a[k] = x; in bitrv1()
1611 x = a[j1]; in bitrv1()
1612 a[j1] = a[k1]; in bitrv1()
1613 a[k1] = x; in bitrv1()
1616 x = a[j1]; in bitrv1()
1617 a[j1] = a[k1]; in bitrv1()
1618 a[k1] = x; in bitrv1()
1621 x = a[j1]; in bitrv1()
1622 a[j1] = a[k1]; in bitrv1()
1623 a[k1] = x; in bitrv1()
1626 x = a[j1]; in bitrv1()
1627 a[j1] = a[k1]; in bitrv1()
1628 a[k1] = x; in bitrv1()
1631 x = a[j1]; in bitrv1()
1632 a[j1] = a[k1]; in bitrv1()
1633 a[k1] = x; in bitrv1()
1636 x = a[j1]; in bitrv1()
1637 a[j1] = a[k1]; in bitrv1()
1638 a[k1] = x; in bitrv1()
1641 x = a[j1]; in bitrv1()
1642 a[j1] = a[k1]; in bitrv1()
1643 a[k1] = x; in bitrv1()
1651 x = a[j1]; in bitrv1()
1652 a[j1] = a[k1]; in bitrv1()
1653 a[k1] = x; in bitrv1()
1656 x = a[j1]; in bitrv1()
1657 a[j1] = a[k1]; in bitrv1()
1658 a[k1] = x; in bitrv1()
1668 x = a[j]; in bitrv1()
1669 a[j] = a[k]; in bitrv1()
1670 a[k] = x; in bitrv1()
1673 x = a[j1]; in bitrv1()
1674 a[j1] = a[k1]; in bitrv1()
1675 a[k1] = x; in bitrv1()
1678 x = a[j1]; in bitrv1()
1679 a[j1] = a[k1]; in bitrv1()
1680 a[k1] = x; in bitrv1()
1683 x = a[j1]; in bitrv1()
1684 a[j1] = a[k1]; in bitrv1()
1685 a[k1] = x; in bitrv1()
1693 x = a[j1]; in bitrv1()
1694 a[j1] = a[k1]; in bitrv1()
1695 a[k1] = x; in bitrv1()
1705 void cftb1st(picoos_int32 n, PICOFFTSG_FFTTYPE *a) in cftb1st() argument
1717 x0r = a[0] + a[j2]; in cftb1st()
1718 x0i = -a[1] - a[j2 + 1]; in cftb1st()
1719 x1r = a[0] - a[j2]; in cftb1st()
1720 x1i = -a[1] + a[j2 + 1]; in cftb1st()
1721 x2r = a[j1] + a[j3]; in cftb1st()
1722 x2i = a[j1 + 1] + a[j3 + 1]; in cftb1st()
1723 x3r = a[j1] - a[j3]; in cftb1st()
1724 x3i = a[j1 + 1] - a[j3 + 1]; in cftb1st()
1725 a[0] = x0r + x2r; in cftb1st()
1726 a[1] = x0i - x2i; in cftb1st()
1727 a[j1] = x0r - x2r; in cftb1st()
1728 a[j1 + 1] = x0i + x2i; in cftb1st()
1729 a[j2] = x1r + x3i; in cftb1st()
1730 a[j2 + 1] = x1i + x3r; in cftb1st()
1731 a[j3] = x1r - x3i; in cftb1st()
1732 a[j3 + 1] = x1i - x3r; in cftb1st()
1761 x0r = a[j] + a[j2]; in cftb1st()
1762 x0i = -a[j + 1] - a[j2 + 1]; in cftb1st()
1763 x1r = a[j] - a[j2]; in cftb1st()
1764 x1i = -a[j + 1] + a[j2 + 1]; in cftb1st()
1765 x2r = a[j1] + a[j3]; in cftb1st()
1766 x2i = a[j1 + 1] + a[j3 + 1]; in cftb1st()
1767 x3r = a[j1] - a[j3]; in cftb1st()
1768 x3i = a[j1 + 1] - a[j3 + 1]; in cftb1st()
1769 a[j] = x0r + x2r; in cftb1st()
1770 a[j + 1] = x0i - x2i; in cftb1st()
1771 a[j1] = x0r - x2r; in cftb1st()
1772 a[j1 + 1] = x0i + x2i; in cftb1st()
1775 a[j2] = Mult_W_W(wk1r, x0r) - Mult_W_W(wk1i, x0i); in cftb1st()
1776 a[j2 + 1] = Mult_W_W(wk1r, x0i) + Mult_W_W(wk1i, x0r); in cftb1st()
1779 a[j3] = Mult_W_W(wk3r, x0r) + Mult_W_W(wk3i, x0i); in cftb1st()
1780 a[j3 + 1] = Mult_W_W(wk3r, x0i) - Mult_W_W(wk3i, x0r); in cftb1st()
1781 x0r = a[j + 2] + a[j2 + 2]; in cftb1st()
1782 x0i = -a[j + 3] - a[j2 + 3]; in cftb1st()
1783 x1r = a[j + 2] - a[j2 + 2]; in cftb1st()
1784 x1i = -a[j + 3] + a[j2 + 3]; in cftb1st()
1785 x2r = a[j1 + 2] + a[j3 + 2]; in cftb1st()
1786 x2i = a[j1 + 3] + a[j3 + 3]; in cftb1st()
1787 x3r = a[j1 + 2] - a[j3 + 2]; in cftb1st()
1788 x3i = a[j1 + 3] - a[j3 + 3]; in cftb1st()
1789 a[j + 2] = x0r + x2r; in cftb1st()
1790 a[j + 3] = x0i - x2i; in cftb1st()
1791 a[j1 + 2] = x0r - x2r; in cftb1st()
1792 a[j1 + 3] = x0i + x2i; in cftb1st()
1795 a[j2 + 2] = Mult_W_W(wd1r, x0r) - Mult_W_W(wd1i, x0i); in cftb1st()
1796 a[j2 + 3] = Mult_W_W(wd1r, x0i) + Mult_W_W(wd1i, x0r); in cftb1st()
1799 a[j3 + 2] = Mult_W_W(wd3r, x0r) + Mult_W_W(wd3i, x0i); in cftb1st()
1800 a[j3 + 3] = Mult_W_W(wd3r, x0i) - Mult_W_W(wd3i, x0r); in cftb1st()
1805 x0r = a[j0] + a[j2]; in cftb1st()
1806 x0i = -a[j0 + 1] - a[j2 + 1]; in cftb1st()
1807 x1r = a[j0] - a[j2]; in cftb1st()
1808 x1i = -a[j0 + 1] + a[j2 + 1]; in cftb1st()
1809 x2r = a[j1] + a[j3]; in cftb1st()
1810 x2i = a[j1 + 1] + a[j3 + 1]; in cftb1st()
1811 x3r = a[j1] - a[j3]; in cftb1st()
1812 x3i = a[j1 + 1] - a[j3 + 1]; in cftb1st()
1813 a[j0] = x0r + x2r; in cftb1st()
1814 a[j0 + 1] = x0i - x2i; in cftb1st()
1815 a[j1] = x0r - x2r; in cftb1st()
1816 a[j1 + 1] = x0i + x2i; in cftb1st()
1819 a[j2] = Mult_W_W(wk1i, x0r) - Mult_W_W(wk1r, x0i); in cftb1st()
1820 a[j2 + 1] = Mult_W_W(wk1i, x0i) + Mult_W_W(wk1r, x0r); in cftb1st()
1823 a[j3] = Mult_W_W(wk3i, x0r) + Mult_W_W(wk3r, x0i); in cftb1st()
1824 a[j3 + 1] = Mult_W_W(wk3i, x0i) - Mult_W_W(wk3r, x0r); in cftb1st()
1825 x0r = a[j0 - 2] + a[j2 - 2]; in cftb1st()
1826 x0i = -a[j0 - 1] - a[j2 - 1]; in cftb1st()
1827 x1r = a[j0 - 2] - a[j2 - 2]; in cftb1st()
1828 x1i = -a[j0 - 1] + a[j2 - 1]; in cftb1st()
1829 x2r = a[j1 - 2] + a[j3 - 2]; in cftb1st()
1830 x2i = a[j1 - 1] + a[j3 - 1]; in cftb1st()
1831 x3r = a[j1 - 2] - a[j3 - 2]; in cftb1st()
1832 x3i = a[j1 - 1] - a[j3 - 1]; in cftb1st()
1833 a[j0 - 2] = x0r + x2r; in cftb1st()
1834 a[j0 - 1] = x0i - x2i; in cftb1st()
1835 a[j1 - 2] = x0r - x2r; in cftb1st()
1836 a[j1 - 1] = x0i + x2i; in cftb1st()
1839 a[j2 - 2] = Mult_W_W(wd1i, x0r) - Mult_W_W(wd1r, x0i); in cftb1st()
1840 a[j2 - 1] = Mult_W_W(wd1i, x0i) + Mult_W_W(wd1r, x0r); in cftb1st()
1843 a[j3 - 2] = Mult_W_W(wd3i, x0r) + Mult_W_W(wd3r, x0i); in cftb1st()
1844 a[j3 - 1] = Mult_W_W(wd3i, x0i) - Mult_W_W(wd3r, x0r); in cftb1st()
1859 x0r = a[j0 - 2] + a[j2 - 2]; in cftb1st()
1860 x0i = -a[j0 - 1] - a[j2 - 1]; in cftb1st()
1861 x1r = a[j0 - 2] - a[j2 - 2]; in cftb1st()
1862 x1i = -a[j0 - 1] + a[j2 - 1]; in cftb1st()
1863 x2r = a[j1 - 2] + a[j3 - 2]; in cftb1st()
1864 x2i = a[j1 - 1] + a[j3 - 1]; in cftb1st()
1865 x3r = a[j1 - 2] - a[j3 - 2]; in cftb1st()
1866 x3i = a[j1 - 1] - a[j3 - 1]; in cftb1st()
1867 a[j0 - 2] = x0r + x2r; in cftb1st()
1868 a[j0 - 1] = x0i - x2i; in cftb1st()
1869 a[j1 - 2] = x0r - x2r; in cftb1st()
1870 a[j1 - 1] = x0i + x2i; in cftb1st()
1873 a[j2 - 2] = Mult_W_W(wk1r, x0r) - Mult_W_W(wk1i, x0i); in cftb1st()
1874 a[j2 - 1] = Mult_W_W(wk1r, x0i) + Mult_W_W(wk1i, x0r); in cftb1st()
1877 a[j3 - 2] = Mult_W_W(wk3r, x0r) + Mult_W_W(wk3i, x0i); in cftb1st()
1878 a[j3 - 1] = Mult_W_W(wk3r, x0i) - Mult_W_W(wk3i, x0r); in cftb1st()
1879 x0r = a[j0] + a[j2]; in cftb1st()
1880 x0i = -a[j0 + 1] - a[j2 + 1]; in cftb1st()
1881 x1r = a[j0] - a[j2]; in cftb1st()
1882 x1i = -a[j0 + 1] + a[j2 + 1]; in cftb1st()
1883 x2r = a[j1] + a[j3]; in cftb1st()
1884 x2i = a[j1 + 1] + a[j3 + 1]; in cftb1st()
1885 x3r = a[j1] - a[j3]; in cftb1st()
1886 x3i = a[j1 + 1] - a[j3 + 1]; in cftb1st()
1887 a[j0] = x0r + x2r; in cftb1st()
1888 a[j0 + 1] = x0i - x2i; in cftb1st()
1889 a[j1] = x0r - x2r; in cftb1st()
1890 a[j1 + 1] = x0i + x2i; in cftb1st()
1893 a[j2] = picofftsg_mult_w_a(wd1r, (x0r - x0i)); in cftb1st()
1894 a[j2 + 1] = picofftsg_mult_w_a(wd1r, (x0i + x0r)); in cftb1st()
1897 a[j3] = -picofftsg_mult_w_a(wd1r, (x0r + x0i)); in cftb1st()
1898 a[j3 + 1] = -picofftsg_mult_w_a(wd1r, (x0i - x0r)); in cftb1st()
1899 x0r = a[j0 + 2] + a[j2 + 2]; in cftb1st()
1900 x0i = -a[j0 + 3] - a[j2 + 3]; in cftb1st()
1901 x1r = a[j0 + 2] - a[j2 + 2]; in cftb1st()
1902 x1i = -a[j0 + 3] + a[j2 + 3]; in cftb1st()
1903 x2r = a[j1 + 2] + a[j3 + 2]; in cftb1st()
1904 x2i = a[j1 + 3] + a[j3 + 3]; in cftb1st()
1905 x3r = a[j1 + 2] - a[j3 + 2]; in cftb1st()
1906 x3i = a[j1 + 3] - a[j3 + 3]; in cftb1st()
1907 a[j0 + 2] = x0r + x2r; in cftb1st()
1908 a[j0 + 3] = x0i - x2i; in cftb1st()
1909 a[j1 + 2] = x0r - x2r; in cftb1st()
1910 a[j1 + 3] = x0i + x2i; in cftb1st()
1913 a[j2 + 2] = Mult_W_W(wk1i, x0r) - Mult_W_W(wk1r, x0i); in cftb1st()
1914 a[j2 + 3] = Mult_W_W(wk1i, x0i) + Mult_W_W(wk1r, x0r); in cftb1st()
1917 a[j3 + 2] = Mult_W_W(wk3i, x0r) + Mult_W_W(wk3r, x0i); in cftb1st()
1918 a[j3 + 3] = Mult_W_W(wk3i, x0i) - Mult_W_W(wk3r, x0r); in cftb1st()
1921 void cftrec4(picoos_int32 n, PICOFFTSG_FFTTYPE *a) in cftrec4() argument
1928 cftmdl1(m, &a[n - m]); in cftrec4()
1930 cftleaf(m, 1, &a[n - m]); in cftrec4()
1934 isplt = cfttree(m, j, k, a); in cftrec4()
1935 cftleaf(m, isplt, &a[j - m]); in cftrec4()
1940 picoos_int32 cfttree(picoos_int32 n, picoos_int32 j, picoos_int32 k, PICOFFTSG_FFTTYPE *a) in cfttree() argument
1947 cftmdl1(n, &a[j - n]); in cfttree()
1949 cftmdl2(n, &a[j - n]); in cfttree()
1959 cftmdl1(m, &a[j - m]); in cfttree()
1964 cftmdl2(m, &a[j - m]); in cfttree()
1973 void cftleaf(picoos_int32 n, picoos_int32 isplt, PICOFFTSG_FFTTYPE *a) in cftleaf() argument
1977 cftmdl1(128, a); in cftleaf()
1978 cftf161(a); in cftleaf()
1979 cftf162(&a[32]); in cftleaf()
1980 cftf161(&a[64]); in cftleaf()
1981 cftf161(&a[96]); in cftleaf()
1982 cftmdl2(128, &a[128]); in cftleaf()
1983 cftf161(&a[128]); in cftleaf()
1984 cftf162(&a[160]); in cftleaf()
1985 cftf161(&a[192]); in cftleaf()
1986 cftf162(&a[224]); in cftleaf()
1987 cftmdl1(128, &a[256]); in cftleaf()
1988 cftf161(&a[256]); in cftleaf()
1989 cftf162(&a[288]); in cftleaf()
1990 cftf161(&a[320]); in cftleaf()
1991 cftf161(&a[352]); in cftleaf()
1993 cftmdl1(128, &a[384]); in cftleaf()
1994 cftf161(&a[480]); in cftleaf()
1996 cftmdl2(128, &a[384]); in cftleaf()
1997 cftf162(&a[480]); in cftleaf()
1999 cftf161(&a[384]); in cftleaf()
2000 cftf162(&a[416]); in cftleaf()
2001 cftf161(&a[448]); in cftleaf()
2003 cftmdl1(64, a); in cftleaf()
2004 cftf081(a); in cftleaf()
2005 cftf082(&a[16]); in cftleaf()
2006 cftf081(&a[32]); in cftleaf()
2007 cftf081(&a[48]); in cftleaf()
2008 cftmdl2(64, &a[64]); in cftleaf()
2009 cftf081(&a[64]); in cftleaf()
2010 cftf082(&a[80]); in cftleaf()
2011 cftf081(&a[96]); in cftleaf()
2012 cftf082(&a[112]); in cftleaf()
2013 cftmdl1(64, &a[128]); in cftleaf()
2014 cftf081(&a[128]); in cftleaf()
2015 cftf082(&a[144]); in cftleaf()
2016 cftf081(&a[160]); in cftleaf()
2017 cftf081(&a[176]); in cftleaf()
2019 cftmdl1(64, &a[192]); in cftleaf()
2020 cftf081(&a[240]); in cftleaf()
2022 cftmdl2(64, &a[192]); in cftleaf()
2023 cftf082(&a[240]); in cftleaf()
2025 cftf081(&a[192]); in cftleaf()
2026 cftf082(&a[208]); in cftleaf()
2027 cftf081(&a[224]); in cftleaf()
2032 void cftmdl1(picoos_int32 n, PICOFFTSG_FFTTYPE *a) in cftmdl1() argument
2044 x0r = a[0] + a[j2]; in cftmdl1()
2045 x0i = a[1] + a[j2 + 1]; in cftmdl1()
2046 x1r = a[0] - a[j2]; in cftmdl1()
2047 x1i = a[1] - a[j2 + 1]; in cftmdl1()
2048 x2r = a[j1] + a[j3]; in cftmdl1()
2049 x2i = a[j1 + 1] + a[j3 + 1]; in cftmdl1()
2050 x3r = a[j1] - a[j3]; in cftmdl1()
2051 x3i = a[j1 + 1] - a[j3 + 1]; in cftmdl1()
2052 a[0] = x0r + x2r; in cftmdl1()
2053 a[1] = x0i + x2i; in cftmdl1()
2054 a[j1] = x0r - x2r; in cftmdl1()
2055 a[j1 + 1] = x0i - x2i; in cftmdl1()
2056 a[j2] = x1r - x3i; in cftmdl1()
2057 a[j2 + 1] = x1i + x3r; in cftmdl1()
2058 a[j3] = x1r + x3i; in cftmdl1()
2059 a[j3 + 1] = x1i - x3r; in cftmdl1()
2085 x0r = a[j] + a[j2]; in cftmdl1()
2086 x0i = a[j + 1] + a[j2 + 1]; in cftmdl1()
2087 x1r = a[j] - a[j2]; in cftmdl1()
2088 x1i = a[j + 1] - a[j2 + 1]; in cftmdl1()
2089 x2r = a[j1] + a[j3]; in cftmdl1()
2090 x2i = a[j1 + 1] + a[j3 + 1]; in cftmdl1()
2091 x3r = a[j1] - a[j3]; in cftmdl1()
2092 x3i = a[j1 + 1] - a[j3 + 1]; in cftmdl1()
2093 a[j] = x0r + x2r; in cftmdl1()
2094 a[j + 1] = x0i + x2i; in cftmdl1()
2095 a[j1] = x0r - x2r; in cftmdl1()
2096 a[j1 + 1] = x0i - x2i; in cftmdl1()
2099 a[j2] = Mult_W_W(wk1r, x0r) - Mult_W_W(wk1i, x0i); in cftmdl1()
2100 a[j2 + 1] = Mult_W_W(wk1r, x0i) + Mult_W_W(wk1i, x0r); in cftmdl1()
2103 a[j3] = Mult_W_W(wk3r, x0r) + Mult_W_W(wk3i, x0i); in cftmdl1()
2104 a[j3 + 1] = Mult_W_W(wk3r, x0i) - Mult_W_W(wk3i, x0r); in cftmdl1()
2105 x0r = a[j + 2] + a[j2 + 2]; in cftmdl1()
2106 x0i = a[j + 3] + a[j2 + 3]; in cftmdl1()
2107 x1r = a[j + 2] - a[j2 + 2]; in cftmdl1()
2108 x1i = a[j + 3] - a[j2 + 3]; in cftmdl1()
2109 x2r = a[j1 + 2] + a[j3 + 2]; in cftmdl1()
2110 x2i = a[j1 + 3] + a[j3 + 3]; in cftmdl1()
2111 x3r = a[j1 + 2] - a[j3 + 2]; in cftmdl1()
2112 x3i = a[j1 + 3] - a[j3 + 3]; in cftmdl1()
2113 a[j + 2] = x0r + x2r; in cftmdl1()
2114 a[j + 3] = x0i + x2i; in cftmdl1()
2115 a[j1 + 2] = x0r - x2r; in cftmdl1()
2116 a[j1 + 3] = x0i - x2i; in cftmdl1()
2119 a[j2 + 2] = Mult_W_W(wd1r, x0r) - Mult_W_W(wd1i, x0i); in cftmdl1()
2120 a[j2 + 3] = Mult_W_W(wd1r, x0i) + Mult_W_W(wd1i, x0r); in cftmdl1()
2123 a[j3 + 2] = Mult_W_W(wd3r, x0r) + Mult_W_W(wd3i, x0i); in cftmdl1()
2124 a[j3 + 3] = Mult_W_W(wd3r, x0i) - Mult_W_W(wd3i, x0r); in cftmdl1()
2129 x0r = a[j0] + a[j2]; in cftmdl1()
2130 x0i = a[j0 + 1] + a[j2 + 1]; in cftmdl1()
2131 x1r = a[j0] - a[j2]; in cftmdl1()
2132 x1i = a[j0 + 1] - a[j2 + 1]; in cftmdl1()
2133 x2r = a[j1] + a[j3]; in cftmdl1()
2134 x2i = a[j1 + 1] + a[j3 + 1]; in cftmdl1()
2135 x3r = a[j1] - a[j3]; in cftmdl1()
2136 x3i = a[j1 + 1] - a[j3 + 1]; in cftmdl1()
2137 a[j0] = x0r + x2r; in cftmdl1()
2138 a[j0 + 1] = x0i + x2i; in cftmdl1()
2139 a[j1] = x0r - x2r; in cftmdl1()
2140 a[j1 + 1] = x0i - x2i; in cftmdl1()
2143 a[j2] = Mult_W_W(wk1i, x0r) - Mult_W_W(wk1r, x0i); in cftmdl1()
2144 a[j2 + 1] = Mult_W_W(wk1i, x0i) + Mult_W_W(wk1r, x0r); in cftmdl1()
2147 a[j3] = Mult_W_W(wk3i, x0r) + Mult_W_W(wk3r, x0i); in cftmdl1()
2148 a[j3 + 1] = Mult_W_W(wk3i, x0i) - Mult_W_W(wk3r, x0r); in cftmdl1()
2149 x0r = a[j0 - 2] + a[j2 - 2]; in cftmdl1()
2150 x0i = a[j0 - 1] + a[j2 - 1]; in cftmdl1()
2151 x1r = a[j0 - 2] - a[j2 - 2]; in cftmdl1()
2152 x1i = a[j0 - 1] - a[j2 - 1]; in cftmdl1()
2153 x2r = a[j1 - 2] + a[j3 - 2]; in cftmdl1()
2154 x2i = a[j1 - 1] + a[j3 - 1]; in cftmdl1()
2155 x3r = a[j1 - 2] - a[j3 - 2]; in cftmdl1()
2156 x3i = a[j1 - 1] - a[j3 - 1]; in cftmdl1()
2157 a[j0 - 2] = x0r + x2r; in cftmdl1()
2158 a[j0 - 1] = x0i + x2i; in cftmdl1()
2159 a[j1 - 2] = x0r - x2r; in cftmdl1()
2160 a[j1 - 1] = x0i - x2i; in cftmdl1()
2163 a[j2 - 2] = Mult_W_W(wd1i, x0r) - Mult_W_W(wd1r, x0i); in cftmdl1()
2164 a[j2 - 1] = Mult_W_W(wd1i, x0i) + Mult_W_W(wd1r, x0r); in cftmdl1()
2167 a[j3 - 2] = Mult_W_W(wd3i, x0r) + Mult_W_W(wd3r, x0i); in cftmdl1()
2168 a[j3 - 1] = Mult_W_W(wd3i, x0i) - Mult_W_W(wd3r, x0r); in cftmdl1()
2183 x0r = a[j0 - 2] + a[j2 - 2]; in cftmdl1()
2184 x0i = a[j0 - 1] + a[j2 - 1]; in cftmdl1()
2185 x1r = a[j0 - 2] - a[j2 - 2]; in cftmdl1()
2186 x1i = a[j0 - 1] - a[j2 - 1]; in cftmdl1()
2187 x2r = a[j1 - 2] + a[j3 - 2]; in cftmdl1()
2188 x2i = a[j1 - 1] + a[j3 - 1]; in cftmdl1()
2189 x3r = a[j1 - 2] - a[j3 - 2]; in cftmdl1()
2190 x3i = a[j1 - 1] - a[j3 - 1]; in cftmdl1()
2191 a[j0 - 2] = x0r + x2r; in cftmdl1()
2192 a[j0 - 1] = x0i + x2i; in cftmdl1()
2193 a[j1 - 2] = x0r - x2r; in cftmdl1()
2194 a[j1 - 1] = x0i - x2i; in cftmdl1()
2197 a[j2 - 2] = Mult_W_W(wk1r, x0r) - Mult_W_W(wk1i, x0i); in cftmdl1()
2198 a[j2 - 1] = Mult_W_W(wk1r, x0i) + Mult_W_W(wk1i, x0r); in cftmdl1()
2201 a[j3 - 2] = Mult_W_W(wk3r, x0r) + Mult_W_W(wk3i, x0i); in cftmdl1()
2202 a[j3 - 1] = Mult_W_W(wk3r, x0i) - Mult_W_W(wk3i, x0r); in cftmdl1()
2203 x0r = a[j0] + a[j2]; in cftmdl1()
2204 x0i = a[j0 + 1] + a[j2 + 1]; in cftmdl1()
2205 x1r = a[j0] - a[j2]; in cftmdl1()
2206 x1i = a[j0 + 1] - a[j2 + 1]; in cftmdl1()
2207 x2r = a[j1] + a[j3]; in cftmdl1()
2208 x2i = a[j1 + 1] + a[j3 + 1]; in cftmdl1()
2209 x3r = a[j1] - a[j3]; in cftmdl1()
2210 x3i = a[j1 + 1] - a[j3 + 1]; in cftmdl1()
2211 a[j0] = x0r + x2r; in cftmdl1()
2212 a[j0 + 1] = x0i + x2i; in cftmdl1()
2213 a[j1] = x0r - x2r; in cftmdl1()
2214 a[j1 + 1] = x0i - x2i; in cftmdl1()
2217 a[j2] = picofftsg_mult_w_a(wd1r, (x0r - x0i)); in cftmdl1()
2218 a[j2 + 1] = picofftsg_mult_w_a(wd1r, (x0i + x0r)); in cftmdl1()
2221 a[j3] = -picofftsg_mult_w_a(wd1r, (x0r + x0i)); in cftmdl1()
2222 a[j3 + 1] = -picofftsg_mult_w_a(wd1r, (x0i - x0r)); in cftmdl1()
2223 x0r = a[j0 + 2] + a[j2 + 2]; in cftmdl1()
2224 x0i = a[j0 + 3] + a[j2 + 3]; in cftmdl1()
2225 x1r = a[j0 + 2] - a[j2 + 2]; in cftmdl1()
2226 x1i = a[j0 + 3] - a[j2 + 3]; in cftmdl1()
2227 x2r = a[j1 + 2] + a[j3 + 2]; in cftmdl1()
2228 x2i = a[j1 + 3] + a[j3 + 3]; in cftmdl1()
2229 x3r = a[j1 + 2] - a[j3 + 2]; in cftmdl1()
2230 x3i = a[j1 + 3] - a[j3 + 3]; in cftmdl1()
2231 a[j0 + 2] = x0r + x2r; in cftmdl1()
2232 a[j0 + 3] = x0i + x2i; in cftmdl1()
2233 a[j1 + 2] = x0r - x2r; in cftmdl1()
2234 a[j1 + 3] = x0i - x2i; in cftmdl1()
2237 a[j2 + 2] = Mult_W_W(wk1i, x0r) - Mult_W_W(wk1r, x0i); in cftmdl1()
2238 a[j2 + 3] = Mult_W_W(wk1i, x0i) + Mult_W_W(wk1r, x0r); in cftmdl1()
2241 a[j3 + 2] = Mult_W_W(wk3i, x0r) + Mult_W_W(wk3r, x0i); in cftmdl1()
2242 a[j3 + 3] = Mult_W_W(wk3i, x0i) - Mult_W_W(wk3r, x0r); in cftmdl1()
2246 void cftmdl2(picoos_int32 n, PICOFFTSG_FFTTYPE *a) in cftmdl2() argument
2260 x0r = a[0] - a[j2 + 1]; in cftmdl2()
2261 x0i = a[1] + a[j2]; in cftmdl2()
2262 x1r = a[0] + a[j2 + 1]; in cftmdl2()
2263 x1i = a[1] - a[j2]; in cftmdl2()
2264 x2r = a[j1] - a[j3 + 1]; in cftmdl2()
2265 x2i = a[j1 + 1] + a[j3]; in cftmdl2()
2266 x3r = a[j1] + a[j3 + 1]; in cftmdl2()
2267 x3i = a[j1 + 1] - a[j3]; in cftmdl2()
2270 a[0] = x0r + y0r; in cftmdl2()
2271 a[1] = x0i + y0i; in cftmdl2()
2272 a[j1] = x0r - y0r; in cftmdl2()
2273 a[j1 + 1] = x0i - y0i; in cftmdl2()
2276 a[j2] = x1r - y0i; in cftmdl2()
2277 a[j2 + 1] = x1i + y0r; in cftmdl2()
2278 a[j3] = x1r + y0i; in cftmdl2()
2279 a[j3 + 1] = x1i - y0r; in cftmdl2()
2318 x0r = a[j] - a[j2 + 1]; in cftmdl2()
2319 x0i = a[j + 1] + a[j2]; in cftmdl2()
2320 x1r = a[j] + a[j2 + 1]; in cftmdl2()
2321 x1i = a[j + 1] - a[j2]; in cftmdl2()
2322 x2r = a[j1] - a[j3 + 1]; in cftmdl2()
2323 x2i = a[j1 + 1] + a[j3]; in cftmdl2()
2324 x3r = a[j1] + a[j3 + 1]; in cftmdl2()
2325 x3i = a[j1 + 1] - a[j3]; in cftmdl2()
2330 a[j] = y0r + y2r; in cftmdl2()
2331 a[j + 1] = y0i + y2i; in cftmdl2()
2332 a[j1] = y0r - y2r; in cftmdl2()
2333 a[j1 + 1] = y0i - y2i; in cftmdl2()
2338 a[j2] = y0r + y2r; in cftmdl2()
2339 a[j2 + 1] = y0i + y2i; in cftmdl2()
2340 a[j3] = y0r - y2r; in cftmdl2()
2341 a[j3 + 1] = y0i - y2i; in cftmdl2()
2342 x0r = a[j + 2] - a[j2 + 3]; in cftmdl2()
2343 x0i = a[j + 3] + a[j2 + 2]; in cftmdl2()
2344 x1r = a[j + 2] + a[j2 + 3]; in cftmdl2()
2345 x1i = a[j + 3] - a[j2 + 2]; in cftmdl2()
2346 x2r = a[j1 + 2] - a[j3 + 3]; in cftmdl2()
2347 x2i = a[j1 + 3] + a[j3 + 2]; in cftmdl2()
2348 x3r = a[j1 + 2] + a[j3 + 3]; in cftmdl2()
2349 x3i = a[j1 + 3] - a[j3 + 2]; in cftmdl2()
2354 a[j + 2] = y0r + y2r; in cftmdl2()
2355 a[j + 3] = y0i + y2i; in cftmdl2()
2356 a[j1 + 2] = y0r - y2r; in cftmdl2()
2357 a[j1 + 3] = y0i - y2i; in cftmdl2()
2362 a[j2 + 2] = y0r + y2r; in cftmdl2()
2363 a[j2 + 3] = y0i + y2i; in cftmdl2()
2364 a[j3 + 2] = y0r - y2r; in cftmdl2()
2365 a[j3 + 3] = y0i - y2i; in cftmdl2()
2370 x0r = a[j0] - a[j2 + 1]; in cftmdl2()
2371 x0i = a[j0 + 1] + a[j2]; in cftmdl2()
2372 x1r = a[j0] + a[j2 + 1]; in cftmdl2()
2373 x1i = a[j0 + 1] - a[j2]; in cftmdl2()
2374 x2r = a[j1] - a[j3 + 1]; in cftmdl2()
2375 x2i = a[j1 + 1] + a[j3]; in cftmdl2()
2376 x3r = a[j1] + a[j3 + 1]; in cftmdl2()
2377 x3i = a[j1 + 1] - a[j3]; in cftmdl2()
2382 a[j0] = y0r + y2r; in cftmdl2()
2383 a[j0 + 1] = y0i + y2i; in cftmdl2()
2384 a[j1] = y0r - y2r; in cftmdl2()
2385 a[j1 + 1] = y0i - y2i; in cftmdl2()
2390 a[j2] = y0r + y2r; in cftmdl2()
2391 a[j2 + 1] = y0i + y2i; in cftmdl2()
2392 a[j3] = y0r - y2r; in cftmdl2()
2393 a[j3 + 1] = y0i - y2i; in cftmdl2()
2394 x0r = a[j0 - 2] - a[j2 - 1]; in cftmdl2()
2395 x0i = a[j0 - 1] + a[j2 - 2]; in cftmdl2()
2396 x1r = a[j0 - 2] + a[j2 - 1]; in cftmdl2()
2397 x1i = a[j0 - 1] - a[j2 - 2]; in cftmdl2()
2398 x2r = a[j1 - 2] - a[j3 - 1]; in cftmdl2()
2399 x2i = a[j1 - 1] + a[j3 - 2]; in cftmdl2()
2400 x3r = a[j1 - 2] + a[j3 - 1]; in cftmdl2()
2401 x3i = a[j1 - 1] - a[j3 - 2]; in cftmdl2()
2406 a[j0 - 2] = y0r + y2r; in cftmdl2()
2407 a[j0 - 1] = y0i + y2i; in cftmdl2()
2408 a[j1 - 2] = y0r - y2r; in cftmdl2()
2409 a[j1 - 1] = y0i - y2i; in cftmdl2()
2414 a[j2 - 2] = y0r + y2r; in cftmdl2()
2415 a[j2 - 1] = y0i + y2i; in cftmdl2()
2416 a[j3 - 2] = y0r - y2r; in cftmdl2()
2417 a[j3 - 1] = y0i - y2i; in cftmdl2()
2437 x0r = a[j0 - 2] - a[j2 - 1]; in cftmdl2()
2438 x0i = a[j0 - 1] + a[j2 - 2]; in cftmdl2()
2439 x1r = a[j0 - 2] + a[j2 - 1]; in cftmdl2()
2440 x1i = a[j0 - 1] - a[j2 - 2]; in cftmdl2()
2441 x2r = a[j1 - 2] - a[j3 - 1]; in cftmdl2()
2442 x2i = a[j1 - 1] + a[j3 - 2]; in cftmdl2()
2443 x3r = a[j1 - 2] + a[j3 - 1]; in cftmdl2()
2444 x3i = a[j1 - 1] - a[j3 - 2]; in cftmdl2()
2449 a[j0 - 2] = y0r + y2r; in cftmdl2()
2450 a[j0 - 1] = y0i + y2i; in cftmdl2()
2451 a[j1 - 2] = y0r - y2r; in cftmdl2()
2452 a[j1 - 1] = y0i - y2i; in cftmdl2()
2457 a[j2 - 2] = y0r + y2r; in cftmdl2()
2458 a[j2 - 1] = y0i + y2i; in cftmdl2()
2459 a[j3 - 2] = y0r - y2r; in cftmdl2()
2460 a[j3 - 1] = y0i - y2i; in cftmdl2()
2461 x0r = a[j0] - a[j2 + 1]; in cftmdl2()
2462 x0i = a[j0 + 1] + a[j2]; in cftmdl2()
2463 x1r = a[j0] + a[j2 + 1]; in cftmdl2()
2464 x1i = a[j0 + 1] - a[j2]; in cftmdl2()
2465 x2r = a[j1] - a[j3 + 1]; in cftmdl2()
2466 x2i = a[j1 + 1] + a[j3]; in cftmdl2()
2467 x3r = a[j1] + a[j3 + 1]; in cftmdl2()
2468 x3i = a[j1 + 1] - a[j3]; in cftmdl2()
2473 a[j0] = y0r + y2r; in cftmdl2()
2474 a[j0 + 1] = y0i + y2i; in cftmdl2()
2475 a[j1] = y0r - y2r; in cftmdl2()
2476 a[j1 + 1] = y0i - y2i; in cftmdl2()
2481 a[j2] = y0r - y2r; in cftmdl2()
2482 a[j2 + 1] = y0i - y2i; in cftmdl2()
2483 a[j3] = y0r + y2r; in cftmdl2()
2484 a[j3 + 1] = y0i + y2i; in cftmdl2()
2485 x0r = a[j0 + 2] - a[j2 + 3]; in cftmdl2()
2486 x0i = a[j0 + 3] + a[j2 + 2]; in cftmdl2()
2487 x1r = a[j0 + 2] + a[j2 + 3]; in cftmdl2()
2488 x1i = a[j0 + 3] - a[j2 + 2]; in cftmdl2()
2489 x2r = a[j1 + 2] - a[j3 + 3]; in cftmdl2()
2490 x2i = a[j1 + 3] + a[j3 + 2]; in cftmdl2()
2491 x3r = a[j1 + 2] + a[j3 + 3]; in cftmdl2()
2492 x3i = a[j1 + 3] - a[j3 + 2]; in cftmdl2()
2497 a[j0 + 2] = y0r + y2r; in cftmdl2()
2498 a[j0 + 3] = y0i + y2i; in cftmdl2()
2499 a[j1 + 2] = y0r - y2r; in cftmdl2()
2500 a[j1 + 3] = y0i - y2i; in cftmdl2()
2505 a[j2 + 2] = y0r + y2r; in cftmdl2()
2506 a[j2 + 3] = y0i + y2i; in cftmdl2()
2507 a[j3 + 2] = y0r - y2r; in cftmdl2()
2508 a[j3 + 3] = y0i - y2i; in cftmdl2()
2512 void cftfx41(picoos_int32 n, PICOFFTSG_FFTTYPE *a) in cftfx41() argument
2516 cftf161(a); in cftfx41()
2517 cftf162(&a[32]); in cftfx41()
2518 cftf161(&a[64]); in cftfx41()
2519 cftf161(&a[96]); in cftfx41()
2521 cftf081(a); in cftfx41()
2522 cftf082(&a[16]); in cftfx41()
2523 cftf081(&a[32]); in cftfx41()
2524 cftf081(&a[48]); in cftfx41()
2529 void cftf161(PICOFFTSG_FFTTYPE *a) in cftf161() argument
2541 x0r = a[0] + a[16]; in cftf161()
2542 x0i = a[1] + a[17]; in cftf161()
2543 x1r = a[0] - a[16]; in cftf161()
2544 x1i = a[1] - a[17]; in cftf161()
2545 x2r = a[8] + a[24]; in cftf161()
2546 x2i = a[9] + a[25]; in cftf161()
2547 x3r = a[8] - a[24]; in cftf161()
2548 x3i = a[9] - a[25]; in cftf161()
2557 x0r = a[2] + a[18]; in cftf161()
2558 x0i = a[3] + a[19]; in cftf161()
2559 x1r = a[2] - a[18]; in cftf161()
2560 x1i = a[3] - a[19]; in cftf161()
2561 x2r = a[10] + a[26]; in cftf161()
2562 x2i = a[11] + a[27]; in cftf161()
2563 x3r = a[10] - a[26]; in cftf161()
2564 x3i = a[11] - a[27]; in cftf161()
2577 x0r = a[4] + a[20]; in cftf161()
2578 x0i = a[5] + a[21]; in cftf161()
2579 x1r = a[4] - a[20]; in cftf161()
2580 x1i = a[5] - a[21]; in cftf161()
2581 x2r = a[12] + a[28]; in cftf161()
2582 x2i = a[13] + a[29]; in cftf161()
2583 x3r = a[12] - a[28]; in cftf161()
2584 x3i = a[13] - a[29]; in cftf161()
2597 x0r = a[6] + a[22]; in cftf161()
2598 x0i = a[7] + a[23]; in cftf161()
2599 x1r = a[6] - a[22]; in cftf161()
2600 x1i = a[7] - a[23]; in cftf161()
2601 x2r = a[14] + a[30]; in cftf161()
2602 x2i = a[15] + a[31]; in cftf161()
2603 x3r = a[14] - a[30]; in cftf161()
2604 x3i = a[15] - a[31]; in cftf161()
2625 a[24] = x0r + x2r; in cftf161()
2626 a[25] = x0i + x2i; in cftf161()
2627 a[26] = x0r - x2r; in cftf161()
2628 a[27] = x0i - x2i; in cftf161()
2629 a[28] = x1r - x3i; in cftf161()
2630 a[29] = x1i + x3r; in cftf161()
2631 a[30] = x1r + x3i; in cftf161()
2632 a[31] = x1i - x3r; in cftf161()
2641 a[16] = x0r + x2r; in cftf161()
2642 a[17] = x0i + x2i; in cftf161()
2643 a[18] = x0r - x2r; in cftf161()
2644 a[19] = x0i - x2i; in cftf161()
2645 a[20] = x1r - x3i; in cftf161()
2646 a[21] = x1i + x3r; in cftf161()
2647 a[22] = x1r + x3i; in cftf161()
2648 a[23] = x1i - x3r; in cftf161()
2661 a[8] = x0r + x2r; in cftf161()
2662 a[9] = x0i + x2i; in cftf161()
2663 a[10] = x0r - x2r; in cftf161()
2664 a[11] = x0i - x2i; in cftf161()
2665 a[12] = x1r - x3i; in cftf161()
2666 a[13] = x1i + x3r; in cftf161()
2667 a[14] = x1r + x3i; in cftf161()
2668 a[15] = x1i - x3r; in cftf161()
2677 a[0] = x0r + x2r; in cftf161()
2678 a[1] = x0i + x2i; in cftf161()
2679 a[2] = x0r - x2r; in cftf161()
2680 a[3] = x0i - x2i; in cftf161()
2681 a[4] = x1r - x3i; in cftf161()
2682 a[5] = x1i + x3r; in cftf161()
2683 a[6] = x1r + x3i; in cftf161()
2684 a[7] = x1i - x3r; in cftf161()
2688 void cftf162(PICOFFTSG_FFTTYPE *a) in cftf162() argument
2704 x1r = a[0] - a[17]; in cftf162()
2705 x1i = a[1] + a[16]; in cftf162()
2706 x0r = a[8] - a[25]; in cftf162()
2707 x0i = a[9] + a[24]; in cftf162()
2714 x1r = a[0] + a[17]; in cftf162()
2715 x1i = a[1] - a[16]; in cftf162()
2716 x0r = a[8] + a[25]; in cftf162()
2717 x0i = a[9] - a[24]; in cftf162()
2724 x0r = a[2] - a[19]; in cftf162()
2725 x0i = a[3] + a[18]; in cftf162()
2728 x0r = a[10] - a[27]; in cftf162()
2729 x0i = a[11] + a[26]; in cftf162()
2736 x0r = a[2] + a[19]; in cftf162()
2737 x0i = a[3] - a[18]; in cftf162()
2740 x0r = a[10] + a[27]; in cftf162()
2741 x0i = a[11] - a[26]; in cftf162()
2748 x0r = a[4] - a[21]; in cftf162()
2749 x0i = a[5] + a[20]; in cftf162()
2752 x0r = a[12] - a[29]; in cftf162()
2753 x0i = a[13] + a[28]; in cftf162()
2760 x0r = a[4] + a[21]; in cftf162()
2761 x0i = a[5] - a[20]; in cftf162()
2764 x0r = a[12] + a[29]; in cftf162()
2765 x0i = a[13] - a[28]; in cftf162()
2772 x0r = a[6] - a[23]; in cftf162()
2773 x0i = a[7] + a[22]; in cftf162()
2776 x0r = a[14] - a[31]; in cftf162()
2777 x0i = a[15] + a[30]; in cftf162()
2784 x0r = a[6] + a[23]; in cftf162()
2785 x0i = a[7] - a[22]; in cftf162()
2788 x0r = a[14] + a[31]; in cftf162()
2789 x0i = a[15] - a[30]; in cftf162()
2800 a[0] = x1r + x2r; in cftf162()
2801 a[1] = x1i + x2i; in cftf162()
2802 a[2] = x1r - x2r; in cftf162()
2803 a[3] = x1i - x2i; in cftf162()
2808 a[4] = x1r - x2i; in cftf162()
2809 a[5] = x1i + x2r; in cftf162()
2810 a[6] = x1r + x2i; in cftf162()
2811 a[7] = x1i - x2r; in cftf162()
2818 a[8] = x1r + x2r; in cftf162()
2819 a[9] = x1i + x2i; in cftf162()
2820 a[10] = x1r - x2r; in cftf162()
2821 a[11] = x1i - x2i; in cftf162()
2828 a[12] = x1r - x2i; in cftf162()
2829 a[13] = x1i + x2r; in cftf162()
2830 a[14] = x1r + x2i; in cftf162()
2831 a[15] = x1i - x2r; in cftf162()
2836 a[16] = x1r + x2r; in cftf162()
2837 a[17] = x1i + x2i; in cftf162()
2838 a[18] = x1r - x2r; in cftf162()
2839 a[19] = x1i - x2i; in cftf162()
2844 a[20] = x1r - x2i; in cftf162()
2845 a[21] = x1i + x2r; in cftf162()
2846 a[22] = x1r + x2i; in cftf162()
2847 a[23] = x1i - x2r; in cftf162()
2854 a[24] = x1r + x2r; in cftf162()
2855 a[25] = x1i + x2i; in cftf162()
2856 a[26] = x1r - x2r; in cftf162()
2857 a[27] = x1i - x2i; in cftf162()
2864 a[28] = x1r - x2i; in cftf162()
2865 a[29] = x1i + x2r; in cftf162()
2866 a[30] = x1r + x2i; in cftf162()
2867 a[31] = x1i - x2r; in cftf162()
2871 void cftf081(PICOFFTSG_FFTTYPE *a) in cftf081() argument
2878 x0r = a[0] + a[8]; in cftf081()
2879 x0i = a[1] + a[9]; in cftf081()
2880 x1r = a[0] - a[8]; in cftf081()
2881 x1i = a[1] - a[9]; in cftf081()
2882 x2r = a[4] + a[12]; in cftf081()
2883 x2i = a[5] + a[13]; in cftf081()
2884 x3r = a[4] - a[12]; in cftf081()
2885 x3i = a[5] - a[13]; in cftf081()
2894 x0r = a[2] + a[10]; in cftf081()
2895 x0i = a[3] + a[11]; in cftf081()
2896 x1r = a[2] - a[10]; in cftf081()
2897 x1i = a[3] - a[11]; in cftf081()
2898 x2r = a[6] + a[14]; in cftf081()
2899 x2i = a[7] + a[15]; in cftf081()
2900 x3r = a[6] - a[14]; in cftf081()
2901 x3i = a[7] - a[15]; in cftf081()
2914 a[8] = y1r + y5r; in cftf081()
2915 a[9] = y1i + y5i; in cftf081()
2916 a[10] = y1r - y5r; in cftf081()
2917 a[11] = y1i - y5i; in cftf081()
2918 a[12] = y3r - y7i; in cftf081()
2919 a[13] = y3i + y7r; in cftf081()
2920 a[14] = y3r + y7i; in cftf081()
2921 a[15] = y3i - y7r; in cftf081()
2922 a[0] = y0r + y4r; in cftf081()
2923 a[1] = y0i + y4i; in cftf081()
2924 a[2] = y0r - y4r; in cftf081()
2925 a[3] = y0i - y4i; in cftf081()
2926 a[4] = y2r - y6i; in cftf081()
2927 a[5] = y2i + y6r; in cftf081()
2928 a[6] = y2r + y6i; in cftf081()
2929 a[7] = y2i - y6r; in cftf081()
2933 void cftf082(PICOFFTSG_FFTTYPE *a) in cftf082() argument
2942 y0r = a[0] - a[9]; in cftf082()
2943 y0i = a[1] + a[8]; in cftf082()
2944 y1r = a[0] + a[9]; in cftf082()
2945 y1i = a[1] - a[8]; in cftf082()
2946 x0r = a[4] - a[13]; in cftf082()
2947 x0i = a[5] + a[12]; in cftf082()
2950 x0r = a[4] + a[13]; in cftf082()
2951 x0i = a[5] - a[12]; in cftf082()
2954 x0r = a[2] - a[11]; in cftf082()
2955 x0i = a[3] + a[10]; in cftf082()
2958 x0r = a[2] + a[11]; in cftf082()
2959 x0i = a[3] - a[10]; in cftf082()
2962 x0r = a[6] - a[15]; in cftf082()
2963 x0i = a[7] + a[14]; in cftf082()
2966 x0r = a[6] + a[15]; in cftf082()
2967 x0i = a[7] - a[14]; in cftf082()
2974 a[0] = x0r + x1r; in cftf082()
2975 a[1] = x0i + x1i; in cftf082()
2976 a[2] = x0r - x1r; in cftf082()
2977 a[3] = x0i - x1i; in cftf082()
2982 a[4] = x0r - x1i; in cftf082()
2983 a[5] = x0i + x1r; in cftf082()
2984 a[6] = x0r + x1i; in cftf082()
2985 a[7] = x0i - x1r; in cftf082()
2990 a[8] = x0r + x1r; in cftf082()
2991 a[9] = x0i + x1i; in cftf082()
2992 a[10] = x0r - x1r; in cftf082()
2993 a[11] = x0i - x1i; in cftf082()
2998 a[12] = x0r - x1i; in cftf082()
2999 a[13] = x0i + x1r; in cftf082()
3000 a[14] = x0r + x1i; in cftf082()
3001 a[15] = x0i - x1r; in cftf082()
3005 void cftf040(PICOFFTSG_FFTTYPE *a) in cftf040() argument
3009 x0r = a[0] + a[4]; in cftf040()
3010 x0i = a[1] + a[5]; in cftf040()
3011 x1r = a[0] - a[4]; in cftf040()
3012 x1i = a[1] - a[5]; in cftf040()
3013 x2r = a[2] + a[6]; in cftf040()
3014 x2i = a[3] + a[7]; in cftf040()
3015 x3r = a[2] - a[6]; in cftf040()
3016 x3i = a[3] - a[7]; in cftf040()
3017 a[0] = x0r + x2r; in cftf040()
3018 a[1] = x0i + x2i; in cftf040()
3019 a[2] = x1r - x3i; in cftf040()
3020 a[3] = x1i + x3r; in cftf040()
3021 a[4] = x0r - x2r; in cftf040()
3022 a[5] = x0i - x2i; in cftf040()
3023 a[6] = x1r + x3i; in cftf040()
3024 a[7] = x1i - x3r; in cftf040()
3028 void cftb040(PICOFFTSG_FFTTYPE *a) in cftb040() argument
3032 x0r = a[0] + a[4]; in cftb040()
3033 x0i = a[1] + a[5]; in cftb040()
3034 x1r = a[0] - a[4]; in cftb040()
3035 x1i = a[1] - a[5]; in cftb040()
3036 x2r = a[2] + a[6]; in cftb040()
3037 x2i = a[3] + a[7]; in cftb040()
3038 x3r = a[2] - a[6]; in cftb040()
3039 x3i = a[3] - a[7]; in cftb040()
3040 a[0] = x0r + x2r; in cftb040()
3041 a[1] = x0i + x2i; in cftb040()
3042 a[2] = x1r + x3i; in cftb040()
3043 a[3] = x1i - x3r; in cftb040()
3044 a[4] = x0r - x2r; in cftb040()
3045 a[5] = x0i - x2i; in cftb040()
3046 a[6] = x1r - x3i; in cftb040()
3047 a[7] = x1i + x3r; in cftb040()
3051 void cftx020(PICOFFTSG_FFTTYPE *a) in cftx020() argument
3055 x0r = a[0] - a[2]; in cftx020()
3056 x0i = a[1] - a[3]; in cftx020()
3057 a[0] += a[2]; in cftx020()
3058 a[1] += a[3]; in cftx020()
3059 a[2] = x0r; in cftx020()
3060 a[3] = x0i; in cftx020()
3064 void rftfsub(picoos_int32 n, PICOFFTSG_FFTTYPE *a) in rftfsub() argument
3102 xr = a[j + 2] - a[k - 2]; in rftfsub()
3103 xi = a[j + 3] + a[k - 1]; in rftfsub()
3106 a[j + 2] -= yr; in rftfsub()
3107 a[j + 3] -= yi; in rftfsub()
3108 a[k - 2] += yr; in rftfsub()
3109 a[k - 1] -= yi; in rftfsub()
3112 xr = a[j] - a[k]; in rftfsub()
3113 xi = a[j + 1] + a[k + 1]; in rftfsub()
3116 a[j] -= yr; in rftfsub()
3117 a[j + 1] -= yi; in rftfsub()
3118 a[k] += yr; in rftfsub()
3119 a[k + 1] -= yi; in rftfsub()
3128 xr = a[2] - a[n - 2]; in rftfsub()
3129 xi = a[3] + a[n - 1]; in rftfsub()
3133 a[2] -= yr; in rftfsub()
3134 a[3] -= yi; in rftfsub()
3135 a[n - 2] += yr; in rftfsub()
3136 a[n - 1] -= yi; in rftfsub()
3142 void rftbsub(picoos_int32 n, PICOFFTSG_FFTTYPE *a) in rftbsub() argument
3162 xr = a[j + 2] - a[k - 2]; in rftbsub()
3163 xi = a[j + 3] + a[k - 1]; in rftbsub()
3166 a[j + 2] -= yr; in rftbsub()
3167 a[j + 3] -= yi; in rftbsub()
3168 a[k - 2] += yr; in rftbsub()
3169 a[k - 1] -= yi; in rftbsub()
3172 xr = a[j] - a[k]; in rftbsub()
3173 xi = a[j + 1] + a[k + 1]; in rftbsub()
3176 a[j] -= yr; in rftbsub()
3177 a[j + 1] -= yi; in rftbsub()
3178 a[k] += yr; in rftbsub()
3179 a[k + 1] -= yi; in rftbsub()
3187 xr = a[2] - a[n - 2]; in rftbsub()
3188 xi = a[3] + a[n - 1]; in rftbsub()
3191 a[2] -= yr; in rftbsub()
3192 a[3] -= yi; in rftbsub()
3193 a[n - 2] += yr; in rftbsub()
3194 a[n - 1] -= yi; in rftbsub()
3198 void dctsub(picoos_int32 n, PICOFFTSG_FFTTYPE *a) in dctsub() argument
3231 xr = picofftsg_mult_w_a(wdi, a[j - 1]) - picofftsg_mult_w_a(wdr, a[k + 1]); in dctsub()
3232 xi = picofftsg_mult_w_a(wdr, a[j - 1]) + picofftsg_mult_w_a(wdi, a[k + 1]); in dctsub()
3235 yr = Mult_W_W(wki, a[j]) - Mult_W_W(wkr, a[k]); in dctsub()
3236 yi = Mult_W_W(wkr, a[j]) + Mult_W_W(wki, a[k]); in dctsub()
3239 a[k + 1] = xr; in dctsub()
3240 a[k] = yr; in dctsub()
3241 a[j - 1] = xi; in dctsub()
3242 a[j] = yi; in dctsub()
3248 xr = picofftsg_mult_w_a(wdi, a[m - 1]) - picofftsg_mult_w_a(wdr, a[m + 1]); in dctsub()
3249 a[m - 1] = picofftsg_mult_w_a(wdr, a[m - 1]) + picofftsg_mult_w_a(wdi, a[m + 1]); in dctsub()
3250 a[m + 1] = xr; in dctsub()
3251 a[m] = Mult_W_W(WR5000, a[m]); in dctsub()
3255 void dctsub4(picoos_int32 n, PICOFFTSG_FFTTYPE *a) in dctsub4() argument
3265 xr = Mult_W_W(wdi, a[1]) - Mult_W_W(wdr, a[3]); in dctsub4()
3266 a[1] = Mult_W_W(wdr, a[1]) + Mult_W_W(wdi, a[3]); in dctsub4()
3267 a[3] = xr; in dctsub4()
3269 a[m] = Mult_W_W(wki, a[m]); in dctsub4()