frc_128.c

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

     {
	  fftw_real tmp283;
	  fftw_real tmp671;
	  fftw_real tmp347;
	  fftw_real tmp672;
	  fftw_real tmp290;
	  fftw_real tmp344;
	  fftw_real tmp674;
	  fftw_real tmp675;
	  fftw_real tmp697;
	  fftw_real tmp698;
	  fftw_real tmp298;
	  fftw_real tmp350;
	  fftw_real tmp838;
	  fftw_real tmp694;
	  fftw_real tmp695;
	  fftw_real tmp305;
	  fftw_real tmp349;
	  fftw_real tmp839;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp281;
	       fftw_real tmp282;
	       fftw_real tmp345;
	       fftw_real tmp346;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp281 = input[127 * istride];
	       tmp282 = input[63 * istride];
	       tmp283 = tmp281 - tmp282;
	       tmp671 = tmp281 + tmp282;
	       tmp345 = input[31 * istride];
	       tmp346 = input[95 * istride];
	       tmp347 = tmp345 - tmp346;
	       tmp672 = tmp345 + tmp346;
	  }
	  {
	       fftw_real tmp284;
	       fftw_real tmp285;
	       fftw_real tmp286;
	       fftw_real tmp287;
	       fftw_real tmp288;
	       fftw_real tmp289;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp284 = input[15 * istride];
	       tmp285 = input[79 * istride];
	       tmp286 = tmp284 - tmp285;
	       tmp287 = input[111 * istride];
	       tmp288 = input[47 * istride];
	       tmp289 = tmp287 - tmp288;
	       tmp290 = K707106781 * (tmp286 + tmp289);
	       tmp344 = K707106781 * (tmp289 - tmp286);
	       tmp674 = tmp284 + tmp285;
	       tmp675 = tmp287 + tmp288;
	  }
	  {
	       fftw_real tmp294;
	       fftw_real tmp297;
	       fftw_real tmp301;
	       fftw_real tmp304;
	       ASSERT_ALIGNED_DOUBLE;
	       {
		    fftw_real tmp292;
		    fftw_real tmp293;
		    fftw_real tmp295;
		    fftw_real tmp296;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp292 = input[7 * istride];
		    tmp293 = input[71 * istride];
		    tmp294 = tmp292 - tmp293;
		    tmp697 = tmp292 + tmp293;
		    tmp295 = input[39 * istride];
		    tmp296 = input[103 * istride];
		    tmp297 = tmp295 - tmp296;
		    tmp698 = tmp295 + tmp296;
	       }
	       tmp298 = (K923879532 * tmp294) - (K382683432 * tmp297);
	       tmp350 = (K382683432 * tmp294) + (K923879532 * tmp297);
	       tmp838 = tmp697 - tmp698;
	       {
		    fftw_real tmp299;
		    fftw_real tmp300;
		    fftw_real tmp302;
		    fftw_real tmp303;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp299 = input[119 * istride];
		    tmp300 = input[55 * istride];
		    tmp301 = tmp299 - tmp300;
		    tmp694 = tmp299 + tmp300;
		    tmp302 = input[23 * istride];
		    tmp303 = input[87 * istride];
		    tmp304 = tmp302 - tmp303;
		    tmp695 = tmp302 + tmp303;
	       }
	       tmp305 = (K923879532 * tmp301) + (K382683432 * tmp304);
	       tmp349 = (K382683432 * tmp301) - (K923879532 * tmp304);
	       tmp839 = tmp694 - tmp695;
	  }
	  {
	       fftw_real tmp291;
	       fftw_real tmp306;
	       fftw_real tmp507;
	       fftw_real tmp508;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp291 = tmp283 + tmp290;
	       tmp306 = tmp298 + tmp305;
	       tmp307 = tmp291 + tmp306;
	       tmp419 = tmp291 - tmp306;
	       tmp507 = tmp283 - tmp290;
	       tmp508 = tmp350 + tmp349;
	       tmp509 = tmp507 + tmp508;
	       tmp583 = tmp507 - tmp508;
	  }
	  {
	       fftw_real tmp348;
	       fftw_real tmp351;
	       fftw_real tmp518;
	       fftw_real tmp519;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp348 = tmp344 - tmp347;
	       tmp351 = tmp349 - tmp350;
	       tmp352 = tmp348 + tmp351;
	       tmp423 = tmp351 - tmp348;
	       tmp518 = tmp347 + tmp344;
	       tmp519 = tmp305 - tmp298;
	       tmp520 = tmp518 + tmp519;
	       tmp587 = tmp519 - tmp518;
	  }
	  {
	       fftw_real tmp673;
	       fftw_real tmp676;
	       fftw_real tmp837;
	       fftw_real tmp840;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp673 = tmp671 + tmp672;
	       tmp676 = tmp674 + tmp675;
	       tmp677 = tmp673 - tmp676;
	       tmp747 = tmp673 + tmp676;
	       tmp837 = tmp671 - tmp672;
	       tmp840 = K707106781 * (tmp838 + tmp839);
	       tmp841 = tmp837 + tmp840;
	       tmp915 = tmp837 - tmp840;
	  }
	  {
	       fftw_real tmp850;
	       fftw_real tmp851;
	       fftw_real tmp696;
	       fftw_real tmp699;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp850 = tmp675 - tmp674;
	       tmp851 = K707106781 * (tmp839 - tmp838);
	       tmp852 = tmp850 + tmp851;
	       tmp919 = tmp851 - tmp850;
	       tmp696 = tmp694 + tmp695;
	       tmp699 = tmp697 + tmp698;
	       tmp700 = tmp696 - tmp699;
	       tmp748 = tmp699 + tmp696;
	  }
     }
     {
	  fftw_real tmp310;
	  fftw_real tmp842;
	  fftw_real tmp680;
	  fftw_real tmp322;
	  fftw_real tmp334;
	  fftw_real tmp336;
	  fftw_real tmp690;
	  fftw_real tmp846;
	  fftw_real tmp327;
	  fftw_real tmp845;
	  fftw_real tmp687;
	  fftw_real tmp339;
	  fftw_real tmp317;
	  fftw_real tmp319;
	  fftw_real tmp683;
	  fftw_real tmp843;
	  fftw_real tmp510;
	  fftw_real tmp511;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp308;
	       fftw_real tmp309;
	       fftw_real tmp678;
	       fftw_real tmp320;
	       fftw_real tmp321;
	       fftw_real tmp679;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp308 = input[3 * istride];
	       tmp309 = input[67 * istride];
	       tmp678 = tmp308 + tmp309;
	       tmp320 = input[35 * istride];
	       tmp321 = input[99 * istride];
	       tmp679 = tmp320 + tmp321;
	       tmp310 = tmp308 - tmp309;
	       tmp842 = tmp678 - tmp679;
	       tmp680 = tmp678 + tmp679;
	       tmp322 = tmp320 - tmp321;
	  }
	  {
	       fftw_real tmp330;
	       fftw_real tmp688;
	       fftw_real tmp333;
	       fftw_real tmp689;
	       ASSERT_ALIGNED_DOUBLE;
	       {
		    fftw_real tmp328;
		    fftw_real tmp329;
		    fftw_real tmp331;
		    fftw_real tmp332;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp328 = input[11 * istride];
		    tmp329 = input[75 * istride];
		    tmp330 = tmp328 - tmp329;
		    tmp688 = tmp328 + tmp329;
		    tmp331 = input[107 * istride];
		    tmp332 = input[43 * istride];
		    tmp333 = tmp331 - tmp332;
		    tmp689 = tmp331 + tmp332;
	       }
	       tmp334 = K707106781 * (tmp330 + tmp333);
	       tmp336 = K707106781 * (tmp333 - tmp330);
	       tmp690 = tmp688 + tmp689;
	       tmp846 = tmp689 - tmp688;
	  }
	  {
	       fftw_real tmp325;
	       fftw_real tmp326;
	       fftw_real tmp685;
	       fftw_real tmp337;
	       fftw_real tmp338;
	       fftw_real tmp686;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp325 = input[123 * istride];
	       tmp326 = input[59 * istride];
	       tmp685 = tmp325 + tmp326;
	       tmp337 = input[27 * istride];
	       tmp338 = input[91 * istride];
	       tmp686 = tmp337 + tmp338;
	       tmp327 = tmp325 - tmp326;
	       tmp845 = tmp685 - tmp686;
	       tmp687 = tmp685 + tmp686;
	       tmp339 = tmp337 - tmp338;
	  }
	  {
	       fftw_real tmp313;
	       fftw_real tmp681;
	       fftw_real tmp316;
	       fftw_real tmp682;
	       ASSERT_ALIGNED_DOUBLE;
	       {
		    fftw_real tmp311;
		    fftw_real tmp312;
		    fftw_real tmp314;
		    fftw_real tmp315;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp311 = input[19 * istride];
		    tmp312 = input[83 * istride];
		    tmp313 = tmp311 - tmp312;
		    tmp681 = tmp311 + tmp312;
		    tmp314 = input[115 * istride];
		    tmp315 = input[51 * istride];
		    tmp316 = tmp314 - tmp315;
		    tmp682 = tmp314 + tmp315;
	       }
	       tmp317 = K707106781 * (tmp313 + tmp316);
	       tmp319 = K707106781 * (tmp316 - tmp313);
	       tmp683 = tmp681 + tmp682;
	       tmp843 = tmp682 - tmp681;
	  }
	  {
	       fftw_real tmp684;
	       fftw_real tmp691;
	       fftw_real tmp853;
	       fftw_real tmp854;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp750 = tmp680 + tmp683;
	       tmp684 = tmp680 - tmp683;
	       tmp691 = tmp687 - tmp690;
	       tmp692 = K707106781 * (tmp684 + tmp691);
	       tmp701 = K707106781 * (tmp691 - tmp684);
	       tmp751 = tmp687 + tmp690;
	       tmp853 = (K923879532 * tmp843) - (K382683432 * tmp842);
	       tmp854 = (K382683432 * tmp845) + (K923879532 * tmp846);
	       tmp855 = tmp853 + tmp854;
	       tmp916 = tmp854 - tmp853;
	  }
	  {
	       fftw_real tmp844;
	       fftw_real tmp847;
	       fftw_real tmp318;
	       fftw_real tmp323;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp844 = (K923879532 * tmp842) + (K382683432 * tmp843);
	       tmp847 = (K923879532 * tmp845) - (K382683432 * tmp846);
	       tmp848 = tmp844 + tmp847;
	       tmp918 = tmp847 - tmp844;
	       tmp318 = tmp310 + tmp317;
	       tmp323 = tmp319 - tmp322;
	       tmp324 = (K980785280 * tmp318) + (K195090322 * tmp323);
	       tmp353 = (K980785280 * tmp323) - (K195090322 * tmp318);
	  }
	  tmp510 = tmp310 - tmp317;
	  tmp511 = tmp322 + tmp319;
	  tmp512 = (K831469612 * tmp510) + (K555570233 * tmp511);
	  tmp521 = (K831469612 * tmp511) - (K555570233 * tmp510);
	  {
	       fftw_real tmp513;
	       fftw_real tmp514;
	       fftw_real tmp335;
	       fftw_real tmp340;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp513 = tmp327 - tmp334;
	       tmp514 = tmp339 + tmp336;
	       tmp515 = (K831469612 * tmp513) - (K555570233 * tmp514);
	       tmp522 = (K555570233 * tmp513) + (K831469612 * tmp514);
	       tmp335 = tmp327 + tmp334;
	       tmp340 = tmp336 - tmp339;
	       tmp341 = (K980785280 * tmp335) - (K195090322 * tmp340);
	       tmp354 = (K195090322 * tmp335) + (K980785280 * tmp340);
	  }
     }
     {
	  fftw_real tmp49;
	  fftw_real tmp149;
	  fftw_real tmp52;
	  fftw_real tmp155;
	  fftw_real tmp56;
	  fftw_real tmp154;
	  fftw_real tmp152;
	  fftw_real tmp59;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp47;
	       fftw_real tmp48;
	       fftw_real tmp50;
	       fftw_real tmp51;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp47 = input[124 * istride];
	       tmp48 = input[60 * istride];
	       tmp49 = tmp47 + tmp48;
	       tmp149 = tmp47 - tmp48;
	       tmp50 = input[28 * istride];
	       tmp51 = input[92 * istride];
	       tmp52 = tmp50 + tmp51;
	       tmp155 = tmp50 - tmp51;
	       {
		    fftw_real tmp54;
		    fftw_real tmp55;
		    fftw_real tmp150;
		    fftw_real tmp57;
		    fftw_real tmp58;
		    fftw_real tmp151;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp54 = input[12 * istride];
		    tmp55 = input[76 * istride];
		    tmp150 = tmp54 - tmp55;
		    tmp57 = input[108 * istride];
		    tmp58 = input[44 * istride];
		    tmp151 = tmp57 - tmp58;
		    tmp56 = tmp54 + tmp55;
		    tmp154 = K707106781 * (tmp151 - tmp150);
		    tmp152 = K707106781 * (tmp150 + tmp151);
		    tmp59 = tmp57 + tmp58;
	       }
	  }
	  {
	       fftw_real tmp53;
	       fftw_real tmp60;
	       fftw_real tmp791;
	       fftw_real tmp792;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp53 = tmp49 + tmp52;
	       tmp60 = tmp56 + tmp59;
	       tmp61 = tmp53 + tmp60;
	       tmp627 = tmp53 - tmp60;
	       tmp791 = tmp49 - tmp52;
	       tmp792 = tmp59 - tmp56;
	       tmp793 = (K923879532 * tmp791) - (K382683432 * tmp792);
	       tmp866 = (K382683432 * tmp791) + (K923879532 * tmp792);
	  }
	  {
	       fftw_real tmp153;
	       fftw_real tmp156;
	       fftw_real tmp467;
	       fftw_real tmp468;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp153 = tmp149 + tmp152;
	       tmp156 = tmp154 - tmp155;
	       tmp157 = (K980785280 * tmp153) - (K195090322 * tmp156);
	       tmp370 = (K195090322 * tmp153) + (K980785280 * tmp156);
	       tmp467 = tmp149 - tmp152;
	       tmp468 = tmp155 + tmp154;
	       tmp469 = (K831469612 * tmp467) - (K555570233 * tmp468);
	       tmp534 = (K555570233 * tmp467) + (K831469612 * tmp468);
	  }
     }
     {
	  fftw_real tmp97;
	  fftw_real tmp181;
	  fftw_real tmp119;
	  fftw_real tmp122;
	  fftw_real tmp191;
	  fftw_real tmp197;
	  fftw_real tmp807;
	  fftw_real tmp100;
	  fftw_real tmp195;
	  fftw_real tmp104;
	  fftw_real tmp194;
	  fftw_real tmp184;
	  fftw_real tmp107;
	  fftw_real tmp112;
	  fftw_real tmp115;
	  fftw_real tmp188;
	  fftw_real tmp198;
	  fftw_real tmp806;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp95;
	       fftw_real tmp96;
	       fftw_real tmp189;
	       fftw_real tmp190;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp95 = input[126 * istride];
	       tmp96 = input[62 * istride];
	       tmp97 = tmp95 + tmp96;
	       tmp181 = tmp95 - tmp96;
	       {
		    fftw_real tmp117;
		    fftw_real tmp118;
		    fftw_real tmp120;
		    fftw_real tmp121;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp117 = input[118 * istride];
		    tmp118 = input[54 * istride];
		    tmp119 = tmp117 + tmp118;
		    tmp189 = tmp117 - tmp118;
		    tmp120 = input[22 * istride];
		    tmp121 = input[86 * istride];
		    tmp122 = tmp120 + tmp121;
		    tmp190 = tmp120 - tmp121;
	       }
	       tmp191 = (K923879532 * tmp189) + (K382683432 * tmp190);