fhb_32.c

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

		    ASSERT_ALIGNED_DOUBLE;
		    tmp200 = X[6 * iostride];
		    tmp201 = Y[-22 * iostride];
		    tmp202 = tmp200 + tmp201;
		    tmp455 = tmp200 - tmp201;
		    tmp242 = Y[-6 * iostride];
		    tmp243 = X[22 * iostride];
		    tmp244 = tmp242 - tmp243;
		    tmp459 = tmp242 + tmp243;
	       }
	       {
		    fftw_real tmp196;
		    fftw_real tmp203;
		    fftw_real tmp355;
		    fftw_real tmp356;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp196 = tmp192 + tmp195;
		    tmp203 = tmp199 + tmp202;
		    tmp204 = tmp196 + tmp203;
		    tmp292 = tmp196 - tmp203;
		    tmp355 = tmp241 - tmp244;
		    tmp356 = tmp199 - tmp202;
		    tmp357 = tmp355 - tmp356;
		    tmp384 = tmp356 + tmp355;
	       }
	       {
		    fftw_real tmp245;
		    fftw_real tmp252;
		    fftw_real tmp358;
		    fftw_real tmp359;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp245 = tmp241 + tmp244;
		    tmp252 = tmp248 + tmp251;
		    tmp253 = tmp245 - tmp252;
		    tmp333 = tmp252 + tmp245;
		    tmp358 = tmp192 - tmp195;
		    tmp359 = tmp248 - tmp251;
		    tmp360 = tmp358 + tmp359;
		    tmp383 = tmp358 - tmp359;
	       }
	       {
		    fftw_real tmp457;
		    fftw_real tmp460;
		    fftw_real tmp586;
		    fftw_real tmp587;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp457 = tmp455 - tmp456;
		    tmp460 = tmp458 - tmp459;
		    tmp461 = (K923879532 * tmp457) - (K382683432 * tmp460);
		    tmp520 = (K382683432 * tmp457) + (K923879532 * tmp460);
		    tmp586 = tmp458 + tmp459;
		    tmp587 = tmp455 + tmp456;
		    tmp588 = (K382683432 * tmp586) - (K923879532 * tmp587);
		    tmp612 = (K382683432 * tmp587) + (K923879532 * tmp586);
	       }
	       {
		    fftw_real tmp464;
		    fftw_real tmp467;
		    fftw_real tmp583;
		    fftw_real tmp584;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp464 = tmp462 + tmp463;
		    tmp467 = tmp465 - tmp466;
		    tmp468 = (K923879532 * tmp464) + (K382683432 * tmp467);
		    tmp519 = (K923879532 * tmp467) - (K382683432 * tmp464);
		    tmp583 = tmp465 + tmp466;
		    tmp584 = tmp463 - tmp462;
		    tmp585 = (K382683432 * tmp583) - (K923879532 * tmp584);
		    tmp611 = (K382683432 * tmp584) + (K923879532 * tmp583);
	       }
	  }
	  {
	       fftw_real tmp208;
	       fftw_real tmp471;
	       fftw_real tmp258;
	       fftw_real tmp483;
	       fftw_real tmp211;
	       fftw_real tmp482;
	       fftw_real tmp261;
	       fftw_real tmp472;
	       fftw_real tmp218;
	       fftw_real tmp486;
	       fftw_real tmp479;
	       fftw_real tmp268;
	       fftw_real tmp215;
	       fftw_real tmp485;
	       fftw_real tmp476;
	       fftw_real tmp265;
	       ASSERT_ALIGNED_DOUBLE;
	       {
		    fftw_real tmp206;
		    fftw_real tmp207;
		    fftw_real tmp259;
		    fftw_real tmp260;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp206 = X[iostride];
		    tmp207 = Y[-17 * iostride];
		    tmp208 = tmp206 + tmp207;
		    tmp471 = tmp206 - tmp207;
		    {
			 fftw_real tmp256;
			 fftw_real tmp257;
			 fftw_real tmp209;
			 fftw_real tmp210;
			 ASSERT_ALIGNED_DOUBLE;
			 tmp256 = Y[-iostride];
			 tmp257 = X[17 * iostride];
			 tmp258 = tmp256 - tmp257;
			 tmp483 = tmp256 + tmp257;
			 tmp209 = X[9 * iostride];
			 tmp210 = Y[-25 * iostride];
			 tmp211 = tmp209 + tmp210;
			 tmp482 = tmp209 - tmp210;
		    }
		    tmp259 = Y[-9 * iostride];
		    tmp260 = X[25 * iostride];
		    tmp261 = tmp259 - tmp260;
		    tmp472 = tmp259 + tmp260;
		    {
			 fftw_real tmp216;
			 fftw_real tmp217;
			 fftw_real tmp477;
			 fftw_real tmp266;
			 fftw_real tmp267;
			 fftw_real tmp478;
			 ASSERT_ALIGNED_DOUBLE;
			 tmp216 = Y[-29 * iostride];
			 tmp217 = X[13 * iostride];
			 tmp477 = tmp216 - tmp217;
			 tmp266 = Y[-13 * iostride];
			 tmp267 = X[29 * iostride];
			 tmp478 = tmp266 + tmp267;
			 tmp218 = tmp216 + tmp217;
			 tmp486 = tmp477 + tmp478;
			 tmp479 = tmp477 - tmp478;
			 tmp268 = tmp266 - tmp267;
		    }
		    {
			 fftw_real tmp213;
			 fftw_real tmp214;
			 fftw_real tmp474;
			 fftw_real tmp263;
			 fftw_real tmp264;
			 fftw_real tmp475;
			 ASSERT_ALIGNED_DOUBLE;
			 tmp213 = X[5 * iostride];
			 tmp214 = Y[-21 * iostride];
			 tmp474 = tmp213 - tmp214;
			 tmp263 = Y[-5 * iostride];
			 tmp264 = X[21 * iostride];
			 tmp475 = tmp263 + tmp264;
			 tmp215 = tmp213 + tmp214;
			 tmp485 = tmp474 + tmp475;
			 tmp476 = tmp474 - tmp475;
			 tmp265 = tmp263 - tmp264;
		    }
	       }
	       {
		    fftw_real tmp212;
		    fftw_real tmp219;
		    fftw_real tmp370;
		    fftw_real tmp371;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp212 = tmp208 + tmp211;
		    tmp219 = tmp215 + tmp218;
		    tmp220 = tmp212 + tmp219;
		    tmp255 = tmp212 - tmp219;
		    tmp370 = tmp258 - tmp261;
		    tmp371 = tmp215 - tmp218;
		    tmp372 = tmp370 - tmp371;
		    tmp414 = tmp371 + tmp370;
	       }
	       {
		    fftw_real tmp262;
		    fftw_real tmp269;
		    fftw_real tmp373;
		    fftw_real tmp374;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp262 = tmp258 + tmp261;
		    tmp269 = tmp265 + tmp268;
		    tmp270 = tmp262 - tmp269;
		    tmp335 = tmp262 + tmp269;
		    tmp373 = tmp208 - tmp211;
		    tmp374 = tmp268 - tmp265;
		    tmp375 = tmp373 - tmp374;
		    tmp415 = tmp373 + tmp374;
	       }
	       {
		    fftw_real tmp473;
		    fftw_real tmp480;
		    fftw_real tmp598;
		    fftw_real tmp599;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp473 = tmp471 - tmp472;
		    tmp480 = K707106781 * (tmp476 + tmp479);
		    tmp481 = tmp473 - tmp480;
		    tmp547 = tmp473 + tmp480;
		    tmp598 = tmp483 - tmp482;
		    tmp599 = K707106781 * (tmp476 - tmp479);
		    tmp600 = tmp598 + tmp599;
		    tmp640 = tmp598 - tmp599;
	       }
	       {
		    fftw_real tmp484;
		    fftw_real tmp487;
		    fftw_real tmp601;
		    fftw_real tmp602;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp484 = tmp482 + tmp483;
		    tmp487 = K707106781 * (tmp485 - tmp486);
		    tmp488 = tmp484 - tmp487;
		    tmp548 = tmp484 + tmp487;
		    tmp601 = tmp471 + tmp472;
		    tmp602 = K707106781 * (tmp485 + tmp486);
		    tmp603 = tmp601 - tmp602;
		    tmp639 = tmp601 + tmp602;
	       }
	  }
	  {
	       fftw_real tmp205;
	       fftw_real tmp236;
	       fftw_real tmp330;
	       fftw_real tmp334;
	       fftw_real tmp337;
	       fftw_real tmp338;
	       fftw_real tmp329;
	       fftw_real tmp331;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp205 = tmp189 + tmp204;
	       tmp236 = tmp220 + tmp235;
	       tmp330 = tmp205 - tmp236;
	       tmp334 = tmp332 + tmp333;
	       tmp337 = tmp335 + tmp336;
	       tmp338 = tmp334 - tmp337;
	       X[0] = tmp205 + tmp236;
	       Y[-31 * iostride] = tmp334 + tmp337;
	       tmp329 = c_re(W[15]);
	       tmp331 = c_im(W[15]);
	       X[16 * iostride] = (tmp329 * tmp330) + (tmp331 * tmp338);
	       Y[-15 * iostride] = (tmp329 * tmp338) - (tmp331 * tmp330);
	  }
	  {
	       fftw_real tmp342;
	       fftw_real tmp348;
	       fftw_real tmp346;
	       fftw_real tmp350;
	       ASSERT_ALIGNED_DOUBLE;
	       {
		    fftw_real tmp340;
		    fftw_real tmp341;
		    fftw_real tmp344;
		    fftw_real tmp345;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp340 = tmp332 - tmp333;
		    tmp341 = tmp220 - tmp235;
		    tmp342 = tmp340 - tmp341;
		    tmp348 = tmp341 + tmp340;
		    tmp344 = tmp189 - tmp204;
		    tmp345 = tmp336 - tmp335;
		    tmp346 = tmp344 - tmp345;
		    tmp350 = tmp344 + tmp345;
	       }
	       {
		    fftw_real tmp339;
		    fftw_real tmp343;
		    fftw_real tmp347;
		    fftw_real tmp349;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp339 = c_re(W[23]);
		    tmp343 = c_im(W[23]);
		    Y[-7 * iostride] =
			(tmp339 * tmp342) - (tmp343 * tmp346);
		    X[24 * iostride] =
			(tmp343 * tmp342) + (tmp339 * tmp346);
		    tmp347 = c_re(W[7]);
		    tmp349 = c_im(W[7]);
		    Y[-23 * iostride] =
			(tmp347 * tmp348) - (tmp349 * tmp350);
		    X[8 * iostride] =
			(tmp349 * tmp348) + (tmp347 * tmp350);
	       }
	  }
	  {
	       fftw_real tmp254;
	       fftw_real tmp322;
	       fftw_real tmp289;
	       fftw_real tmp319;
	       fftw_real tmp308;
	       fftw_real tmp318;
	       fftw_real tmp311;
	       fftw_real tmp323;
	       ASSERT_ALIGNED_DOUBLE;
	       {
		    fftw_real tmp271;
		    fftw_real tmp288;
		    fftw_real tmp309;
		    fftw_real tmp310;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp254 = tmp238 + tmp253;
		    tmp322 = tmp238 - tmp253;
		    tmp271 = tmp255 - tmp270;
		    tmp288 = tmp272 + tmp287;
		    tmp289 = K707106781 * (tmp271 + tmp288);
		    tmp319 = K707106781 * (tmp271 - tmp288);
		    tmp308 = tmp292 + tmp307;
		    tmp318 = tmp307 - tmp292;
		    tmp309 = tmp255 + tmp270;
		    tmp310 = tmp287 - tmp272;
		    tmp311 = K707106781 * (tmp309 + tmp310);
		    tmp323 = K707106781 * (tmp310 - tmp309);
	       }
	       {
		    fftw_real tmp290;
		    fftw_real tmp312;
		    fftw_real tmp237;
		    fftw_real tmp291;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp290 = tmp254 - tmp289;
		    tmp312 = tmp308 - tmp311;
		    tmp237 = c_re(W[19]);
		    tmp291 = c_im(W[19]);
		    X[20 * iostride] =
			(tmp237 * tmp290) + (tmp291 * tmp312);
		    Y[-11 * iostride] =
			(tmp237 * tmp312) - (tmp291 * tmp290);
	       }
	       {
		    fftw_real tmp314;
		    fftw_real tmp316;
		    fftw_real tmp313;
		    fftw_real tmp315;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp314 = tmp254 + tmp289;
		    tmp316 = tmp308 + tmp311;
		    tmp313 = c_re(W[3]);
		    tmp315 = c_im(W[3]);
		    X[4 * iostride] =
			(tmp313 * tmp314) + (tmp315 * tmp316);
		    Y[-27 * iostride] =
			(tmp313 * tmp316) - (tmp315 * tmp314);
	       }
	       {
		    fftw_real tmp320;
		    fftw_real tmp324;
		    fftw_real tmp317;
		    fftw_real tmp321;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp320 = tmp318 - tmp319;
		    tmp324 = tmp322 - tmp323;
		    tmp317 = c_re(W[27]);
		    tmp321 = c_im(W[27]);
		    Y[-3 * iostride] =
			(tmp317 * tmp320) - (tmp321 * tmp324);
		    X[28 * iostride] =
			(tmp321 * tmp320) + (tmp317 * tmp324);
	       }
	       {
		    fftw_real tmp326;
		    fftw_real tmp328;
		    fftw_real tmp325;
		    fftw_real tmp327;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp326 = tmp318 + tmp319;
		    tmp328 = tmp322 + tmp323;
		    tmp325 = c_re(W[11]);
		    tmp327 = c_im(W[11]);
		    Y[-19 * iostride] =
			(tmp325 * tmp326) - (tmp327 * tmp328);
		    X[12 * iostride] =
			(tmp327 * tmp326) + (tmp325 * tmp328);
	       }
	  }
	  {
	       fftw_real tmp638;
	       fftw_real tmp664;
	       fftw_real tmp653;
	       fftw_real tmp665;
	       fftw_real tmp645;
	       fftw_real tmp661;
	       fftw_real tmp650;
	       fftw_real tmp660;
	       ASSERT_ALIGNED_DOUBLE;
	       {
		    fftw_real tmp636;
		    fftw_real tmp637;
		    fftw_real tmp651;
		    fftw_real tmp652;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp636 = tmp580 + tmp581;
		    tmp637 = tmp611 + tmp612;
		    tmp638 = tmp636 - tmp637;
		    tmp664 = tmp636 + tmp637;
		    tmp651 = (K195090322 * tmp640) + (K980785280 * tmp639);
		    tmp652 = (K195090322 * tmp643) + (K980785280 * tmp642);
		    tmp653 = tmp651 - tmp652;
		    tmp665 = tmp651 + tmp652;
	       }
	       {
		    fftw_real tmp641;
		    fftw_real tmp644;
		    fftw_real tmp648;
		    fftw_real tmp649;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp641 = (K195090322 * tmp639) - (K980785280 * tmp640);
		    tmp644 = (K195090322 * tmp642) - (K980785280 * tmp643);
		    tmp645 = tmp641 + tmp644;
		    tmp661 = tmp641 - tmp644;
		    tmp648 = tmp608 - tmp609;
		    tmp649 = tmp585 - tmp588;
		    tmp650 = tmp648 + tmp649;
		    tmp660 = tmp648 - tmp649;
	       }
	       {
		    fftw_real tmp646;
		    fftw_real tmp654;
		    fftw_real tmp635;
		    fftw_real tmp647;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp646 = tmp638 - tmp645;
		    tmp654 = tmp650 - tmp653;
		    tmp635 = c_re(W[22]);
		    tmp647 = c_im(W[22]);
		    X[23 * iostride] =
			(tmp635 * tmp646) + (tmp647 * tmp654);
		    Y[-8 * iostride] =
			(tmp635 * tmp654) - (tmp647 * tmp646);
	       }
	       {
		    fftw_real tmp656;
		    fftw_real tmp658;
		    fftw_real tmp655;
		    fftw_real tmp657;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp656 = tmp638 + tmp645;
		    tmp658 = tmp650 + tmp653;
		    tmp655 = c_re(W[6]);
		    tmp657 = c_im(W[6]);
		    X[7 * iostride] =
			(tmp655 * tmp656) + (tmp657 * tmp658);
		    Y[-24 * iostride] =
			(tmp655 * tmp658) - (tmp657 * tmp656);
	       }
	       {
		    fftw_real tmp662;
		    fftw_real tmp666;
		    fftw_real tmp659;
		    fftw_real tmp663;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp662 = tmp660 + tmp661;
		    tmp666 = tmp664 - tmp665;
		    tmp659 = c_re(W[14]);
		    tmp663 = c_im(W[14]);
		    Y[-16 * iostride] =
			(tmp659 * tmp662) - (tmp663 * tmp666);
		    X[15 * iostride] =
			(tmp663 * tmp662) + (tmp659 * tmp666);
	       }
	       {
		    fftw_real tmp668;
		    fftw_real tmp670;
		    fftw_real tmp667;
		    fftw_real tmp669;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp668 = tmp660 - tmp661;
		    tmp670 = tmp664 + tmp665;
		    tmp667 = c_re(W[30]);
		    tmp669 = c_im(W[30]);
		    Y[0] = (tmp667 * tmp668) - (tmp669 * tmp670);
		    X[31 * iostride] =
			(tmp669 * tmp668) + (tmp667 * tmp670);