ftwi_32.c

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

	       {
		    fftw_real tmp351;
		    fftw_real tmp358;
		    fftw_real tmp448;
		    fftw_real tmp449;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp351 = tmp19 - tmp42;
		    tmp358 = tmp354 - tmp357;
		    tmp359 = tmp351 + tmp358;
		    tmp379 = tmp351 - tmp358;
		    tmp448 = tmp66 - tmp89;
		    tmp449 = tmp442 - tmp437;
		    tmp450 = tmp448 + tmp449;
		    tmp452 = tmp449 - tmp448;
	       }
	       {
		    fftw_real tmp360;
		    fftw_real tmp367;
		    fftw_real tmp369;
		    fftw_real tmp376;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp360 = tmp114 - tmp137;
		    tmp367 = tmp363 - tmp366;
		    tmp368 = tmp360 - tmp367;
		    tmp381 = tmp360 + tmp367;
		    tmp369 = tmp161 - tmp184;
		    tmp376 = tmp372 - tmp375;
		    tmp377 = tmp369 + tmp376;
		    tmp380 = tmp376 - tmp369;
	       }
	       {
		    fftw_real tmp378;
		    fftw_real tmp451;
		    fftw_real tmp382;
		    fftw_real tmp447;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp378 = K707106781 * (tmp368 + tmp377);
		    c_re(inout[20 * iostride]) = tmp359 - tmp378;
		    c_re(inout[4 * iostride]) = tmp359 + tmp378;
		    tmp451 = K707106781 * (tmp368 - tmp377);
		    c_im(inout[12 * iostride]) = tmp451 + tmp452;
		    c_im(inout[28 * iostride]) = tmp452 - tmp451;
		    tmp382 = K707106781 * (tmp380 - tmp381);
		    c_re(inout[28 * iostride]) = tmp379 - tmp382;
		    c_re(inout[12 * iostride]) = tmp379 + tmp382;
		    tmp447 = K707106781 * (tmp381 + tmp380);
		    c_im(inout[4 * iostride]) = tmp447 + tmp450;
		    c_im(inout[20 * iostride]) = tmp450 - tmp447;
	       }
	  }
	  {
	       fftw_real tmp391;
	       fftw_real tmp419;
	       fftw_real tmp398;
	       fftw_real tmp454;
	       fftw_real tmp422;
	       fftw_real tmp462;
	       fftw_real tmp406;
	       fftw_real tmp417;
	       fftw_real tmp457;
	       fftw_real tmp463;
	       fftw_real tmp426;
	       fftw_real tmp433;
	       fftw_real tmp413;
	       fftw_real tmp416;
	       fftw_real tmp429;
	       fftw_real tmp432;
	       ASSERT_ALIGNED_DOUBLE;
	       {
		    fftw_real tmp394;
		    fftw_real tmp397;
		    fftw_real tmp424;
		    fftw_real tmp425;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp391 = tmp387 - tmp390;
		    tmp419 = tmp387 + tmp390;
		    tmp394 = tmp392 - tmp393;
		    tmp397 = tmp395 + tmp396;
		    tmp398 = K707106781 * (tmp394 - tmp397);
		    tmp454 = K707106781 * (tmp397 + tmp394);
		    {
			 fftw_real tmp420;
			 fftw_real tmp421;
			 fftw_real tmp402;
			 fftw_real tmp405;
			 ASSERT_ALIGNED_DOUBLE;
			 tmp420 = tmp395 - tmp396;
			 tmp421 = tmp393 + tmp392;
			 tmp422 = K707106781 * (tmp420 + tmp421);
			 tmp462 = K707106781 * (tmp420 - tmp421);
			 tmp402 = tmp400 - tmp401;
			 tmp405 = tmp403 - tmp404;
			 tmp406 =
			     (K382683432 * tmp402) - (K923879532 * tmp405);
			 tmp417 =
			     (K923879532 * tmp402) + (K382683432 * tmp405);
		    }
		    tmp457 = tmp455 + tmp456;
		    tmp463 = tmp456 - tmp455;
		    tmp424 = tmp400 + tmp401;
		    tmp425 = tmp403 + tmp404;
		    tmp426 = (K923879532 * tmp424) - (K382683432 * tmp425);
		    tmp433 = (K382683432 * tmp424) + (K923879532 * tmp425);
		    {
			 fftw_real tmp409;
			 fftw_real tmp412;
			 fftw_real tmp427;
			 fftw_real tmp428;
			 ASSERT_ALIGNED_DOUBLE;
			 tmp409 = tmp407 - tmp408;
			 tmp412 = tmp410 - tmp411;
			 tmp413 =
			     (K382683432 * tmp409) + (K923879532 * tmp412);
			 tmp416 =
			     (K382683432 * tmp412) - (K923879532 * tmp409);
			 tmp427 = tmp407 + tmp408;
			 tmp428 = tmp410 + tmp411;
			 tmp429 =
			     (K923879532 * tmp427) + (K382683432 * tmp428);
			 tmp432 =
			     (K923879532 * tmp428) - (K382683432 * tmp427);
		    }
	       }
	       {
		    fftw_real tmp399;
		    fftw_real tmp414;
		    fftw_real tmp415;
		    fftw_real tmp418;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp399 = tmp391 - tmp398;
		    tmp414 = tmp406 - tmp413;
		    c_re(inout[30 * iostride]) = tmp399 - tmp414;
		    c_re(inout[14 * iostride]) = tmp399 + tmp414;
		    tmp415 = tmp391 + tmp398;
		    tmp418 = tmp416 + tmp417;
		    c_re(inout[22 * iostride]) = tmp415 - tmp418;
		    c_re(inout[6 * iostride]) = tmp415 + tmp418;
	       }
	       {
		    fftw_real tmp465;
		    fftw_real tmp466;
		    fftw_real tmp461;
		    fftw_real tmp464;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp465 = tmp416 - tmp417;
		    tmp466 = tmp463 - tmp462;
		    c_im(inout[14 * iostride]) = tmp465 + tmp466;
		    c_im(inout[30 * iostride]) = tmp466 - tmp465;
		    tmp461 = tmp413 + tmp406;
		    tmp464 = tmp462 + tmp463;
		    c_im(inout[6 * iostride]) = tmp461 + tmp464;
		    c_im(inout[22 * iostride]) = tmp464 - tmp461;
	       }
	       {
		    fftw_real tmp423;
		    fftw_real tmp430;
		    fftw_real tmp431;
		    fftw_real tmp434;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp423 = tmp419 - tmp422;
		    tmp430 = tmp426 - tmp429;
		    c_re(inout[26 * iostride]) = tmp423 - tmp430;
		    c_re(inout[10 * iostride]) = tmp423 + tmp430;
		    tmp431 = tmp419 + tmp422;
		    tmp434 = tmp432 + tmp433;
		    c_re(inout[18 * iostride]) = tmp431 - tmp434;
		    c_re(inout[2 * iostride]) = tmp431 + tmp434;
	       }
	       {
		    fftw_real tmp459;
		    fftw_real tmp460;
		    fftw_real tmp453;
		    fftw_real tmp458;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp459 = tmp432 - tmp433;
		    tmp460 = tmp457 - tmp454;
		    c_im(inout[10 * iostride]) = tmp459 + tmp460;
		    c_im(inout[26 * iostride]) = tmp460 - tmp459;
		    tmp453 = tmp429 + tmp426;
		    tmp458 = tmp454 + tmp457;
		    c_im(inout[2 * iostride]) = tmp453 + tmp458;
		    c_im(inout[18 * iostride]) = tmp458 - tmp453;
	       }
	  }
	  {
	       fftw_real tmp307;
	       fftw_real tmp335;
	       fftw_real tmp338;
	       fftw_real tmp492;
	       fftw_real tmp487;
	       fftw_real tmp493;
	       fftw_real tmp314;
	       fftw_real tmp484;
	       fftw_real tmp322;
	       fftw_real tmp333;
	       fftw_real tmp342;
	       fftw_real tmp349;
	       fftw_real tmp329;
	       fftw_real tmp332;
	       fftw_real tmp345;
	       fftw_real tmp348;
	       ASSERT_ALIGNED_DOUBLE;
	       {
		    fftw_real tmp306;
		    fftw_real tmp336;
		    fftw_real tmp337;
		    fftw_real tmp485;
		    fftw_real tmp310;
		    fftw_real tmp313;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp306 = K707106781 * (tmp304 - tmp305);
		    tmp307 = tmp303 - tmp306;
		    tmp335 = tmp303 + tmp306;
		    tmp336 = (K382683432 * tmp312) - (K923879532 * tmp311);
		    tmp337 = (K923879532 * tmp308) + (K382683432 * tmp309);
		    tmp338 = tmp336 + tmp337;
		    tmp492 = tmp336 - tmp337;
		    tmp485 = K707106781 * (tmp196 - tmp201);
		    tmp487 = tmp485 + tmp486;
		    tmp493 = tmp486 - tmp485;
		    tmp310 = (K382683432 * tmp308) - (K923879532 * tmp309);
		    tmp313 = (K382683432 * tmp311) + (K923879532 * tmp312);
		    tmp314 = tmp310 - tmp313;
		    tmp484 = tmp313 + tmp310;
	       }
	       {
		    fftw_real tmp318;
		    fftw_real tmp321;
		    fftw_real tmp340;
		    fftw_real tmp341;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp318 = tmp316 - tmp317;
		    tmp321 = tmp319 - tmp320;
		    tmp322 = (K195090322 * tmp318) - (K980785280 * tmp321);
		    tmp333 = (K980785280 * tmp318) + (K195090322 * tmp321);
		    tmp340 = tmp316 + tmp317;
		    tmp341 = tmp319 + tmp320;
		    tmp342 = (K831469612 * tmp340) - (K555570233 * tmp341);
		    tmp349 = (K555570233 * tmp340) + (K831469612 * tmp341);
	       }
	       {
		    fftw_real tmp325;
		    fftw_real tmp328;
		    fftw_real tmp343;
		    fftw_real tmp344;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp325 = tmp323 - tmp324;
		    tmp328 = tmp326 - tmp327;
		    tmp329 = (K195090322 * tmp325) + (K980785280 * tmp328);
		    tmp332 = (K195090322 * tmp328) - (K980785280 * tmp325);
		    tmp343 = tmp323 + tmp324;
		    tmp344 = tmp326 + tmp327;
		    tmp345 = (K831469612 * tmp343) + (K555570233 * tmp344);
		    tmp348 = (K831469612 * tmp344) - (K555570233 * tmp343);
	       }
	       {
		    fftw_real tmp315;
		    fftw_real tmp330;
		    fftw_real tmp331;
		    fftw_real tmp334;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp315 = tmp307 - tmp314;
		    tmp330 = tmp322 - tmp329;
		    c_re(inout[31 * iostride]) = tmp315 - tmp330;
		    c_re(inout[15 * iostride]) = tmp315 + tmp330;
		    tmp331 = tmp307 + tmp314;
		    tmp334 = tmp332 + tmp333;
		    c_re(inout[23 * iostride]) = tmp331 - tmp334;
		    c_re(inout[7 * iostride]) = tmp331 + tmp334;
	       }
	       {
		    fftw_real tmp495;
		    fftw_real tmp496;
		    fftw_real tmp491;
		    fftw_real tmp494;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp495 = tmp332 - tmp333;
		    tmp496 = tmp493 - tmp492;
		    c_im(inout[15 * iostride]) = tmp495 + tmp496;
		    c_im(inout[31 * iostride]) = tmp496 - tmp495;
		    tmp491 = tmp329 + tmp322;
		    tmp494 = tmp492 + tmp493;
		    c_im(inout[7 * iostride]) = tmp491 + tmp494;
		    c_im(inout[23 * iostride]) = tmp494 - tmp491;
	       }
	       {
		    fftw_real tmp339;
		    fftw_real tmp346;
		    fftw_real tmp347;
		    fftw_real tmp350;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp339 = tmp335 - tmp338;
		    tmp346 = tmp342 - tmp345;
		    c_re(inout[27 * iostride]) = tmp339 - tmp346;
		    c_re(inout[11 * iostride]) = tmp339 + tmp346;
		    tmp347 = tmp335 + tmp338;
		    tmp350 = tmp348 + tmp349;
		    c_re(inout[19 * iostride]) = tmp347 - tmp350;
		    c_re(inout[3 * iostride]) = tmp347 + tmp350;
	       }
	       {
		    fftw_real tmp489;
		    fftw_real tmp490;
		    fftw_real tmp483;
		    fftw_real tmp488;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp489 = tmp348 - tmp349;
		    tmp490 = tmp487 - tmp484;
		    c_im(inout[11 * iostride]) = tmp489 + tmp490;
		    c_im(inout[27 * iostride]) = tmp490 - tmp489;
		    tmp483 = tmp345 + tmp342;
		    tmp488 = tmp484 + tmp487;
		    c_im(inout[3 * iostride]) = tmp483 + tmp488;
		    c_im(inout[19 * iostride]) = tmp488 - tmp483;
	       }
	  }
	  {
	       fftw_real tmp203;
	       fftw_real tmp287;
	       fftw_real tmp290;
	       fftw_real tmp478;
	       fftw_real tmp473;
	       fftw_real tmp479;
	       fftw_real tmp226;
	       fftw_real tmp468;
	       fftw_real tmp254;
	       fftw_real tmp285;
	       fftw_real tmp294;
	       fftw_real tmp301;
	       fftw_real tmp281;
	       fftw_real tmp284;
	       fftw_real tmp297;
	       fftw_real tmp300;
	       ASSERT_ALIGNED_DOUBLE;
	       {
		    fftw_real tmp202;
		    fftw_real tmp288;
		    fftw_real tmp289;
		    fftw_real tmp469;
		    fftw_real tmp214;
		    fftw_real tmp225;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp202 = K707106781 * (tmp196 + tmp201);
		    tmp203 = tmp191 - tmp202;
		    tmp287 = tmp191 + tmp202;
		    tmp288 = (K923879532 * tmp224) - (K382683432 * tmp219);
		    tmp289 = (K382683432 * tmp208) + (K923879532 * tmp213);
		    tmp290 = tmp288 + tmp289;
		    tmp478 = tmp288 - tmp289;
		    tmp469 = K707106781 * (tmp305 + tmp304);
		    tmp473 = tmp469 + tmp472;
		    tmp479 = tmp472 - tmp469;
		    tmp214 = (K923879532 * tmp208) - (K382683432 * tmp213);
		    tmp225 = (K923879532 * tmp219) + (K382683432 * tmp224);
		    tmp226 = tmp214 - tmp225;
		    tmp468 = tmp225 + tmp214;
	       }
	       {
		    fftw_real tmp244;
		    fftw_real tmp253;
		    fftw_real tmp292;
		    fftw_real tmp293;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp244 = tmp232 - tmp243;
		    tmp253 = tmp249 - tmp252;
		    tmp254 = (K555570233 * tmp244) - (K831469612 * tmp253);
		    tmp285 = (K831469612 * tmp244) + (K555570233 * tmp253);
		    tmp292 = tmp232 + tmp243;
		    tmp293 = tmp249 + tmp252;
		    tmp294 = (K980785280 * tmp292) - (K195090322 * tmp293);
		    tmp301 = (K195090322 * tmp292) + (K980785280 * tmp293);
	       }
	       {
		    fftw_real tmp271;
		    fftw_real tmp280;
		    fftw_real tmp295;
		    fftw_real tmp296;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp271 = tmp259 - tmp270;
		    tmp280 = tmp276 - tmp279;
		    tmp281 = (K555570233 * tmp271) + (K831469612 * tmp280);
		    tmp284 = (K555570233 * tmp280) - (K831469612 * tmp271);
		    tmp295 = tmp259 + tmp270;
		    tmp296 = tmp276 + tmp279;
		    tmp297 = (K980785280 * tmp295) + (K195090322 * tmp296);
		    tmp300 = (K980785280 * tmp296) - (K195090322 * tmp295);
	       }
	       {
		    fftw_real tmp227;
		    fftw_real tmp282;
		    fftw_real tmp283;
		    fftw_real tmp286;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp227 = tmp203 + tmp226;
		    tmp282 = tmp254 + tmp281;
		    c_re(inout[21 * iostride]) = tmp227 - tmp282;
		    c_re(inout[5 * iostride]) = tmp227 + tmp282;
		    tmp283 = tmp203 - tmp226;
		    tmp286 = tmp284 - tmp285;
		    c_re(inout[29 * iostride]) = tmp283 - tmp286;
		    c_re(inout[13 * iostride]) = tmp283 + tmp286;
	       }
	       {
		    fftw_real tmp477;
		    fftw_real tmp480;
		    fftw_real tmp481;
		    fftw_real tmp482;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp477 = tmp285 + tmp284;
		    tmp480 = tmp478 + tmp479;
		    c_im(inout[5 * iostride]) = tmp477 + tmp480;
		    c_im(inout[21 * iostride]) = tmp480 - tmp477;
		    tmp481 = tmp254 - tmp281;
		    tmp482 = tmp479 - tmp478;
		    c_im(inout[13 * iostride]) = tmp481 + tmp482;
		    c_im(inout[29 * iostride]) = tmp482 - tmp481;
	       }
	       {
		    fftw_real tmp291;
		    fftw_real tmp298;
		    fftw_real tmp299;
		    fftw_real tmp302;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp291 = tmp287 + tmp290;
		    tmp298 = tmp294 + tmp297;
		    c_re(inout[17 * iostride]) = tmp291 - tmp298;
		    c_re(inout[iostride]) = tmp291 + tmp298;
		    tmp299 = tmp287 - tmp290;
		    tmp302 = tmp300 - tmp301;
		    c_re(inout[25 * iostride]) = tmp299 - tmp302;
		    c_re(inout[9 * iostride]) = tmp299 + tmp302;
	       }
	       {
		    fftw_real tmp467;
		    fftw_real tmp474;
		    fftw_real tmp475;
		    fftw_real tmp476;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp467 = tmp301 + tmp300;
		    tmp474 = tmp468 + tmp473;
		    c_im(inout[iostride]) = tmp467 + tmp474;
		    c_im(inout[17 * iostride]) = tmp474 - tmp467;
		    tmp475 = tmp294 - tmp297;
		    tmp476 = tmp473 - tmp468;
		    c_im(inout[9 * iostride]) = tmp475 + tmp476;
		    c_im(inout[25 * iostride]) = tmp476 - tmp475;
	       }
	  }
     }
}

static const int twiddle_order[] =
    { 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19,
20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31 };
fftw_codelet_desc fftwi_twiddle_32_desc = {
     "fftwi_twiddle_32",
     (void (*)()) fftwi_twiddle_32,
     32,
     FFTW_BACKWARD,
     FFTW_TWIDDLE,
     715,
     31,
     twiddle_order,
};