fcr_128.c

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

	       tmp233 = tmp94 - tmp91;
	       tmp231 = tmp227 + tmp230;
	       tmp403 = tmp230 - tmp227;
	       tmp603 = tmp601 - tmp602;
	       tmp814 = tmp601 + tmp602;
	       tmp606 = tmp604 + tmp605;
	       tmp815 = tmp605 - tmp604;
	  }
	  {
	       fftw_real tmp76;
	       fftw_real tmp586;
	       fftw_real tmp211;
	       fftw_real tmp590;
	       fftw_real tmp79;
	       fftw_real tmp589;
	       fftw_real tmp208;
	       fftw_real tmp587;
	       ASSERT_ALIGNED_DOUBLE;
	       {
		    fftw_real tmp74;
		    fftw_real tmp75;
		    fftw_real tmp209;
		    fftw_real tmp210;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp74 = real_input[15 * real_istride];
		    tmp75 = real_input[49 * real_istride];
		    tmp76 = tmp74 + tmp75;
		    tmp586 = tmp74 - tmp75;
		    tmp209 = imag_input[15 * imag_istride];
		    tmp210 = imag_input[49 * imag_istride];
		    tmp211 = tmp209 - tmp210;
		    tmp590 = tmp209 + tmp210;
	       }
	       {
		    fftw_real tmp77;
		    fftw_real tmp78;
		    fftw_real tmp206;
		    fftw_real tmp207;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp77 = real_input[17 * real_istride];
		    tmp78 = real_input[47 * real_istride];
		    tmp79 = tmp77 + tmp78;
		    tmp589 = tmp78 - tmp77;
		    tmp206 = imag_input[17 * imag_istride];
		    tmp207 = imag_input[47 * imag_istride];
		    tmp208 = tmp206 - tmp207;
		    tmp587 = tmp206 + tmp207;
	       }
	       tmp80 = tmp76 + tmp79;
	       tmp214 = tmp79 - tmp76;
	       tmp212 = tmp208 + tmp211;
	       tmp400 = tmp211 - tmp208;
	       tmp588 = tmp586 - tmp587;
	       tmp807 = tmp586 + tmp587;
	       tmp591 = tmp589 + tmp590;
	       tmp808 = tmp590 - tmp589;
	  }
	  {
	       fftw_real tmp84;
	       fftw_real tmp608;
	       fftw_real tmp239;
	       fftw_real tmp612;
	       fftw_real tmp87;
	       fftw_real tmp611;
	       fftw_real tmp236;
	       fftw_real tmp609;
	       ASSERT_ALIGNED_DOUBLE;
	       {
		    fftw_real tmp82;
		    fftw_real tmp83;
		    fftw_real tmp237;
		    fftw_real tmp238;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp82 = real_input[7 * real_istride];
		    tmp83 = real_input[57 * real_istride];
		    tmp84 = tmp82 + tmp83;
		    tmp608 = tmp82 - tmp83;
		    tmp237 = imag_input[7 * imag_istride];
		    tmp238 = imag_input[57 * imag_istride];
		    tmp239 = tmp237 - tmp238;
		    tmp612 = tmp237 + tmp238;
	       }
	       {
		    fftw_real tmp85;
		    fftw_real tmp86;
		    fftw_real tmp234;
		    fftw_real tmp235;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp85 = real_input[25 * real_istride];
		    tmp86 = real_input[39 * real_istride];
		    tmp87 = tmp85 + tmp86;
		    tmp611 = tmp86 - tmp85;
		    tmp234 = imag_input[25 * imag_istride];
		    tmp235 = imag_input[39 * imag_istride];
		    tmp236 = tmp234 - tmp235;
		    tmp609 = tmp234 + tmp235;
	       }
	       tmp88 = tmp84 + tmp87;
	       tmp224 = tmp84 - tmp87;
	       tmp240 = tmp236 + tmp239;
	       tmp402 = tmp239 - tmp236;
	       tmp610 = tmp608 - tmp609;
	       tmp811 = tmp608 + tmp609;
	       tmp613 = tmp611 + tmp612;
	       tmp812 = tmp612 - tmp611;
	  }
	  {
	       fftw_real tmp813;
	       fftw_real tmp816;
	       fftw_real tmp592;
	       fftw_real tmp599;
	       ASSERT_ALIGNED_DOUBLE;
	       {
		    fftw_real tmp469;
		    fftw_real tmp470;
		    fftw_real tmp607;
		    fftw_real tmp614;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp81 = tmp73 + tmp80;
		    tmp469 = tmp73 - tmp80;
		    tmp470 = tmp403 + tmp402;
		    tmp471 = tmp469 - tmp470;
		    tmp503 = tmp469 + tmp470;
		    tmp404 = tmp402 - tmp403;
		    tmp607 = (K903989293 * tmp603) - (K427555093 * tmp606);
		    tmp614 = (K941544065 * tmp610) - (K336889853 * tmp613);
		    tmp615 = tmp607 + tmp614;
		    tmp724 = tmp607 - tmp614;
	       }
	       {
		    fftw_real tmp653;
		    fftw_real tmp654;
		    fftw_real tmp839;
		    fftw_real tmp840;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp653 = (K336889853 * tmp610) + (K941544065 * tmp613);
		    tmp654 = (K427555093 * tmp603) + (K903989293 * tmp606);
		    tmp655 = tmp653 - tmp654;
		    tmp721 = tmp654 + tmp653;
		    tmp839 = (K242980179 * tmp814) - (K970031253 * tmp815);
		    tmp840 = (K514102744 * tmp811) - (K857728610 * tmp812);
		    tmp841 = tmp839 + tmp840;
		    tmp910 = tmp839 - tmp840;
	       }
	       tmp813 = (K857728610 * tmp811) + (K514102744 * tmp812);
	       tmp816 = (K970031253 * tmp814) + (K242980179 * tmp815);
	       tmp817 = tmp813 - tmp816;
	       tmp907 = tmp816 + tmp813;
	       {
		    fftw_real tmp806;
		    fftw_real tmp809;
		    fftw_real tmp836;
		    fftw_real tmp837;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp806 = (K146730474 * tmp804) + (K989176509 * tmp805);
		    tmp809 = (K803207531 * tmp807) - (K595699304 * tmp808);
		    tmp810 = tmp806 - tmp809;
		    tmp909 = tmp809 + tmp806;
		    tmp836 = (K989176509 * tmp804) - (K146730474 * tmp805);
		    tmp837 = (K595699304 * tmp807) + (K803207531 * tmp808);
		    tmp838 = tmp836 - tmp837;
		    tmp906 = tmp837 + tmp836;
	       }
	       {
		    fftw_real tmp472;
		    fftw_real tmp473;
		    fftw_real tmp213;
		    fftw_real tmp222;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp96 = tmp88 + tmp95;
		    tmp472 = tmp95 - tmp88;
		    tmp473 = tmp400 + tmp399;
		    tmp474 = tmp472 + tmp473;
		    tmp504 = tmp473 - tmp472;
		    tmp401 = tmp399 - tmp400;
		    tmp213 = tmp205 - tmp212;
		    tmp222 = tmp214 + tmp221;
		    tmp223 = (K995184726 * tmp213) - (K098017140 * tmp222);
		    tmp285 = (K098017140 * tmp213) + (K995184726 * tmp222);
	       }
	       {
		    fftw_real tmp341;
		    fftw_real tmp342;
		    fftw_real tmp650;
		    fftw_real tmp651;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp341 = tmp205 + tmp212;
		    tmp342 = tmp221 - tmp214;
		    tmp343 = (K956940335 * tmp341) - (K290284677 * tmp342);
		    tmp357 = (K290284677 * tmp341) + (K956940335 * tmp342);
		    tmp650 = (K049067674 * tmp595) + (K998795456 * tmp598);
		    tmp651 = (K671558954 * tmp588) + (K740951125 * tmp591);
		    tmp652 = tmp650 - tmp651;
		    tmp723 = tmp651 + tmp650;
	       }
	       tmp592 = (K740951125 * tmp588) - (K671558954 * tmp591);
	       tmp599 = (K998795456 * tmp595) - (K049067674 * tmp598);
	       tmp600 = tmp592 + tmp599;
	       tmp720 = tmp599 - tmp592;
	       {
		    fftw_real tmp232;
		    fftw_real tmp241;
		    fftw_real tmp344;
		    fftw_real tmp345;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp232 = tmp224 - tmp231;
		    tmp241 = tmp233 + tmp240;
		    tmp242 = (K773010453 * tmp232) - (K634393284 * tmp241);
		    tmp286 = (K634393284 * tmp232) + (K773010453 * tmp241);
		    tmp344 = tmp224 + tmp231;
		    tmp345 = tmp240 - tmp233;
		    tmp346 = (K471396736 * tmp344) + (K881921264 * tmp345);
		    tmp358 = (K881921264 * tmp344) - (K471396736 * tmp345);
	       }
	  }
     }
     {
	  fftw_real tmp65;
	  fftw_real tmp421;
	  fftw_real tmp34;
	  fftw_real tmp414;
	  fftw_real tmp129;
	  fftw_real tmp423;
	  fftw_real tmp413;
	  fftw_real tmp424;
	  fftw_real tmp18;
	  fftw_real tmp66;
	  fftw_real tmp398;
	  ASSERT_ALIGNED_DOUBLE;
	  tmp65 = K2_000000000 * (tmp49 + tmp64);
	  tmp421 = K2_000000000 * (tmp417 - tmp420);
	  tmp18 = tmp10 + tmp17;
	  tmp34 = tmp18 + tmp33;
	  tmp414 = tmp18 - tmp33;
	  {
	       fftw_real tmp97;
	       fftw_real tmp128;
	       fftw_real tmp405;
	       fftw_real tmp412;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp97 = tmp81 + tmp96;
	       tmp128 = tmp112 + tmp127;
	       tmp129 = K2_000000000 * (tmp97 + tmp128);
	       tmp423 = tmp97 - tmp128;
	       tmp405 = tmp401 - tmp404;
	       tmp412 = tmp408 - tmp411;
	       tmp413 = K2_000000000 * (tmp405 - tmp412);
	       tmp424 = tmp405 + tmp412;
	  }
	  tmp66 = tmp34 + tmp65;
	  output[64 * ostride] = tmp66 - tmp129;
	  output[0] = tmp66 + tmp129;
	  tmp398 = tmp34 - tmp65;
	  output[32 * ostride] = tmp398 - tmp413;
	  output[96 * ostride] = tmp398 + tmp413;
	  {
	       fftw_real tmp422;
	       fftw_real tmp425;
	       fftw_real tmp426;
	       fftw_real tmp427;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp422 = tmp414 - tmp421;
	       tmp425 = K1_414213562 * (tmp423 - tmp424);
	       output[80 * ostride] = tmp422 - tmp425;
	       output[16 * ostride] = tmp422 + tmp425;
	       tmp426 = tmp414 + tmp421;
	       tmp427 = K1_414213562 * (tmp423 + tmp424);
	       output[48 * ostride] = tmp426 - tmp427;
	       output[112 * ostride] = tmp426 + tmp427;
	  }
     }
     {
	  fftw_real tmp432;
	  fftw_real tmp446;
	  fftw_real tmp442;
	  fftw_real tmp450;
	  fftw_real tmp435;
	  fftw_real tmp447;
	  fftw_real tmp439;
	  fftw_real tmp449;
	  fftw_real tmp428;
	  fftw_real tmp440;
	  fftw_real tmp441;
	  ASSERT_ALIGNED_DOUBLE;
	  tmp428 = tmp10 - tmp17;
	  tmp432 = tmp428 - tmp431;
	  tmp446 = tmp428 + tmp431;
	  tmp440 = tmp127 - tmp112;
	  tmp441 = tmp401 + tmp404;
	  tmp442 = tmp440 + tmp441;
	  tmp450 = tmp441 - tmp440;
	  {
	       fftw_real tmp433;
	       fftw_real tmp434;
	       fftw_real tmp437;
	       fftw_real tmp438;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp433 = tmp49 - tmp64;
	       tmp434 = tmp417 + tmp420;
	       tmp435 = K1_414213562 * (tmp433 - tmp434);
	       tmp447 = K1_414213562 * (tmp433 + tmp434);
	       tmp437 = tmp81 - tmp96;
	       tmp438 = tmp408 + tmp411;
	       tmp439 = tmp437 - tmp438;
	       tmp449 = tmp437 + tmp438;
	  }
	  {
	       fftw_real tmp436;
	       fftw_real tmp443;
	       fftw_real tmp444;
	       fftw_real tmp445;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp436 = tmp432 + tmp435;
	       tmp443 = (K1_847759065 * tmp439) - (K765366864 * tmp442);
	       output[72 * ostride] = tmp436 - tmp443;
	       output[8 * ostride] = tmp436 + tmp443;
	       tmp444 = tmp432 - tmp435;
	       tmp445 = (K1_847759065 * tmp442) + (K765366864 * tmp439);
	       output[40 * ostride] = tmp444 - tmp445;
	       output[104 * ostride] = tmp444 + tmp445;
	  }
	  {
	       fftw_real tmp448;
	       fftw_real tmp451;
	       fftw_real tmp452;
	       fftw_real tmp453;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp448 = tmp446 - tmp447;
	       tmp451 = (K765366864 * tmp449) - (K1_847759065 * tmp450);
	       output[88 * ostride] = tmp448 - tmp451;
	       output[24 * ostride] = tmp448 + tmp451;
	       tmp452 = tmp446 + tmp447;
	       tmp453 = (K765366864 * tmp450) + (K1_847759065 * tmp449);
	       output[56 * ostride] = tmp452 - tmp453;
	       output[120 * ostride] = tmp452 + tmp453;
	  }
     }
     {
	  fftw_real tmp203;
	  fftw_real tmp295;
	  fftw_real tmp164;
	  fftw_real tmp292;
	  fftw_real tmp283;
	  fftw_real tmp297;
	  fftw_real tmp291;
	  fftw_real tmp298;
	  fftw_real tmp144;
	  fftw_real tmp163;
	  fftw_real tmp204;
	  fftw_real tmp284;
	  ASSERT_ALIGNED_DOUBLE;
	  tmp203 = K2_000000000 * (tmp183 + tmp202);
	  tmp295 = K2_000000000 * (tmp293 - tmp294);
	  tmp144 = tmp134 + tmp143;
	  tmp163 = (K1_847759065 * tmp153) - (K765366864 * tmp162);
	  tmp164 = tmp144 + tmp163;
	  tmp292 = tmp144 - tmp163;
	  {
	       fftw_real tmp243;
	       fftw_real tmp282;
	       fftw_real tmp287;
	       fftw_real tmp290;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp243 = tmp223 + tmp242;
	       tmp282 = tmp262 + tmp281;
	       tmp283 = K2_000000000 * (tmp243 + tmp282);
	       tmp297 = tmp243 - tmp282;
	       tmp287 = tmp285 - tmp286;
	       tmp290 = tmp288 - tmp289;
	       tmp291 = K2_000000000 * (tmp287 - tmp290);
	       tmp298 = tmp287 + tmp290;
	  }
	  tmp204 = tmp164 + tmp203;
	  output[66 * ostride] = tmp204 - tmp283;
	  output[2 * ostride] = tmp204 + tmp283;
	  tmp284 = tmp164 - tmp203;
	  output[34 * ostride] = tmp284 - tmp291;
	  output[98 * ostride] = tmp284 + tmp291;
	  {
	       fftw_real tmp296;
	       fftw_real tmp299;
	       fftw_real tmp300;
	       fftw_real tmp301;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp296 = tmp292 - tmp295;
	       tmp299 = K1_414213562 * (tmp297 - tmp298);
	       output[82 * ostride] = tmp296 - tmp299;
	       output[18 * ostride] = tmp296 + tmp299;
	       tmp300 = tmp292 + tmp295;
	       tmp301 = K1_414213562 * (tmp297 + tmp298);
	       output[50 * ostride] = tmp300 - tmp301;
	       output[114 * ostride] = tmp300 + tmp301;
	  }
     }
     {
	  fftw_real tmp304;
	  fftw_real tmp318;
	  fftw_real tmp314;
	  fftw_real tmp322;
	  fftw_real tmp307;
	  fftw_real tmp319;
	  fftw_real tmp311;
	  fftw_real tmp321;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp302;
	       fftw_real tmp303;
	       fftw_real tmp312;
	       fftw_real tmp313;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp302 = tmp134 - tmp143;
	       tmp303 = (K765366864 * tmp153) + (K1_847759065 * tmp162);
	       tmp304 = tmp302 - tmp303;
	       tmp318 = tmp302 + tmp303;
	       tmp312 = tmp285 + tmp286;
	       tmp313 = tmp281 - tmp262;
	       tmp314 = tmp312 + tmp313;
	       tmp322 = tmp312 - tmp313;
	  }
	  {
	       fftw_real tmp305;
	       fftw_real tmp306;
	       fftw_real tmp309;
	       fftw_real tmp310;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp305 = tmp183 - tmp202;
	       tmp306 = tmp293 + tmp294;
	       tmp307 = K1_414213562 * (tmp305 - tmp306);
	       tmp319 = K1_414213562 * (tmp305 + tmp306);
	       tmp309 = tmp223 - tmp242;
	       tmp310 = tmp288 + tmp289;
	       tmp311 = tmp309 - tmp310;
	       tmp321 = tmp309 + tmp310;
	  }
	  {
	       fftw_real tmp308;