fni_64.c

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

	       fftw_real tmp309;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp85 = tmp81 + tmp84;
	       tmp92 = tmp88 + tmp91;
	       tmp93 = tmp85 + tmp92;
	       tmp865 = tmp85 - tmp92;
	       tmp308 = (K923879532 * tmp295) - (K382683432 * tmp290);
	       tmp309 = (K382683432 * tmp279) + (K923879532 * tmp284);
	       tmp310 = tmp308 - tmp309;
	       tmp426 = tmp308 + tmp309;
	  }
	  {
	       fftw_real tmp285;
	       fftw_real tmp296;
	       fftw_real tmp822;
	       fftw_real tmp823;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp285 = (K923879532 * tmp279) - (K382683432 * tmp284);
	       tmp296 = (K923879532 * tmp290) + (K382683432 * tmp295);
	       tmp297 = tmp285 - tmp296;
	       tmp423 = tmp296 + tmp285;
	       tmp822 = tmp668 + tmp669;
	       tmp823 = tmp673 + tmp674;
	       tmp824 = tmp822 + tmp823;
	       tmp868 = tmp823 - tmp822;
	  }
	  {
	       fftw_real tmp503;
	       fftw_real tmp506;
	       fftw_real tmp512;
	       fftw_real tmp513;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp503 = (K382683432 * tmp501) - (K923879532 * tmp502);
	       tmp506 = (K382683432 * tmp504) + (K923879532 * tmp505);
	       tmp507 = tmp503 - tmp506;
	       tmp587 = tmp506 + tmp503;
	       tmp512 = (K382683432 * tmp505) - (K923879532 * tmp504);
	       tmp513 = (K923879532 * tmp501) + (K382683432 * tmp502);
	       tmp514 = tmp512 - tmp513;
	       tmp590 = tmp512 + tmp513;
	  }
	  {
	       fftw_real tmp671;
	       fftw_real tmp676;
	       fftw_real tmp684;
	       fftw_real tmp685;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp671 = tmp667 - tmp670;
	       tmp676 = tmp672 + tmp675;
	       tmp677 = K707106781 * (tmp671 + tmp676);
	       tmp767 = K707106781 * (tmp671 - tmp676);
	       tmp684 = tmp667 + tmp670;
	       tmp685 = tmp675 - tmp672;
	       tmp686 = K707106781 * (tmp684 + tmp685);
	       tmp770 = K707106781 * (tmp685 - tmp684);
	  }
     }
     {
	  fftw_real tmp63;
	  fftw_real tmp851;
	  fftw_real tmp850;
	  fftw_real tmp852;
	  fftw_real tmp126;
	  fftw_real tmp127;
	  fftw_real tmp190;
	  fftw_real tmp847;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp31;
	       fftw_real tmp62;
	       fftw_real tmp848;
	       fftw_real tmp849;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp31 = tmp15 + tmp30;
	       tmp62 = tmp46 + tmp61;
	       tmp63 = tmp31 + tmp62;
	       tmp851 = tmp31 - tmp62;
	       tmp848 = tmp821 + tmp824;
	       tmp849 = tmp830 + tmp833;
	       tmp850 = tmp848 + tmp849;
	       tmp852 = tmp849 - tmp848;
	  }
	  {
	       fftw_real tmp94;
	       fftw_real tmp125;
	       fftw_real tmp158;
	       fftw_real tmp189;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp94 = tmp78 + tmp93;
	       tmp125 = tmp109 + tmp124;
	       tmp126 = tmp94 + tmp125;
	       tmp127 = tmp94 - tmp125;
	       tmp158 = tmp142 + tmp157;
	       tmp189 = tmp173 + tmp188;
	       tmp190 = tmp158 - tmp189;
	       tmp847 = tmp158 + tmp189;
	  }
	  c_re(output[32 * ostride]) = tmp63 - tmp126;
	  c_re(output[0]) = tmp63 + tmp126;
	  c_im(output[16 * ostride]) = tmp127 + tmp190;
	  c_im(output[48 * ostride]) = tmp190 - tmp127;
	  c_im(output[32 * ostride]) = tmp847 - tmp850;
	  c_im(output[0]) = tmp847 + tmp850;
	  c_re(output[48 * ostride]) = tmp851 - tmp852;
	  c_re(output[16 * ostride]) = tmp851 + tmp852;
     }
     {
	  fftw_real tmp817;
	  fftw_real tmp841;
	  fftw_real tmp839;
	  fftw_real tmp845;
	  fftw_real tmp826;
	  fftw_real tmp842;
	  fftw_real tmp835;
	  fftw_real tmp843;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp815;
	       fftw_real tmp816;
	       fftw_real tmp837;
	       fftw_real tmp838;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp815 = tmp142 - tmp157;
	       tmp816 = tmp46 - tmp61;
	       tmp817 = tmp815 - tmp816;
	       tmp841 = tmp816 + tmp815;
	       tmp837 = tmp15 - tmp30;
	       tmp838 = tmp188 - tmp173;
	       tmp839 = tmp837 + tmp838;
	       tmp845 = tmp837 - tmp838;
	  }
	  {
	       fftw_real tmp818;
	       fftw_real tmp825;
	       fftw_real tmp827;
	       fftw_real tmp834;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp818 = tmp78 - tmp93;
	       tmp825 = tmp821 - tmp824;
	       tmp826 = tmp818 - tmp825;
	       tmp842 = tmp818 + tmp825;
	       tmp827 = tmp109 - tmp124;
	       tmp834 = tmp830 - tmp833;
	       tmp835 = tmp827 + tmp834;
	       tmp843 = tmp834 - tmp827;
	  }
	  {
	       fftw_real tmp836;
	       fftw_real tmp840;
	       fftw_real tmp844;
	       fftw_real tmp846;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp836 = K707106781 * (tmp826 - tmp835);
	       c_im(output[56 * ostride]) = tmp817 - tmp836;
	       c_im(output[24 * ostride]) = tmp817 + tmp836;
	       tmp840 = K707106781 * (tmp826 + tmp835);
	       c_re(output[40 * ostride]) = tmp839 - tmp840;
	       c_re(output[8 * ostride]) = tmp839 + tmp840;
	       tmp844 = K707106781 * (tmp842 + tmp843);
	       c_im(output[40 * ostride]) = tmp841 - tmp844;
	       c_im(output[8 * ostride]) = tmp841 + tmp844;
	       tmp846 = K707106781 * (tmp843 - tmp842);
	       c_re(output[56 * ostride]) = tmp845 - tmp846;
	       c_re(output[24 * ostride]) = tmp845 + tmp846;
	  }
     }
     {
	  fftw_real tmp217;
	  fftw_real tmp391;
	  fftw_real tmp396;
	  fftw_real tmp406;
	  fftw_real tmp399;
	  fftw_real tmp407;
	  fftw_real tmp367;
	  fftw_real tmp387;
	  fftw_real tmp312;
	  fftw_real tmp386;
	  fftw_real tmp379;
	  fftw_real tmp401;
	  fftw_real tmp382;
	  fftw_real tmp392;
	  fftw_real tmp256;
	  fftw_real tmp402;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp201;
	       fftw_real tmp216;
	       fftw_real tmp394;
	       fftw_real tmp395;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp201 = tmp193 - tmp200;
	       tmp216 = tmp208 - tmp215;
	       tmp217 = tmp201 - tmp216;
	       tmp391 = tmp201 + tmp216;
	       tmp394 = tmp274 + tmp297;
	       tmp395 = tmp307 + tmp310;
	       tmp396 = (K881921264 * tmp394) - (K471396736 * tmp395);
	       tmp406 = (K471396736 * tmp394) + (K881921264 * tmp395);
	  }
	  {
	       fftw_real tmp397;
	       fftw_real tmp398;
	       fftw_real tmp353;
	       fftw_real tmp366;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp397 = tmp329 + tmp352;
	       tmp398 = tmp362 + tmp365;
	       tmp399 = (K881921264 * tmp397) + (K471396736 * tmp398);
	       tmp407 = (K881921264 * tmp398) - (K471396736 * tmp397);
	       tmp353 = tmp329 - tmp352;
	       tmp366 = tmp362 - tmp365;
	       tmp367 = (K290284677 * tmp353) + (K956940335 * tmp366);
	       tmp387 = (K290284677 * tmp366) - (K956940335 * tmp353);
	  }
	  {
	       fftw_real tmp298;
	       fftw_real tmp311;
	       fftw_real tmp375;
	       fftw_real tmp378;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp298 = tmp274 - tmp297;
	       tmp311 = tmp307 - tmp310;
	       tmp312 = (K290284677 * tmp298) - (K956940335 * tmp311);
	       tmp386 = (K956940335 * tmp298) + (K290284677 * tmp311);
	       tmp375 = tmp371 - tmp374;
	       tmp378 = tmp376 - tmp377;
	       tmp379 = tmp375 - tmp378;
	       tmp401 = tmp375 + tmp378;
	  }
	  {
	       fftw_real tmp380;
	       fftw_real tmp381;
	       fftw_real tmp236;
	       fftw_real tmp255;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp380 = (K555570233 * tmp254) - (K831469612 * tmp247);
	       tmp381 = (K831469612 * tmp228) + (K555570233 * tmp235);
	       tmp382 = tmp380 - tmp381;
	       tmp392 = tmp381 + tmp380;
	       tmp236 = (K555570233 * tmp228) - (K831469612 * tmp235);
	       tmp255 = (K555570233 * tmp247) + (K831469612 * tmp254);
	       tmp256 = tmp236 - tmp255;
	       tmp402 = tmp236 + tmp255;
	  }
	  {
	       fftw_real tmp257;
	       fftw_real tmp368;
	       fftw_real tmp383;
	       fftw_real tmp384;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp257 = tmp217 - tmp256;
	       tmp368 = tmp312 - tmp367;
	       c_im(output[61 * ostride]) = tmp257 - tmp368;
	       c_im(output[29 * ostride]) = tmp257 + tmp368;
	       tmp383 = tmp379 + tmp382;
	       tmp384 = tmp312 + tmp367;
	       c_re(output[45 * ostride]) = tmp383 - tmp384;
	       c_re(output[13 * ostride]) = tmp383 + tmp384;
	  }
	  {
	       fftw_real tmp389;
	       fftw_real tmp390;
	       fftw_real tmp385;
	       fftw_real tmp388;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp389 = tmp379 - tmp382;
	       tmp390 = tmp387 - tmp386;
	       c_re(output[61 * ostride]) = tmp389 - tmp390;
	       c_re(output[29 * ostride]) = tmp389 + tmp390;
	       tmp385 = tmp217 + tmp256;
	       tmp388 = tmp386 + tmp387;
	       c_im(output[45 * ostride]) = tmp385 - tmp388;
	       c_im(output[13 * ostride]) = tmp385 + tmp388;
	  }
	  {
	       fftw_real tmp393;
	       fftw_real tmp400;
	       fftw_real tmp403;
	       fftw_real tmp404;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp393 = tmp391 - tmp392;
	       tmp400 = tmp396 - tmp399;
	       c_im(output[53 * ostride]) = tmp393 - tmp400;
	       c_im(output[21 * ostride]) = tmp393 + tmp400;
	       tmp403 = tmp401 + tmp402;
	       tmp404 = tmp396 + tmp399;
	       c_re(output[37 * ostride]) = tmp403 - tmp404;
	       c_re(output[5 * ostride]) = tmp403 + tmp404;
	  }
	  {
	       fftw_real tmp409;
	       fftw_real tmp410;
	       fftw_real tmp405;
	       fftw_real tmp408;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp409 = tmp401 - tmp402;
	       tmp410 = tmp407 - tmp406;
	       c_re(output[53 * ostride]) = tmp409 - tmp410;
	       c_re(output[21 * ostride]) = tmp409 + tmp410;
	       tmp405 = tmp391 + tmp392;
	       tmp408 = tmp406 + tmp407;
	       c_im(output[37 * ostride]) = tmp405 - tmp408;
	       c_im(output[5 * ostride]) = tmp405 + tmp408;
	  }
     }
     {
	  fftw_real tmp413;
	  fftw_real tmp451;
	  fftw_real tmp456;
	  fftw_real tmp466;
	  fftw_real tmp459;
	  fftw_real tmp467;
	  fftw_real tmp435;
	  fftw_real tmp447;
	  fftw_real tmp428;
	  fftw_real tmp446;
	  fftw_real tmp439;
	  fftw_real tmp461;
	  fftw_real tmp442;
	  fftw_real tmp452;
	  fftw_real tmp420;
	  fftw_real tmp462;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp411;
	       fftw_real tmp412;
	       fftw_real tmp454;
	       fftw_real tmp455;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp411 = tmp193 + tmp200;
	       tmp412 = tmp377 + tmp376;
	       tmp413 = tmp411 - tmp412;
	       tmp451 = tmp411 + tmp412;
	       tmp454 = tmp422 + tmp423;
	       tmp455 = tmp425 + tmp426;
	       tmp456 = (K995184726 * tmp454) + (K098017140 * tmp455);
	       tmp466 = (K995184726 * tmp455) - (K098017140 * tmp454);
	  }
	  {
	       fftw_real tmp457;
	       fftw_real tmp458;
	       fftw_real tmp431;
	       fftw_real tmp434;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp457 = tmp429 + tmp430;
	       tmp458 = tmp432 + tmp433;
	       tmp459 = (K995184726 * tmp457) - (K098017140 * tmp458);
	       tmp467 = (K098017140 * tmp457) + (K995184726 * tmp458);
	       tmp431 = tmp429 - tmp430;
	       tmp434 = tmp432 - tmp433;
	       tmp435 = (K634393284 * tmp431) - (K773010453 * tmp434);
	       tmp447 = (K773010453 * tmp431) + (K634393284 * tmp434);
	  }
	  {
	       fftw_real tmp424;
	       fftw_real tmp427;
	       fftw_real tmp437;
	       fftw_real tmp438;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp424 = tmp422 - tmp423;
	       tmp427 = tmp425 - tmp426;
	       tmp428 = (K634393284 * tmp424) + (K773010453 * tmp427);
	       tmp446 = (K634393284 * tmp427) - (K773010453 * tmp424);
	       tmp437 = tmp371 + tmp374;
	       tmp438 = tmp208 + tmp215;
	       tmp439 = tmp437 - tmp438;
	       tmp461 = tmp437 + tmp438;
	  }
	  {
	       fftw_real tmp440;
	       fftw_real tmp441;
	       fftw_real tmp416;
	       fftw_real tmp419;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp440 = (K980785280 * tmp418) - (K195090322 * tmp417);
	       tmp441 = (K195090322 * tmp414) + (K980785280 * tmp415);
	       tmp442 = tmp440 - tmp441;
	       tmp452 = tmp441 + tmp440;
	       tmp416 = (K980785280 * tmp414) - (K195090322 * tmp415);
	       tmp419 = (K980785280 * tmp417) + (K195090322 * tmp418);
	       tmp420 = tmp416 - tmp419;
	       tmp462 = tmp416 + tmp419;
	  }
	  {
	       fftw_real tmp421;
	       fftw_real tmp436;
	       fftw_real tmp443;
	       fftw_real tmp444;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp421 = tmp413 + tmp420;
	       tmp436 = tmp428 + tmp435;
	       c_im(output[41 * ostride]) = tmp421 - tmp436;
	       c_im(output[9 * ostride]) = tmp421 + tmp436;
	       tmp443 = tmp439 - tmp442;
	       tmp444 = tmp435 - tmp428;
	       c_re(output[57 * ostride]) = tmp443 - tmp444;
	       c_re(output[25 * ostride]) = tmp443 + tmp444;
	  }
	  {
	       fftw_real tmp449;
	       fftw_real tmp450;
	       fftw_real tmp445;
	       fftw_real tmp448;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp449 = tmp439 + tmp442;
	       tmp450 = tmp446 + tmp447;
	       c_re(output[41 * ostride]) = tmp449 - tmp450;
	       c_re(output[9 * ostride]) = tmp449 + tmp450;
	       tmp445 = tmp413 - tmp420;
	       tmp448 = tmp446 - tmp447;
	       c_im(output[57 * ostride]) = tmp445 - tmp448;
	       c_im(output[25 * ostride]) = tmp445 + tmp448;
	  }
	  {
	       fftw_real tmp453;
	       fftw_real tmp460;
	       fftw_real tmp463;
	       fftw_real tmp464;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp453 = tmp451 + tmp452;
	       tmp460 = tmp456 + tmp459;
	       c_im(output[33 * ostride]) = tmp453 - tmp460;
	       c_im(output[ostride]) = tmp453 + tmp460;
	       tmp463 = tmp461 - tmp462;
	       tmp464 = tmp459 - tmp456;
	       c_re(output[49 * ostride]) = tmp463 - tmp464;
	       c_re(output[17 * ostride]) = tmp463 + tmp464;
	  }
	  {
	       fftw_real tmp469;
	       fftw_real tmp470;
	       fftw_real tmp465;
	       fftw_real tmp468;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp469 = tmp461 + tmp462;
	       tmp470 = tmp466 + tmp467;
	       c_re(output[33 * ostride]) = tmp469 - tmp470;
	       c_re(output[ostride]) = tmp469 + tmp470;
	       tmp465 = tmp451 - tmp452;
	       tmp468 = tmp466 - tmp467;
	       c_im(output[49 * ostride]) = tmp465 - tmp468;
	       c_im(output[17 * ostride]) = tmp465 + tmp468;
	  }
     }
     {
	  fftw_real tmp855;