ftw_32.c

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

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

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 fftw_twiddle_32_desc = {
     "fftw_twiddle_32",
     (void (*)()) fftw_twiddle_32,
     32,
     FFTW_FORWARD,
     FFTW_TWIDDLE,
     704,
     31,
     twiddle_order,
};