frc_128.c

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

	  fftw_real tmp240;
	  fftw_real tmp242;
	  fftw_real tmp650;
	  fftw_real tmp822;
	  fftw_real tmp491;
	  fftw_real tmp492;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp231;
	       fftw_real tmp232;
	       fftw_real tmp645;
	       fftw_real tmp243;
	       fftw_real tmp244;
	       fftw_real tmp646;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp231 = input[5 * istride];
	       tmp232 = input[69 * istride];
	       tmp645 = tmp231 + tmp232;
	       tmp243 = input[37 * istride];
	       tmp244 = input[101 * istride];
	       tmp646 = tmp243 + tmp244;
	       tmp233 = tmp231 - tmp232;
	       tmp821 = tmp645 - tmp646;
	       tmp647 = tmp645 + tmp646;
	       tmp245 = tmp243 - tmp244;
	  }
	  {
	       fftw_real tmp253;
	       fftw_real tmp655;
	       fftw_real tmp256;
	       fftw_real tmp656;
	       ASSERT_ALIGNED_DOUBLE;
	       {
		    fftw_real tmp251;
		    fftw_real tmp252;
		    fftw_real tmp254;
		    fftw_real tmp255;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp251 = input[13 * istride];
		    tmp252 = input[77 * istride];
		    tmp253 = tmp251 - tmp252;
		    tmp655 = tmp251 + tmp252;
		    tmp254 = input[109 * istride];
		    tmp255 = input[45 * istride];
		    tmp256 = tmp254 - tmp255;
		    tmp656 = tmp254 + tmp255;
	       }
	       tmp257 = K707106781 * (tmp253 + tmp256);
	       tmp259 = K707106781 * (tmp256 - tmp253);
	       tmp657 = tmp655 + tmp656;
	       tmp825 = tmp656 - tmp655;
	  }
	  {
	       fftw_real tmp248;
	       fftw_real tmp249;
	       fftw_real tmp652;
	       fftw_real tmp260;
	       fftw_real tmp261;
	       fftw_real tmp653;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp248 = input[125 * istride];
	       tmp249 = input[61 * istride];
	       tmp652 = tmp248 + tmp249;
	       tmp260 = input[29 * istride];
	       tmp261 = input[93 * istride];
	       tmp653 = tmp260 + tmp261;
	       tmp250 = tmp248 - tmp249;
	       tmp824 = tmp652 - tmp653;
	       tmp654 = tmp652 + tmp653;
	       tmp262 = tmp260 - tmp261;
	  }
	  {
	       fftw_real tmp236;
	       fftw_real tmp648;
	       fftw_real tmp239;
	       fftw_real tmp649;
	       ASSERT_ALIGNED_DOUBLE;
	       {
		    fftw_real tmp234;
		    fftw_real tmp235;
		    fftw_real tmp237;
		    fftw_real tmp238;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp234 = input[21 * istride];
		    tmp235 = input[85 * istride];
		    tmp236 = tmp234 - tmp235;
		    tmp648 = tmp234 + tmp235;
		    tmp237 = input[117 * istride];
		    tmp238 = input[53 * istride];
		    tmp239 = tmp237 - tmp238;
		    tmp649 = tmp237 + tmp238;
	       }
	       tmp240 = K707106781 * (tmp236 + tmp239);
	       tmp242 = K707106781 * (tmp239 - tmp236);
	       tmp650 = tmp648 + tmp649;
	       tmp822 = tmp649 - tmp648;
	  }
	  {
	       fftw_real tmp651;
	       fftw_real tmp658;
	       fftw_real tmp832;
	       fftw_real tmp833;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp743 = tmp647 + tmp650;
	       tmp651 = tmp647 - tmp650;
	       tmp658 = tmp654 - tmp657;
	       tmp659 = K707106781 * (tmp651 + tmp658);
	       tmp668 = K707106781 * (tmp658 - tmp651);
	       tmp744 = tmp654 + tmp657;
	       tmp832 = (K923879532 * tmp822) - (K382683432 * tmp821);
	       tmp833 = (K382683432 * tmp824) + (K923879532 * tmp825);
	       tmp834 = tmp832 + tmp833;
	       tmp912 = tmp833 - tmp832;
	  }
	  {
	       fftw_real tmp823;
	       fftw_real tmp826;
	       fftw_real tmp241;
	       fftw_real tmp246;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp823 = (K923879532 * tmp821) + (K382683432 * tmp822);
	       tmp826 = (K923879532 * tmp824) - (K382683432 * tmp825);
	       tmp827 = tmp823 + tmp826;
	       tmp908 = tmp826 - tmp823;
	       tmp241 = tmp233 + tmp240;
	       tmp246 = tmp242 - tmp245;
	       tmp247 = (K980785280 * tmp241) + (K195090322 * tmp246);
	       tmp276 = (K980785280 * tmp246) - (K195090322 * tmp241);
	  }
	  tmp491 = tmp233 - tmp240;
	  tmp492 = tmp245 + tmp242;
	  tmp493 = (K831469612 * tmp491) + (K555570233 * tmp492);
	  tmp502 = (K831469612 * tmp492) - (K555570233 * tmp491);
	  {
	       fftw_real tmp494;
	       fftw_real tmp495;
	       fftw_real tmp258;
	       fftw_real tmp263;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp494 = tmp250 - tmp257;
	       tmp495 = tmp262 + tmp259;
	       tmp496 = (K831469612 * tmp494) - (K555570233 * tmp495);
	       tmp503 = (K555570233 * tmp494) + (K831469612 * tmp495);
	       tmp258 = tmp250 + tmp257;
	       tmp263 = tmp259 - tmp262;
	       tmp264 = (K980785280 * tmp258) - (K195090322 * tmp263);
	       tmp277 = (K195090322 * tmp258) + (K980785280 * tmp263);
	  }
     }
     {
	  fftw_real tmp63;
	  fftw_real tmp739;
	  fftw_real tmp753;
	  fftw_real tmp759;
	  fftw_real tmp126;
	  fftw_real tmp755;
	  fftw_real tmp746;
	  fftw_real tmp758;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp31;
	       fftw_real tmp62;
	       fftw_real tmp749;
	       fftw_real tmp752;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp31 = tmp15 + tmp30;
	       tmp62 = tmp46 + tmp61;
	       tmp63 = tmp31 + tmp62;
	       tmp739 = tmp31 - tmp62;
	       tmp749 = tmp747 + tmp748;
	       tmp752 = tmp750 + tmp751;
	       tmp753 = tmp749 - tmp752;
	       tmp759 = tmp749 + tmp752;
	  }
	  {
	       fftw_real tmp94;
	       fftw_real tmp125;
	       fftw_real tmp742;
	       fftw_real tmp745;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp94 = tmp78 + tmp93;
	       tmp125 = tmp109 + tmp124;
	       tmp126 = tmp94 + tmp125;
	       tmp755 = tmp125 - tmp94;
	       tmp742 = tmp740 + tmp741;
	       tmp745 = tmp743 + tmp744;
	       tmp746 = tmp742 - tmp745;
	       tmp758 = tmp742 + tmp745;
	  }
	  {
	       fftw_real tmp757;
	       fftw_real tmp760;
	       fftw_real tmp754;
	       fftw_real tmp756;
	       ASSERT_ALIGNED_DOUBLE;
	       real_output[32 * real_ostride] = tmp63 - tmp126;
	       tmp757 = tmp63 + tmp126;
	       tmp760 = tmp758 + tmp759;
	       real_output[64 * real_ostride] = tmp757 - tmp760;
	       real_output[0] = tmp757 + tmp760;
	       imag_output[32 * imag_ostride] = tmp759 - tmp758;
	       tmp754 = K707106781 * (tmp746 + tmp753);
	       real_output[48 * real_ostride] = tmp739 - tmp754;
	       real_output[16 * real_ostride] = tmp739 + tmp754;
	       tmp756 = K707106781 * (tmp753 - tmp746);
	       imag_output[16 * imag_ostride] = tmp755 + tmp756;
	       imag_output[48 * imag_ostride] = tmp756 - tmp755;
	  }
     }
     {
	  fftw_real tmp761;
	  fftw_real tmp774;
	  fftw_real tmp764;
	  fftw_real tmp773;
	  fftw_real tmp768;
	  fftw_real tmp778;
	  fftw_real tmp771;
	  fftw_real tmp779;
	  fftw_real tmp762;
	  fftw_real tmp763;
	  ASSERT_ALIGNED_DOUBLE;
	  tmp761 = tmp15 - tmp30;
	  tmp774 = tmp61 - tmp46;
	  tmp762 = tmp78 - tmp93;
	  tmp763 = tmp109 - tmp124;
	  tmp764 = K707106781 * (tmp762 + tmp763);
	  tmp773 = K707106781 * (tmp763 - tmp762);
	  {
	       fftw_real tmp766;
	       fftw_real tmp767;
	       fftw_real tmp769;
	       fftw_real tmp770;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp766 = tmp740 - tmp741;
	       tmp767 = tmp744 - tmp743;
	       tmp768 = (K923879532 * tmp766) + (K382683432 * tmp767);
	       tmp778 = (K923879532 * tmp767) - (K382683432 * tmp766);
	       tmp769 = tmp747 - tmp748;
	       tmp770 = tmp751 - tmp750;
	       tmp771 = (K923879532 * tmp769) - (K382683432 * tmp770);
	       tmp779 = (K382683432 * tmp769) + (K923879532 * tmp770);
	  }
	  {
	       fftw_real tmp765;
	       fftw_real tmp772;
	       fftw_real tmp781;
	       fftw_real tmp782;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp765 = tmp761 + tmp764;
	       tmp772 = tmp768 + tmp771;
	       real_output[56 * real_ostride] = tmp765 - tmp772;
	       real_output[8 * real_ostride] = tmp765 + tmp772;
	       tmp781 = tmp761 - tmp764;
	       tmp782 = tmp779 - tmp778;
	       real_output[40 * real_ostride] = tmp781 - tmp782;
	       real_output[24 * real_ostride] = tmp781 + tmp782;
	  }
	  {
	       fftw_real tmp777;
	       fftw_real tmp780;
	       fftw_real tmp775;
	       fftw_real tmp776;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp777 = tmp774 + tmp773;
	       tmp780 = tmp778 + tmp779;
	       imag_output[8 * imag_ostride] = tmp777 + tmp780;
	       imag_output[56 * imag_ostride] = tmp780 - tmp777;
	       tmp775 = tmp773 - tmp774;
	       tmp776 = tmp771 - tmp768;
	       imag_output[24 * imag_ostride] = tmp775 + tmp776;
	       imag_output[40 * imag_ostride] = tmp776 - tmp775;
	  }
     }
     {
	  fftw_real tmp159;
	  fftw_real tmp381;
	  fftw_real tmp202;
	  fftw_real tmp391;
	  fftw_real tmp372;
	  fftw_real tmp392;
	  fftw_real tmp361;
	  fftw_real tmp382;
	  fftw_real tmp357;
	  fftw_real tmp377;
	  fftw_real tmp389;
	  fftw_real tmp397;
	  fftw_real tmp280;
	  fftw_real tmp376;
	  fftw_real tmp386;
	  fftw_real tmp396;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp139;
	       fftw_real tmp158;
	       fftw_real tmp359;
	       fftw_real tmp360;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp139 = tmp131 + tmp138;
	       tmp158 = tmp148 + tmp157;
	       tmp159 = tmp139 + tmp158;
	       tmp381 = tmp139 - tmp158;
	       {
		    fftw_real tmp180;
		    fftw_real tmp201;
		    fftw_real tmp368;
		    fftw_real tmp371;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp180 = (K098017140 * tmp172) + (K995184726 * tmp179);
		    tmp201 = (K995184726 * tmp193) - (K098017140 * tmp200);
		    tmp202 = tmp180 + tmp201;
		    tmp391 = tmp201 - tmp180;
		    tmp368 = tmp364 + tmp367;
		    tmp371 = tmp369 + tmp370;
		    tmp372 = tmp368 + tmp371;
		    tmp392 = tmp371 - tmp368;
	       }
	       tmp359 = (K995184726 * tmp172) - (K098017140 * tmp179);
	       tmp360 = (K995184726 * tmp200) + (K098017140 * tmp193);
	       tmp361 = tmp359 + tmp360;
	       tmp382 = tmp360 - tmp359;
	       {
		    fftw_real tmp343;
		    fftw_real tmp387;
		    fftw_real tmp356;
		    fftw_real tmp388;
		    fftw_real tmp342;
		    fftw_real tmp355;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp342 = tmp324 + tmp341;
		    tmp343 = tmp307 + tmp342;
		    tmp387 = tmp307 - tmp342;
		    tmp355 = tmp353 + tmp354;
		    tmp356 = tmp352 + tmp355;
		    tmp388 = tmp355 - tmp352;
		    tmp357 = (K998795456 * tmp343) - (K049067674 * tmp356);
		    tmp377 = (K049067674 * tmp343) + (K998795456 * tmp356);
		    tmp389 = (K740951125 * tmp387) - (K671558954 * tmp388);
		    tmp397 = (K671558954 * tmp387) + (K740951125 * tmp388);
	       }
	       {
		    fftw_real tmp266;
		    fftw_real tmp384;
		    fftw_real tmp279;
		    fftw_real tmp385;
		    fftw_real tmp265;
		    fftw_real tmp278;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp265 = tmp247 + tmp264;
		    tmp266 = tmp230 + tmp265;
		    tmp384 = tmp230 - tmp265;
		    tmp278 = tmp276 + tmp277;
		    tmp279 = tmp275 + tmp278;
		    tmp385 = tmp278 - tmp275;
		    tmp280 = (K998795456 * tmp266) + (K049067674 * tmp279);
		    tmp376 = (K998795456 * tmp279) - (K049067674 * tmp266);
		    tmp386 = (K740951125 * tmp384) + (K671558954 * tmp385);
		    tmp396 = (K740951125 * tmp385) - (K671558954 * tmp384);
	       }
	  }
	  {
	       fftw_real tmp203;
	       fftw_real tmp358;
	       fftw_real tmp373;
	       fftw_real tmp374;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp203 = tmp159 + tmp202;
	       tmp358 = tmp280 + tmp357;
	       real_output[63 * real_ostride] = tmp203 - tmp358;
	       real_output[real_ostride] = tmp203 + tmp358;
	       tmp373 = tmp361 - tmp372;
	       tmp374 = tmp357 - tmp280;
	       imag_output[31 * imag_ostride] = tmp373 + tmp374;
	       imag_output[33 * imag_ostride] = tmp374 - tmp373;
	  }
	  {
	       fftw_real tmp375;
	       fftw_real tmp378;
	       fftw_real tmp379;
	       fftw_real tmp380;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp375 = tmp372 + tmp361;
	       tmp378 = tmp376 + tmp377;
	       imag_output[imag_ostride] = tmp375 + tmp378;
	       imag_output[63 * imag_ostride] = tmp378 - tmp375;
	       tmp379 = tmp159 - tmp202;
	       tmp380 = tmp377 - tmp376;
	       real_output[33 * real_ostride] = tmp379 - tmp380;
	       real_output[31 * real_ostride] = tmp379 + tmp380;
	  }
	  {
	       fftw_real tmp383;
	       fftw_real tmp390;
	       fftw_real tmp393;
	       fftw_real tmp394;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp383 = tmp381 + tmp382;
	       tmp390 = tmp386 + tmp389;
	       real_output[49 * real_ostride] = tmp383 - tmp390;
	       real_output[15 * real_ostride] = tmp383 + tmp390;
	       tmp393 = tmp391 - tmp392;
	       tmp394 = tmp389 - tmp386;
	       imag_output[17 * imag_ostride] = tmp393 + tmp394;
	       imag_output[47 * imag_ostride] = tmp394 - tmp393;
	  }
	  {
	       fftw_real tmp395;
	       fftw_real tmp398;
	       fftw_real tmp399;
	       fftw_real tmp400;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp395 = tmp392 + tmp391;
	       tmp398 = tmp396 + tmp397;
	       imag_output[15 * imag_ostride] = tmp395 + tmp398;
	       imag_output[49 * imag_ostride] = tmp398 - tmp395;
	       tmp399 = tmp381 - tmp382;
	       tmp400 = tmp397 - tmp396;
	       real_output[47 * real_ostride] = tmp399 - tmp400;
	       real_output[17 * real_ostride] = tmp399 + tmp400;
	  }
     }
     {
	  fftw_real tmp403;
	  fftw_real tmp441;
	  fftw_real tmp410;
	  fftw_real tmp451;
	  fftw_real tmp432;
	  fftw_real tmp452;
	  fftw_real tmp429;
	  fftw_real tmp442;
	  fftw_real tmp425;
	  fftw_real tmp437;
	  fftw_real tmp449;
	  fftw_real tmp457;
	  fftw_real tmp418;
	  fftw_real tmp436;
	  fftw_real tmp446;
	  fftw_real tmp456;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp401;
	       fftw_real tmp402;