fn_64.c

FFTW, a collection of fast C routines to compute the Discrete Fourier Transform in one or more dime

	  fftw_real tmp863;
	  fftw_real tmp824;
	  fftw_real tmp874;
	  fftw_real tmp854;
	  fftw_real tmp864;
	  fftw_real tmp836;
	  fftw_real tmp858;
	  fftw_real tmp851;
	  fftw_real tmp873;
	  fftw_real tmp868;
	  fftw_real tmp878;
	  fftw_real tmp847;
	  fftw_real tmp859;
	  fftw_real tmp871;
	  fftw_real tmp879;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp820;
	       fftw_real tmp823;
	       fftw_real tmp866;
	       fftw_real tmp867;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp817 = tmp815 - tmp816;
	       tmp863 = tmp815 + tmp816;
	       tmp820 = tmp818 - tmp819;
	       tmp823 = tmp821 + tmp822;
	       tmp824 = K707106781 * (tmp820 - tmp823);
	       tmp874 = K707106781 * (tmp820 + tmp823);
	       {
		    fftw_real tmp852;
		    fftw_real tmp853;
		    fftw_real tmp830;
		    fftw_real tmp835;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp852 = tmp821 - tmp822;
		    tmp853 = tmp819 + tmp818;
		    tmp854 = K707106781 * (tmp852 - tmp853);
		    tmp864 = K707106781 * (tmp853 + tmp852);
		    tmp830 = tmp828 - tmp829;
		    tmp835 = tmp831 - tmp834;
		    tmp836 = (K923879532 * tmp830) + (K382683432 * tmp835);
		    tmp858 = (K382683432 * tmp830) - (K923879532 * tmp835);
	       }
	       tmp851 = tmp849 - tmp850;
	       tmp873 = tmp850 + tmp849;
	       tmp866 = tmp829 + tmp828;
	       tmp867 = tmp831 + tmp834;
	       tmp868 = (K382683432 * tmp866) + (K923879532 * tmp867);
	       tmp878 = (K923879532 * tmp866) - (K382683432 * tmp867);
	       {
		    fftw_real tmp841;
		    fftw_real tmp846;
		    fftw_real tmp869;
		    fftw_real tmp870;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp841 = tmp837 - tmp840;
		    tmp846 = tmp844 - tmp845;
		    tmp847 = (K382683432 * tmp841) - (K923879532 * tmp846);
		    tmp859 = (K382683432 * tmp846) + (K923879532 * tmp841);
		    tmp869 = tmp837 + tmp840;
		    tmp870 = tmp845 + tmp844;
		    tmp871 = (K923879532 * tmp869) - (K382683432 * tmp870);
		    tmp879 = (K923879532 * tmp870) + (K382683432 * tmp869);
	       }
	  }
	  {
	       fftw_real tmp825;
	       fftw_real tmp848;
	       fftw_real tmp855;
	       fftw_real tmp856;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp825 = tmp817 + tmp824;
	       tmp848 = tmp836 + tmp847;
	       c_re(output[44 * ostride]) = tmp825 - tmp848;
	       c_re(output[12 * ostride]) = tmp825 + tmp848;
	       tmp855 = tmp851 - tmp854;
	       tmp856 = tmp847 - tmp836;
	       c_im(output[60 * ostride]) = tmp855 - tmp856;
	       c_im(output[28 * ostride]) = tmp855 + tmp856;
	  }
	  {
	       fftw_real tmp861;
	       fftw_real tmp862;
	       fftw_real tmp857;
	       fftw_real tmp860;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp861 = tmp851 + tmp854;
	       tmp862 = tmp858 + tmp859;
	       c_im(output[44 * ostride]) = tmp861 - tmp862;
	       c_im(output[12 * ostride]) = tmp861 + tmp862;
	       tmp857 = tmp817 - tmp824;
	       tmp860 = tmp858 - tmp859;
	       c_re(output[60 * ostride]) = tmp857 - tmp860;
	       c_re(output[28 * ostride]) = tmp857 + tmp860;
	  }
	  {
	       fftw_real tmp865;
	       fftw_real tmp872;
	       fftw_real tmp875;
	       fftw_real tmp876;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp865 = tmp863 + tmp864;
	       tmp872 = tmp868 + tmp871;
	       c_re(output[36 * ostride]) = tmp865 - tmp872;
	       c_re(output[4 * ostride]) = tmp865 + tmp872;
	       tmp875 = tmp873 - tmp874;
	       tmp876 = tmp871 - tmp868;
	       c_im(output[52 * ostride]) = tmp875 - tmp876;
	       c_im(output[20 * ostride]) = tmp875 + tmp876;
	  }
	  {
	       fftw_real tmp881;
	       fftw_real tmp882;
	       fftw_real tmp877;
	       fftw_real tmp880;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp881 = tmp873 + tmp874;
	       tmp882 = tmp878 + tmp879;
	       c_im(output[36 * ostride]) = tmp881 - tmp882;
	       c_im(output[4 * ostride]) = tmp881 + tmp882;
	       tmp877 = tmp863 - tmp864;
	       tmp880 = tmp878 - tmp879;
	       c_re(output[52 * ostride]) = tmp877 - tmp880;
	       c_re(output[20 * ostride]) = tmp877 + tmp880;
	  }
     }
     {
	  fftw_real tmp757;
	  fftw_real tmp795;
	  fftw_real tmp800;
	  fftw_real tmp810;
	  fftw_real tmp803;
	  fftw_real tmp811;
	  fftw_real tmp779;
	  fftw_real tmp791;
	  fftw_real tmp783;
	  fftw_real tmp805;
	  fftw_real tmp764;
	  fftw_real tmp806;
	  fftw_real tmp786;
	  fftw_real tmp796;
	  fftw_real tmp772;
	  fftw_real tmp790;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp756;
	       fftw_real tmp798;
	       fftw_real tmp799;
	       fftw_real tmp782;
	       fftw_real tmp760;
	       fftw_real tmp763;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp756 = K707106781 * (tmp721 + tmp720);
	       tmp757 = tmp755 - tmp756;
	       tmp795 = tmp755 + tmp756;
	       tmp798 = tmp766 + tmp767;
	       tmp799 = tmp769 + tmp770;
	       tmp800 = (K980785280 * tmp798) + (K195090322 * tmp799);
	       tmp810 = (K980785280 * tmp799) - (K195090322 * tmp798);
	       {
		    fftw_real tmp801;
		    fftw_real tmp802;
		    fftw_real tmp775;
		    fftw_real tmp778;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp801 = tmp773 + tmp774;
		    tmp802 = tmp776 + tmp777;
		    tmp803 = (K980785280 * tmp801) - (K195090322 * tmp802);
		    tmp811 = (K195090322 * tmp801) + (K980785280 * tmp802);
		    tmp775 = tmp773 - tmp774;
		    tmp778 = tmp776 - tmp777;
		    tmp779 = (K555570233 * tmp775) - (K831469612 * tmp778);
		    tmp791 = (K831469612 * tmp775) + (K555570233 * tmp778);
	       }
	       tmp782 = K707106781 * (tmp640 + tmp643);
	       tmp783 = tmp781 - tmp782;
	       tmp805 = tmp781 + tmp782;
	       tmp760 = (K923879532 * tmp758) - (K382683432 * tmp759);
	       tmp763 = (K923879532 * tmp761) + (K382683432 * tmp762);
	       tmp764 = tmp760 - tmp763;
	       tmp806 = tmp760 + tmp763;
	       {
		    fftw_real tmp784;
		    fftw_real tmp785;
		    fftw_real tmp768;
		    fftw_real tmp771;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp784 = (K923879532 * tmp762) - (K382683432 * tmp761);
		    tmp785 = (K382683432 * tmp758) + (K923879532 * tmp759);
		    tmp786 = tmp784 - tmp785;
		    tmp796 = tmp785 + tmp784;
		    tmp768 = tmp766 - tmp767;
		    tmp771 = tmp769 - tmp770;
		    tmp772 = (K555570233 * tmp768) + (K831469612 * tmp771);
		    tmp790 = (K555570233 * tmp771) - (K831469612 * tmp768);
	       }
	  }
	  {
	       fftw_real tmp765;
	       fftw_real tmp780;
	       fftw_real tmp787;
	       fftw_real tmp788;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp765 = tmp757 + tmp764;
	       tmp780 = tmp772 + tmp779;
	       c_re(output[42 * ostride]) = tmp765 - tmp780;
	       c_re(output[10 * ostride]) = tmp765 + tmp780;
	       tmp787 = tmp783 - tmp786;
	       tmp788 = tmp779 - tmp772;
	       c_im(output[58 * ostride]) = tmp787 - tmp788;
	       c_im(output[26 * ostride]) = tmp787 + tmp788;
	  }
	  {
	       fftw_real tmp793;
	       fftw_real tmp794;
	       fftw_real tmp789;
	       fftw_real tmp792;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp793 = tmp783 + tmp786;
	       tmp794 = tmp790 + tmp791;
	       c_im(output[42 * ostride]) = tmp793 - tmp794;
	       c_im(output[10 * ostride]) = tmp793 + tmp794;
	       tmp789 = tmp757 - tmp764;
	       tmp792 = tmp790 - tmp791;
	       c_re(output[58 * ostride]) = tmp789 - tmp792;
	       c_re(output[26 * ostride]) = tmp789 + tmp792;
	  }
	  {
	       fftw_real tmp797;
	       fftw_real tmp804;
	       fftw_real tmp807;
	       fftw_real tmp808;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp797 = tmp795 + tmp796;
	       tmp804 = tmp800 + tmp803;
	       c_re(output[34 * ostride]) = tmp797 - tmp804;
	       c_re(output[2 * ostride]) = tmp797 + tmp804;
	       tmp807 = tmp805 - tmp806;
	       tmp808 = tmp803 - tmp800;
	       c_im(output[50 * ostride]) = tmp807 - tmp808;
	       c_im(output[18 * ostride]) = tmp807 + tmp808;
	  }
	  {
	       fftw_real tmp813;
	       fftw_real tmp814;
	       fftw_real tmp809;
	       fftw_real tmp812;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp813 = tmp805 + tmp806;
	       tmp814 = tmp810 + tmp811;
	       c_im(output[34 * ostride]) = tmp813 - tmp814;
	       c_im(output[2 * ostride]) = tmp813 + tmp814;
	       tmp809 = tmp795 - tmp796;
	       tmp812 = tmp810 - tmp811;
	       c_re(output[50 * ostride]) = tmp809 - tmp812;
	       c_re(output[18 * ostride]) = tmp809 + tmp812;
	  }
     }
     {
	  fftw_real tmp645;
	  fftw_real tmp735;
	  fftw_real tmp740;
	  fftw_real tmp750;
	  fftw_real tmp743;
	  fftw_real tmp751;
	  fftw_real tmp715;
	  fftw_real tmp731;
	  fftw_real tmp723;
	  fftw_real tmp745;
	  fftw_real tmp660;
	  fftw_real tmp746;
	  fftw_real tmp726;
	  fftw_real tmp736;
	  fftw_real tmp688;
	  fftw_real tmp730;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp644;
	       fftw_real tmp738;
	       fftw_real tmp739;
	       fftw_real tmp722;
	       fftw_real tmp652;
	       fftw_real tmp659;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp644 = K707106781 * (tmp640 - tmp643);
	       tmp645 = tmp637 - tmp644;
	       tmp735 = tmp637 + tmp644;
	       tmp738 = tmp666 + tmp677;
	       tmp739 = tmp683 + tmp686;
	       tmp740 = (K555570233 * tmp738) + (K831469612 * tmp739);
	       tmp750 = (K831469612 * tmp738) - (K555570233 * tmp739);
	       {
		    fftw_real tmp741;
		    fftw_real tmp742;
		    fftw_real tmp705;
		    fftw_real tmp714;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp741 = tmp693 + tmp704;
		    tmp742 = tmp710 + tmp713;
		    tmp743 = (K831469612 * tmp741) - (K555570233 * tmp742);
		    tmp751 = (K831469612 * tmp742) + (K555570233 * tmp741);
		    tmp705 = tmp693 - tmp704;
		    tmp714 = tmp710 - tmp713;
		    tmp715 = (K195090322 * tmp705) - (K980785280 * tmp714);
		    tmp731 = (K195090322 * tmp714) + (K980785280 * tmp705);
	       }
	       tmp722 = K707106781 * (tmp720 - tmp721);
	       tmp723 = tmp719 - tmp722;
	       tmp745 = tmp719 + tmp722;
	       tmp652 = (K382683432 * tmp648) - (K923879532 * tmp651);
	       tmp659 = (K382683432 * tmp655) + (K923879532 * tmp658);
	       tmp660 = tmp652 - tmp659;
	       tmp746 = tmp652 + tmp659;
	       {
		    fftw_real tmp724;
		    fftw_real tmp725;
		    fftw_real tmp678;
		    fftw_real tmp687;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp724 = (K382683432 * tmp658) - (K923879532 * tmp655);
		    tmp725 = (K923879532 * tmp648) + (K382683432 * tmp651);
		    tmp726 = tmp724 - tmp725;
		    tmp736 = tmp725 + tmp724;
		    tmp678 = tmp666 - tmp677;
		    tmp687 = tmp683 - tmp686;
		    tmp688 = (K980785280 * tmp678) + (K195090322 * tmp687);
		    tmp730 = (K195090322 * tmp678) - (K980785280 * tmp687);
	       }
	  }
	  {
	       fftw_real tmp661;
	       fftw_real tmp716;
	       fftw_real tmp727;
	       fftw_real tmp728;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp661 = tmp645 + tmp660;
	       tmp716 = tmp688 + tmp715;
	       c_re(output[46 * ostride]) = tmp661 - tmp716;
	       c_re(output[14 * ostride]) = tmp661 + tmp716;
	       tmp727 = tmp723 - tmp726;
	       tmp728 = tmp715 - tmp688;
	       c_im(output[62 * ostride]) = tmp727 - tmp728;
	       c_im(output[30 * ostride]) = tmp727 + tmp728;
	  }
	  {
	       fftw_real tmp733;
	       fftw_real tmp734;
	       fftw_real tmp729;
	       fftw_real tmp732;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp733 = tmp723 + tmp726;
	       tmp734 = tmp730 + tmp731;
	       c_im(output[46 * ostride]) = tmp733 - tmp734;
	       c_im(output[14 * ostride]) = tmp733 + tmp734;
	       tmp729 = tmp645 - tmp660;
	       tmp732 = tmp730 - tmp731;
	       c_re(output[62 * ostride]) = tmp729 - tmp732;
	       c_re(output[30 * ostride]) = tmp729 + tmp732;
	  }
	  {
	       fftw_real tmp737;
	       fftw_real tmp744;
	       fftw_real tmp747;
	       fftw_real tmp748;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp737 = tmp735 + tmp736;
	       tmp744 = tmp740 + tmp743;
	       c_re(output[38 * ostride]) = tmp737 - tmp744;
	       c_re(output[6 * ostride]) = tmp737 + tmp744;
	       tmp747 = tmp745 - tmp746;
	       tmp748 = tmp743 - tmp740;
	       c_im(output[54 * ostride]) = tmp747 - tmp748;
	       c_im(output[22 * ostride]) = tmp747 + tmp748;
	  }
	  {
	       fftw_real tmp753;
	       fftw_real tmp754;
	       fftw_real tmp749;
	       fftw_real tmp752;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp753 = tmp745 + tmp746;
	       tmp754 = tmp750 + tmp751;
	       c_im(output[38 * ostride]) = tmp753 - tmp754;
	       c_im(output[6 * ostride]) = tmp753 + tmp754;
	       tmp749 = tmp735 - tmp736;
	       tmp752 = tmp750 - tmp751;
	       c_re(output[54 * ostride]) = tmp749 - tmp752;
	       c_re(output[22 * ostride]) = tmp749 + tmp752;
	  }
     }
     {
	  fftw_real tmp481;
	  fftw_real tmp555;
	  fftw_real tmp560;
	  fftw_real tmp570;
	  fftw_real tmp563;
	  fftw_real tmp571;
	  fftw_real tmp535;
	  fftw_real tmp551;
	  fftw_real tmp516;
	  fftw_real tmp550;
	  fftw_real tmp543;
	  fftw_real tmp565;
	  fftw_real tmp546;
	  fftw_real tmp556;
	  fftw_real tmp496;
	  fftw_real tmp566;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp473;
	       fftw_real tmp480;
	       fftw_real tmp558;
	       fftw_real tmp559;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp473 = tmp471 - tmp472;
	       tmp480 = tmp476 - tmp479;
	       tmp481 = tmp473 - tmp480;
	       tmp555 = tmp473 + tmp480;
	       tmp558 = tmp500 + tmp507;
	       tmp559 = tmp511 + tmp514;
	       tmp560 = (K471396736 * tmp558) + (K881921264 * tmp559);
	       tmp570 = (K881921264 * tmp558) - (K471396736 * tmp559);
	  }
	  {
	       fftw_real tmp561;
	       fftw_real tmp562;
	       fftw_real tmp527;
	       fftw_real tmp534;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp561 = tmp519 + tmp526;
	       tmp562 = tmp530 + tmp533;
	       tmp563 = (K881921264 * tmp561) - (K471396736 * tmp562);
	       tmp571 = (K881921264 * tmp562) + (K471396736 * tmp561);
	       tmp527 = tmp519 - tmp526;
	       tmp534 = tmp530 - tmp533;
	       tmp535 = (K290284677 * tmp527) - (K956940335 * tmp534);
	       tmp551 = (K290284677 * tmp534) + (K956940335 * tmp527);
	  }
	  {
	       fftw_real tmp508;
	       fftw_real tmp515;
	       fftw_real tmp539;
	       fftw_real tmp542;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp508 = tmp500 - tmp507