ftw_64.c

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

		    ASSERT_ALIGNED_DOUBLE;
		    tmp978 = tmp843 + tmp844;
		    tmp979 = tmp838 + tmp839;
		    tmp980 = tmp978 - tmp979;
		    tmp1022 = tmp978 + tmp979;
		    tmp842 = tmp245 - tmp256;
		    tmp845 = tmp843 - tmp844;
		    tmp846 = tmp842 + tmp845;
		    tmp854 = tmp845 - tmp842;
	       }
	       {
		    fftw_real tmp497;
		    fftw_real tmp502;
		    fftw_real tmp687;
		    fftw_real tmp688;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp497 = tmp493 - tmp496;
		    tmp502 = tmp500 + tmp501;
		    tmp503 = (K382683432 * tmp497) - (K923879532 * tmp502);
		    tmp527 = (K382683432 * tmp502) + (K923879532 * tmp497);
		    tmp687 = tmp493 + tmp496;
		    tmp688 = tmp500 - tmp501;
		    tmp689 = (K923879532 * tmp687) - (K382683432 * tmp688);
		    tmp699 = (K923879532 * tmp688) + (K382683432 * tmp687);
	       }
	       {
		    fftw_real tmp508;
		    fftw_real tmp513;
		    fftw_real tmp690;
		    fftw_real tmp691;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp508 = tmp506 + tmp507;
		    tmp513 = tmp509 - tmp512;
		    tmp514 = (K923879532 * tmp508) + (K382683432 * tmp513);
		    tmp526 = (K382683432 * tmp508) - (K923879532 * tmp513);
		    tmp690 = tmp506 - tmp507;
		    tmp691 = tmp509 + tmp512;
		    tmp692 = (K382683432 * tmp690) + (K923879532 * tmp691);
		    tmp698 = (K923879532 * tmp690) - (K382683432 * tmp691);
	       }
	  }
	  {
	       fftw_real tmp91;
	       fftw_real tmp1015;
	       fftw_real tmp1038;
	       fftw_real tmp1039;
	       fftw_real tmp1049;
	       fftw_real tmp1055;
	       fftw_real tmp186;
	       fftw_real tmp1054;
	       fftw_real tmp1024;
	       fftw_real tmp1032;
	       fftw_real tmp378;
	       fftw_real tmp1051;
	       fftw_real tmp1029;
	       fftw_real tmp1033;
	       fftw_real tmp1018;
	       fftw_real tmp1040;
	       ASSERT_ALIGNED_DOUBLE;
	       {
		    fftw_real tmp43;
		    fftw_real tmp90;
		    fftw_real tmp1036;
		    fftw_real tmp1037;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp43 = tmp19 + tmp42;
		    tmp90 = tmp66 + tmp89;
		    tmp91 = tmp43 + tmp90;
		    tmp1015 = tmp43 - tmp90;
		    tmp1036 = tmp1021 + tmp1022;
		    tmp1037 = tmp1026 + tmp1027;
		    tmp1038 = tmp1036 - tmp1037;
		    tmp1039 = tmp1036 + tmp1037;
	       }
	       {
		    fftw_real tmp1041;
		    fftw_real tmp1048;
		    fftw_real tmp138;
		    fftw_real tmp185;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp1041 = tmp956 + tmp957;
		    tmp1048 = tmp1042 + tmp1047;
		    tmp1049 = tmp1041 + tmp1048;
		    tmp1055 = tmp1048 - tmp1041;
		    tmp138 = tmp114 + tmp137;
		    tmp185 = tmp161 + tmp184;
		    tmp186 = tmp138 + tmp185;
		    tmp1054 = tmp185 - tmp138;
	       }
	       {
		    fftw_real tmp1020;
		    fftw_real tmp1023;
		    fftw_real tmp282;
		    fftw_real tmp377;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp1020 = tmp234 - tmp281;
		    tmp1023 = tmp1021 - tmp1022;
		    tmp1024 = tmp1020 + tmp1023;
		    tmp1032 = tmp1023 - tmp1020;
		    tmp282 = tmp234 + tmp281;
		    tmp377 = tmp329 + tmp376;
		    tmp378 = tmp282 + tmp377;
		    tmp1051 = tmp377 - tmp282;
	       }
	       {
		    fftw_real tmp1025;
		    fftw_real tmp1028;
		    fftw_real tmp1016;
		    fftw_real tmp1017;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp1025 = tmp329 - tmp376;
		    tmp1028 = tmp1026 - tmp1027;
		    tmp1029 = tmp1025 - tmp1028;
		    tmp1033 = tmp1025 + tmp1028;
		    tmp1016 = tmp960 + tmp961;
		    tmp1017 = tmp966 + tmp967;
		    tmp1018 = tmp1016 - tmp1017;
		    tmp1040 = tmp1016 + tmp1017;
	       }
	       {
		    fftw_real tmp187;
		    fftw_real tmp1035;
		    fftw_real tmp1050;
		    fftw_real tmp1052;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp187 = tmp91 + tmp186;
		    c_re(inout[32 * iostride]) = tmp187 - tmp378;
		    c_re(inout[0]) = tmp187 + tmp378;
		    tmp1035 = tmp91 - tmp186;
		    c_re(inout[48 * iostride]) = tmp1035 - tmp1038;
		    c_re(inout[16 * iostride]) = tmp1035 + tmp1038;
		    {
			 fftw_real tmp1019;
			 fftw_real tmp1030;
			 fftw_real tmp1057;
			 fftw_real tmp1058;
			 ASSERT_ALIGNED_DOUBLE;
			 tmp1019 = tmp1015 + tmp1018;
			 tmp1030 = K707106781 * (tmp1024 + tmp1029);
			 c_re(inout[40 * iostride]) = tmp1019 - tmp1030;
			 c_re(inout[8 * iostride]) = tmp1019 + tmp1030;
			 tmp1057 = K707106781 * (tmp1029 - tmp1024);
			 tmp1058 = tmp1055 - tmp1054;
			 c_im(inout[24 * iostride]) = tmp1057 + tmp1058;
			 c_im(inout[56 * iostride]) = tmp1058 - tmp1057;
		    }
		    tmp1050 = tmp1040 + tmp1049;
		    c_im(inout[0]) = tmp1039 + tmp1050;
		    c_im(inout[32 * iostride]) = tmp1050 - tmp1039;
		    tmp1052 = tmp1049 - tmp1040;
		    c_im(inout[16 * iostride]) = tmp1051 + tmp1052;
		    c_im(inout[48 * iostride]) = tmp1052 - tmp1051;
		    {
			 fftw_real tmp1053;
			 fftw_real tmp1056;
			 fftw_real tmp1031;
			 fftw_real tmp1034;
			 ASSERT_ALIGNED_DOUBLE;
			 tmp1053 = K707106781 * (tmp1032 + tmp1033);
			 tmp1056 = tmp1054 + tmp1055;
			 c_im(inout[8 * iostride]) = tmp1053 + tmp1056;
			 c_im(inout[40 * iostride]) = tmp1056 - tmp1053;
			 tmp1031 = tmp1015 - tmp1018;
			 tmp1034 = K707106781 * (tmp1032 - tmp1033);
			 c_re(inout[56 * iostride]) = tmp1031 - tmp1034;
			 c_re(inout[24 * iostride]) = tmp1031 + tmp1034;
		    }
	       }
	  }
	  {
	       fftw_real tmp959;
	       fftw_real tmp999;
	       fftw_real tmp1002;
	       fftw_real tmp1068;
	       fftw_real tmp970;
	       fftw_real tmp1060;
	       fftw_real tmp1063;
	       fftw_real tmp1069;
	       fftw_real tmp982;
	       fftw_real tmp996;
	       fftw_real tmp1006;
	       fftw_real tmp1012;
	       fftw_real tmp993;
	       fftw_real tmp997;
	       fftw_real tmp1009;
	       fftw_real tmp1013;
	       ASSERT_ALIGNED_DOUBLE;
	       {
		    fftw_real tmp955;
		    fftw_real tmp958;
		    fftw_real tmp1000;
		    fftw_real tmp1001;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp955 = tmp19 - tmp42;
		    tmp958 = tmp956 - tmp957;
		    tmp959 = tmp955 - tmp958;
		    tmp999 = tmp955 + tmp958;
		    tmp1000 = tmp963 + tmp962;
		    tmp1001 = tmp965 - tmp968;
		    tmp1002 = K707106781 * (tmp1000 + tmp1001);
		    tmp1068 = K707106781 * (tmp1001 - tmp1000);
	       }
	       {
		    fftw_real tmp964;
		    fftw_real tmp969;
		    fftw_real tmp1061;
		    fftw_real tmp1062;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp964 = tmp962 - tmp963;
		    tmp969 = tmp965 + tmp968;
		    tmp970 = K707106781 * (tmp964 - tmp969);
		    tmp1060 = K707106781 * (tmp964 + tmp969);
		    tmp1061 = tmp89 - tmp66;
		    tmp1062 = tmp1047 - tmp1042;
		    tmp1063 = tmp1061 + tmp1062;
		    tmp1069 = tmp1062 - tmp1061;
	       }
	       {
		    fftw_real tmp976;
		    fftw_real tmp981;
		    fftw_real tmp1004;
		    fftw_real tmp1005;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp976 = tmp974 - tmp975;
		    tmp981 = tmp977 - tmp980;
		    tmp982 = (K923879532 * tmp976) + (K382683432 * tmp981);
		    tmp996 = (K382683432 * tmp976) - (K923879532 * tmp981);
		    tmp1004 = tmp974 + tmp975;
		    tmp1005 = tmp977 + tmp980;
		    tmp1006 =
			(K382683432 * tmp1004) + (K923879532 * tmp1005);
		    tmp1012 =
			(K923879532 * tmp1004) - (K382683432 * tmp1005);
	       }
	       {
		    fftw_real tmp987;
		    fftw_real tmp992;
		    fftw_real tmp1007;
		    fftw_real tmp1008;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp987 = tmp983 - tmp986;
		    tmp992 = tmp990 - tmp991;
		    tmp993 = (K382683432 * tmp987) - (K923879532 * tmp992);
		    tmp997 = (K382683432 * tmp992) + (K923879532 * tmp987);
		    tmp1007 = tmp983 + tmp986;
		    tmp1008 = tmp990 + tmp991;
		    tmp1009 =
			(K923879532 * tmp1007) - (K382683432 * tmp1008);
		    tmp1013 =
			(K923879532 * tmp1008) + (K382683432 * tmp1007);
	       }
	       {
		    fftw_real tmp971;
		    fftw_real tmp994;
		    fftw_real tmp995;
		    fftw_real tmp998;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp971 = tmp959 + tmp970;
		    tmp994 = tmp982 + tmp993;
		    c_re(inout[44 * iostride]) = tmp971 - tmp994;
		    c_re(inout[12 * iostride]) = tmp971 + tmp994;
		    tmp995 = tmp959 - tmp970;
		    tmp998 = tmp996 - tmp997;
		    c_re(inout[60 * iostride]) = tmp995 - tmp998;
		    c_re(inout[28 * iostride]) = tmp995 + tmp998;
	       }
	       {
		    fftw_real tmp1067;
		    fftw_real tmp1070;
		    fftw_real tmp1071;
		    fftw_real tmp1072;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp1067 = tmp996 + tmp997;
		    tmp1070 = tmp1068 + tmp1069;
		    c_im(inout[12 * iostride]) = tmp1067 + tmp1070;
		    c_im(inout[44 * iostride]) = tmp1070 - tmp1067;
		    tmp1071 = tmp993 - tmp982;
		    tmp1072 = tmp1069 - tmp1068;
		    c_im(inout[28 * iostride]) = tmp1071 + tmp1072;
		    c_im(inout[60 * iostride]) = tmp1072 - tmp1071;
	       }
	       {
		    fftw_real tmp1003;
		    fftw_real tmp1010;
		    fftw_real tmp1011;
		    fftw_real tmp1014;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp1003 = tmp999 + tmp1002;
		    tmp1010 = tmp1006 + tmp1009;
		    c_re(inout[36 * iostride]) = tmp1003 - tmp1010;
		    c_re(inout[4 * iostride]) = tmp1003 + tmp1010;
		    tmp1011 = tmp999 - tmp1002;
		    tmp1014 = tmp1012 - tmp1013;
		    c_re(inout[52 * iostride]) = tmp1011 - tmp1014;
		    c_re(inout[20 * iostride]) = tmp1011 + tmp1014;
	       }
	       {
		    fftw_real tmp1059;
		    fftw_real tmp1064;
		    fftw_real tmp1065;
		    fftw_real tmp1066;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp1059 = tmp1012 + tmp1013;
		    tmp1064 = tmp1060 + tmp1063;
		    c_im(inout[4 * iostride]) = tmp1059 + tmp1064;
		    c_im(inout[36 * iostride]) = tmp1064 - tmp1059;
		    tmp1065 = tmp1009 - tmp1006;
		    tmp1066 = tmp1063 - tmp1060;
		    c_im(inout[20 * iostride]) = tmp1065 + tmp1066;
		    c_im(inout[52 * iostride]) = tmp1066 - tmp1065;
	       }
	  }
	  {
	       fftw_real tmp419;
	       fftw_real tmp591;
	       fftw_real tmp1155;
	       fftw_real tmp1161;
	       fftw_real tmp474;
	       fftw_real tmp1152;
	       fftw_real tmp594;
	       fftw_real tmp1160;
	       fftw_real tmp530;
	       fftw_real tmp588;
	       fftw_real tmp598;
	       fftw_real tmp604;
	       fftw_real tmp585;
	       fftw_real tmp589;
	       fftw_real tmp601;
	       fftw_real tmp605;
	       ASSERT_ALIGNED_DOUBLE;
	       {
		    fftw_real tmp395;
		    fftw_real tmp418;
		    fftw_real tmp1153;
		    fftw_real tmp1154;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp395 = tmp383 - tmp394;
		    tmp418 = tmp406 - tmp417;
		    tmp419 = tmp395 - tmp418;
		    tmp591 = tmp395 + tmp418;
		    tmp1153 = tmp609 - tmp608;
		    tmp1154 = tmp1139 - tmp1138;
		    tmp1155 = tmp1153 + tmp1154;
		    tmp1161 = tmp1154 - tmp1153;
	       }
	       {
		    fftw_real tmp446;
		    fftw_real tmp592;
		    fftw_real tmp473;
		    fftw_real tmp593;
		    ASSERT_ALIGNED_DOUBLE;
		    {
			 fftw_real tmp436;
			 fftw_real tmp445;
			 fftw_real tmp463;
			 fftw_real tmp472;
			 ASSERT_ALIGNED_DOUBLE;
			 tmp436 = tmp424 - tmp435;
			 tmp445 = tmp441 - tmp444;
			 tmp446 =
			     (K195090322 * tmp436) - (K980785280 * tmp445);
			 tmp592 =
			     (K980785280 * tmp436) + (K195090322 * tmp445);
			 tmp463 = tmp451 - tmp462;
			 tmp472 = tmp468 - tmp471;
			 tmp473 =
			     (K195090322 * tmp463) + (K980785280 * tmp472);
			 tmp593 =
			     (K195090322 * tmp472) - (K980785280 * tmp463);
		    }
		    tmp474 = tmp446 - tmp473;
		    tmp1152 = tmp446 + tmp473;
		    tmp594 = tmp592 + tmp593;
		    tmp1160 = tmp593 - tmp592;
	       }
	       {
		    fftw_real tmp516;
		    fftw_real tmp596;
		    fftw_real tmp529;
		    fftw_real tmp597;
		    ASSERT_ALIGNED_DOUBLE;
		    {
			 fftw_real tmp492;
			 fftw_real tmp515;
			 fftw_real tmp525;
			 fftw_real tmp528;
			 ASSERT_ALIGNED_DOUBLE;
			 tmp492 = tmp480 - tmp491;
			 tmp515 = tmp503 - tmp514;
			 tmp516 = tmp492 - tmp515;
			 tmp596 = tmp492 + tmp515;
			 tmp525 = tmp521 - tmp524;
			 tmp528 = tmp526 - tmp527;
			 tmp529 = tmp525 - tmp528;
			 tmp597 = tmp525 + tmp528;
		    }
		    tmp530 = (K995184726 * tmp516) + (K098017140 * tmp529);
		    tmp588 = (K098017140 * tmp516) - (K995184726 * tmp529);
		    tmp598 = (K634393284 * tmp596) + (K773010453 * tmp597);
		    tmp604 = (K773010453 * tmp596) - (K634393284 * tmp597);
	       }
	       {
		    fftw_real tmp571;
		    fftw_real tmp599;
		    fftw_real tmp584;
		    fftw_real tmp600;
		    ASSERT_ALIGNED_DOUBLE;
		    {
			 fftw_real tmp547;
			 fftw_real tmp570;
			 fftw_real tmp580;
			 fftw_real tmp583;
			 ASSERT_ALIGNED_DOUBLE;
			 tmp547 = tmp535 - tmp546;
			 tmp570 = tmp558 - tmp569;
			 tmp571 = tmp547 - tmp570;
			 tmp599 = tmp547 + tmp570;
			 tmp580 = tmp576 - tmp579;
			 tmp583 = tmp581 - tmp582;
			 tmp584 = tmp580 - tmp583;
			 tmp600 = tmp580 + tmp583;
		    }
		    tmp585 = (K098017140 * tmp571) - (K995184726 * tmp584);
		    tmp589 = (K098017140 * tmp584) + (K995184726 * tmp571);
		    tmp601 = (K773010453 * tmp599) - (K634393284 * tmp600);
		    tmp605 = (K773010453 * tmp600) + (K634393284 * tmp599);
	       }
	       {
		    fftw_real tmp475;
		    fftw_real tmp586;
		    fftw_real tmp587;
		    fftw_real tmp590;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp475 = tmp419 + tmp474;
		    tmp586 = tmp530 + tmp585;
		    c_re(inout[47 * iostride]) = tmp475 - tmp586;
		    c_re(inout[15 * iostride]) = tmp475 + tmp586;
		    tmp587 = tmp419 - tmp474;
		    tmp590 = tmp588 - tmp589;
		    c_re(inout[63 * iostride]) = tmp587 - tmp590;
		    c_re(inout[31 * iostride]) = tmp587 + tmp590;
	       }
	       {
		    fftw_real tmp1159;
		    fftw_real tmp1162;
		    fftw_real tmp1163;
		    fftw_real tmp1164;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp1159 = tmp588 + tmp589;
		    tmp1162 = tmp1160 + tmp1161;
		    c_im(inout[15 * iostride]) = tmp1159 + tmp1162;
		    c_im(inout[47 * iostride]) = tmp1162 - tmp1159;
		    tmp1163 = tmp585 - tmp530;
		    tmp1164 = tmp1161 - tmp1160;
		    c_im(inout[31 * iostride]) = tmp1163 + tmp1164;
		    c_im(inout[63 * iostride]) = tmp1164 - tmp1163;
	       }
	       {
		    fftw_real tmp595;
		    fftw_real tmp602;
		    fftw_real tmp603;
		    fftw_real tmp6