ftwi_64.c

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

		    tmp1015 = tmp969 - tmp968;
		    tmp842 = tmp268 - tmp279;
		    tmp845 = tmp843 - tmp844;
		    tmp846 = tmp842 + tmp845;
		    tmp855 = tmp845 - tmp842;
	       }
	       {
		    fftw_real tmp497;
		    fftw_real tmp502;
		    fftw_real tmp687;
		    fftw_real tmp688;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp497 = tmp495 + tmp496;
		    tmp502 = tmp498 - tmp501;
		    tmp503 = (K923879532 * tmp497) - (K382683432 * tmp502);
		    tmp527 = (K382683432 * tmp497) + (K923879532 * tmp502);
		    tmp687 = tmp495 - tmp496;
		    tmp688 = tmp498 + tmp501;
		    tmp689 = (K382683432 * tmp687) - (K923879532 * tmp688);
		    tmp699 = (K923879532 * tmp687) + (K382683432 * tmp688);
	       }
	       {
		    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 = (K923879532 * tmp513) - (K382683432 * tmp508);
		    tmp690 = tmp506 - tmp507;
		    tmp691 = tmp509 + tmp512;
		    tmp692 = (K382683432 * tmp690) + (K923879532 * tmp691);
		    tmp698 = (K382683432 * tmp691) - (K923879532 * tmp690);
	       }
	  }
	  {
	       fftw_real tmp91;
	       fftw_real tmp955;
	       fftw_real tmp990;
	       fftw_real tmp1039;
	       fftw_real tmp1049;
	       fftw_real tmp1055;
	       fftw_real tmp186;
	       fftw_real tmp1054;
	       fftw_real tmp972;
	       fftw_real tmp985;
	       fftw_real tmp378;
	       fftw_real tmp1051;
	       fftw_real tmp981;
	       fftw_real tmp984;
	       fftw_real tmp962;
	       fftw_real tmp1040;
	       ASSERT_ALIGNED_DOUBLE;
	       {
		    fftw_real tmp43;
		    fftw_real tmp90;
		    fftw_real tmp988;
		    fftw_real tmp989;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp43 = tmp19 + tmp42;
		    tmp90 = tmp66 + tmp89;
		    tmp91 = tmp43 + tmp90;
		    tmp955 = tmp43 - tmp90;
		    tmp988 = tmp976 + tmp979;
		    tmp989 = tmp967 + tmp970;
		    tmp990 = tmp988 - tmp989;
		    tmp1039 = tmp989 + tmp988;
	       }
	       {
		    fftw_real tmp1041;
		    fftw_real tmp1048;
		    fftw_real tmp138;
		    fftw_real tmp185;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp1041 = tmp993 + tmp992;
		    tmp1048 = tmp1042 + tmp1047;
		    tmp1049 = tmp1041 + tmp1048;
		    tmp1055 = tmp1048 - tmp1041;
		    tmp138 = tmp114 + tmp137;
		    tmp185 = tmp161 + tmp184;
		    tmp186 = tmp138 + tmp185;
		    tmp1054 = tmp138 - tmp185;
	       }
	       {
		    fftw_real tmp964;
		    fftw_real tmp971;
		    fftw_real tmp282;
		    fftw_real tmp377;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp964 = tmp234 - tmp281;
		    tmp971 = tmp967 - tmp970;
		    tmp972 = tmp964 - tmp971;
		    tmp985 = tmp964 + tmp971;
		    tmp282 = tmp234 + tmp281;
		    tmp377 = tmp329 + tmp376;
		    tmp378 = tmp282 + tmp377;
		    tmp1051 = tmp282 - tmp377;
	       }
	       {
		    fftw_real tmp973;
		    fftw_real tmp980;
		    fftw_real tmp958;
		    fftw_real tmp961;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp973 = tmp329 - tmp376;
		    tmp980 = tmp976 - tmp979;
		    tmp981 = tmp973 + tmp980;
		    tmp984 = tmp980 - tmp973;
		    tmp958 = tmp956 + tmp957;
		    tmp961 = tmp959 + tmp960;
		    tmp962 = tmp958 - tmp961;
		    tmp1040 = tmp961 + tmp958;
	       }
	       {
		    fftw_real tmp187;
		    fftw_real tmp987;
		    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;
		    tmp987 = tmp91 - tmp186;
		    c_re(inout[48 * iostride]) = tmp987 - tmp990;
		    c_re(inout[16 * iostride]) = tmp987 + tmp990;
		    {
			 fftw_real tmp963;
			 fftw_real tmp982;
			 fftw_real tmp1057;
			 fftw_real tmp1058;
			 ASSERT_ALIGNED_DOUBLE;
			 tmp963 = tmp955 + tmp962;
			 tmp982 = K707106781 * (tmp972 + tmp981);
			 c_re(inout[40 * iostride]) = tmp963 - tmp982;
			 c_re(inout[8 * iostride]) = tmp963 + tmp982;
			 tmp1057 = K707106781 * (tmp972 - tmp981);
			 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 tmp983;
			 fftw_real tmp986;
			 ASSERT_ALIGNED_DOUBLE;
			 tmp1053 = K707106781 * (tmp985 + tmp984);
			 tmp1056 = tmp1054 + tmp1055;
			 c_im(inout[8 * iostride]) = tmp1053 + tmp1056;
			 c_im(inout[40 * iostride]) = tmp1056 - tmp1053;
			 tmp983 = tmp955 - tmp962;
			 tmp986 = K707106781 * (tmp984 - tmp985);
			 c_re(inout[56 * iostride]) = tmp983 - tmp986;
			 c_re(inout[24 * iostride]) = tmp983 + tmp986;
		    }
	       }
	  }
	  {
	       fftw_real tmp995;
	       fftw_real tmp1023;
	       fftw_real tmp1026;
	       fftw_real tmp1068;
	       fftw_real tmp1002;
	       fftw_real tmp1060;
	       fftw_real tmp1063;
	       fftw_real tmp1069;
	       fftw_real tmp1010;
	       fftw_real tmp1021;
	       fftw_real tmp1030;
	       fftw_real tmp1037;
	       fftw_real tmp1017;
	       fftw_real tmp1020;
	       fftw_real tmp1033;
	       fftw_real tmp1036;
	       ASSERT_ALIGNED_DOUBLE;
	       {
		    fftw_real tmp991;
		    fftw_real tmp994;
		    fftw_real tmp1024;
		    fftw_real tmp1025;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp991 = tmp19 - tmp42;
		    tmp994 = tmp992 - tmp993;
		    tmp995 = tmp991 - tmp994;
		    tmp1023 = tmp991 + tmp994;
		    tmp1024 = tmp999 - tmp1000;
		    tmp1025 = tmp997 + tmp996;
		    tmp1026 = K707106781 * (tmp1024 + tmp1025);
		    tmp1068 = K707106781 * (tmp1024 - tmp1025);
	       }
	       {
		    fftw_real tmp998;
		    fftw_real tmp1001;
		    fftw_real tmp1061;
		    fftw_real tmp1062;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp998 = tmp996 - tmp997;
		    tmp1001 = tmp999 + tmp1000;
		    tmp1002 = K707106781 * (tmp998 - tmp1001);
		    tmp1060 = K707106781 * (tmp1001 + tmp998);
		    tmp1061 = tmp66 - tmp89;
		    tmp1062 = tmp1047 - tmp1042;
		    tmp1063 = tmp1061 + tmp1062;
		    tmp1069 = tmp1062 - tmp1061;
	       }
	       {
		    fftw_real tmp1006;
		    fftw_real tmp1009;
		    fftw_real tmp1028;
		    fftw_real tmp1029;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp1006 = tmp1004 - tmp1005;
		    tmp1009 = tmp1007 - tmp1008;
		    tmp1010 =
			(K382683432 * tmp1006) - (K923879532 * tmp1009);
		    tmp1021 =
			(K923879532 * tmp1006) + (K382683432 * tmp1009);
		    tmp1028 = tmp1004 + tmp1005;
		    tmp1029 = tmp1007 + tmp1008;
		    tmp1030 =
			(K923879532 * tmp1028) - (K382683432 * tmp1029);
		    tmp1037 =
			(K382683432 * tmp1028) + (K923879532 * tmp1029);
	       }
	       {
		    fftw_real tmp1013;
		    fftw_real tmp1016;
		    fftw_real tmp1031;
		    fftw_real tmp1032;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp1013 = tmp1011 - tmp1012;
		    tmp1016 = tmp1014 - tmp1015;
		    tmp1017 =
			(K382683432 * tmp1013) + (K923879532 * tmp1016);
		    tmp1020 =
			(K382683432 * tmp1016) - (K923879532 * tmp1013);
		    tmp1031 = tmp1011 + tmp1012;
		    tmp1032 = tmp1014 + tmp1015;
		    tmp1033 =
			(K923879532 * tmp1031) + (K382683432 * tmp1032);
		    tmp1036 =
			(K923879532 * tmp1032) - (K382683432 * tmp1031);
	       }
	       {
		    fftw_real tmp1003;
		    fftw_real tmp1018;
		    fftw_real tmp1019;
		    fftw_real tmp1022;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp1003 = tmp995 - tmp1002;
		    tmp1018 = tmp1010 - tmp1017;
		    c_re(inout[60 * iostride]) = tmp1003 - tmp1018;
		    c_re(inout[28 * iostride]) = tmp1003 + tmp1018;
		    tmp1019 = tmp995 + tmp1002;
		    tmp1022 = tmp1020 + tmp1021;
		    c_re(inout[44 * iostride]) = tmp1019 - tmp1022;
		    c_re(inout[12 * iostride]) = tmp1019 + tmp1022;
	       }
	       {
		    fftw_real tmp1071;
		    fftw_real tmp1072;
		    fftw_real tmp1067;
		    fftw_real tmp1070;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp1071 = tmp1020 - tmp1021;
		    tmp1072 = tmp1069 - tmp1068;
		    c_im(inout[28 * iostride]) = tmp1071 + tmp1072;
		    c_im(inout[60 * iostride]) = tmp1072 - tmp1071;
		    tmp1067 = tmp1017 + tmp1010;
		    tmp1070 = tmp1068 + tmp1069;
		    c_im(inout[12 * iostride]) = tmp1067 + tmp1070;
		    c_im(inout[44 * iostride]) = tmp1070 - tmp1067;
	       }
	       {
		    fftw_real tmp1027;
		    fftw_real tmp1034;
		    fftw_real tmp1035;
		    fftw_real tmp1038;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp1027 = tmp1023 - tmp1026;
		    tmp1034 = tmp1030 - tmp1033;
		    c_re(inout[52 * iostride]) = tmp1027 - tmp1034;
		    c_re(inout[20 * iostride]) = tmp1027 + tmp1034;
		    tmp1035 = tmp1023 + tmp1026;
		    tmp1038 = tmp1036 + tmp1037;
		    c_re(inout[36 * iostride]) = tmp1035 - tmp1038;
		    c_re(inout[4 * iostride]) = tmp1035 + tmp1038;
	       }
	       {
		    fftw_real tmp1065;
		    fftw_real tmp1066;
		    fftw_real tmp1059;
		    fftw_real tmp1064;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp1065 = tmp1036 - tmp1037;
		    tmp1066 = tmp1063 - tmp1060;
		    c_im(inout[20 * iostride]) = tmp1065 + tmp1066;
		    c_im(inout[52 * iostride]) = tmp1066 - tmp1065;
		    tmp1059 = tmp1033 + tmp1030;
		    tmp1064 = tmp1060 + tmp1063;
		    c_im(inout[4 * iostride]) = tmp1059 + tmp1064;
		    c_im(inout[36 * iostride]) = tmp1064 - tmp1059;
	       }
	  }
	  {
	       fftw_real tmp419;
	       fftw_real tmp591;
	       fftw_real tmp1125;
	       fftw_real tmp1131;
	       fftw_real tmp474;
	       fftw_real tmp1122;
	       fftw_real tmp594;
	       fftw_real tmp1130;
	       fftw_real tmp530;
	       fftw_real tmp589;
	       fftw_real tmp598;
	       fftw_real tmp605;
	       fftw_real tmp585;
	       fftw_real tmp588;
	       fftw_real tmp601;
	       fftw_real tmp604;
	       ASSERT_ALIGNED_DOUBLE;
	       {
		    fftw_real tmp395;
		    fftw_real tmp418;
		    fftw_real tmp1123;
		    fftw_real tmp1124;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp395 = tmp383 - tmp394;
		    tmp418 = tmp406 - tmp417;
		    tmp419 = tmp395 - tmp418;
		    tmp591 = tmp395 + tmp418;
		    tmp1123 = tmp608 - tmp609;
		    tmp1124 = tmp1109 - tmp1106;
		    tmp1125 = tmp1123 + tmp1124;
		    tmp1131 = tmp1124 - tmp1123;
	       }
	       {
		    fftw_real tmp446;
		    fftw_real tmp593;
		    fftw_real tmp473;
		    fftw_real tmp592;
		    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 =
			     (K555570233 * tmp436) - (K831469612 * tmp445);
			 tmp593 =
			     (K555570233 * tmp445) + (K831469612 * tmp436);
			 tmp463 = tmp451 - tmp462;
			 tmp472 = tmp468 - tmp471;
			 tmp473 =
			     (K831469612 * tmp463) + (K555570233 * tmp472);
			 tmp592 =
			     (K555570233 * tmp463) - (K831469612 * tmp472);
		    }
		    tmp474 = tmp446 - tmp473;
		    tmp1122 = tmp473 + tmp446;
		    tmp594 = tmp592 + tmp593;
		    tmp1130 = tmp592 - tmp593;
	       }
	       {
		    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 = (K290284677 * tmp516) - (K956940335 * tmp529);
		    tmp589 = (K956940335 * tmp516) + (K290284677 * tmp529);
		    tmp598 = (K881921264 * tmp596) - (K471396736 * tmp597);
		    tmp605 = (K471396736 * tmp596) + (K881921264 * 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 = (K290284677 * tmp571) + (K956940335 * tmp584);
		    tmp588 = (K290284677 * tmp584) - (K956940335 * tmp571);
		    tmp601 = (K881921264 * tmp599) + (K471396736 * tmp600);
		    tmp604 = (K881921264 * tmp600) - (K471396736 * 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[45 * iostride]) = tmp475 - tmp586;
		    c_re(inout[13 * iostride]) = tmp475 + tmp586;
		    tmp587 = tmp419 - tmp474;
		    tmp590 = tmp588 - tmp589;
		    c_re(inout[61 * iostride]) = tmp587 - tmp590;
		    c_re(inout[29 * iostride]) = tmp587 + tmp590;
	       }
	       {
		    fftw_real tmp1129;
		    fftw_real tmp1132;
		    fftw_real tmp1133;
		    fftw_real tmp1134;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp1129 = tmp589 + tmp588;
		    tmp1132 = tmp1130 + tmp1131;
		    c_im(inout[13 * iostride]) = tmp1129 + tmp1132;
		    c_im(inout[45 * iostride]) = tmp1132 - tmp1129;
		    tmp1133 = tmp530 - tmp585;
		    tmp1134 = tmp1131 - tmp1130;
		    c_im(inout[29 * iostride]) = tmp1133 + tmp1134;
		    c_im(inout[61 * iostride]) = tmp1134 - tmp1133;
	       }
	       {
		    fftw_real tmp595;
		    fftw_real tmp602;
		    fftw_real tmp603;