fn_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;
	       tmp829 = tmp92 - tmp85;
	       tmp308 = (K382683432 * tmp290) - (K923879532 * tmp295);
	       tmp309 = (K382683432 * tmp284) + (K923879532 * tmp279);
	       tmp310 = tmp308 - tmp309;
	       tmp426 = tmp308 + tmp309;
	  }
	  {
	       fftw_real tmp285;
	       fftw_real tmp296;
	       fftw_real tmp832;
	       fftw_real tmp833;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp285 = (K382683432 * tmp279) - (K923879532 * tmp284);
	       tmp296 = (K923879532 * tmp290) + (K382683432 * tmp295);
	       tmp297 = tmp285 - tmp296;
	       tmp423 = tmp296 + tmp285;
	       tmp832 = tmp673 + tmp674;
	       tmp833 = tmp668 + tmp669;
	       tmp834 = tmp832 - tmp833;
	       tmp888 = tmp832 + tmp833;
	  }
	  {
	       fftw_real tmp503;
	       fftw_real tmp506;
	       fftw_real tmp512;
	       fftw_real tmp513;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp503 = (K923879532 * tmp501) - (K382683432 * tmp502);
	       tmp506 = (K382683432 * tmp504) + (K923879532 * tmp505);
	       tmp507 = tmp503 - tmp506;
	       tmp587 = tmp506 + tmp503;
	       tmp512 = (K923879532 * tmp504) - (K382683432 * tmp505);
	       tmp513 = (K923879532 * tmp502) + (K382683432 * tmp501);
	       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 * (tmp676 + tmp671);
	       tmp684 = tmp675 - tmp672;
	       tmp685 = tmp667 + tmp670;
	       tmp686 = K707106781 * (tmp684 - tmp685);
	       tmp770 = K707106781 * (tmp684 + tmp685);
	  }
     }
     {
	  fftw_real tmp63;
	  fftw_real tmp907;
	  fftw_real tmp910;
	  fftw_real tmp912;
	  fftw_real tmp126;
	  fftw_real tmp127;
	  fftw_real tmp190;
	  fftw_real tmp911;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp31;
	       fftw_real tmp62;
	       fftw_real tmp908;
	       fftw_real tmp909;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp31 = tmp15 + tmp30;
	       tmp62 = tmp46 + tmp61;
	       tmp63 = tmp31 + tmp62;
	       tmp907 = tmp31 - tmp62;
	       tmp908 = tmp887 + tmp888;
	       tmp909 = tmp892 + tmp893;
	       tmp910 = tmp908 - tmp909;
	       tmp912 = tmp908 + tmp909;
	  }
	  {
	       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 = tmp125 - tmp94;
	       tmp158 = tmp142 + tmp157;
	       tmp189 = tmp173 + tmp188;
	       tmp190 = tmp158 - tmp189;
	       tmp911 = 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_re(output[48 * ostride]) = tmp907 - tmp910;
	  c_re(output[16 * ostride]) = tmp907 + tmp910;
	  c_im(output[32 * ostride]) = tmp911 - tmp912;
	  c_im(output[0]) = tmp911 + tmp912;
     }
     {
	  fftw_real tmp885;
	  fftw_real tmp901;
	  fftw_real tmp899;
	  fftw_real tmp905;
	  fftw_real tmp890;
	  fftw_real tmp902;
	  fftw_real tmp895;
	  fftw_real tmp903;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp883;
	       fftw_real tmp884;
	       fftw_real tmp897;
	       fftw_real tmp898;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp883 = tmp15 - tmp30;
	       tmp884 = tmp173 - tmp188;
	       tmp885 = tmp883 + tmp884;
	       tmp901 = tmp883 - tmp884;
	       tmp897 = tmp142 - tmp157;
	       tmp898 = tmp61 - tmp46;
	       tmp899 = tmp897 - tmp898;
	       tmp905 = tmp898 + tmp897;
	  }
	  {
	       fftw_real tmp886;
	       fftw_real tmp889;
	       fftw_real tmp891;
	       fftw_real tmp894;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp886 = tmp78 - tmp93;
	       tmp889 = tmp887 - tmp888;
	       tmp890 = tmp886 + tmp889;
	       tmp902 = tmp889 - tmp886;
	       tmp891 = tmp109 - tmp124;
	       tmp894 = tmp892 - tmp893;
	       tmp895 = tmp891 - tmp894;
	       tmp903 = tmp891 + tmp894;
	  }
	  {
	       fftw_real tmp896;
	       fftw_real tmp900;
	       fftw_real tmp904;
	       fftw_real tmp906;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp896 = K707106781 * (tmp890 + tmp895);
	       c_re(output[40 * ostride]) = tmp885 - tmp896;
	       c_re(output[8 * ostride]) = tmp885 + tmp896;
	       tmp900 = K707106781 * (tmp895 - tmp890);
	       c_im(output[56 * ostride]) = tmp899 - tmp900;
	       c_im(output[24 * ostride]) = tmp899 + tmp900;
	       tmp904 = K707106781 * (tmp902 - tmp903);
	       c_re(output[56 * ostride]) = tmp901 - tmp904;
	       c_re(output[24 * ostride]) = tmp901 + tmp904;
	       tmp906 = K707106781 * (tmp902 + tmp903);
	       c_im(output[40 * ostride]) = tmp905 - tmp906;
	       c_im(output[8 * ostride]) = tmp905 + tmp906;
	  }
     }
     {
	  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 = (K634393284 * tmp394) + (K773010453 * tmp395);
	       tmp406 = (K773010453 * tmp394) - (K634393284 * tmp395);
	  }
	  {
	       fftw_real tmp397;
	       fftw_real tmp398;
	       fftw_real tmp353;
	       fftw_real tmp366;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp397 = tmp329 + tmp352;
	       tmp398 = tmp362 + tmp365;
	       tmp399 = (K773010453 * tmp397) - (K634393284 * tmp398);
	       tmp407 = (K773010453 * tmp398) + (K634393284 * tmp397);
	       tmp353 = tmp329 - tmp352;
	       tmp366 = tmp362 - tmp365;
	       tmp367 = (K098017140 * tmp353) - (K995184726 * tmp366);
	       tmp387 = (K098017140 * tmp366) + (K995184726 * tmp353);
	  }
	  {
	       fftw_real tmp298;
	       fftw_real tmp311;
	       fftw_real tmp375;
	       fftw_real tmp378;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp298 = tmp274 - tmp297;
	       tmp311 = tmp307 - tmp310;
	       tmp312 = (K995184726 * tmp298) + (K098017140 * tmp311);
	       tmp386 = (K098017140 * tmp298) - (K995184726 * 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 = (K195090322 * tmp254) - (K980785280 * tmp247);
	       tmp381 = (K980785280 * tmp228) + (K195090322 * tmp235);
	       tmp382 = tmp380 - tmp381;
	       tmp392 = tmp381 + tmp380;
	       tmp236 = (K195090322 * tmp228) - (K980785280 * tmp235);
	       tmp255 = (K195090322 * tmp247) + (K980785280 * 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_re(output[47 * ostride]) = tmp257 - tmp368;
	       c_re(output[15 * ostride]) = tmp257 + tmp368;
	       tmp383 = tmp379 - tmp382;
	       tmp384 = tmp367 - tmp312;
	       c_im(output[63 * ostride]) = tmp383 - tmp384;
	       c_im(output[31 * ostride]) = tmp383 + tmp384;
	  }
	  {
	       fftw_real tmp389;
	       fftw_real tmp390;
	       fftw_real tmp385;
	       fftw_real tmp388;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp389 = tmp379 + tmp382;
	       tmp390 = tmp386 + tmp387;
	       c_im(output[47 * ostride]) = tmp389 - tmp390;
	       c_im(output[15 * ostride]) = tmp389 + tmp390;
	       tmp385 = tmp217 - tmp256;
	       tmp388 = tmp386 - tmp387;
	       c_re(output[63 * ostride]) = tmp385 - tmp388;
	       c_re(output[31 * 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_re(output[39 * ostride]) = tmp393 - tmp400;
	       c_re(output[7 * ostride]) = tmp393 + tmp400;
	       tmp403 = tmp401 - tmp402;
	       tmp404 = tmp399 - tmp396;
	       c_im(output[55 * ostride]) = tmp403 - tmp404;
	       c_im(output[23 * ostride]) = tmp403 + tmp404;
	  }
	  {
	       fftw_real tmp409;
	       fftw_real tmp410;
	       fftw_real tmp405;
	       fftw_real tmp408;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp409 = tmp401 + tmp402;
	       tmp410 = tmp406 + tmp407;
	       c_im(output[39 * ostride]) = tmp409 - tmp410;
	       c_im(output[7 * ostride]) = tmp409 + tmp410;
	       tmp405 = tmp391 - tmp392;
	       tmp408 = tmp406 - tmp407;
	       c_re(output[55 * ostride]) = tmp405 - tmp408;
	       c_re(output[23 * 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 = (K956940335 * tmp454) + (K290284677 * tmp455);
	       tmp466 = (K956940335 * tmp455) - (K290284677 * tmp454);
	  }
	  {
	       fftw_real tmp457;
	       fftw_real tmp458;
	       fftw_real tmp431;
	       fftw_real tmp434;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp457 = tmp429 + tmp430;
	       tmp458 = tmp432 + tmp433;
	       tmp459 = (K956940335 * tmp457) - (K290284677 * tmp458);
	       tmp467 = (K290284677 * tmp457) + (K956940335 * tmp458);
	       tmp431 = tmp429 - tmp430;
	       tmp434 = tmp432 - tmp433;
	       tmp435 = (K471396736 * tmp431) - (K881921264 * tmp434);
	       tmp447 = (K881921264 * tmp431) + (K471396736 * tmp434);
	  }
	  {
	       fftw_real tmp424;
	       fftw_real tmp427;
	       fftw_real tmp437;
	       fftw_real tmp438;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp424 = tmp422 - tmp423;
	       tmp427 = tmp425 - tmp426;
	       tmp428 = (K471396736 * tmp424) + (K881921264 * tmp427);
	       tmp446 = (K471396736 * tmp427) - (K881921264 * 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 = (K831469612 * tmp418) - (K555570233 * tmp417);
	       tmp441 = (K555570233 * tmp414) + (K831469612 * tmp415);
	       tmp442 = tmp440 - tmp441;
	       tmp452 = tmp441 + tmp440;
	       tmp416 = (K831469612 * tmp414) - (K555570233 * tmp415);
	       tmp419 = (K831469612 * tmp417) + (K555570233 * 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_re(output[43 * ostride]) = tmp421 - tmp436;
	       c_re(output[11 * ostride]) = tmp421 + tmp436;
	       tmp443 = tmp439 - tmp442;
	       tmp444 = tmp435 - tmp428;
	       c_im(output[59 * ostride]) = tmp443 - tmp444;
	       c_im(output[27 * 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_im(output[43 * ostride]) = tmp449 - tmp450;
	       c_im(output[11 * ostride]) = tmp449 + tmp450;
	       tmp445 = tmp413 - tmp420;
	       tmp448 = tmp446 - tmp447;
	       c_re(output[59 * ostride]) = tmp445 - tmp448;
	       c_re(output[27 * 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_re(output[35 * ostride]) = tmp453 - tmp460;
	       c_re(output[3 * ostride]) = tmp453 + tmp460;
	       tmp463 = tmp461 - tmp462;
	       tmp464 = tmp459 - tmp456;
	       c_im(output[51 * ostride]) = tmp463 - tmp464;
	       c_im(output[19 * 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_im(output[35 * ostride]) = tmp469 - tmp470;
	       c_im(output[3 * ostride]) = tmp469 + tmp470;
	       tmp465 = tmp451 - tmp452;
	       tmp468 = tmp466 - tmp467;
	       c_re(output[51 * ostride]) = tmp465 - tmp468;
	       c_re(output[19 * ostride]) = tmp465 + tmp468;
	  }
     }
     {
	  fftw_real tmp817;