fcr_128.c

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

	       fftw_real tmp315;
	       fftw_real tmp316;
	       fftw_real tmp317;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp308 = tmp304 + tmp307;
	       tmp315 = (K1_847759065 * tmp311) - (K765366864 * tmp314);
	       output[74 * ostride] = tmp308 - tmp315;
	       output[10 * ostride] = tmp308 + tmp315;
	       tmp316 = tmp304 - tmp307;
	       tmp317 = (K1_847759065 * tmp314) + (K765366864 * tmp311);
	       output[42 * ostride] = tmp316 - tmp317;
	       output[106 * ostride] = tmp316 + tmp317;
	  }
	  {
	       fftw_real tmp320;
	       fftw_real tmp323;
	       fftw_real tmp324;
	       fftw_real tmp325;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp320 = tmp318 - tmp319;
	       tmp323 = (K765366864 * tmp321) - (K1_847759065 * tmp322);
	       output[90 * ostride] = tmp320 - tmp323;
	       output[26 * ostride] = tmp320 + tmp323;
	       tmp324 = tmp318 + tmp319;
	       tmp325 = (K765366864 * tmp322) + (K1_847759065 * tmp321);
	       output[58 * ostride] = tmp324 - tmp325;
	       output[122 * ostride] = tmp324 + tmp325;
	  }
     }
     {
	  fftw_real tmp460;
	  fftw_real tmp488;
	  fftw_real tmp487;
	  fftw_real tmp492;
	  fftw_real tmp483;
	  fftw_real tmp491;
	  fftw_real tmp467;
	  fftw_real tmp489;
	  fftw_real tmp468;
	  fftw_real tmp484;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp456;
	       fftw_real tmp459;
	       fftw_real tmp485;
	       fftw_real tmp486;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp456 = tmp454 - tmp455;
	       tmp459 = K1_414213562 * (tmp457 - tmp458);
	       tmp460 = tmp456 + tmp459;
	       tmp488 = tmp456 - tmp459;
	       tmp485 = (K980785280 * tmp474) + (K195090322 * tmp471);
	       tmp486 = (K831469612 * tmp481) + (K555570233 * tmp478);
	       tmp487 = K2_000000000 * (tmp485 - tmp486);
	       tmp492 = tmp485 + tmp486;
	  }
	  {
	       fftw_real tmp475;
	       fftw_real tmp482;
	       fftw_real tmp463;
	       fftw_real tmp466;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp475 = (K980785280 * tmp471) - (K195090322 * tmp474);
	       tmp482 = (K831469612 * tmp478) - (K555570233 * tmp481);
	       tmp483 = K2_000000000 * (tmp475 + tmp482);
	       tmp491 = tmp475 - tmp482;
	       tmp463 = tmp461 - tmp462;
	       tmp466 = tmp464 + tmp465;
	       tmp467 = (K1_847759065 * tmp463) - (K765366864 * tmp466);
	       tmp489 = (K1_847759065 * tmp466) + (K765366864 * tmp463);
	  }
	  tmp468 = tmp460 + tmp467;
	  output[68 * ostride] = tmp468 - tmp483;
	  output[4 * ostride] = tmp468 + tmp483;
	  tmp484 = tmp460 - tmp467;
	  output[36 * ostride] = tmp484 - tmp487;
	  output[100 * ostride] = tmp484 + tmp487;
	  {
	       fftw_real tmp490;
	       fftw_real tmp493;
	       fftw_real tmp494;
	       fftw_real tmp495;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp490 = tmp488 - tmp489;
	       tmp493 = K1_414213562 * (tmp491 - tmp492);
	       output[84 * ostride] = tmp490 - tmp493;
	       output[20 * ostride] = tmp490 + tmp493;
	       tmp494 = tmp488 + tmp489;
	       tmp495 = K1_414213562 * (tmp491 + tmp492);
	       output[52 * ostride] = tmp494 - tmp495;
	       output[116 * ostride] = tmp494 + tmp495;
	  }
     }
     {
	  fftw_real tmp498;
	  fftw_real tmp514;
	  fftw_real tmp513;
	  fftw_real tmp518;
	  fftw_real tmp509;
	  fftw_real tmp517;
	  fftw_real tmp501;
	  fftw_real tmp515;
	  fftw_real tmp502;
	  fftw_real tmp510;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp496;
	       fftw_real tmp497;
	       fftw_real tmp511;
	       fftw_real tmp512;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp496 = tmp454 + tmp455;
	       tmp497 = K1_414213562 * (tmp457 + tmp458);
	       tmp498 = tmp496 - tmp497;
	       tmp514 = tmp496 + tmp497;
	       tmp511 = (K831469612 * tmp504) + (K555570233 * tmp503);
	       tmp512 = (K980785280 * tmp507) - (K195090322 * tmp506);
	       tmp513 = K2_000000000 * (tmp511 - tmp512);
	       tmp518 = tmp511 + tmp512;
	  }
	  {
	       fftw_real tmp505;
	       fftw_real tmp508;
	       fftw_real tmp499;
	       fftw_real tmp500;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp505 = (K831469612 * tmp503) - (K555570233 * tmp504);
	       tmp508 = (K980785280 * tmp506) + (K195090322 * tmp507);
	       tmp509 = K2_000000000 * (tmp505 - tmp508);
	       tmp517 = tmp505 + tmp508;
	       tmp499 = tmp461 + tmp462;
	       tmp500 = tmp465 - tmp464;
	       tmp501 = (K765366864 * tmp499) - (K1_847759065 * tmp500);
	       tmp515 = (K765366864 * tmp500) + (K1_847759065 * tmp499);
	  }
	  tmp502 = tmp498 + tmp501;
	  output[76 * ostride] = tmp502 - tmp509;
	  output[12 * ostride] = tmp502 + tmp509;
	  tmp510 = tmp498 - tmp501;
	  output[44 * ostride] = tmp510 - tmp513;
	  output[108 * ostride] = tmp510 + tmp513;
	  {
	       fftw_real tmp516;
	       fftw_real tmp519;
	       fftw_real tmp520;
	       fftw_real tmp521;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp516 = tmp514 - tmp515;
	       tmp519 = K1_414213562 * (tmp517 - tmp518);
	       output[92 * ostride] = tmp516 - tmp519;
	       output[28 * ostride] = tmp516 + tmp519;
	       tmp520 = tmp514 + tmp515;
	       tmp521 = K1_414213562 * (tmp517 + tmp518);
	       output[60 * ostride] = tmp520 - tmp521;
	       output[124 * ostride] = tmp520 + tmp521;
	  }
     }
     {
	  fftw_real tmp339;
	  fftw_real tmp367;
	  fftw_real tmp332;
	  fftw_real tmp364;
	  fftw_real tmp355;
	  fftw_real tmp369;
	  fftw_real tmp363;
	  fftw_real tmp370;
	  fftw_real tmp328;
	  fftw_real tmp331;
	  fftw_real tmp340;
	  fftw_real tmp356;
	  ASSERT_ALIGNED_DOUBLE;
	  tmp339 = K2_000000000 * (tmp335 - tmp338);
	  tmp367 = K2_000000000 * (tmp365 - tmp366);
	  tmp328 = tmp326 - tmp327;
	  tmp331 = (K765366864 * tmp329) - (K1_847759065 * tmp330);
	  tmp332 = tmp328 + tmp331;
	  tmp364 = tmp328 - tmp331;
	  {
	       fftw_real tmp347;
	       fftw_real tmp354;
	       fftw_real tmp359;
	       fftw_real tmp362;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp347 = tmp343 - tmp346;
	       tmp354 = tmp350 + tmp353;
	       tmp355 = K2_000000000 * (tmp347 + tmp354);
	       tmp369 = tmp347 - tmp354;
	       tmp359 = tmp357 - tmp358;
	       tmp362 = tmp360 - tmp361;
	       tmp363 = K2_000000000 * (tmp359 - tmp362);
	       tmp370 = tmp359 + tmp362;
	  }
	  tmp340 = tmp332 + tmp339;
	  output[70 * ostride] = tmp340 - tmp355;
	  output[6 * ostride] = tmp340 + tmp355;
	  tmp356 = tmp332 - tmp339;
	  output[38 * ostride] = tmp356 - tmp363;
	  output[102 * ostride] = tmp356 + tmp363;
	  {
	       fftw_real tmp368;
	       fftw_real tmp371;
	       fftw_real tmp372;
	       fftw_real tmp373;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp368 = tmp364 - tmp367;
	       tmp371 = K1_414213562 * (tmp369 - tmp370);
	       output[86 * ostride] = tmp368 - tmp371;
	       output[22 * ostride] = tmp368 + tmp371;
	       tmp372 = tmp364 + tmp367;
	       tmp373 = K1_414213562 * (tmp369 + tmp370);
	       output[54 * ostride] = tmp372 - tmp373;
	       output[118 * ostride] = tmp372 + tmp373;
	  }
     }
     {
	  fftw_real tmp376;
	  fftw_real tmp390;
	  fftw_real tmp386;
	  fftw_real tmp394;
	  fftw_real tmp379;
	  fftw_real tmp391;
	  fftw_real tmp383;
	  fftw_real tmp393;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp374;
	       fftw_real tmp375;
	       fftw_real tmp384;
	       fftw_real tmp385;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp374 = tmp326 + tmp327;
	       tmp375 = (K1_847759065 * tmp329) + (K765366864 * tmp330);
	       tmp376 = tmp374 - tmp375;
	       tmp390 = tmp374 + tmp375;
	       tmp384 = tmp357 + tmp358;
	       tmp385 = tmp353 - tmp350;
	       tmp386 = tmp384 + tmp385;
	       tmp394 = tmp384 - tmp385;
	  }
	  {
	       fftw_real tmp377;
	       fftw_real tmp378;
	       fftw_real tmp381;
	       fftw_real tmp382;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp377 = tmp335 + tmp338;
	       tmp378 = tmp365 + tmp366;
	       tmp379 = K1_414213562 * (tmp377 - tmp378);
	       tmp391 = K1_414213562 * (tmp378 + tmp377);
	       tmp381 = tmp343 + tmp346;
	       tmp382 = tmp360 + tmp361;
	       tmp383 = tmp381 - tmp382;
	       tmp393 = tmp381 + tmp382;
	  }
	  {
	       fftw_real tmp380;
	       fftw_real tmp387;
	       fftw_real tmp388;
	       fftw_real tmp389;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp380 = tmp376 + tmp379;
	       tmp387 = (K1_847759065 * tmp383) - (K765366864 * tmp386);
	       output[78 * ostride] = tmp380 - tmp387;
	       output[14 * ostride] = tmp380 + tmp387;
	       tmp388 = tmp376 - tmp379;
	       tmp389 = (K1_847759065 * tmp386) + (K765366864 * tmp383);
	       output[46 * ostride] = tmp388 - tmp389;
	       output[110 * ostride] = tmp388 + tmp389;
	  }
	  {
	       fftw_real tmp392;
	       fftw_real tmp395;
	       fftw_real tmp396;
	       fftw_real tmp397;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp392 = tmp390 - tmp391;
	       tmp395 = (K765366864 * tmp393) - (K1_847759065 * tmp394);
	       output[94 * ostride] = tmp392 - tmp395;
	       output[30 * ostride] = tmp392 + tmp395;
	       tmp396 = tmp390 + tmp391;
	       tmp397 = (K765366864 * tmp394) + (K1_847759065 * tmp393);
	       output[62 * ostride] = tmp396 - tmp397;
	       output[126 * ostride] = tmp396 + tmp397;
	  }
     }
     {
	  fftw_real tmp552;
	  fftw_real tmp672;
	  fftw_real tmp584;
	  fftw_real tmp665;
	  fftw_real tmp648;
	  fftw_real tmp674;
	  fftw_real tmp664;
	  fftw_real tmp675;
	  fftw_real tmp567;
	  fftw_real tmp583;
	  fftw_real tmp585;
	  fftw_real tmp649;
	  ASSERT_ALIGNED_DOUBLE;
	  tmp552 = K2_000000000 * (tmp536 + tmp551);
	  tmp672 = K2_000000000 * (tmp668 - tmp671);
	  tmp567 = K2_000000000 * (tmp559 + tmp566);
	  tmp583 = tmp574 + tmp582;
	  tmp584 = tmp567 + tmp583;
	  tmp665 = tmp583 - tmp567;
	  {
	       fftw_real tmp616;
	       fftw_real tmp647;
	       fftw_real tmp656;
	       fftw_real tmp663;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp616 = tmp600 + tmp615;
	       tmp647 = tmp631 + tmp646;
	       tmp648 = K2_000000000 * (tmp616 + tmp647);
	       tmp674 = tmp616 - tmp647;
	       tmp656 = tmp652 - tmp655;
	       tmp663 = tmp659 - tmp662;
	       tmp664 = K2_000000000 * (tmp656 - tmp663);
	       tmp675 = tmp656 + tmp663;
	  }
	  tmp585 = tmp552 + tmp584;
	  output[65 * ostride] = tmp585 - tmp648;
	  output[ostride] = tmp648 + tmp585;
	  tmp649 = tmp584 - tmp552;
	  output[33 * ostride] = tmp649 - tmp664;
	  output[97 * ostride] = tmp664 + tmp649;
	  {
	       fftw_real tmp673;
	       fftw_real tmp676;
	       fftw_real tmp677;
	       fftw_real tmp678;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp673 = tmp665 - tmp672;
	       tmp676 = K1_414213562 * (tmp674 - tmp675);
	       output[81 * ostride] = tmp673 - tmp676;
	       output[17 * ostride] = tmp676 + tmp673;
	       tmp677 = tmp672 + tmp665;
	       tmp678 = K1_414213562 * (tmp674 + tmp675);
	       output[49 * ostride] = tmp677 - tmp678;
	       output[113 * ostride] = tmp678 + tmp677;
	  }
     }
     {
	  fftw_real tmp681;
	  fftw_real tmp698;
	  fftw_real tmp693;
	  fftw_real tmp701;
	  fftw_real tmp686;
	  fftw_real tmp697;
	  fftw_real tmp690;
	  fftw_real tmp700;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp679;
	       fftw_real tmp680;
	       fftw_real tmp691;
	       fftw_real tmp692;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp679 = tmp536 - tmp551;
	       tmp680 = tmp668 + tmp671;
	       tmp681 = K1_414213562 * (tmp679 - tmp680);
	       tmp698 = K1_414213562 * (tmp679 + tmp680);
	       tmp691 = tmp652 + tmp655;
	       tmp692 = tmp646 - tmp631;
	       tmp693 = tmp691 + tmp692;
	       tmp701 = tmp691 - tmp692;
	  }
	  {
	       fftw_real tmp682;
	       fftw_real tmp685;
	       fftw_real tmp688;
	       fftw_real tmp689;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp682 = tmp582 - tmp574;
	       tmp685 = K2_000000000 * (tmp683 - tmp684);
	       tmp686 = tmp682 - tmp685;
	       tmp697 = tmp685 + tmp682;
	       tmp688 = tmp600 - tmp615;
	       tmp689 = tmp659 + tmp662;
	       tmp690 = tmp688 - tmp689;
	       tmp700 = tmp688 + tmp689;
	  }
	  {
	       fftw_real tmp687;
	       fftw_real tmp694;
	       fftw_real tmp695;
	       fftw_real tmp696;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp687 = tmp681 + tmp686;
	       tmp694 = (K1_847759065 * tmp690) - (K765366864 * tmp693);
	       output[73 * ostride] = tmp687 - tmp694;
	       output[9 * ostride] = tmp694 + tmp687;
	       tmp695 = tmp686 - tmp681;
	       tmp696 = (K1_847759065 * tmp693) + (K765366864 * tmp690);
	       output[41 * ostride] = tmp695 - tmp696;
	       output[105 * ostride] = tmp696 + tmp695;
	  }
	  {
	       fftw_real tmp699;
	       fftw_real tmp702;
	       fftw_real tmp703;
	       fftw_real tmp704;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp699 = tmp697 - tmp698;
	       tmp702 = (K765366864 * tmp700) - (K1_847759065 * tmp701);
	       output[89 * ostride] = tmp699 - tmp702;
	       output[25 * ostride] = tmp702 + tmp699;
	       tmp703 = tmp698 + tmp697;
	       tmp704 = (K765366864 * tmp701) + (K1_847759065 * tmp700);
	       output[57 * ostride] = tmp703 - tmp704;
	       output[121 * ostride] = tmp704 + tmp703;
	  }
     }
     {
	  fftw_real tmp711;
	  fftw_real tmp740;
	  fftw_real tmp718;
	  fftw_real tmp739;
	  fftw_real tmp726;
	  fftw_real tmp736;
	  fftw_real tmp733;
	  fftw_real tmp737;
	  fftw_real tmp742;
	  fftw_real tmp743;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp707;
	       fftw_real tmp710;
	       fftw_real tmp714;
	       fftw_real tmp717;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp707 = tmp705 -