fhb_32.c

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

		    tmp797 = tmp795 - tmp796;
		    tmp801 = tmp796 + tmp795;
		    tmp798 = tmp789 - tmp782;
		    tmp799 = tmp793 + tmp792;
		    tmp800 = K1_414213562 * (tmp798 - tmp799);
		    tmp802 = K1_414213562 * (tmp798 + tmp799);
	       }
	       X[21 * iostride] = tmp797 - tmp800;
	       X[5 * iostride] = tmp800 + tmp797;
	       X[13 * iostride] = tmp801 - tmp802;
	       X[29 * iostride] = tmp802 + tmp801;
	  }
	  {
	       fftw_real tmp816;
	       fftw_real tmp820;
	       fftw_real tmp809;
	       fftw_real tmp817;
	       fftw_real tmp805;
	       fftw_real tmp808;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp816 = K2_000000000 * (tmp812 - tmp815);
	       tmp820 = K2_000000000 * (tmp818 - tmp819);
	       tmp805 = tmp803 - tmp804;
	       tmp808 = (K765366864 * tmp806) - (K1_847759065 * tmp807);
	       tmp809 = tmp805 - tmp808;
	       tmp817 = tmp808 + tmp805;
	       X[11 * iostride] = tmp809 - tmp816;
	       X[27 * iostride] = tmp816 + tmp809;
	       X[19 * iostride] = tmp817 - tmp820;
	       X[3 * iostride] = tmp820 + tmp817;
	  }
	  {
	       fftw_real tmp823;
	       fftw_real tmp827;
	       fftw_real tmp826;
	       fftw_real tmp828;
	       ASSERT_ALIGNED_DOUBLE;
	       {
		    fftw_real tmp821;
		    fftw_real tmp822;
		    fftw_real tmp824;
		    fftw_real tmp825;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp821 = tmp804 + tmp803;
		    tmp822 =
			(K1_847759065 * tmp806) + (K765366864 * tmp807);
		    tmp823 = tmp821 - tmp822;
		    tmp827 = tmp822 + tmp821;
		    tmp824 = tmp819 + tmp818;
		    tmp825 = tmp815 + tmp812;
		    tmp826 = K1_414213562 * (tmp824 - tmp825);
		    tmp828 = K1_414213562 * (tmp824 + tmp825);
	       }
	       X[23 * iostride] = tmp823 - tmp826;
	       X[7 * iostride] = tmp826 + tmp823;
	       X[15 * iostride] = tmp827 - tmp828;
	       X[31 * iostride] = tmp828 + tmp827;
	  }
     }
     X = X + dist;
     Y = Y - dist;
     for (i = 2; i < m; i = i + 2, X = X + dist, Y = Y - dist, W = W + 31) {
	  fftw_real tmp446;
	  fftw_real tmp580;
	  fftw_real tmp189;
	  fftw_real tmp238;
	  fftw_real tmp382;
	  fftw_real tmp420;
	  fftw_real tmp307;
	  fftw_real tmp332;
	  fftw_real tmp514;
	  fftw_real tmp608;
	  fftw_real tmp354;
	  fftw_real tmp408;
	  fftw_real tmp453;
	  fftw_real tmp609;
	  fftw_real tmp517;
	  fftw_real tmp581;
	  fftw_real tmp235;
	  fftw_real tmp272;
	  fftw_real tmp365;
	  fftw_real tmp411;
	  fftw_real tmp287;
	  fftw_real tmp336;
	  fftw_real tmp368;
	  fftw_real tmp412;
	  fftw_real tmp500;
	  fftw_real tmp550;
	  fftw_real tmp593;
	  fftw_real tmp643;
	  fftw_real tmp507;
	  fftw_real tmp551;
	  fftw_real tmp596;
	  fftw_real tmp642;
	  fftw_real tmp204;
	  fftw_real tmp292;
	  fftw_real tmp357;
	  fftw_real tmp384;
	  fftw_real tmp253;
	  fftw_real tmp333;
	  fftw_real tmp360;
	  fftw_real tmp383;
	  fftw_real tmp461;
	  fftw_real tmp520;
	  fftw_real tmp588;
	  fftw_real tmp612;
	  fftw_real tmp468;
	  fftw_real tmp519;
	  fftw_real tmp585;
	  fftw_real tmp611;
	  fftw_real tmp220;
	  fftw_real tmp255;
	  fftw_real tmp372;
	  fftw_real tmp414;
	  fftw_real tmp270;
	  fftw_real tmp335;
	  fftw_real tmp375;
	  fftw_real tmp415;
	  fftw_real tmp481;
	  fftw_real tmp547;
	  fftw_real tmp600;
	  fftw_real tmp640;
	  fftw_real tmp488;
	  fftw_real tmp548;
	  fftw_real tmp603;
	  fftw_real tmp639;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp177;
	       fftw_real tmp444;
	       fftw_real tmp295;
	       fftw_real tmp513;
	       fftw_real tmp180;
	       fftw_real tmp512;
	       fftw_real tmp298;
	       fftw_real tmp445;
	       fftw_real tmp184;
	       fftw_real tmp447;
	       fftw_real tmp302;
	       fftw_real tmp448;
	       fftw_real tmp187;
	       fftw_real tmp450;
	       fftw_real tmp305;
	       fftw_real tmp451;
	       ASSERT_ALIGNED_DOUBLE;
	       {
		    fftw_real tmp175;
		    fftw_real tmp176;
		    fftw_real tmp293;
		    fftw_real tmp294;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp175 = X[0];
		    tmp176 = Y[-16 * iostride];
		    tmp177 = tmp175 + tmp176;
		    tmp444 = tmp175 - tmp176;
		    tmp293 = Y[0];
		    tmp294 = X[16 * iostride];
		    tmp295 = tmp293 - tmp294;
		    tmp513 = tmp293 + tmp294;
	       }
	       {
		    fftw_real tmp178;
		    fftw_real tmp179;
		    fftw_real tmp296;
		    fftw_real tmp297;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp178 = X[8 * iostride];
		    tmp179 = Y[-24 * iostride];
		    tmp180 = tmp178 + tmp179;
		    tmp512 = tmp178 - tmp179;
		    tmp296 = Y[-8 * iostride];
		    tmp297 = X[24 * iostride];
		    tmp298 = tmp296 - tmp297;
		    tmp445 = tmp296 + tmp297;
	       }
	       {
		    fftw_real tmp182;
		    fftw_real tmp183;
		    fftw_real tmp300;
		    fftw_real tmp301;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp182 = X[4 * iostride];
		    tmp183 = Y[-20 * iostride];
		    tmp184 = tmp182 + tmp183;
		    tmp447 = tmp182 - tmp183;
		    tmp300 = Y[-4 * iostride];
		    tmp301 = X[20 * iostride];
		    tmp302 = tmp300 - tmp301;
		    tmp448 = tmp300 + tmp301;
	       }
	       {
		    fftw_real tmp185;
		    fftw_real tmp186;
		    fftw_real tmp303;
		    fftw_real tmp304;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp185 = Y[-28 * iostride];
		    tmp186 = X[12 * iostride];
		    tmp187 = tmp185 + tmp186;
		    tmp450 = tmp185 - tmp186;
		    tmp303 = Y[-12 * iostride];
		    tmp304 = X[28 * iostride];
		    tmp305 = tmp303 - tmp304;
		    tmp451 = tmp303 + tmp304;
	       }
	       {
		    fftw_real tmp181;
		    fftw_real tmp188;
		    fftw_real tmp352;
		    fftw_real tmp353;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp446 = tmp444 - tmp445;
		    tmp580 = tmp444 + tmp445;
		    tmp181 = tmp177 + tmp180;
		    tmp188 = tmp184 + tmp187;
		    tmp189 = tmp181 + tmp188;
		    tmp238 = tmp181 - tmp188;
		    {
			 fftw_real tmp380;
			 fftw_real tmp381;
			 fftw_real tmp299;
			 fftw_real tmp306;
			 ASSERT_ALIGNED_DOUBLE;
			 tmp380 = tmp295 - tmp298;
			 tmp381 = tmp184 - tmp187;
			 tmp382 = tmp380 - tmp381;
			 tmp420 = tmp381 + tmp380;
			 tmp299 = tmp295 + tmp298;
			 tmp306 = tmp302 + tmp305;
			 tmp307 = tmp299 - tmp306;
			 tmp332 = tmp299 + tmp306;
		    }
		    tmp514 = tmp512 + tmp513;
		    tmp608 = tmp513 - tmp512;
		    tmp352 = tmp177 - tmp180;
		    tmp353 = tmp305 - tmp302;
		    tmp354 = tmp352 - tmp353;
		    tmp408 = tmp352 + tmp353;
		    {
			 fftw_real tmp449;
			 fftw_real tmp452;
			 fftw_real tmp515;
			 fftw_real tmp516;
			 ASSERT_ALIGNED_DOUBLE;
			 tmp449 = tmp447 - tmp448;
			 tmp452 = tmp450 - tmp451;
			 tmp453 = K707106781 * (tmp449 + tmp452);
			 tmp609 = K707106781 * (tmp449 - tmp452);
			 tmp515 = tmp447 + tmp448;
			 tmp516 = tmp450 + tmp451;
			 tmp517 = K707106781 * (tmp515 - tmp516);
			 tmp581 = K707106781 * (tmp515 + tmp516);
		    }
	       }
	  }
	  {
	       fftw_real tmp223;
	       fftw_real tmp490;
	       fftw_real tmp275;
	       fftw_real tmp502;
	       fftw_real tmp226;
	       fftw_real tmp501;
	       fftw_real tmp278;
	       fftw_real tmp491;
	       fftw_real tmp233;
	       fftw_real tmp505;
	       fftw_real tmp498;
	       fftw_real tmp285;
	       fftw_real tmp230;
	       fftw_real tmp504;
	       fftw_real tmp495;
	       fftw_real tmp282;
	       ASSERT_ALIGNED_DOUBLE;
	       {
		    fftw_real tmp221;
		    fftw_real tmp222;
		    fftw_real tmp276;
		    fftw_real tmp277;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp221 = Y[-31 * iostride];
		    tmp222 = X[15 * iostride];
		    tmp223 = tmp221 + tmp222;
		    tmp490 = tmp221 - tmp222;
		    {
			 fftw_real tmp273;
			 fftw_real tmp274;
			 fftw_real tmp224;
			 fftw_real tmp225;
			 ASSERT_ALIGNED_DOUBLE;
			 tmp273 = Y[-15 * iostride];
			 tmp274 = X[31 * iostride];
			 tmp275 = tmp273 - tmp274;
			 tmp502 = tmp273 + tmp274;
			 tmp224 = X[7 * iostride];
			 tmp225 = Y[-23 * iostride];
			 tmp226 = tmp224 + tmp225;
			 tmp501 = tmp224 - tmp225;
		    }
		    tmp276 = Y[-7 * iostride];
		    tmp277 = X[23 * iostride];
		    tmp278 = tmp276 - tmp277;
		    tmp491 = tmp276 + tmp277;
		    {
			 fftw_real tmp231;
			 fftw_real tmp232;
			 fftw_real tmp496;
			 fftw_real tmp283;
			 fftw_real tmp284;
			 fftw_real tmp497;
			 ASSERT_ALIGNED_DOUBLE;
			 tmp231 = Y[-27 * iostride];
			 tmp232 = X[11 * iostride];
			 tmp496 = tmp231 - tmp232;
			 tmp283 = Y[-11 * iostride];
			 tmp284 = X[27 * iostride];
			 tmp497 = tmp283 + tmp284;
			 tmp233 = tmp231 + tmp232;
			 tmp505 = tmp496 + tmp497;
			 tmp498 = tmp496 - tmp497;
			 tmp285 = tmp283 - tmp284;
		    }
		    {
			 fftw_real tmp228;
			 fftw_real tmp229;
			 fftw_real tmp493;
			 fftw_real tmp280;
			 fftw_real tmp281;
			 fftw_real tmp494;
			 ASSERT_ALIGNED_DOUBLE;
			 tmp228 = X[3 * iostride];
			 tmp229 = Y[-19 * iostride];
			 tmp493 = tmp228 - tmp229;
			 tmp280 = Y[-3 * iostride];
			 tmp281 = X[19 * iostride];
			 tmp494 = tmp280 + tmp281;
			 tmp230 = tmp228 + tmp229;
			 tmp504 = tmp493 + tmp494;
			 tmp495 = tmp493 - tmp494;
			 tmp282 = tmp280 - tmp281;
		    }
	       }
	       {
		    fftw_real tmp227;
		    fftw_real tmp234;
		    fftw_real tmp363;
		    fftw_real tmp364;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp227 = tmp223 + tmp226;
		    tmp234 = tmp230 + tmp233;
		    tmp235 = tmp227 + tmp234;
		    tmp272 = tmp227 - tmp234;
		    tmp363 = tmp275 - tmp278;
		    tmp364 = tmp230 - tmp233;
		    tmp365 = tmp363 - tmp364;
		    tmp411 = tmp364 + tmp363;
	       }
	       {
		    fftw_real tmp279;
		    fftw_real tmp286;
		    fftw_real tmp366;
		    fftw_real tmp367;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp279 = tmp275 + tmp278;
		    tmp286 = tmp282 + tmp285;
		    tmp287 = tmp279 - tmp286;
		    tmp336 = tmp279 + tmp286;
		    tmp366 = tmp223 - tmp226;
		    tmp367 = tmp285 - tmp282;
		    tmp368 = tmp366 - tmp367;
		    tmp412 = tmp366 + tmp367;
	       }
	       {
		    fftw_real tmp492;
		    fftw_real tmp499;
		    fftw_real tmp591;
		    fftw_real tmp592;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp492 = tmp490 - tmp491;
		    tmp499 = K707106781 * (tmp495 + tmp498);
		    tmp500 = tmp492 - tmp499;
		    tmp550 = tmp492 + tmp499;
		    tmp591 = K707106781 * (tmp495 - tmp498);
		    tmp592 = tmp501 + tmp502;
		    tmp593 = tmp591 - tmp592;
		    tmp643 = tmp592 + tmp591;
	       }
	       {
		    fftw_real tmp503;
		    fftw_real tmp506;
		    fftw_real tmp594;
		    fftw_real tmp595;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp503 = tmp501 - tmp502;
		    tmp506 = K707106781 * (tmp504 - tmp505);
		    tmp507 = tmp503 - tmp506;
		    tmp551 = tmp503 + tmp506;
		    tmp594 = tmp490 + tmp491;
		    tmp595 = K707106781 * (tmp504 + tmp505);
		    tmp596 = tmp594 - tmp595;
		    tmp642 = tmp594 + tmp595;
	       }
	  }
	  {
	       fftw_real tmp192;
	       fftw_real tmp465;
	       fftw_real tmp248;
	       fftw_real tmp463;
	       fftw_real tmp195;
	       fftw_real tmp462;
	       fftw_real tmp251;
	       fftw_real tmp466;
	       fftw_real tmp199;
	       fftw_real tmp458;
	       fftw_real tmp241;
	       fftw_real tmp456;
	       fftw_real tmp202;
	       fftw_real tmp455;
	       fftw_real tmp244;
	       fftw_real tmp459;
	       ASSERT_ALIGNED_DOUBLE;
	       {
		    fftw_real tmp190;
		    fftw_real tmp191;
		    fftw_real tmp246;
		    fftw_real tmp247;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp190 = X[2 * iostride];
		    tmp191 = Y[-18 * iostride];
		    tmp192 = tmp190 + tmp191;
		    tmp465 = tmp190 - tmp191;
		    tmp246 = Y[-2 * iostride];
		    tmp247 = X[18 * iostride];
		    tmp248 = tmp246 - tmp247;
		    tmp463 = tmp246 + tmp247;
	       }
	       {
		    fftw_real tmp193;
		    fftw_real tmp194;
		    fftw_real tmp249;
		    fftw_real tmp250;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp193 = X[10 * iostride];
		    tmp194 = Y[-26 * iostride];
		    tmp195 = tmp193 + tmp194;
		    tmp462 = tmp193 - tmp194;
		    tmp249 = Y[-10 * iostride];
		    tmp250 = X[26 * iostride];
		    tmp251 = tmp249 - tmp250;
		    tmp466 = tmp249 + tmp250;
	       }
	       {
		    fftw_real tmp197;
		    fftw_real tmp198;
		    fftw_real tmp239;
		    fftw_real tmp240;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp197 = Y[-30 * iostride];
		    tmp198 = X[14 * iostride];
		    tmp199 = tmp197 + tmp198;
		    tmp458 = tmp197 - tmp198;
		    tmp239 = Y[-14 * iostride];
		    tmp240 = X[30 * iostride];
		    tmp241 = tmp239 - tmp240;
		    tmp456 = tmp239 + tmp240;
	       }
	       {
		    fftw_real tmp200;
		    fftw_real tmp201;
		    fftw_real tmp242;
		    fftw_real tmp243;