fhb_32.c

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

	       }
	  }
	  {
	       fftw_real tmp410;
	       fftw_real tmp436;
	       fftw_real tmp425;
	       fftw_real tmp437;
	       fftw_real tmp422;
	       fftw_real tmp432;
	       fftw_real tmp417;
	       fftw_real tmp433;
	       ASSERT_ALIGNED_DOUBLE;
	       {
		    fftw_real tmp409;
		    fftw_real tmp423;
		    fftw_real tmp424;
		    fftw_real tmp421;
		    fftw_real tmp413;
		    fftw_real tmp416;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp409 = K707106781 * (tmp383 + tmp384);
		    tmp410 = tmp408 - tmp409;
		    tmp436 = tmp408 + tmp409;
		    tmp423 = (K923879532 * tmp415) - (K382683432 * tmp414);
		    tmp424 = (K382683432 * tmp411) + (K923879532 * tmp412);
		    tmp425 = tmp423 - tmp424;
		    tmp437 = tmp423 + tmp424;
		    tmp421 = K707106781 * (tmp360 + tmp357);
		    tmp422 = tmp420 - tmp421;
		    tmp432 = tmp420 + tmp421;
		    tmp413 = (K923879532 * tmp411) - (K382683432 * tmp412);
		    tmp416 = (K923879532 * tmp414) + (K382683432 * tmp415);
		    tmp417 = tmp413 - tmp416;
		    tmp433 = tmp416 + tmp413;
	       }
	       {
		    fftw_real tmp418;
		    fftw_real tmp426;
		    fftw_real tmp407;
		    fftw_real tmp419;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp418 = tmp410 - tmp417;
		    tmp426 = tmp422 - tmp425;
		    tmp407 = c_re(W[25]);
		    tmp419 = c_im(W[25]);
		    X[26 * iostride] =
			(tmp407 * tmp418) + (tmp419 * tmp426);
		    Y[-5 * iostride] =
			(tmp407 * tmp426) - (tmp419 * tmp418);
	       }
	       {
		    fftw_real tmp428;
		    fftw_real tmp430;
		    fftw_real tmp427;
		    fftw_real tmp429;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp428 = tmp410 + tmp417;
		    tmp430 = tmp422 + tmp425;
		    tmp427 = c_re(W[9]);
		    tmp429 = c_im(W[9]);
		    X[10 * iostride] =
			(tmp427 * tmp428) + (tmp429 * tmp430);
		    Y[-21 * iostride] =
			(tmp427 * tmp430) - (tmp429 * tmp428);
	       }
	       {
		    fftw_real tmp434;
		    fftw_real tmp438;
		    fftw_real tmp431;
		    fftw_real tmp435;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp434 = tmp432 - tmp433;
		    tmp438 = tmp436 - tmp437;
		    tmp431 = c_re(W[17]);
		    tmp435 = c_im(W[17]);
		    Y[-13 * iostride] =
			(tmp431 * tmp434) - (tmp435 * tmp438);
		    X[18 * iostride] =
			(tmp435 * tmp434) + (tmp431 * tmp438);
	       }
	       {
		    fftw_real tmp440;
		    fftw_real tmp442;
		    fftw_real tmp439;
		    fftw_real tmp441;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp440 = tmp432 + tmp433;
		    tmp442 = tmp436 + tmp437;
		    tmp439 = c_re(W[1]);
		    tmp441 = c_im(W[1]);
		    Y[-29 * iostride] =
			(tmp439 * tmp440) - (tmp441 * tmp442);
		    X[2 * iostride] =
			(tmp441 * tmp440) + (tmp439 * tmp442);
	       }
	  }
	  {
	       fftw_real tmp362;
	       fftw_real tmp400;
	       fftw_real tmp389;
	       fftw_real tmp401;
	       fftw_real tmp386;
	       fftw_real tmp396;
	       fftw_real tmp377;
	       fftw_real tmp397;
	       ASSERT_ALIGNED_DOUBLE;
	       {
		    fftw_real tmp361;
		    fftw_real tmp387;
		    fftw_real tmp388;
		    fftw_real tmp385;
		    fftw_real tmp369;
		    fftw_real tmp376;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp361 = K707106781 * (tmp357 - tmp360);
		    tmp362 = tmp354 - tmp361;
		    tmp400 = tmp354 + tmp361;
		    tmp387 = (K382683432 * tmp375) - (K923879532 * tmp372);
		    tmp388 = (K923879532 * tmp365) + (K382683432 * tmp368);
		    tmp389 = tmp387 - tmp388;
		    tmp401 = tmp387 + tmp388;
		    tmp385 = K707106781 * (tmp383 - tmp384);
		    tmp386 = tmp382 - tmp385;
		    tmp396 = tmp382 + tmp385;
		    tmp369 = (K382683432 * tmp365) - (K923879532 * tmp368);
		    tmp376 = (K382683432 * tmp372) + (K923879532 * tmp375);
		    tmp377 = tmp369 - tmp376;
		    tmp397 = tmp376 + tmp369;
	       }
	       {
		    fftw_real tmp378;
		    fftw_real tmp390;
		    fftw_real tmp351;
		    fftw_real tmp379;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp378 = tmp362 - tmp377;
		    tmp390 = tmp386 - tmp389;
		    tmp351 = c_re(W[29]);
		    tmp379 = c_im(W[29]);
		    X[30 * iostride] =
			(tmp351 * tmp378) + (tmp379 * tmp390);
		    Y[-iostride] = (tmp351 * tmp390) - (tmp379 * tmp378);
	       }
	       {
		    fftw_real tmp392;
		    fftw_real tmp394;
		    fftw_real tmp391;
		    fftw_real tmp393;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp392 = tmp362 + tmp377;
		    tmp394 = tmp386 + tmp389;
		    tmp391 = c_re(W[13]);
		    tmp393 = c_im(W[13]);
		    X[14 * iostride] =
			(tmp391 * tmp392) + (tmp393 * tmp394);
		    Y[-17 * iostride] =
			(tmp391 * tmp394) - (tmp393 * tmp392);
	       }
	       {
		    fftw_real tmp398;
		    fftw_real tmp402;
		    fftw_real tmp395;
		    fftw_real tmp399;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp398 = tmp396 - tmp397;
		    tmp402 = tmp400 - tmp401;
		    tmp395 = c_re(W[21]);
		    tmp399 = c_im(W[21]);
		    Y[-9 * iostride] =
			(tmp395 * tmp398) - (tmp399 * tmp402);
		    X[22 * iostride] =
			(tmp399 * tmp398) + (tmp395 * tmp402);
	       }
	       {
		    fftw_real tmp404;
		    fftw_real tmp406;
		    fftw_real tmp403;
		    fftw_real tmp405;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp404 = tmp396 + tmp397;
		    tmp406 = tmp400 + tmp401;
		    tmp403 = c_re(W[5]);
		    tmp405 = c_im(W[5]);
		    Y[-25 * iostride] =
			(tmp403 * tmp404) - (tmp405 * tmp406);
		    X[6 * iostride] =
			(tmp405 * tmp404) + (tmp403 * tmp406);
	       }
	  }
	  {
	       fftw_real tmp590;
	       fftw_real tmp628;
	       fftw_real tmp617;
	       fftw_real tmp629;
	       fftw_real tmp605;
	       fftw_real tmp625;
	       fftw_real tmp614;
	       fftw_real tmp624;
	       ASSERT_ALIGNED_DOUBLE;
	       {
		    fftw_real tmp582;
		    fftw_real tmp589;
		    fftw_real tmp615;
		    fftw_real tmp616;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp582 = tmp580 - tmp581;
		    tmp589 = tmp585 + tmp588;
		    tmp590 = tmp582 - tmp589;
		    tmp628 = tmp582 + tmp589;
		    tmp615 = (K831469612 * tmp603) - (K555570233 * tmp600);
		    tmp616 = (K555570233 * tmp593) + (K831469612 * tmp596);
		    tmp617 = tmp615 - tmp616;
		    tmp629 = tmp615 + tmp616;
	       }
	       {
		    fftw_real tmp597;
		    fftw_real tmp604;
		    fftw_real tmp610;
		    fftw_real tmp613;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp597 = (K831469612 * tmp593) - (K555570233 * tmp596);
		    tmp604 = (K831469612 * tmp600) + (K555570233 * tmp603);
		    tmp605 = tmp597 - tmp604;
		    tmp625 = tmp604 + tmp597;
		    tmp610 = tmp608 + tmp609;
		    tmp613 = tmp611 - tmp612;
		    tmp614 = tmp610 - tmp613;
		    tmp624 = tmp610 + tmp613;
	       }
	       {
		    fftw_real tmp606;
		    fftw_real tmp618;
		    fftw_real tmp579;
		    fftw_real tmp607;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp606 = tmp590 - tmp605;
		    tmp618 = tmp614 - tmp617;
		    tmp579 = c_re(W[26]);
		    tmp607 = c_im(W[26]);
		    X[27 * iostride] =
			(tmp579 * tmp606) + (tmp607 * tmp618);
		    Y[-4 * iostride] =
			(tmp579 * tmp618) - (tmp607 * tmp606);
	       }
	       {
		    fftw_real tmp620;
		    fftw_real tmp622;
		    fftw_real tmp619;
		    fftw_real tmp621;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp620 = tmp590 + tmp605;
		    tmp622 = tmp614 + tmp617;
		    tmp619 = c_re(W[10]);
		    tmp621 = c_im(W[10]);
		    X[11 * iostride] =
			(tmp619 * tmp620) + (tmp621 * tmp622);
		    Y[-20 * iostride] =
			(tmp619 * tmp622) - (tmp621 * tmp620);
	       }
	       {
		    fftw_real tmp626;
		    fftw_real tmp630;
		    fftw_real tmp623;
		    fftw_real tmp627;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp626 = tmp624 - tmp625;
		    tmp630 = tmp628 - tmp629;
		    tmp623 = c_re(W[18]);
		    tmp627 = c_im(W[18]);
		    Y[-12 * iostride] =
			(tmp623 * tmp626) - (tmp627 * tmp630);
		    X[19 * iostride] =
			(tmp627 * tmp626) + (tmp623 * tmp630);
	       }
	       {
		    fftw_real tmp632;
		    fftw_real tmp634;
		    fftw_real tmp631;
		    fftw_real tmp633;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp632 = tmp624 + tmp625;
		    tmp634 = tmp628 + tmp629;
		    tmp631 = c_re(W[2]);
		    tmp633 = c_im(W[2]);
		    Y[-28 * iostride] =
			(tmp631 * tmp632) - (tmp633 * tmp634);
		    X[3 * iostride] =
			(tmp633 * tmp632) + (tmp631 * tmp634);
	       }
	  }
	  {
	       fftw_real tmp546;
	       fftw_real tmp572;
	       fftw_real tmp561;
	       fftw_real tmp573;
	       fftw_real tmp553;
	       fftw_real tmp569;
	       fftw_real tmp558;
	       fftw_real tmp568;
	       ASSERT_ALIGNED_DOUBLE;
	       {
		    fftw_real tmp544;
		    fftw_real tmp545;
		    fftw_real tmp559;
		    fftw_real tmp560;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp544 = tmp446 + tmp453;
		    tmp545 = tmp519 + tmp520;
		    tmp546 = tmp544 + tmp545;
		    tmp572 = tmp544 - tmp545;
		    tmp559 = (K195090322 * tmp547) + (K980785280 * tmp548);
		    tmp560 = (K980785280 * tmp551) - (K195090322 * tmp550);
		    tmp561 = tmp559 + tmp560;
		    tmp573 = tmp560 - tmp559;
	       }
	       {
		    fftw_real tmp549;
		    fftw_real tmp552;
		    fftw_real tmp556;
		    fftw_real tmp557;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp549 = (K980785280 * tmp547) - (K195090322 * tmp548);
		    tmp552 = (K980785280 * tmp550) + (K195090322 * tmp551);
		    tmp553 = tmp549 + tmp552;
		    tmp569 = tmp549 - tmp552;
		    tmp556 = tmp514 + tmp517;
		    tmp557 = tmp468 + tmp461;
		    tmp558 = tmp556 + tmp557;
		    tmp568 = tmp556 - tmp557;
	       }
	       {
		    fftw_real tmp554;
		    fftw_real tmp562;
		    fftw_real tmp543;
		    fftw_real tmp555;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp554 = tmp546 - tmp553;
		    tmp562 = tmp558 - tmp561;
		    tmp543 = c_re(W[16]);
		    tmp555 = c_im(W[16]);
		    X[17 * iostride] =
			(tmp543 * tmp554) + (tmp555 * tmp562);
		    Y[-14 * iostride] =
			(tmp543 * tmp562) - (tmp555 * tmp554);
	       }
	       {
		    fftw_real tmp564;
		    fftw_real tmp566;
		    fftw_real tmp563;
		    fftw_real tmp565;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp564 = tmp546 + tmp553;
		    tmp566 = tmp558 + tmp561;
		    tmp563 = c_re(W[0]);
		    tmp565 = c_im(W[0]);
		    X[iostride] = (tmp563 * tmp564) + (tmp565 * tmp566);
		    Y[-30 * iostride] =
			(tmp563 * tmp566) - (tmp565 * tmp564);
	       }
	       {
		    fftw_real tmp570;
		    fftw_real tmp574;
		    fftw_real tmp567;
		    fftw_real tmp571;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp570 = tmp568 - tmp569;
		    tmp574 = tmp572 - tmp573;
		    tmp567 = c_re(W[24]);
		    tmp571 = c_im(W[24]);
		    Y[-6 * iostride] =
			(tmp567 * tmp570) - (tmp571 * tmp574);
		    X[25 * iostride] =
			(tmp571 * tmp570) + (tmp567 * tmp574);
	       }
	       {
		    fftw_real tmp576;
		    fftw_real tmp578;
		    fftw_real tmp575;
		    fftw_real tmp577;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp576 = tmp568 + tmp569;
		    tmp578 = tmp572 + tmp573;
		    tmp575 = c_re(W[8]);
		    tmp577 = c_im(W[8]);
		    Y[-22 * iostride] =
			(tmp575 * tmp576) - (tmp577 * tmp578);
		    X[9 * iostride] =
			(tmp577 * tmp576) + (tmp575 * tmp578);
	       }
	  }
	  {
	       fftw_real tmp470;
	       fftw_real tmp536;
	       fftw_real tmp525;
	       fftw_real tmp537;
	       fftw_real tmp509;
	       fftw_real tmp533;
	       fftw_real tmp522;
	       fftw_real tmp532;
	       ASSERT_ALIGNED_DOUBLE;
	       {
		    fftw_real tmp454;
		    fftw_real tmp469;
		    fftw_real tmp523;
		    fftw_real tmp524;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp454 = tmp446 - tmp453;
		    tmp469 = tmp461 - tmp468;
		    tmp470 = tmp454 + tmp469;
		    tmp536 = tmp454 - tmp469;
		    tmp523 = (K831469612 * tmp481) + (K555570233 * tmp488);
		    tmp524 = (K555570233 * tmp507) - (K831469612 * tmp500);
		    tmp525 = tmp523 + tmp524;
		    tmp537 = tmp524 - tmp523;
	       }
	       {
		    fftw_real tmp489;
		    fftw_real tmp508;
		    fftw_real tmp518;
		    fftw_real tmp521;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp489 = (K555570233 * tmp481) - (K831469612 * tmp488);
		    tmp508 = (K555570233 * tmp500) + (K831469612 * tmp507);
		    tmp509 = tmp489 + tmp508;
		    tmp533 = tmp489 - tmp508;
		    tmp518 = tmp514 - tmp517;
		    tmp521 = tmp519 - tmp520;
		    tmp522 = tmp518 + tmp521;
		    tmp532 = tmp518 - tmp521;
	       }
	       {
		    fftw_real tmp510;
		    fftw_real tmp526;
		    fftw_real tmp443;
		    fftw_real tmp511;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp510 = tmp470 - tmp509;
		    tmp526 = tmp522 - tmp525;
		    tmp443 = c_re(W[20]);
		    tmp511 = c_im(W[20]);
		    X[21 * iostride] =
			(tmp443 * tmp510) + (tmp511 * tmp526);
		    Y[-10 * iostride] =
			(tmp443 * tmp526) - (tmp511 * tmp510);
	       }
	       {
		    fftw_real tmp528;
		    fftw_real tmp530;
		    fftw_real tmp527;
		    fftw_real tmp529;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp528 = tmp470 + tmp509;
		    tmp530 = tmp522 + tmp525;
		    tmp527 = c_re(W[4]);
		    tmp529 = c_im(W[4]);
		    X[5 * iostride] =
			(tmp527 * tmp528) + (tmp529 * tmp530);
		    Y[-26 * iostride] =
			(tmp527 * tmp530) - (tmp529 * tmp528);
	       }
	       {
		    fftw_real tmp534;
		    fftw_real tmp538;
		    fftw_real tmp531;
		    fftw_real tmp535;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp534 = tmp532 - tmp533;
		    tmp538 = tmp536 - tmp537;
		    tmp531 = c_re(W[28]);
		    tmp535 = c_im(W[28]);