frc_64.c

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

		    tmp205 = (K995184726 * tmp196) + (K098017140 * tmp195);
	       }
	  }
	  {
	       fftw_real tmp95;
	       fftw_real tmp170;
	       fftw_real tmp181;
	       fftw_real tmp182;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp95 = tmp75 + tmp94;
	       tmp170 = tmp132 + tmp169;
	       real_output[25 * real_ostride] = tmp95 - tmp170;
	       real_output[7 * real_ostride] = tmp95 + tmp170;
	       tmp181 = tmp173 - tmp180;
	       tmp182 = tmp169 - tmp132;
	       imag_output[9 * imag_ostride] = tmp181 + tmp182;
	       imag_output[23 * imag_ostride] = tmp182 - tmp181;
	  }
	  {
	       fftw_real tmp183;
	       fftw_real tmp186;
	       fftw_real tmp187;
	       fftw_real tmp188;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp183 = tmp180 + tmp173;
	       tmp186 = tmp184 + tmp185;
	       imag_output[7 * imag_ostride] = tmp183 + tmp186;
	       imag_output[25 * imag_ostride] = tmp186 - tmp183;
	       tmp187 = tmp75 - tmp94;
	       tmp188 = tmp185 - tmp184;
	       real_output[23 * real_ostride] = tmp187 - tmp188;
	       real_output[9 * real_ostride] = tmp187 + tmp188;
	  }
	  {
	       fftw_real tmp191;
	       fftw_real tmp198;
	       fftw_real tmp201;
	       fftw_real tmp202;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp191 = tmp189 + tmp190;
	       tmp198 = tmp194 + tmp197;
	       real_output[31 * real_ostride] = tmp191 - tmp198;
	       real_output[real_ostride] = tmp191 + tmp198;
	       tmp201 = tmp199 - tmp200;
	       tmp202 = tmp197 - tmp194;
	       imag_output[15 * imag_ostride] = tmp201 + tmp202;
	       imag_output[17 * imag_ostride] = tmp202 - tmp201;
	  }
	  {
	       fftw_real tmp203;
	       fftw_real tmp206;
	       fftw_real tmp207;
	       fftw_real tmp208;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp203 = tmp200 + tmp199;
	       tmp206 = tmp204 + tmp205;
	       imag_output[imag_ostride] = tmp203 + tmp206;
	       imag_output[31 * imag_ostride] = tmp206 - tmp203;
	       tmp207 = tmp189 - tmp190;
	       tmp208 = tmp205 - tmp204;
	       real_output[17 * real_ostride] = tmp207 - tmp208;
	       real_output[15 * real_ostride] = tmp207 + tmp208;
	  }
     }
     {
	  fftw_real tmp377;
	  fftw_real tmp393;
	  fftw_real tmp387;
	  fftw_real tmp389;
	  fftw_real tmp384;
	  fftw_real tmp388;
	  fftw_real tmp392;
	  fftw_real tmp394;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp373;
	       fftw_real tmp376;
	       fftw_real tmp385;
	       fftw_real tmp386;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp373 = tmp7 - tmp14;
	       tmp376 = K707106781 * (tmp374 + tmp375);
	       tmp377 = tmp373 + tmp376;
	       tmp393 = tmp373 - tmp376;
	       tmp385 = K707106781 * (tmp375 - tmp374);
	       tmp386 = tmp29 - tmp22;
	       tmp387 = tmp385 - tmp386;
	       tmp389 = tmp386 + tmp385;
	  }
	  {
	       fftw_real tmp380;
	       fftw_real tmp383;
	       fftw_real tmp390;
	       fftw_real tmp391;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp380 = (K923879532 * tmp378) + (K382683432 * tmp379);
	       tmp383 = (K923879532 * tmp381) - (K382683432 * tmp382);
	       tmp384 = tmp380 + tmp383;
	       tmp388 = tmp383 - tmp380;
	       tmp390 = (K923879532 * tmp379) - (K382683432 * tmp378);
	       tmp391 = (K382683432 * tmp381) + (K923879532 * tmp382);
	       tmp392 = tmp390 + tmp391;
	       tmp394 = tmp391 - tmp390;
	  }
	  real_output[28 * real_ostride] = tmp377 - tmp384;
	  real_output[4 * real_ostride] = tmp377 + tmp384;
	  imag_output[12 * imag_ostride] = tmp387 + tmp388;
	  imag_output[20 * imag_ostride] = tmp388 - tmp387;
	  imag_output[4 * imag_ostride] = tmp389 + tmp392;
	  imag_output[28 * imag_ostride] = tmp392 - tmp389;
	  real_output[20 * real_ostride] = tmp393 - tmp394;
	  real_output[12 * real_ostride] = tmp393 + tmp394;
     }
     {
	  fftw_real tmp281;
	  fftw_real tmp329;
	  fftw_real tmp323;
	  fftw_real tmp325;
	  fftw_real tmp298;
	  fftw_real tmp326;
	  fftw_real tmp315;
	  fftw_real tmp327;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp273;
	       fftw_real tmp280;
	       fftw_real tmp319;
	       fftw_real tmp322;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp273 = tmp269 + tmp272;
	       tmp280 = tmp276 + tmp279;
	       tmp281 = tmp273 + tmp280;
	       tmp329 = tmp273 - tmp280;
	       tmp319 = tmp317 + tmp318;
	       tmp322 = tmp320 + tmp321;
	       tmp323 = tmp319 - tmp322;
	       tmp325 = tmp322 + tmp319;
	  }
	  {
	       fftw_real tmp292;
	       fftw_real tmp297;
	       fftw_real tmp309;
	       fftw_real tmp314;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp292 = tmp284 + tmp291;
	       tmp297 = tmp295 + tmp296;
	       tmp298 = (K980785280 * tmp292) + (K195090322 * tmp297);
	       tmp326 = (K980785280 * tmp297) - (K195090322 * tmp292);
	       tmp309 = tmp301 + tmp308;
	       tmp314 = tmp312 + tmp313;
	       tmp315 = (K980785280 * tmp309) - (K195090322 * tmp314);
	       tmp327 = (K195090322 * tmp309) + (K980785280 * tmp314);
	  }
	  {
	       fftw_real tmp316;
	       fftw_real tmp324;
	       fftw_real tmp328;
	       fftw_real tmp330;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp316 = tmp298 + tmp315;
	       real_output[30 * real_ostride] = tmp281 - tmp316;
	       real_output[2 * real_ostride] = tmp281 + tmp316;
	       tmp324 = tmp315 - tmp298;
	       imag_output[14 * imag_ostride] = tmp323 + tmp324;
	       imag_output[18 * imag_ostride] = tmp324 - tmp323;
	       tmp328 = tmp326 + tmp327;
	       imag_output[2 * imag_ostride] = tmp325 + tmp328;
	       imag_output[30 * imag_ostride] = tmp328 - tmp325;
	       tmp330 = tmp327 - tmp326;
	       real_output[18 * real_ostride] = tmp329 - tmp330;
	       real_output[14 * real_ostride] = tmp329 + tmp330;
	  }
     }
     {
	  fftw_real tmp333;
	  fftw_real tmp349;
	  fftw_real tmp343;
	  fftw_real tmp345;
	  fftw_real tmp336;
	  fftw_real tmp346;
	  fftw_real tmp339;
	  fftw_real tmp347;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp331;
	       fftw_real tmp332;
	       fftw_real tmp341;
	       fftw_real tmp342;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp331 = tmp269 - tmp272;
	       tmp332 = tmp318 - tmp317;
	       tmp333 = tmp331 + tmp332;
	       tmp349 = tmp331 - tmp332;
	       tmp341 = tmp279 - tmp276;
	       tmp342 = tmp321 - tmp320;
	       tmp343 = tmp341 - tmp342;
	       tmp345 = tmp342 + tmp341;
	  }
	  {
	       fftw_real tmp334;
	       fftw_real tmp335;
	       fftw_real tmp337;
	       fftw_real tmp338;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp334 = tmp284 - tmp291;
	       tmp335 = tmp296 - tmp295;
	       tmp336 = (K831469612 * tmp334) + (K555570233 * tmp335);
	       tmp346 = (K831469612 * tmp335) - (K555570233 * tmp334);
	       tmp337 = tmp301 - tmp308;
	       tmp338 = tmp313 - tmp312;
	       tmp339 = (K831469612 * tmp337) - (K555570233 * tmp338);
	       tmp347 = (K555570233 * tmp337) + (K831469612 * tmp338);
	  }
	  {
	       fftw_real tmp340;
	       fftw_real tmp344;
	       fftw_real tmp348;
	       fftw_real tmp350;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp340 = tmp336 + tmp339;
	       real_output[26 * real_ostride] = tmp333 - tmp340;
	       real_output[6 * real_ostride] = tmp333 + tmp340;
	       tmp344 = tmp339 - tmp336;
	       imag_output[10 * imag_ostride] = tmp343 + tmp344;
	       imag_output[22 * imag_ostride] = tmp344 - tmp343;
	       tmp348 = tmp346 + tmp347;
	       imag_output[6 * imag_ostride] = tmp345 + tmp348;
	       imag_output[26 * imag_ostride] = tmp348 - tmp345;
	       tmp350 = tmp347 - tmp346;
	       real_output[22 * real_ostride] = tmp349 - tmp350;
	       real_output[10 * real_ostride] = tmp349 + tmp350;
	  }
     }
     {
	  fftw_real tmp211;
	  fftw_real tmp249;
	  fftw_real tmp218;
	  fftw_real tmp259;
	  fftw_real tmp237;
	  fftw_real tmp250;
	  fftw_real tmp226;
	  fftw_real tmp244;
	  fftw_real tmp240;
	  fftw_real tmp260;
	  fftw_real tmp254;
	  fftw_real tmp264;
	  fftw_real tmp233;
	  fftw_real tmp245;
	  fftw_real tmp257;
	  fftw_real tmp265;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp214;
	       fftw_real tmp217;
	       fftw_real tmp252;
	       fftw_real tmp253;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp211 = tmp209 + tmp210;
	       tmp249 = tmp209 - tmp210;
	       tmp214 = (K831469612 * tmp212) + (K555570233 * tmp213);
	       tmp217 = (K831469612 * tmp215) - (K555570233 * tmp216);
	       tmp218 = tmp214 + tmp217;
	       tmp259 = tmp217 - tmp214;
	       {
		    fftw_real tmp235;
		    fftw_real tmp236;
		    fftw_real tmp222;
		    fftw_real tmp225;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp235 = (K831469612 * tmp213) - (K555570233 * tmp212);
		    tmp236 = (K555570233 * tmp215) + (K831469612 * tmp216);
		    tmp237 = tmp235 + tmp236;
		    tmp250 = tmp236 - tmp235;
		    tmp222 = tmp220 + tmp221;
		    tmp225 = tmp223 + tmp224;
		    tmp226 = (K956940335 * tmp222) + (K290284677 * tmp225);
		    tmp244 = (K956940335 * tmp225) - (K290284677 * tmp222);
	       }
	       tmp240 = tmp238 + tmp239;
	       tmp260 = tmp239 - tmp238;
	       tmp252 = tmp220 - tmp221;
	       tmp253 = tmp224 - tmp223;
	       tmp254 = (K881921264 * tmp252) + (K471396736 * tmp253);
	       tmp264 = (K881921264 * tmp253) - (K471396736 * tmp252);
	       {
		    fftw_real tmp229;
		    fftw_real tmp232;
		    fftw_real tmp255;
		    fftw_real tmp256;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp229 = tmp227 + tmp228;
		    tmp232 = tmp230 + tmp231;
		    tmp233 = (K956940335 * tmp229) - (K290284677 * tmp232);
		    tmp245 = (K290284677 * tmp229) + (K956940335 * tmp232);
		    tmp255 = tmp227 - tmp228;
		    tmp256 = tmp231 - tmp230;
		    tmp257 = (K881921264 * tmp255) - (K471396736 * tmp256);
		    tmp265 = (K471396736 * tmp255) + (K881921264 * tmp256);
	       }
	  }
	  {
	       fftw_real tmp219;
	       fftw_real tmp234;
	       fftw_real tmp241;
	       fftw_real tmp242;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp219 = tmp211 + tmp218;
	       tmp234 = tmp226 + tmp233;
	       real_output[29 * real_ostride] = tmp219 - tmp234;
	       real_output[3 * real_ostride] = tmp219 + tmp234;
	       tmp241 = tmp237 - tmp240;
	       tmp242 = tmp233 - tmp226;
	       imag_output[13 * imag_ostride] = tmp241 + tmp242;
	       imag_output[19 * imag_ostride] = tmp242 - tmp241;
	  }
	  {
	       fftw_real tmp243;
	       fftw_real tmp246;
	       fftw_real tmp247;
	       fftw_real tmp248;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp243 = tmp240 + tmp237;
	       tmp246 = tmp244 + tmp245;
	       imag_output[3 * imag_ostride] = tmp243 + tmp246;
	       imag_output[29 * imag_ostride] = tmp246 - tmp243;
	       tmp247 = tmp211 - tmp218;
	       tmp248 = tmp245 - tmp244;
	       real_output[19 * real_ostride] = tmp247 - tmp248;
	       real_output[13 * real_ostride] = tmp247 + tmp248;
	  }
	  {
	       fftw_real tmp251;
	       fftw_real tmp258;
	       fftw_real tmp261;
	       fftw_real tmp262;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp251 = tmp249 + tmp250;
	       tmp258 = tmp254 + tmp257;
	       real_output[27 * real_ostride] = tmp251 - tmp258;
	       real_output[5 * real_ostride] = tmp251 + tmp258;
	       tmp261 = tmp259 - tmp260;
	       tmp262 = tmp257 - tmp254;
	       imag_output[11 * imag_ostride] = tmp261 + tmp262;
	       imag_output[21 * imag_ostride] = tmp262 - tmp261;
	  }
	  {
	       fftw_real tmp263;
	       fftw_real tmp266;
	       fftw_real tmp267;
	       fftw_real tmp268;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp263 = tmp260 + tmp259;
	       tmp266 = tmp264 + tmp265;
	       imag_output[5 * imag_ostride] = tmp263 + tmp266;
	       imag_output[27 * imag_ostride] = tmp266 - tmp263;
	       tmp267 = tmp249 - tmp250;
	       tmp268 = tmp265 - tmp264;
	       real_output[21 * real_ostride] = tmp267 - tmp268;
	       real_output[11 * real_ostride] = tmp267 + tmp268;
	  }
     }
}

fftw_codelet_desc fftw_real2hc_64_desc = {
     "fftw_real2hc_64",
     (void (*)()) fftw_real2hc_64,
     64,
     FFTW_FORWARD,
     FFTW_REAL2HC,
     1410,
     0,
     (const int *) 0,
};