fcr_64.c

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

	       fftw_real tmp141;
	       fftw_real tmp142;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp119 = (K980785280 * tmp109) - (K195090322 * tmp118);
	       tmp138 = (K831469612 * tmp128) - (K555570233 * tmp137);
	       tmp139 = K2_000000000 * (tmp119 + tmp138);
	       tmp147 = tmp119 - tmp138;
	       tmp141 = (K195090322 * tmp109) + (K980785280 * tmp118);
	       tmp142 = (K555570233 * tmp128) + (K831469612 * tmp137);
	       tmp143 = K2_000000000 * (tmp141 - tmp142);
	       tmp148 = tmp141 + tmp142;
	  }
	  output[34 * ostride] = tmp100 - tmp139;
	  output[2 * ostride] = tmp100 + tmp139;
	  output[18 * ostride] = tmp140 - tmp143;
	  output[50 * ostride] = tmp140 + tmp143;
	  tmp151 = K1_414213562 * (tmp147 + tmp148);
	  output[26 * ostride] = tmp150 - tmp151;
	  output[58 * ostride] = tmp150 + tmp151;
	  tmp149 = K1_414213562 * (tmp147 - tmp148);
	  output[42 * ostride] = tmp146 - tmp149;
	  output[10 * ostride] = tmp146 + tmp149;
     }
     {
	  fftw_real tmp339;
	  fftw_real tmp355;
	  fftw_real tmp367;
	  fftw_real tmp371;
	  fftw_real tmp354;
	  fftw_real tmp368;
	  fftw_real tmp362;
	  fftw_real tmp369;
	  fftw_real tmp372;
	  fftw_real tmp370;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp331;
	       fftw_real tmp338;
	       fftw_real tmp363;
	       fftw_real tmp366;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp331 = K2_000000000 * (tmp327 - tmp330);
	       tmp338 = tmp334 + tmp337;
	       tmp339 = tmp331 + tmp338;
	       tmp355 = tmp338 - tmp331;
	       tmp363 = tmp337 - tmp334;
	       tmp366 = K2_000000000 * (tmp364 - tmp365);
	       tmp367 = tmp363 - tmp366;
	       tmp371 = tmp366 + tmp363;
	  }
	  {
	       fftw_real tmp346;
	       fftw_real tmp353;
	       fftw_real tmp358;
	       fftw_real tmp361;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp346 = tmp342 - tmp345;
	       tmp353 = tmp349 + tmp352;
	       tmp354 = K2_000000000 * (tmp346 + tmp353);
	       tmp368 = tmp346 - tmp353;
	       tmp358 = tmp356 - tmp357;
	       tmp361 = tmp359 - tmp360;
	       tmp362 = K2_000000000 * (tmp358 - tmp361);
	       tmp369 = tmp358 + tmp361;
	  }
	  output[35 * ostride] = tmp339 - tmp354;
	  output[3 * ostride] = tmp354 + tmp339;
	  output[19 * ostride] = tmp355 - tmp362;
	  output[51 * ostride] = tmp362 + tmp355;
	  tmp372 = K1_414213562 * (tmp368 + tmp369);
	  output[27 * ostride] = tmp371 - tmp372;
	  output[59 * ostride] = tmp372 + tmp371;
	  tmp370 = K1_414213562 * (tmp368 - tmp369);
	  output[43 * ostride] = tmp367 - tmp370;
	  output[11 * ostride] = tmp370 + tmp367;
     }
     {
	  fftw_real tmp378;
	  fftw_real tmp389;
	  fftw_real tmp375;
	  fftw_real tmp390;
	  fftw_real tmp382;
	  fftw_real tmp392;
	  fftw_real tmp385;
	  fftw_real tmp393;
	  fftw_real tmp373;
	  fftw_real tmp374;
	  ASSERT_ALIGNED_DOUBLE;
	  tmp378 = tmp376 - tmp377;
	  tmp389 = tmp377 + tmp376;
	  tmp373 = tmp330 + tmp327;
	  tmp374 = tmp365 + tmp364;
	  tmp375 = K1_414213562 * (tmp373 - tmp374);
	  tmp390 = K1_414213562 * (tmp373 + tmp374);
	  {
	       fftw_real tmp380;
	       fftw_real tmp381;
	       fftw_real tmp383;
	       fftw_real tmp384;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp380 = tmp345 + tmp342;
	       tmp381 = tmp360 + tmp359;
	       tmp382 = tmp380 - tmp381;
	       tmp392 = tmp380 + tmp381;
	       tmp383 = tmp357 + tmp356;
	       tmp384 = tmp349 - tmp352;
	       tmp385 = tmp383 + tmp384;
	       tmp393 = tmp383 - tmp384;
	  }
	  {
	       fftw_real tmp379;
	       fftw_real tmp386;
	       fftw_real tmp387;
	       fftw_real tmp388;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp379 = tmp375 + tmp378;
	       tmp386 = (K1_847759065 * tmp382) - (K765366864 * tmp385);
	       output[39 * ostride] = tmp379 - tmp386;
	       output[7 * ostride] = tmp386 + tmp379;
	       tmp387 = tmp378 - tmp375;
	       tmp388 = (K1_847759065 * tmp385) + (K765366864 * tmp382);
	       output[23 * ostride] = tmp387 - tmp388;
	       output[55 * ostride] = tmp388 + tmp387;
	  }
	  {
	       fftw_real tmp391;
	       fftw_real tmp394;
	       fftw_real tmp395;
	       fftw_real tmp396;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp391 = tmp389 - tmp390;
	       tmp394 = (K765366864 * tmp392) - (K1_847759065 * tmp393);
	       output[47 * ostride] = tmp391 - tmp394;
	       output[15 * ostride] = tmp394 + tmp391;
	       tmp395 = tmp390 + tmp389;
	       tmp396 = (K765366864 * tmp393) + (K1_847759065 * tmp392);
	       output[31 * ostride] = tmp395 - tmp396;
	       output[63 * ostride] = tmp396 + tmp395;
	  }
     }
     {
	  fftw_real tmp214;
	  fftw_real tmp218;
	  fftw_real tmp217;
	  fftw_real tmp219;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp212;
	       fftw_real tmp213;
	       fftw_real tmp215;
	       fftw_real tmp216;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp212 = tmp196 + tmp197;
	       tmp213 = K1_414213562 * (tmp199 + tmp200);
	       tmp214 = tmp212 - tmp213;
	       tmp218 = tmp212 + tmp213;
	       tmp215 = tmp203 + tmp204;
	       tmp216 = tmp207 - tmp206;
	       tmp217 = (K765366864 * tmp215) - (K1_847759065 * tmp216);
	       tmp219 = (K765366864 * tmp216) + (K1_847759065 * tmp215);
	  }
	  output[44 * ostride] = tmp214 - tmp217;
	  output[12 * ostride] = tmp214 + tmp217;
	  output[28 * ostride] = tmp218 - tmp219;
	  output[60 * ostride] = tmp218 + tmp219;
     }
     {
	  fftw_real tmp202;
	  fftw_real tmp210;
	  fftw_real tmp209;
	  fftw_real tmp211;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp198;
	       fftw_real tmp201;
	       fftw_real tmp205;
	       fftw_real tmp208;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp198 = tmp196 - tmp197;
	       tmp201 = K1_414213562 * (tmp199 - tmp200);
	       tmp202 = tmp198 + tmp201;
	       tmp210 = tmp198 - tmp201;
	       tmp205 = tmp203 - tmp204;
	       tmp208 = tmp206 + tmp207;
	       tmp209 = (K1_847759065 * tmp205) - (K765366864 * tmp208);
	       tmp211 = (K1_847759065 * tmp208) + (K765366864 * tmp205);
	  }
	  output[36 * ostride] = tmp202 - tmp209;
	  output[4 * ostride] = tmp202 + tmp209;
	  output[20 * ostride] = tmp210 - tmp211;
	  output[52 * ostride] = tmp210 + tmp211;
     }
     {
	  fftw_real tmp158;
	  fftw_real tmp166;
	  fftw_real tmp172;
	  fftw_real tmp176;
	  fftw_real tmp165;
	  fftw_real tmp174;
	  fftw_real tmp169;
	  fftw_real tmp173;
	  fftw_real tmp177;
	  fftw_real tmp175;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp154;
	       fftw_real tmp157;
	       fftw_real tmp170;
	       fftw_real tmp171;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp154 = tmp152 - tmp153;
	       tmp157 = (K765366864 * tmp155) - (K1_847759065 * tmp156);
	       tmp158 = tmp154 - tmp157;
	       tmp166 = tmp154 + tmp157;
	       tmp170 = tmp152 + tmp153;
	       tmp171 = (K1_847759065 * tmp155) + (K765366864 * tmp156);
	       tmp172 = tmp170 - tmp171;
	       tmp176 = tmp170 + tmp171;
	  }
	  {
	       fftw_real tmp161;
	       fftw_real tmp164;
	       fftw_real tmp167;
	       fftw_real tmp168;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp161 = (K555570233 * tmp159) + (K831469612 * tmp160);
	       tmp164 = (K980785280 * tmp162) - (K195090322 * tmp163);
	       tmp165 = K2_000000000 * (tmp161 - tmp164);
	       tmp174 = tmp161 + tmp164;
	       tmp167 = (K831469612 * tmp159) - (K555570233 * tmp160);
	       tmp168 = (K195090322 * tmp162) + (K980785280 * tmp163);
	       tmp169 = K2_000000000 * (tmp167 - tmp168);
	       tmp173 = tmp167 + tmp168;
	  }
	  output[22 * ostride] = tmp158 - tmp165;
	  output[54 * ostride] = tmp158 + tmp165;
	  output[38 * ostride] = tmp166 - tmp169;
	  output[6 * ostride] = tmp166 + tmp169;
	  tmp177 = K1_414213562 * (tmp174 + tmp173);
	  output[30 * ostride] = tmp176 - tmp177;
	  output[62 * ostride] = tmp176 + tmp177;
	  tmp175 = K1_414213562 * (tmp173 - tmp174);
	  output[46 * ostride] = tmp172 - tmp175;
	  output[14 * ostride] = tmp172 + tmp175;
     }
     {
	  fftw_real tmp251;
	  fftw_real tmp283;
	  fftw_real tmp295;
	  fftw_real tmp299;
	  fftw_real tmp282;
	  fftw_real tmp296;
	  fftw_real tmp290;
	  fftw_real tmp297;
	  fftw_real tmp300;
	  fftw_real tmp298;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp234;
	       fftw_real tmp250;
	       fftw_real tmp291;
	       fftw_real tmp294;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp234 = K2_000000000 * (tmp226 + tmp233);
	       tmp250 = tmp241 + tmp249;
	       tmp251 = tmp234 + tmp250;
	       tmp283 = tmp250 - tmp234;
	       tmp291 = tmp249 - tmp241;
	       tmp294 = K2_000000000 * (tmp292 - tmp293);
	       tmp295 = tmp291 - tmp294;
	       tmp299 = tmp294 + tmp291;
	  }
	  {
	       fftw_real tmp266;
	       fftw_real tmp281;
	       fftw_real tmp286;
	       fftw_real tmp289;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp266 = tmp258 + tmp265;
	       tmp281 = tmp273 + tmp280;
	       tmp282 = K2_000000000 * (tmp266 + tmp281);
	       tmp296 = tmp266 - tmp281;
	       tmp286 = tmp284 - tmp285;
	       tmp289 = tmp287 - tmp288;
	       tmp290 = K2_000000000 * (tmp286 - tmp289);
	       tmp297 = tmp286 + tmp289;
	  }
	  output[33 * ostride] = tmp251 - tmp282;
	  output[ostride] = tmp282 + tmp251;
	  output[17 * ostride] = tmp283 - tmp290;
	  output[49 * ostride] = tmp290 + tmp283;
	  tmp300 = K1_414213562 * (tmp296 + tmp297);
	  output[25 * ostride] = tmp299 - tmp300;
	  output[57 * ostride] = tmp300 + tmp299;
	  tmp298 = K1_414213562 * (tmp296 - tmp297);
	  output[41 * ostride] = tmp295 - tmp298;
	  output[9 * ostride] = tmp298 + tmp295;
     }
     {
	  fftw_real tmp306;
	  fftw_real tmp317;
	  fftw_real tmp303;
	  fftw_real tmp318;
	  fftw_real tmp310;
	  fftw_real tmp320;
	  fftw_real tmp313;
	  fftw_real tmp321;
	  fftw_real tmp301;
	  fftw_real tmp302;
	  ASSERT_ALIGNED_DOUBLE;
	  tmp306 = tmp304 - tmp305;
	  tmp317 = tmp305 + tmp304;
	  tmp301 = tmp233 - tmp226;
	  tmp302 = tmp293 + tmp292;
	  tmp303 = K1_414213562 * (tmp301 - tmp302);
	  tmp318 = K1_414213562 * (tmp301 + tmp302);
	  {
	       fftw_real tmp308;
	       fftw_real tmp309;
	       fftw_real tmp311;
	       fftw_real tmp312;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp308 = tmp265 - tmp258;
	       tmp309 = tmp288 + tmp287;
	       tmp310 = tmp308 - tmp309;
	       tmp320 = tmp308 + tmp309;
	       tmp311 = tmp285 + tmp284;
	       tmp312 = tmp273 - tmp280;
	       tmp313 = tmp311 + tmp312;
	       tmp321 = tmp311 - tmp312;
	  }
	  {
	       fftw_real tmp307;
	       fftw_real tmp314;
	       fftw_real tmp315;
	       fftw_real tmp316;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp307 = tmp303 + tmp306;
	       tmp314 = (K1_847759065 * tmp310) - (K765366864 * tmp313);
	       output[37 * ostride] = tmp307 - tmp314;
	       output[5 * ostride] = tmp314 + tmp307;
	       tmp315 = tmp306 - tmp303;
	       tmp316 = (K1_847759065 * tmp313) + (K765366864 * tmp310);
	       output[21 * ostride] = tmp315 - tmp316;
	       output[53 * ostride] = tmp316 + tmp315;
	  }
	  {
	       fftw_real tmp319;
	       fftw_real tmp322;
	       fftw_real tmp323;
	       fftw_real tmp324;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp319 = tmp317 - tmp318;
	       tmp322 = (K765366864 * tmp320) - (K1_847759065 * tmp321);
	       output[45 * ostride] = tmp319 - tmp322;
	       output[13 * ostride] = tmp322 + tmp319;
	       tmp323 = tmp318 + tmp317;
	       tmp324 = (K765366864 * tmp321) + (K1_847759065 * tmp320);
	       output[29 * ostride] = tmp323 - tmp324;
	       output[61 * ostride] = tmp324 + tmp323;
	  }
     }
}

fftw_codelet_desc fftw_hc2real_64_desc = {
     "fftw_hc2real_64",
     (void (*)()) fftw_hc2real_64,
     64,
     FFTW_BACKWARD,
     FFTW_HC2REAL,
     1423,
     0,
     (const int *) 0,
};