ftw_32.c

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

		    fftw_real tmp380;
		    fftw_real tmp381;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp172 = tmp166 + tmp171;
		    tmp183 = tmp177 + tmp182;
		    tmp184 = tmp172 + tmp183;
		    tmp387 = tmp183 - tmp172;
		    tmp380 = tmp260 + tmp261;
		    tmp381 = tmp266 + tmp267;
		    tmp382 = tmp380 - tmp381;
		    tmp423 = tmp380 + tmp381;
	       }
	       {
		    fftw_real tmp264;
		    fftw_real tmp269;
		    fftw_real tmp277;
		    fftw_real tmp278;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp264 = tmp262 - tmp263;
		    tmp269 = tmp265 + tmp268;
		    tmp270 = K707106781 * (tmp264 - tmp269);
		    tmp327 = K707106781 * (tmp264 + tmp269);
		    tmp277 = tmp265 - tmp268;
		    tmp278 = tmp263 + tmp262;
		    tmp279 = K707106781 * (tmp277 - tmp278);
		    tmp324 = K707106781 * (tmp278 + tmp277);
	       }
	  }
	  {
	       fftw_real tmp48;
	       fftw_real tmp204;
	       fftw_real tmp64;
	       fftw_real tmp211;
	       fftw_real tmp53;
	       fftw_real tmp205;
	       fftw_real tmp59;
	       fftw_real tmp210;
	       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);
		    tmp204 = (tmp46 * tmp45) + (tmp44 * tmp47);
	       }
	       {
		    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);
		    tmp211 = (tmp62 * tmp61) + (tmp60 * tmp63);
	       }
	       {
		    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);
		    tmp205 = (tmp51 * tmp50) + (tmp49 * tmp52);
	       }
	       {
		    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);
		    tmp210 = (tmp57 * tmp56) + (tmp55 * tmp58);
	       }
	       {
		    fftw_real tmp54;
		    fftw_real tmp65;
		    fftw_real tmp209;
		    fftw_real tmp212;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp54 = tmp48 + tmp53;
		    tmp65 = tmp59 + tmp64;
		    tmp66 = tmp54 + tmp65;
		    tmp359 = tmp54 - tmp65;
		    tmp209 = tmp48 - tmp53;
		    tmp212 = tmp210 - tmp211;
		    tmp213 = tmp209 - tmp212;
		    tmp309 = tmp209 + tmp212;
	       }
	       {
		    fftw_real tmp356;
		    fftw_real tmp357;
		    fftw_real tmp206;
		    fftw_real tmp207;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp356 = tmp204 + tmp205;
		    tmp357 = tmp210 + tmp211;
		    tmp358 = tmp356 - tmp357;
		    tmp412 = tmp356 + tmp357;
		    tmp206 = tmp204 - tmp205;
		    tmp207 = tmp59 - tmp64;
		    tmp208 = tmp206 + tmp207;
		    tmp308 = tmp206 - tmp207;
	       }
	  }
	  {
	       fftw_real tmp96;
	       fftw_real tmp228;
	       fftw_real tmp112;
	       fftw_real tmp247;
	       fftw_real tmp101;
	       fftw_real tmp229;
	       fftw_real tmp107;
	       fftw_real tmp246;
	       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);
		    tmp228 = (tmp94 * tmp93) + (tmp92 * tmp95);
	       }
	       {
		    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);
		    tmp247 = (tmp110 * tmp109) + (tmp108 * tmp111);
	       }
	       {
		    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);
		    tmp229 = (tmp99 * tmp98) + (tmp97 * tmp100);
	       }
	       {
		    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);
		    tmp246 = (tmp105 * tmp104) + (tmp103 * tmp106);
	       }
	       {
		    fftw_real tmp102;
		    fftw_real tmp113;
		    fftw_real tmp245;
		    fftw_real tmp248;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp102 = tmp96 + tmp101;
		    tmp113 = tmp107 + tmp112;
		    tmp114 = tmp102 + tmp113;
		    tmp373 = tmp102 - tmp113;
		    tmp245 = tmp96 - tmp101;
		    tmp248 = tmp246 - tmp247;
		    tmp249 = tmp245 - tmp248;
		    tmp316 = tmp245 + tmp248;
	       }
	       {
		    fftw_real tmp368;
		    fftw_real tmp369;
		    fftw_real tmp230;
		    fftw_real tmp231;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp368 = tmp228 + tmp229;
		    tmp369 = tmp246 + tmp247;
		    tmp370 = tmp368 - tmp369;
		    tmp417 = tmp368 + tmp369;
		    tmp230 = tmp228 - tmp229;
		    tmp231 = tmp107 - tmp112;
		    tmp232 = tmp230 + tmp231;
		    tmp319 = tmp230 - tmp231;
	       }
	  }
	  {
	       fftw_real tmp71;
	       fftw_real tmp215;
	       fftw_real tmp87;
	       fftw_real tmp222;
	       fftw_real tmp76;
	       fftw_real tmp216;
	       fftw_real tmp82;
	       fftw_real tmp221;
	       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);
		    tmp215 = (tmp69 * tmp68) + (tmp67 * tmp70);
	       }
	       {
		    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);
		    tmp222 = (tmp85 * tmp84) + (tmp83 * tmp86);
	       }
	       {
		    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);
		    tmp216 = (tmp74 * tmp73) + (tmp72 * tmp75);
	       }
	       {
		    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);
		    tmp221 = (tmp80 * tmp79) + (tmp78 * tmp81);
	       }
	       {
		    fftw_real tmp77;
		    fftw_real tmp88;
		    fftw_real tmp220;
		    fftw_real tmp223;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp77 = tmp71 + tmp76;
		    tmp88 = tmp82 + tmp87;
		    tmp89 = tmp77 + tmp88;
		    tmp361 = tmp77 - tmp88;
		    tmp220 = tmp71 - tmp76;
		    tmp223 = tmp221 - tmp222;
		    tmp224 = tmp220 - tmp223;
		    tmp312 = tmp220 + tmp223;
	       }
	       {
		    fftw_real tmp362;
		    fftw_real tmp363;
		    fftw_real tmp217;
		    fftw_real tmp218;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp362 = tmp215 + tmp216;
		    tmp363 = tmp221 + tmp222;
		    tmp364 = tmp362 - tmp363;
		    tmp413 = tmp362 + tmp363;
		    tmp217 = tmp215 - tmp216;
		    tmp218 = tmp82 - tmp87;
		    tmp219 = tmp217 + tmp218;
		    tmp311 = tmp217 - tmp218;
	       }
	  }
	  {
	       fftw_real tmp119;
	       fftw_real tmp239;
	       fftw_real tmp124;
	       fftw_real tmp240;
	       fftw_real tmp238;
	       fftw_real tmp241;
	       fftw_real tmp130;
	       fftw_real tmp234;
	       fftw_real tmp135;
	       fftw_real tmp235;
	       fftw_real tmp233;
	       fftw_real tmp236;
	       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);
		    tmp239 = (tmp117 * tmp116) + (tmp115 * tmp118);
	       }
	       {
		    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);
		    tmp240 = (tmp122 * tmp121) + (tmp120 * tmp123);
	       }
	       tmp238 = tmp119 - tmp124;
	       tmp241 = tmp239 - tmp240;
	       {
		    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);
		    tmp234 = (tmp128 * tmp127) + (tmp126 * tmp129);
	       }
	       {
		    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);
		    tmp235 = (tmp133 * tmp132) + (tmp131 * tmp134);
	       }
	       tmp233 = tmp130 - tmp135;
	       tmp236 = tmp234 - tmp235;
	       {
		    fftw_real tmp125;
		    fftw_real tmp136;
		    fftw_real tmp374;
		    fftw_real tmp375;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp125 = tmp119 + tmp124;
		    tmp136 = tmp130 + tmp135;
		    tmp137 = tmp125 + tmp136;
		    tmp371 = tmp136 - tmp125;
		    tmp374 = tmp239 + tmp240;
		    tmp375 = tmp234 + tmp235;
		    tmp376 = tmp374 - tmp375;
		    tmp418 = tmp374 + tmp375;
	       }
	       {
		    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);
		    tmp317 = K707106781 * (tmp242 + tmp237);
		    tmp250 = tmp241 - tmp238;
		    tmp251 = tmp233 + tmp236;
		    tmp252 = K707106781 * (tmp250 - tmp251);
		    tmp320 = K707106781 * (tmp250 + tmp251);
	       }
	  }
	  {
	       fftw_real tmp91;
	       fftw_real tmp431;
	       fftw_real tmp444;
	       fftw_real tmp446;
	       fftw_real tmp186;
	       fftw_real tmp445;
	       fftw_real tmp434;
	       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;
		    tmp431 = tmp43 - tmp90;
		    tmp436 = tmp412 + tmp413;
		    tmp443 = tmp437 + tmp442;
		    tmp444 = tmp436 + tmp443;
		    tmp446 = tmp443 - tmp436;
	       }
	       {
		    fftw_real tmp138;
		    fftw_real tmp185;
		    fftw_real tmp432;
		    fftw_real tmp433;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp138 = tmp114 + tmp137;
		    tmp185 = tmp161 + tmp184;
		    tmp186 = tmp138 + tmp185;
		    tmp445 = tmp185 - tmp138;
		    tmp432 = tmp417 + tmp418;
		    tmp433 = tmp422 + tmp423;
		    tmp434 = tmp432 - tmp433;
		    tmp435 = tmp432 + tmp433;
	       }
	       c_re(inout[16 * iostride]) = tmp91 - tmp186;
	       c_re(inout[0]) = tmp91 + tmp186;
	       c_re(inout[24 * iostride]) = tmp431 - tmp434;
	       c_re(inout[8 * iostride]) = tmp431 + tmp434;
	       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 tmp415;
	       fftw_real tmp427;
	       fftw_real tmp450;
	       fftw_real tmp452;
	       fftw_real tmp420;
	       fftw_real tmp428;
	       fftw_real tmp425;
	       fftw_real tmp429;
	       ASSERT_ALIGNED_DOUBLE;