ftwi_32.c

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

		    fftw_real tmp373;
		    fftw_real tmp374;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp172 = tmp166 + tmp171;
		    tmp183 = tmp177 + tmp182;
		    tmp184 = tmp172 + tmp183;
		    tmp401 = tmp172 - tmp183;
		    tmp373 = tmp261 + tmp262;
		    tmp374 = tmp266 + tmp267;
		    tmp375 = tmp373 + tmp374;
		    tmp404 = tmp374 - tmp373;
	       }
	       {
		    fftw_real tmp264;
		    fftw_real tmp269;
		    fftw_real tmp277;
		    fftw_real tmp278;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp264 = tmp260 - tmp263;
		    tmp269 = tmp265 + tmp268;
		    tmp270 = K707106781 * (tmp264 + tmp269);
		    tmp317 = K707106781 * (tmp264 - tmp269);
		    tmp277 = tmp260 + tmp263;
		    tmp278 = tmp268 - tmp265;
		    tmp279 = K707106781 * (tmp277 + tmp278);
		    tmp320 = K707106781 * (tmp278 - tmp277);
	       }
	  }
	  {
	       fftw_real tmp48;
	       fftw_real tmp215;
	       fftw_real tmp64;
	       fftw_real tmp222;
	       fftw_real tmp53;
	       fftw_real tmp216;
	       fftw_real tmp59;
	       fftw_real tmp221;
	       ASSERT_ALIGNED_DOUBLE;
	       {
		    fftw_real tmp45;
		    fftw_real tmp47;
		    fftw_real tmp44;
		    fftw_real tmp46;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp45 = c_re(inout[2 * iostride]);
		    tmp47 = c_im(inout[2 * iostride]);
		    tmp44 = c_re(W[1]);
		    tmp46 = c_im(W[1]);
		    tmp48 = (tmp44 * tmp45) + (tmp46 * tmp47);
		    tmp215 = (tmp44 * tmp47) - (tmp46 * tmp45);
	       }
	       {
		    fftw_real tmp61;
		    fftw_real tmp63;
		    fftw_real tmp60;
		    fftw_real tmp62;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp61 = c_re(inout[26 * iostride]);
		    tmp63 = c_im(inout[26 * iostride]);
		    tmp60 = c_re(W[25]);
		    tmp62 = c_im(W[25]);
		    tmp64 = (tmp60 * tmp61) + (tmp62 * tmp63);
		    tmp222 = (tmp60 * tmp63) - (tmp62 * tmp61);
	       }
	       {
		    fftw_real tmp50;
		    fftw_real tmp52;
		    fftw_real tmp49;
		    fftw_real tmp51;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp50 = c_re(inout[18 * iostride]);
		    tmp52 = c_im(inout[18 * iostride]);
		    tmp49 = c_re(W[17]);
		    tmp51 = c_im(W[17]);
		    tmp53 = (tmp49 * tmp50) + (tmp51 * tmp52);
		    tmp216 = (tmp49 * tmp52) - (tmp51 * tmp50);
	       }
	       {
		    fftw_real tmp56;
		    fftw_real tmp58;
		    fftw_real tmp55;
		    fftw_real tmp57;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp56 = c_re(inout[10 * iostride]);
		    tmp58 = c_im(inout[10 * iostride]);
		    tmp55 = c_re(W[9]);
		    tmp57 = c_im(W[9]);
		    tmp59 = (tmp55 * tmp56) + (tmp57 * tmp58);
		    tmp221 = (tmp55 * tmp58) - (tmp57 * tmp56);
	       }
	       {
		    fftw_real tmp54;
		    fftw_real tmp65;
		    fftw_real tmp220;
		    fftw_real tmp223;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp54 = tmp48 + tmp53;
		    tmp65 = tmp59 + tmp64;
		    tmp66 = tmp54 + tmp65;
		    tmp395 = tmp54 - tmp65;
		    tmp220 = tmp48 - tmp53;
		    tmp223 = tmp221 - tmp222;
		    tmp224 = tmp220 - tmp223;
		    tmp312 = tmp220 + tmp223;
	       }
	       {
		    fftw_real tmp355;
		    fftw_real tmp356;
		    fftw_real tmp217;
		    fftw_real tmp218;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp355 = tmp215 + tmp216;
		    tmp356 = tmp221 + tmp222;
		    tmp357 = tmp355 + tmp356;
		    tmp396 = tmp355 - tmp356;
		    tmp217 = tmp215 - tmp216;
		    tmp218 = tmp59 - tmp64;
		    tmp219 = tmp217 + tmp218;
		    tmp311 = tmp217 - tmp218;
	       }
	  }
	  {
	       fftw_real tmp96;
	       fftw_real tmp245;
	       fftw_real tmp112;
	       fftw_real tmp230;
	       fftw_real tmp101;
	       fftw_real tmp246;
	       fftw_real tmp107;
	       fftw_real tmp229;
	       ASSERT_ALIGNED_DOUBLE;
	       {
		    fftw_real tmp93;
		    fftw_real tmp95;
		    fftw_real tmp92;
		    fftw_real tmp94;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp93 = c_re(inout[iostride]);
		    tmp95 = c_im(inout[iostride]);
		    tmp92 = c_re(W[0]);
		    tmp94 = c_im(W[0]);
		    tmp96 = (tmp92 * tmp93) + (tmp94 * tmp95);
		    tmp245 = (tmp92 * tmp95) - (tmp94 * tmp93);
	       }
	       {
		    fftw_real tmp109;
		    fftw_real tmp111;
		    fftw_real tmp108;
		    fftw_real tmp110;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp109 = c_re(inout[25 * iostride]);
		    tmp111 = c_im(inout[25 * iostride]);
		    tmp108 = c_re(W[24]);
		    tmp110 = c_im(W[24]);
		    tmp112 = (tmp108 * tmp109) + (tmp110 * tmp111);
		    tmp230 = (tmp108 * tmp111) - (tmp110 * tmp109);
	       }
	       {
		    fftw_real tmp98;
		    fftw_real tmp100;
		    fftw_real tmp97;
		    fftw_real tmp99;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp98 = c_re(inout[17 * iostride]);
		    tmp100 = c_im(inout[17 * iostride]);
		    tmp97 = c_re(W[16]);
		    tmp99 = c_im(W[16]);
		    tmp101 = (tmp97 * tmp98) + (tmp99 * tmp100);
		    tmp246 = (tmp97 * tmp100) - (tmp99 * tmp98);
	       }
	       {
		    fftw_real tmp104;
		    fftw_real tmp106;
		    fftw_real tmp103;
		    fftw_real tmp105;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp104 = c_re(inout[9 * iostride]);
		    tmp106 = c_im(inout[9 * iostride]);
		    tmp103 = c_re(W[8]);
		    tmp105 = c_im(W[8]);
		    tmp107 = (tmp103 * tmp104) + (tmp105 * tmp106);
		    tmp229 = (tmp103 * tmp106) - (tmp105 * tmp104);
	       }
	       {
		    fftw_real tmp102;
		    fftw_real tmp113;
		    fftw_real tmp247;
		    fftw_real tmp248;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp102 = tmp96 + tmp101;
		    tmp113 = tmp107 + tmp112;
		    tmp114 = tmp102 + tmp113;
		    tmp410 = tmp102 - tmp113;
		    tmp247 = tmp245 - tmp246;
		    tmp248 = tmp107 - tmp112;
		    tmp249 = tmp247 + tmp248;
		    tmp323 = tmp247 - tmp248;
	       }
	       {
		    fftw_real tmp361;
		    fftw_real tmp362;
		    fftw_real tmp228;
		    fftw_real tmp231;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp361 = tmp245 + tmp246;
		    tmp362 = tmp229 + tmp230;
		    tmp363 = tmp361 + tmp362;
		    tmp407 = tmp361 - tmp362;
		    tmp228 = tmp96 - tmp101;
		    tmp231 = tmp229 - tmp230;
		    tmp232 = tmp228 - tmp231;
		    tmp326 = tmp228 + tmp231;
	       }
	  }
	  {
	       fftw_real tmp71;
	       fftw_real tmp204;
	       fftw_real tmp87;
	       fftw_real tmp211;
	       fftw_real tmp76;
	       fftw_real tmp205;
	       fftw_real tmp82;
	       fftw_real tmp210;
	       ASSERT_ALIGNED_DOUBLE;
	       {
		    fftw_real tmp68;
		    fftw_real tmp70;
		    fftw_real tmp67;
		    fftw_real tmp69;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp68 = c_re(inout[30 * iostride]);
		    tmp70 = c_im(inout[30 * iostride]);
		    tmp67 = c_re(W[29]);
		    tmp69 = c_im(W[29]);
		    tmp71 = (tmp67 * tmp68) + (tmp69 * tmp70);
		    tmp204 = (tmp67 * tmp70) - (tmp69 * tmp68);
	       }
	       {
		    fftw_real tmp84;
		    fftw_real tmp86;
		    fftw_real tmp83;
		    fftw_real tmp85;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp84 = c_re(inout[22 * iostride]);
		    tmp86 = c_im(inout[22 * iostride]);
		    tmp83 = c_re(W[21]);
		    tmp85 = c_im(W[21]);
		    tmp87 = (tmp83 * tmp84) + (tmp85 * tmp86);
		    tmp211 = (tmp83 * tmp86) - (tmp85 * tmp84);
	       }
	       {
		    fftw_real tmp73;
		    fftw_real tmp75;
		    fftw_real tmp72;
		    fftw_real tmp74;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp73 = c_re(inout[14 * iostride]);
		    tmp75 = c_im(inout[14 * iostride]);
		    tmp72 = c_re(W[13]);
		    tmp74 = c_im(W[13]);
		    tmp76 = (tmp72 * tmp73) + (tmp74 * tmp75);
		    tmp205 = (tmp72 * tmp75) - (tmp74 * tmp73);
	       }
	       {
		    fftw_real tmp79;
		    fftw_real tmp81;
		    fftw_real tmp78;
		    fftw_real tmp80;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp79 = c_re(inout[6 * iostride]);
		    tmp81 = c_im(inout[6 * iostride]);
		    tmp78 = c_re(W[5]);
		    tmp80 = c_im(W[5]);
		    tmp82 = (tmp78 * tmp79) + (tmp80 * tmp81);
		    tmp210 = (tmp78 * tmp81) - (tmp80 * tmp79);
	       }
	       {
		    fftw_real tmp77;
		    fftw_real tmp88;
		    fftw_real tmp209;
		    fftw_real tmp212;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp77 = tmp71 + tmp76;
		    tmp88 = tmp82 + tmp87;
		    tmp89 = tmp77 + tmp88;
		    tmp393 = tmp77 - tmp88;
		    tmp209 = tmp71 - tmp76;
		    tmp212 = tmp210 - tmp211;
		    tmp213 = tmp209 - tmp212;
		    tmp309 = tmp209 + tmp212;
	       }
	       {
		    fftw_real tmp352;
		    fftw_real tmp353;
		    fftw_real tmp206;
		    fftw_real tmp207;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp352 = tmp204 + tmp205;
		    tmp353 = tmp210 + tmp211;
		    tmp354 = tmp352 + tmp353;
		    tmp392 = tmp352 - tmp353;
		    tmp206 = tmp204 - tmp205;
		    tmp207 = tmp82 - tmp87;
		    tmp208 = tmp206 + tmp207;
		    tmp308 = tmp206 - tmp207;
	       }
	  }
	  {
	       fftw_real tmp119;
	       fftw_real tmp234;
	       fftw_real tmp124;
	       fftw_real tmp235;
	       fftw_real tmp233;
	       fftw_real tmp236;
	       fftw_real tmp130;
	       fftw_real tmp239;
	       fftw_real tmp135;
	       fftw_real tmp240;
	       fftw_real tmp238;
	       fftw_real tmp241;
	       ASSERT_ALIGNED_DOUBLE;
	       {
		    fftw_real tmp116;
		    fftw_real tmp118;
		    fftw_real tmp115;
		    fftw_real tmp117;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp116 = c_re(inout[5 * iostride]);
		    tmp118 = c_im(inout[5 * iostride]);
		    tmp115 = c_re(W[4]);
		    tmp117 = c_im(W[4]);
		    tmp119 = (tmp115 * tmp116) + (tmp117 * tmp118);
		    tmp234 = (tmp115 * tmp118) - (tmp117 * tmp116);
	       }
	       {
		    fftw_real tmp121;
		    fftw_real tmp123;
		    fftw_real tmp120;
		    fftw_real tmp122;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp121 = c_re(inout[21 * iostride]);
		    tmp123 = c_im(inout[21 * iostride]);
		    tmp120 = c_re(W[20]);
		    tmp122 = c_im(W[20]);
		    tmp124 = (tmp120 * tmp121) + (tmp122 * tmp123);
		    tmp235 = (tmp120 * tmp123) - (tmp122 * tmp121);
	       }
	       tmp233 = tmp119 - tmp124;
	       tmp236 = tmp234 - tmp235;
	       {
		    fftw_real tmp127;
		    fftw_real tmp129;
		    fftw_real tmp126;
		    fftw_real tmp128;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp127 = c_re(inout[29 * iostride]);
		    tmp129 = c_im(inout[29 * iostride]);
		    tmp126 = c_re(W[28]);
		    tmp128 = c_im(W[28]);
		    tmp130 = (tmp126 * tmp127) + (tmp128 * tmp129);
		    tmp239 = (tmp126 * tmp129) - (tmp128 * tmp127);
	       }
	       {
		    fftw_real tmp132;
		    fftw_real tmp134;
		    fftw_real tmp131;
		    fftw_real tmp133;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp132 = c_re(inout[13 * iostride]);
		    tmp134 = c_im(inout[13 * iostride]);
		    tmp131 = c_re(W[12]);
		    tmp133 = c_im(W[12]);
		    tmp135 = (tmp131 * tmp132) + (tmp133 * tmp134);
		    tmp240 = (tmp131 * tmp134) - (tmp133 * tmp132);
	       }
	       tmp238 = tmp130 - tmp135;
	       tmp241 = tmp239 - tmp240;
	       {
		    fftw_real tmp125;
		    fftw_real tmp136;
		    fftw_real tmp364;
		    fftw_real tmp365;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp125 = tmp119 + tmp124;
		    tmp136 = tmp130 + tmp135;
		    tmp137 = tmp125 + tmp136;
		    tmp408 = tmp125 - tmp136;
		    tmp364 = tmp234 + tmp235;
		    tmp365 = tmp239 + tmp240;
		    tmp366 = tmp364 + tmp365;
		    tmp411 = tmp365 - tmp364;
	       }
	       {
		    fftw_real tmp237;
		    fftw_real tmp242;
		    fftw_real tmp250;
		    fftw_real tmp251;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp237 = tmp233 - tmp236;
		    tmp242 = tmp238 + tmp241;
		    tmp243 = K707106781 * (tmp237 + tmp242);
		    tmp324 = K707106781 * (tmp237 - tmp242);
		    tmp250 = tmp233 + tmp236;
		    tmp251 = tmp241 - tmp238;
		    tmp252 = K707106781 * (tmp250 + tmp251);
		    tmp327 = K707106781 * (tmp251 - tmp250);
	       }
	  }
	  {
	       fftw_real tmp91;
	       fftw_real tmp383;
	       fftw_real tmp444;
	       fftw_real tmp446;
	       fftw_real tmp186;
	       fftw_real tmp445;
	       fftw_real tmp386;
	       fftw_real tmp435;
	       ASSERT_ALIGNED_DOUBLE;
	       {
		    fftw_real tmp43;
		    fftw_real tmp90;
		    fftw_real tmp436;
		    fftw_real tmp443;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp43 = tmp19 + tmp42;
		    tmp90 = tmp66 + tmp89;
		    tmp91 = tmp43 + tmp90;
		    tmp383 = tmp43 - tmp90;
		    tmp436 = tmp357 + tmp354;
		    tmp443 = tmp437 + tmp442;
		    tmp444 = tmp436 + tmp443;
		    tmp446 = tmp443 - tmp436;
	       }
	       {
		    fftw_real tmp138;
		    fftw_real tmp185;
		    fftw_real tmp384;
		    fftw_real tmp385;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp138 = tmp114 + tmp137;
		    tmp185 = tmp161 + tmp184;
		    tmp186 = tmp138 + tmp185;
		    tmp445 = tmp138 - tmp185;
		    tmp384 = tmp372 + tmp375;
		    tmp385 = tmp363 + tmp366;
		    tmp386 = tmp384 - tmp385;
		    tmp435 = tmp385 + tmp384;
	       }
	       c_re(inout[16 * iostride]) = tmp91 - tmp186;
	       c_re(inout[0]) = tmp91 + tmp186;
	       c_re(inout[24 * iostride]) = tmp383 - tmp386;
	       c_re(inout[8 * iostride]) = tmp383 + tmp386;
	       c_im(inout[0]) = tmp435 + tmp444;
	       c_im(inout[16 * iostride]) = tmp444 - tmp435;
	       c_im(inout[8 * iostride]) = tmp445 + tmp446;
	       c_im(inout[24 * iostride]) = tmp446 - tmp445;
	  }
	  {
	       fftw_real tmp359;
	       fftw_real tmp379;
	       fftw_real tmp450;
	       fftw_real tmp452;
	       fftw_real tmp368;
	       fftw_real tmp381;
	       fftw_real tmp377;
	       fftw_real tmp380;
	       ASSERT_ALIGNED_DOUBLE;