fcr_128.c

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

		    tmp103 = tmp101 + tmp102;
		    tmp627 = tmp102 - tmp101;
		    tmp254 = imag_input[29 * imag_istride];
		    tmp255 = imag_input[35 * imag_istride];
		    tmp256 = tmp254 - tmp255;
		    tmp625 = tmp254 + tmp255;
	       }
	       tmp104 = tmp100 + tmp103;
	       tmp244 = tmp100 - tmp103;
	       tmp260 = tmp256 + tmp259;
	       tmp406 = tmp259 - tmp256;
	       tmp626 = tmp624 - tmp625;
	       tmp819 = tmp624 + tmp625;
	       tmp629 = tmp627 + tmp628;
	       tmp820 = tmp628 - tmp627;
	  }
	  {
	       fftw_real tmp122;
	       fftw_real tmp632;
	       fftw_real tmp269;
	       fftw_real tmp636;
	       fftw_real tmp125;
	       fftw_real tmp635;
	       fftw_real tmp266;
	       fftw_real tmp633;
	       ASSERT_ALIGNED_DOUBLE;
	       {
		    fftw_real tmp120;
		    fftw_real tmp121;
		    fftw_real tmp267;
		    fftw_real tmp268;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp120 = real_input[11 * real_istride];
		    tmp121 = real_input[53 * real_istride];
		    tmp122 = tmp120 + tmp121;
		    tmp632 = tmp120 - tmp121;
		    tmp267 = imag_input[11 * imag_istride];
		    tmp268 = imag_input[53 * imag_istride];
		    tmp269 = tmp267 - tmp268;
		    tmp636 = tmp267 + tmp268;
	       }
	       {
		    fftw_real tmp123;
		    fftw_real tmp124;
		    fftw_real tmp264;
		    fftw_real tmp265;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp123 = real_input[21 * real_istride];
		    tmp124 = real_input[43 * real_istride];
		    tmp125 = tmp123 + tmp124;
		    tmp635 = tmp124 - tmp123;
		    tmp264 = imag_input[21 * imag_istride];
		    tmp265 = imag_input[43 * imag_istride];
		    tmp266 = tmp264 - tmp265;
		    tmp633 = tmp264 + tmp265;
	       }
	       tmp126 = tmp122 + tmp125;
	       tmp272 = tmp125 - tmp122;
	       tmp270 = tmp266 + tmp269;
	       tmp410 = tmp269 - tmp266;
	       tmp634 = tmp632 - tmp633;
	       tmp829 = tmp632 + tmp633;
	       tmp637 = tmp635 + tmp636;
	       tmp830 = tmp636 - tmp635;
	  }
	  {
	       fftw_real tmp107;
	       fftw_real tmp617;
	       fftw_real tmp250;
	       fftw_real tmp621;
	       fftw_real tmp110;
	       fftw_real tmp620;
	       fftw_real tmp247;
	       fftw_real tmp618;
	       ASSERT_ALIGNED_DOUBLE;
	       {
		    fftw_real tmp105;
		    fftw_real tmp106;
		    fftw_real tmp248;
		    fftw_real tmp249;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp105 = real_input[13 * real_istride];
		    tmp106 = real_input[51 * real_istride];
		    tmp107 = tmp105 + tmp106;
		    tmp617 = tmp105 - tmp106;
		    tmp248 = imag_input[13 * imag_istride];
		    tmp249 = imag_input[51 * imag_istride];
		    tmp250 = tmp248 - tmp249;
		    tmp621 = tmp248 + tmp249;
	       }
	       {
		    fftw_real tmp108;
		    fftw_real tmp109;
		    fftw_real tmp245;
		    fftw_real tmp246;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp108 = real_input[19 * real_istride];
		    tmp109 = real_input[45 * real_istride];
		    tmp110 = tmp108 + tmp109;
		    tmp620 = tmp109 - tmp108;
		    tmp245 = imag_input[19 * imag_istride];
		    tmp246 = imag_input[45 * imag_istride];
		    tmp247 = tmp245 - tmp246;
		    tmp618 = tmp245 + tmp246;
	       }
	       tmp111 = tmp107 + tmp110;
	       tmp253 = tmp110 - tmp107;
	       tmp251 = tmp247 + tmp250;
	       tmp407 = tmp250 - tmp247;
	       tmp619 = tmp617 - tmp618;
	       tmp822 = tmp617 + tmp618;
	       tmp622 = tmp620 + tmp621;
	       tmp823 = tmp621 - tmp620;
	  }
	  {
	       fftw_real tmp115;
	       fftw_real tmp639;
	       fftw_real tmp278;
	       fftw_real tmp643;
	       fftw_real tmp118;
	       fftw_real tmp642;
	       fftw_real tmp275;
	       fftw_real tmp640;
	       ASSERT_ALIGNED_DOUBLE;
	       {
		    fftw_real tmp113;
		    fftw_real tmp114;
		    fftw_real tmp276;
		    fftw_real tmp277;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp113 = real_input[5 * real_istride];
		    tmp114 = real_input[59 * real_istride];
		    tmp115 = tmp113 + tmp114;
		    tmp639 = tmp113 - tmp114;
		    tmp276 = imag_input[5 * imag_istride];
		    tmp277 = imag_input[59 * imag_istride];
		    tmp278 = tmp276 - tmp277;
		    tmp643 = tmp276 + tmp277;
	       }
	       {
		    fftw_real tmp116;
		    fftw_real tmp117;
		    fftw_real tmp273;
		    fftw_real tmp274;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp116 = real_input[27 * real_istride];
		    tmp117 = real_input[37 * real_istride];
		    tmp118 = tmp116 + tmp117;
		    tmp642 = tmp117 - tmp116;
		    tmp273 = imag_input[27 * imag_istride];
		    tmp274 = imag_input[37 * imag_istride];
		    tmp275 = tmp273 - tmp274;
		    tmp640 = tmp273 + tmp274;
	       }
	       tmp119 = tmp115 + tmp118;
	       tmp263 = tmp115 - tmp118;
	       tmp279 = tmp275 + tmp278;
	       tmp409 = tmp278 - tmp275;
	       tmp641 = tmp639 - tmp640;
	       tmp826 = tmp639 + tmp640;
	       tmp644 = tmp642 + tmp643;
	       tmp827 = tmp643 - tmp642;
	  }
	  {
	       fftw_real tmp828;
	       fftw_real tmp831;
	       fftw_real tmp623;
	       fftw_real tmp630;
	       ASSERT_ALIGNED_DOUBLE;
	       {
		    fftw_real tmp476;
		    fftw_real tmp477;
		    fftw_real tmp638;
		    fftw_real tmp645;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp112 = tmp104 + tmp111;
		    tmp476 = tmp104 - tmp111;
		    tmp477 = tmp410 + tmp409;
		    tmp478 = tmp476 - tmp477;
		    tmp507 = tmp476 + tmp477;
		    tmp411 = tmp409 - tmp410;
		    tmp638 = (K857728610 * tmp634) - (K514102744 * tmp637);
		    tmp645 = (K970031253 * tmp641) - (K242980179 * tmp644);
		    tmp646 = tmp638 + tmp645;
		    tmp731 = tmp638 - tmp645;
	       }
	       {
		    fftw_real tmp660;
		    fftw_real tmp661;
		    fftw_real tmp846;
		    fftw_real tmp847;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp660 = (K242980179 * tmp641) + (K970031253 * tmp644);
		    tmp661 = (K514102744 * tmp634) + (K857728610 * tmp637);
		    tmp662 = tmp660 - tmp661;
		    tmp728 = tmp661 + tmp660;
		    tmp846 = (K740951125 * tmp826) - (K671558954 * tmp827);
		    tmp847 = (K049067674 * tmp829) + (K998795456 * tmp830);
		    tmp848 = tmp846 - tmp847;
		    tmp917 = tmp847 + tmp846;
	       }
	       tmp828 = (K671558954 * tmp826) + (K740951125 * tmp827);
	       tmp831 = (K998795456 * tmp829) - (K049067674 * tmp830);
	       tmp832 = tmp828 - tmp831;
	       tmp914 = tmp831 + tmp828;
	       {
		    fftw_real tmp821;
		    fftw_real tmp824;
		    fftw_real tmp843;
		    fftw_real tmp844;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp821 = (K427555093 * tmp819) + (K903989293 * tmp820);
		    tmp824 = (K941544065 * tmp822) - (K336889853 * tmp823);
		    tmp825 = tmp821 - tmp824;
		    tmp916 = tmp824 + tmp821;
		    tmp843 = (K903989293 * tmp819) - (K427555093 * tmp820);
		    tmp844 = (K336889853 * tmp822) + (K941544065 * tmp823);
		    tmp845 = tmp843 - tmp844;
		    tmp913 = tmp844 + tmp843;
	       }
	       {
		    fftw_real tmp479;
		    fftw_real tmp480;
		    fftw_real tmp252;
		    fftw_real tmp261;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp127 = tmp119 + tmp126;
		    tmp479 = tmp126 - tmp119;
		    tmp480 = tmp407 + tmp406;
		    tmp481 = tmp479 + tmp480;
		    tmp506 = tmp480 - tmp479;
		    tmp408 = tmp406 - tmp407;
		    tmp252 = tmp244 - tmp251;
		    tmp261 = tmp253 + tmp260;
		    tmp262 = (K956940335 * tmp252) - (K290284677 * tmp261);
		    tmp288 = (K290284677 * tmp252) + (K956940335 * tmp261);
	       }
	       {
		    fftw_real tmp348;
		    fftw_real tmp349;
		    fftw_real tmp657;
		    fftw_real tmp658;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp348 = tmp244 + tmp251;
		    tmp349 = tmp260 - tmp253;
		    tmp350 = (K634393284 * tmp348) - (K773010453 * tmp349);
		    tmp360 = (K773010453 * tmp348) + (K634393284 * tmp349);
		    tmp657 = (K146730474 * tmp626) + (K989176509 * tmp629);
		    tmp658 = (K595699304 * tmp619) + (K803207531 * tmp622);
		    tmp659 = tmp657 - tmp658;
		    tmp730 = tmp658 + tmp657;
	       }
	       tmp623 = (K803207531 * tmp619) - (K595699304 * tmp622);
	       tmp630 = (K989176509 * tmp626) - (K146730474 * tmp629);
	       tmp631 = tmp623 + tmp630;
	       tmp727 = tmp630 - tmp623;
	       {
		    fftw_real tmp271;
		    fftw_real tmp280;
		    fftw_real tmp351;
		    fftw_real tmp352;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp271 = tmp263 - tmp270;
		    tmp280 = tmp272 + tmp279;
		    tmp281 = (K881921264 * tmp271) - (K471396736 * tmp280);
		    tmp289 = (K471396736 * tmp271) + (K881921264 * tmp280);
		    tmp351 = tmp263 + tmp270;
		    tmp352 = tmp279 - tmp272;
		    tmp353 = (K098017140 * tmp351) - (K995184726 * tmp352);
		    tmp361 = (K995184726 * tmp351) + (K098017140 * tmp352);
	       }
	  }
     }
     {
	  fftw_real tmp21;
	  fftw_real tmp773;
	  fftw_real tmp562;
	  fftw_real tmp157;
	  fftw_real tmp24;
	  fftw_real tmp774;
	  fftw_real tmp565;
	  fftw_real tmp160;
	  fftw_real tmp28;
	  fftw_real tmp776;
	  fftw_real tmp555;
	  fftw_real tmp148;
	  fftw_real tmp31;
	  fftw_real tmp777;
	  fftw_real tmp558;
	  fftw_real tmp151;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp19;
	       fftw_real tmp20;
	       fftw_real tmp560;
	       fftw_real tmp155;
	       fftw_real tmp156;
	       fftw_real tmp561;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp19 = real_input[4 * real_istride];
	       tmp20 = real_input[60 * real_istride];
	       tmp560 = tmp19 - tmp20;
	       tmp155 = imag_input[28 * imag_istride];
	       tmp156 = imag_input[36 * imag_istride];
	       tmp561 = tmp155 + tmp156;
	       tmp21 = tmp19 + tmp20;
	       tmp773 = tmp560 + tmp561;
	       tmp562 = tmp560 - tmp561;
	       tmp157 = tmp155 - tmp156;
	  }
	  {
	       fftw_real tmp22;
	       fftw_real tmp23;
	       fftw_real tmp563;
	       fftw_real tmp158;
	       fftw_real tmp159;
	       fftw_real tmp564;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp22 = real_input[28 * real_istride];
	       tmp23 = real_input[36 * real_istride];
	       tmp563 = tmp23 - tmp22;
	       tmp158 = imag_input[4 * imag_istride];
	       tmp159 = imag_input[60 * imag_istride];
	       tmp564 = tmp158 + tmp159;
	       tmp24 = tmp22 + tmp23;
	       tmp774 = tmp564 - tmp563;
	       tmp565 = tmp563 + tmp564;
	       tmp160 = tmp158 - tmp159;
	  }
	  {
	       fftw_real tmp26;
	       fftw_real tmp27;
	       fftw_real tmp553;
	       fftw_real tmp146;
	       fftw_real tmp147;
	       fftw_real tmp554;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp26 = real_input[12 * real_istride];
	       tmp27 = real_input[52 * real_istride];
	       tmp553 = tmp26 - tmp27;
	       tmp146 = imag_input[20 * imag_istride];
	       tmp147 = imag_input[44 * imag_istride];
	       tmp554 = tmp146 + tmp147;
	       tmp28 = tmp26 + tmp27;
	       tmp776 = tmp553 + tmp554;
	       tmp555 = tmp553 - tmp554;
	       tmp148 = tmp146 - tmp147;
	  }
	  {
	       fftw_real tmp29;
	       fftw_real tmp30;
	       fftw_real tmp556;
	       fftw_real tmp149;
	       fftw_real tmp150;
	       fftw_real tmp557;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp29 = real_input[20 * real_istride];
	       tmp30 = real_input[44 * real_istride];
	       tmp556 = tmp30 - tmp29;
	       tmp149 = imag_input[12 * imag_istride];
	       tmp150 = imag_input[52 * imag_istride];
	       tmp557 = tmp149 + tmp150;
	       tmp31 = tmp29 + tmp30;
	       tmp777 = tmp557 - tmp556;
	       tmp558 = tmp556 + tmp557;
	       tmp151 = tmp149 - tmp150;
	  }
	  {
	       fftw_real tmp25;
	       fftw_real tmp32;
	       fftw_real tmp145;
	       fftw_real tmp152;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp25 = tmp21 + tmp24;
	       tmp32 = tmp28 + tmp31;
	       tmp33 = K2_000000000 * (tmp25 + tmp32);
	       tmp457 = tmp25 - tmp32;
	       tmp145 = tmp21 - tmp24;
	       tmp152 = tmp148 + tmp151;
	       tmp153 = tmp145 - tmp152;
	       tmp329 = tmp145 + tmp152;
	  }
	  {
	       fftw_real tmp429;
	       fftw_real tmp430;
	       fftw_real tmp154;
	       fftw_real tmp161;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp429 = tmp160 - tmp157;
	       tmp430 = tmp151 - tmp148;
	       tmp431 = K2_000000000 * (tmp429 - tmp430);
	       tmp458 = tmp430 + tmp429;
	       tmp154 = tmp31 - tmp28;
	       tmp161 = tmp157 + tmp160;
	       tmp162 = tmp154 + tmp161;
	       tmp330 = tmp161 - tmp154;
	  }
	  tmp559 = (K831469612 * tmp555) - (K555570233 * tmp558);
	  tmp566 = (K980785280 * tmp562) - (K195090322 * tmp565);
	  tmp712 = tmp566 - tmp559;
	  tmp683 = (K195090322 * tmp562) + (K980785280 * tmp565);
	  tmp684 = (K555570233 * tmp555) + (K831469612 * tmp558);
	  tmp713 = tmp684 + tmp683;
	  tmp775 = (K831469612 * tmp773) - (K555570233 * tmp774);
	  tmp778 = (K195090322 * tmp776) + (K980785280 * tmp777);
	  tmp898 = tmp778 + tmp775;
	  tmp865 = (K555570233 * tmp773) + (K831469612 * tmp774);
	  tmp866 = (K980785280 * tmp776) - (K195090322 * tmp777);
	  tmp899 = tmp866 + tmp865;
     }
     {
	  fftw_real tmp37;
	  fftw_real tmp180;
	  fftw_real tmp40;
	  fftw_real tmp177;
	  fftw_real tmp165;
	  fftw_real tmp534;
	  fftw_real tmp789;
	  fftw_real tmp788;
	  fftw_real tmp531;
	  fftw_real tmp181;
	  fftw_real tmp44;
	  fftw_real tmp171;
	  fftw_real tmp47;
	  fftw_real tmp168;
	  fftw_real tmp172;
	  fftw_real tmp527;
	  fftw_real tmp792;
	  fftw_real tmp791;
	  fftw_real tmp524;
	  fftw_real tmp174;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp529;
	       fftw_real tmp533;
	       fftw_real tmp532;