fni_64.c

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

	  fftw_real tmp893;
	  fftw_real tmp862;
	  fftw_real tmp904;
	  fftw_real tmp884;
	  fftw_real tmp894;
	  fftw_real tmp870;
	  fftw_real tmp888;
	  fftw_real tmp881;
	  fftw_real tmp903;
	  fftw_real tmp898;
	  fftw_real tmp908;
	  fftw_real tmp877;
	  fftw_real tmp889;
	  fftw_real tmp901;
	  fftw_real tmp909;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp858;
	       fftw_real tmp861;
	       fftw_real tmp896;
	       fftw_real tmp897;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp855 = tmp853 - tmp854;
	       tmp893 = tmp854 + tmp853;
	       tmp858 = tmp856 - tmp857;
	       tmp861 = tmp859 + tmp860;
	       tmp862 = K707106781 * (tmp858 - tmp861);
	       tmp904 = K707106781 * (tmp858 + tmp861);
	       {
		    fftw_real tmp882;
		    fftw_real tmp883;
		    fftw_real tmp866;
		    fftw_real tmp869;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp882 = tmp860 - tmp859;
		    tmp883 = tmp856 + tmp857;
		    tmp884 = K707106781 * (tmp882 - tmp883);
		    tmp894 = K707106781 * (tmp883 + tmp882);
		    tmp866 = tmp864 - tmp865;
		    tmp869 = tmp867 - tmp868;
		    tmp870 = (K382683432 * tmp866) + (K923879532 * tmp869);
		    tmp888 = (K382683432 * tmp869) - (K923879532 * tmp866);
	       }
	       tmp881 = tmp879 - tmp880;
	       tmp903 = tmp879 + tmp880;
	       tmp896 = tmp865 + tmp864;
	       tmp897 = tmp867 + tmp868;
	       tmp898 = (K923879532 * tmp896) + (K382683432 * tmp897);
	       tmp908 = (K923879532 * tmp897) - (K382683432 * tmp896);
	       {
		    fftw_real tmp873;
		    fftw_real tmp876;
		    fftw_real tmp899;
		    fftw_real tmp900;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp873 = tmp871 - tmp872;
		    tmp876 = tmp874 - tmp875;
		    tmp877 = (K382683432 * tmp873) - (K923879532 * tmp876);
		    tmp889 = (K923879532 * tmp873) + (K382683432 * tmp876);
		    tmp899 = tmp872 + tmp871;
		    tmp900 = tmp874 + tmp875;
		    tmp901 = (K923879532 * tmp899) - (K382683432 * tmp900);
		    tmp909 = (K382683432 * tmp899) + (K923879532 * tmp900);
	       }
	  }
	  {
	       fftw_real tmp863;
	       fftw_real tmp878;
	       fftw_real tmp885;
	       fftw_real tmp886;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp863 = tmp855 + tmp862;
	       tmp878 = tmp870 + tmp877;
	       c_im(output[44 * ostride]) = tmp863 - tmp878;
	       c_im(output[12 * ostride]) = tmp863 + tmp878;
	       tmp885 = tmp881 - tmp884;
	       tmp886 = tmp877 - tmp870;
	       c_re(output[60 * ostride]) = tmp885 - tmp886;
	       c_re(output[28 * ostride]) = tmp885 + tmp886;
	  }
	  {
	       fftw_real tmp891;
	       fftw_real tmp892;
	       fftw_real tmp887;
	       fftw_real tmp890;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp891 = tmp881 + tmp884;
	       tmp892 = tmp888 + tmp889;
	       c_re(output[44 * ostride]) = tmp891 - tmp892;
	       c_re(output[12 * ostride]) = tmp891 + tmp892;
	       tmp887 = tmp855 - tmp862;
	       tmp890 = tmp888 - tmp889;
	       c_im(output[60 * ostride]) = tmp887 - tmp890;
	       c_im(output[28 * ostride]) = tmp887 + tmp890;
	  }
	  {
	       fftw_real tmp895;
	       fftw_real tmp902;
	       fftw_real tmp905;
	       fftw_real tmp906;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp895 = tmp893 + tmp894;
	       tmp902 = tmp898 + tmp901;
	       c_im(output[36 * ostride]) = tmp895 - tmp902;
	       c_im(output[4 * ostride]) = tmp895 + tmp902;
	       tmp905 = tmp903 - tmp904;
	       tmp906 = tmp901 - tmp898;
	       c_re(output[52 * ostride]) = tmp905 - tmp906;
	       c_re(output[20 * ostride]) = tmp905 + tmp906;
	  }
	  {
	       fftw_real tmp911;
	       fftw_real tmp912;
	       fftw_real tmp907;
	       fftw_real tmp910;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp911 = tmp903 + tmp904;
	       tmp912 = tmp908 + tmp909;
	       c_re(output[36 * ostride]) = tmp911 - tmp912;
	       c_re(output[4 * ostride]) = tmp911 + tmp912;
	       tmp907 = tmp893 - tmp894;
	       tmp910 = tmp908 - tmp909;
	       c_im(output[52 * ostride]) = tmp907 - tmp910;
	       c_im(output[20 * ostride]) = tmp907 + tmp910;
	  }
     }
     {
	  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 * (tmp720 - tmp721);
	       tmp757 = tmp755 - tmp756;
	       tmp795 = tmp755 + tmp756;
	       tmp798 = tmp766 + tmp767;
	       tmp799 = tmp769 + tmp770;
	       tmp800 = (K831469612 * tmp798) + (K555570233 * tmp799);
	       tmp810 = (K831469612 * tmp799) - (K555570233 * tmp798);
	       {
		    fftw_real tmp801;
		    fftw_real tmp802;
		    fftw_real tmp775;
		    fftw_real tmp778;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp801 = tmp773 + tmp774;
		    tmp802 = tmp776 + tmp777;
		    tmp803 = (K831469612 * tmp801) - (K555570233 * tmp802);
		    tmp811 = (K555570233 * tmp801) + (K831469612 * tmp802);
		    tmp775 = tmp773 - tmp774;
		    tmp778 = tmp776 - tmp777;
		    tmp779 = (K195090322 * tmp775) - (K980785280 * tmp778);
		    tmp791 = (K980785280 * tmp775) + (K195090322 * tmp778);
	       }
	       tmp782 = K707106781 * (tmp643 - tmp640);
	       tmp783 = tmp781 - tmp782;
	       tmp805 = tmp781 + tmp782;
	       tmp760 = (K382683432 * tmp758) - (K923879532 * 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 = (K382683432 * tmp761) - (K923879532 * tmp762);
		    tmp785 = (K382683432 * tmp759) + (K923879532 * tmp758);
		    tmp786 = tmp784 - tmp785;
		    tmp796 = tmp785 + tmp784;
		    tmp768 = tmp766 - tmp767;
		    tmp771 = tmp769 - tmp770;
		    tmp772 = (K195090322 * tmp768) + (K980785280 * tmp771);
		    tmp790 = (K195090322 * tmp771) - (K980785280 * tmp768);
	       }
	  }
	  {
	       fftw_real tmp765;
	       fftw_real tmp780;
	       fftw_real tmp787;
	       fftw_real tmp788;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp765 = tmp757 + tmp764;
	       tmp780 = tmp772 + tmp779;
	       c_im(output[46 * ostride]) = tmp765 - tmp780;
	       c_im(output[14 * ostride]) = tmp765 + tmp780;
	       tmp787 = tmp783 - tmp786;
	       tmp788 = tmp779 - tmp772;
	       c_re(output[62 * ostride]) = tmp787 - tmp788;
	       c_re(output[30 * 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_re(output[46 * ostride]) = tmp793 - tmp794;
	       c_re(output[14 * ostride]) = tmp793 + tmp794;
	       tmp789 = tmp757 - tmp764;
	       tmp792 = tmp790 - tmp791;
	       c_im(output[62 * ostride]) = tmp789 - tmp792;
	       c_im(output[30 * 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_im(output[38 * ostride]) = tmp797 - tmp804;
	       c_im(output[6 * ostride]) = tmp797 + tmp804;
	       tmp807 = tmp805 - tmp806;
	       tmp808 = tmp803 - tmp800;
	       c_re(output[54 * ostride]) = tmp807 - tmp808;
	       c_re(output[22 * 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_re(output[38 * ostride]) = tmp813 - tmp814;
	       c_re(output[6 * ostride]) = tmp813 + tmp814;
	       tmp809 = tmp795 - tmp796;
	       tmp812 = tmp810 - tmp811;
	       c_im(output[54 * ostride]) = tmp809 - tmp812;
	       c_im(output[22 * 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 = (K980785280 * tmp738) - (K195090322 * tmp739);
	       tmp750 = (K195090322 * tmp738) + (K980785280 * tmp739);
	       {
		    fftw_real tmp741;
		    fftw_real tmp742;
		    fftw_real tmp705;
		    fftw_real tmp714;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp741 = tmp693 + tmp704;
		    tmp742 = tmp710 + tmp713;
		    tmp743 = (K980785280 * tmp741) + (K195090322 * tmp742);
		    tmp751 = (K980785280 * tmp742) - (K195090322 * tmp741);
		    tmp705 = tmp693 - tmp704;
		    tmp714 = tmp710 - tmp713;
		    tmp715 = (K555570233 * tmp705) + (K831469612 * tmp714);
		    tmp731 = (K555570233 * tmp714) - (K831469612 * tmp705);
	       }
	       tmp722 = K707106781 * (tmp720 + tmp721);
	       tmp723 = tmp719 - tmp722;
	       tmp745 = tmp719 + tmp722;
	       tmp652 = (K923879532 * tmp648) - (K382683432 * 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 = (K923879532 * tmp655) - (K382683432 * tmp658);
		    tmp725 = (K923879532 * tmp651) + (K382683432 * tmp648);
		    tmp726 = tmp724 - tmp725;
		    tmp736 = tmp725 + tmp724;
		    tmp678 = tmp666 - tmp677;
		    tmp687 = tmp683 - tmp686;
		    tmp688 = (K555570233 * tmp678) - (K831469612 * tmp687);
		    tmp730 = (K831469612 * tmp678) + (K555570233 * tmp687);
	       }
	  }
	  {
	       fftw_real tmp661;
	       fftw_real tmp716;
	       fftw_real tmp727;
	       fftw_real tmp728;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp661 = tmp645 - tmp660;
	       tmp716 = tmp688 - tmp715;
	       c_im(output[58 * ostride]) = tmp661 - tmp716;
	       c_im(output[26 * ostride]) = tmp661 + tmp716;
	       tmp727 = tmp723 + tmp726;
	       tmp728 = tmp688 + tmp715;
	       c_re(output[42 * ostride]) = tmp727 - tmp728;
	       c_re(output[10 * ostride]) = tmp727 + tmp728;
	  }
	  {
	       fftw_real tmp733;
	       fftw_real tmp734;
	       fftw_real tmp729;
	       fftw_real tmp732;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp733 = tmp723 - tmp726;
	       tmp734 = tmp731 - tmp730;
	       c_re(output[58 * ostride]) = tmp733 - tmp734;
	       c_re(output[26 * ostride]) = tmp733 + tmp734;
	       tmp729 = tmp645 + tmp660;
	       tmp732 = tmp730 + tmp731;
	       c_im(output[42 * ostride]) = tmp729 - tmp732;
	       c_im(output[10 * 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_im(output[50 * ostride]) = tmp737 - tmp744;
	       c_im(output[18 * ostride]) = tmp737 + tmp744;
	       tmp747 = tmp745 + tmp746;
	       tmp748 = tmp740 + tmp743;
	       c_re(output[34 * ostride]) = tmp747 - tmp748;
	       c_re(output[2 * ostride]) = tmp747 + tmp748;
	  }
	  {
	       fftw_real tmp753;
	       fftw_real tmp754;
	       fftw_real tmp749;
	       fftw_real tmp752;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp753 = tmp745 - tmp746;
	       tmp754 = tmp751 - tmp750;
	       c_re(output[50 * ostride]) = tmp753 - tmp754;
	       c_re(output[18 * ostride]) = tmp753 + tmp754;
	       tmp749 = tmp735 + tmp736;
	       tmp752 = tmp750 + tmp751;
	       c_im(output[34 * ostride]) = tmp749 - tmp752;
	       c_im(output[2 * 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 = (K773010453 * tmp558) - (K634393284 * tmp559);
	       tmp570 = (K634393284 * tmp558) + (K773010453 * tmp559);
	  }
	  {
	       fftw_real tmp561;
	       fftw_real tmp562;
	       fftw_real tmp527;
	       fftw_real tmp534;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp561 = tmp519 + tmp526;
	       tmp562 = tmp530 + tmp533;
	       tmp563 = (K773010453 * tmp561) + (K634393284 * tmp562);
	       tmp571 = (K773010453 * tmp562) - (K634393284 * tmp561);
	       tmp527 = tmp519 - tmp526;
	       tmp534 = tmp530 - tmp533;
	       tmp535 = (K098017140 * tmp527) + (K995184726 * tmp534);
	       tmp551 = (K098017140 * tmp534) - (K995184726 * tmp527);
	  }
	  {
	       fftw_real tmp508;
	       fftw_real tmp515;
	       fftw_real tmp539;
	       fftw_real tmp542;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp508 = tmp500 - tmp507;