fcr_64.c

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

	       fftw_real tmp73;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp14 = real_input[12 * real_istride];
	       tmp15 = real_input[20 * real_istride];
	       tmp16 = tmp14 + tmp15;
	       tmp238 = tmp15 - tmp14;
	       tmp72 = imag_input[12 * imag_istride];
	       tmp73 = imag_input[20 * imag_istride];
	       tmp74 = tmp72 - tmp73;
	       tmp236 = tmp72 + tmp73;
	  }
	  tmp17 = K2_000000000 * (tmp13 + tmp16);
	  tmp71 = tmp13 - tmp16;
	  tmp78 = tmp74 + tmp77;
	  tmp79 = K1_414213562 * (tmp71 - tmp78);
	  tmp153 = K1_414213562 * (tmp71 + tmp78);
	  tmp197 = K2_000000000 * (tmp77 - tmp74);
	  {
	       fftw_real tmp237;
	       fftw_real tmp240;
	       fftw_real tmp332;
	       fftw_real tmp333;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp237 = tmp235 - tmp236;
	       tmp240 = tmp238 + tmp239;
	       tmp241 = (K1_847759065 * tmp237) - (K765366864 * tmp240);
	       tmp305 = (K765366864 * tmp237) + (K1_847759065 * tmp240);
	       tmp332 = tmp235 + tmp236;
	       tmp333 = tmp239 - tmp238;
	       tmp334 = (K765366864 * tmp332) - (K1_847759065 * tmp333);
	       tmp377 = (K1_847759065 * tmp332) + (K765366864 * tmp333);
	  }
     }
     {
	  fftw_real tmp21;
	  fftw_real tmp227;
	  fftw_real tmp96;
	  fftw_real tmp231;
	  fftw_real tmp24;
	  fftw_real tmp230;
	  fftw_real tmp93;
	  fftw_real tmp228;
	  fftw_real tmp28;
	  fftw_real tmp220;
	  fftw_real tmp87;
	  fftw_real tmp224;
	  fftw_real tmp31;
	  fftw_real tmp223;
	  fftw_real tmp84;
	  fftw_real tmp221;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp19;
	       fftw_real tmp20;
	       fftw_real tmp94;
	       fftw_real tmp95;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp19 = real_input[2 * real_istride];
	       tmp20 = real_input[30 * real_istride];
	       tmp21 = tmp19 + tmp20;
	       tmp227 = tmp19 - tmp20;
	       tmp94 = imag_input[2 * imag_istride];
	       tmp95 = imag_input[30 * imag_istride];
	       tmp96 = tmp94 - tmp95;
	       tmp231 = tmp94 + tmp95;
	  }
	  {
	       fftw_real tmp22;
	       fftw_real tmp23;
	       fftw_real tmp91;
	       fftw_real tmp92;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp22 = real_input[14 * real_istride];
	       tmp23 = real_input[18 * real_istride];
	       tmp24 = tmp22 + tmp23;
	       tmp230 = tmp23 - tmp22;
	       tmp91 = imag_input[14 * imag_istride];
	       tmp92 = imag_input[18 * imag_istride];
	       tmp93 = tmp91 - tmp92;
	       tmp228 = tmp91 + tmp92;
	  }
	  {
	       fftw_real tmp26;
	       fftw_real tmp27;
	       fftw_real tmp85;
	       fftw_real tmp86;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp26 = real_input[6 * real_istride];
	       tmp27 = real_input[26 * real_istride];
	       tmp28 = tmp26 + tmp27;
	       tmp220 = tmp26 - tmp27;
	       tmp85 = imag_input[6 * imag_istride];
	       tmp86 = imag_input[26 * imag_istride];
	       tmp87 = tmp85 - tmp86;
	       tmp224 = tmp85 + tmp86;
	  }
	  {
	       fftw_real tmp29;
	       fftw_real tmp30;
	       fftw_real tmp82;
	       fftw_real tmp83;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp29 = real_input[10 * real_istride];
	       tmp30 = real_input[22 * real_istride];
	       tmp31 = tmp29 + tmp30;
	       tmp223 = tmp30 - tmp29;
	       tmp82 = imag_input[10 * imag_istride];
	       tmp83 = imag_input[22 * imag_istride];
	       tmp84 = tmp82 - tmp83;
	       tmp221 = tmp82 + tmp83;
	  }
	  {
	       fftw_real tmp81;
	       fftw_real tmp88;
	       fftw_real tmp90;
	       fftw_real tmp97;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp25 = tmp21 + tmp24;
	       tmp32 = tmp28 + tmp31;
	       tmp199 = tmp25 - tmp32;
	       tmp81 = tmp21 - tmp24;
	       tmp88 = tmp84 + tmp87;
	       tmp89 = tmp81 - tmp88;
	       tmp155 = tmp81 + tmp88;
	       tmp187 = tmp96 - tmp93;
	       tmp188 = tmp87 - tmp84;
	       tmp200 = tmp188 + tmp187;
	       tmp90 = tmp31 - tmp28;
	       tmp97 = tmp93 + tmp96;
	       tmp98 = tmp90 + tmp97;
	       tmp156 = tmp97 - tmp90;
	  }
	  {
	       fftw_real tmp222;
	       fftw_real tmp225;
	       fftw_real tmp328;
	       fftw_real tmp329;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp222 = tmp220 - tmp221;
	       tmp225 = tmp223 + tmp224;
	       tmp226 = (K831469612 * tmp222) - (K555570233 * tmp225);
	       tmp293 = (K555570233 * tmp222) + (K831469612 * tmp225);
	       tmp328 = tmp220 + tmp221;
	       tmp329 = tmp224 - tmp223;
	       tmp330 = (K195090322 * tmp328) + (K980785280 * tmp329);
	       tmp365 = (K980785280 * tmp328) - (K195090322 * tmp329);
	  }
	  {
	       fftw_real tmp229;
	       fftw_real tmp232;
	       fftw_real tmp325;
	       fftw_real tmp326;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp229 = tmp227 - tmp228;
	       tmp232 = tmp230 + tmp231;
	       tmp233 = (K980785280 * tmp229) - (K195090322 * tmp232);
	       tmp292 = (K195090322 * tmp229) + (K980785280 * tmp232);
	       tmp325 = tmp227 + tmp228;
	       tmp326 = tmp231 - tmp230;
	       tmp327 = (K831469612 * tmp325) - (K555570233 * tmp326);
	       tmp364 = (K555570233 * tmp325) + (K831469612 * tmp326);
	  }
     }
     {
	  fftw_real tmp37;
	  fftw_real tmp259;
	  fftw_real tmp116;
	  fftw_real tmp263;
	  fftw_real tmp40;
	  fftw_real tmp262;
	  fftw_real tmp113;
	  fftw_real tmp260;
	  fftw_real tmp44;
	  fftw_real tmp252;
	  fftw_real tmp107;
	  fftw_real tmp256;
	  fftw_real tmp47;
	  fftw_real tmp255;
	  fftw_real tmp104;
	  fftw_real tmp253;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp35;
	       fftw_real tmp36;
	       fftw_real tmp114;
	       fftw_real tmp115;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp35 = real_input[real_istride];
	       tmp36 = real_input[31 * real_istride];
	       tmp37 = tmp35 + tmp36;
	       tmp259 = tmp35 - tmp36;
	       tmp114 = imag_input[imag_istride];
	       tmp115 = imag_input[31 * imag_istride];
	       tmp116 = tmp114 - tmp115;
	       tmp263 = tmp114 + tmp115;
	  }
	  {
	       fftw_real tmp38;
	       fftw_real tmp39;
	       fftw_real tmp111;
	       fftw_real tmp112;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp38 = real_input[15 * real_istride];
	       tmp39 = real_input[17 * real_istride];
	       tmp40 = tmp38 + tmp39;
	       tmp262 = tmp39 - tmp38;
	       tmp111 = imag_input[15 * imag_istride];
	       tmp112 = imag_input[17 * imag_istride];
	       tmp113 = tmp111 - tmp112;
	       tmp260 = tmp111 + tmp112;
	  }
	  {
	       fftw_real tmp42;
	       fftw_real tmp43;
	       fftw_real tmp105;
	       fftw_real tmp106;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp42 = real_input[7 * real_istride];
	       tmp43 = real_input[25 * real_istride];
	       tmp44 = tmp42 + tmp43;
	       tmp252 = tmp42 - tmp43;
	       tmp105 = imag_input[7 * imag_istride];
	       tmp106 = imag_input[25 * imag_istride];
	       tmp107 = tmp105 - tmp106;
	       tmp256 = tmp105 + tmp106;
	  }
	  {
	       fftw_real tmp45;
	       fftw_real tmp46;
	       fftw_real tmp102;
	       fftw_real tmp103;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp45 = real_input[9 * real_istride];
	       tmp46 = real_input[23 * real_istride];
	       tmp47 = tmp45 + tmp46;
	       tmp255 = tmp46 - tmp45;
	       tmp102 = imag_input[9 * imag_istride];
	       tmp103 = imag_input[23 * imag_istride];
	       tmp104 = tmp102 - tmp103;
	       tmp253 = tmp102 + tmp103;
	  }
	  {
	       fftw_real tmp41;
	       fftw_real tmp48;
	       fftw_real tmp101;
	       fftw_real tmp108;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp41 = tmp37 + tmp40;
	       tmp48 = tmp44 + tmp47;
	       tmp49 = tmp41 + tmp48;
	       tmp203 = tmp41 - tmp48;
	       tmp101 = tmp37 - tmp40;
	       tmp108 = tmp104 + tmp107;
	       tmp109 = tmp101 - tmp108;
	       tmp159 = tmp101 + tmp108;
	  }
	  {
	       fftw_real tmp179;
	       fftw_real tmp180;
	       fftw_real tmp110;
	       fftw_real tmp117;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp179 = tmp116 - tmp113;
	       tmp180 = tmp107 - tmp104;
	       tmp181 = tmp179 - tmp180;
	       tmp207 = tmp180 + tmp179;
	       tmp110 = tmp47 - tmp44;
	       tmp117 = tmp113 + tmp116;
	       tmp118 = tmp110 + tmp117;
	       tmp160 = tmp117 - tmp110;
	  }
	  {
	       fftw_real tmp254;
	       fftw_real tmp257;
	       fftw_real tmp343;
	       fftw_real tmp344;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp254 = tmp252 - tmp253;
	       tmp257 = tmp255 + tmp256;
	       tmp258 = (K773010453 * tmp254) - (K634393284 * tmp257);
	       tmp285 = (K634393284 * tmp254) + (K773010453 * tmp257);
	       tmp343 = tmp252 + tmp253;
	       tmp344 = tmp256 - tmp255;
	       tmp345 = (K471396736 * tmp343) + (K881921264 * tmp344);
	       tmp357 = (K881921264 * tmp343) - (K471396736 * tmp344);
	  }
	  {
	       fftw_real tmp261;
	       fftw_real tmp264;
	       fftw_real tmp340;
	       fftw_real tmp341;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp261 = tmp259 - tmp260;
	       tmp264 = tmp262 + tmp263;
	       tmp265 = (K995184726 * tmp261) - (K098017140 * tmp264);
	       tmp284 = (K098017140 * tmp261) + (K995184726 * tmp264);
	       tmp340 = tmp259 + tmp260;
	       tmp341 = tmp263 - tmp262;
	       tmp342 = (K956940335 * tmp340) - (K290284677 * tmp341);
	       tmp356 = (K290284677 * tmp340) + (K956940335 * tmp341);
	  }
     }
     {
	  fftw_real tmp65;
	  fftw_real tmp185;
	  fftw_real tmp34;
	  fftw_real tmp178;
	  fftw_real tmp18;
	  fftw_real tmp33;
	  ASSERT_ALIGNED_DOUBLE;
	  tmp65 = K2_000000000 * (tmp49 + tmp64);
	  tmp185 = K2_000000000 * (tmp181 - tmp184);
	  tmp18 = tmp10 + tmp17;
	  tmp33 = K2_000000000 * (tmp25 + tmp32);
	  tmp34 = tmp18 + tmp33;
	  tmp178 = tmp18 - tmp33;
	  output[32 * ostride] = tmp34 - tmp65;
	  output[0] = tmp34 + tmp65;
	  output[16 * ostride] = tmp178 - tmp185;
	  output[48 * ostride] = tmp178 + tmp185;
     }
     {
	  fftw_real tmp190;
	  fftw_real tmp194;
	  fftw_real tmp193;
	  fftw_real tmp195;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp186;
	       fftw_real tmp189;
	       fftw_real tmp191;
	       fftw_real tmp192;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp186 = tmp10 - tmp17;
	       tmp189 = K2_000000000 * (tmp187 - tmp188);
	       tmp190 = tmp186 - tmp189;
	       tmp194 = tmp186 + tmp189;
	       tmp191 = tmp49 - tmp64;
	       tmp192 = tmp181 + tmp184;
	       tmp193 = K1_414213562 * (tmp191 - tmp192);
	       tmp195 = K1_414213562 * (tmp191 + tmp192);
	  }
	  output[40 * ostride] = tmp190 - tmp193;
	  output[8 * ostride] = tmp190 + tmp193;
	  output[24 * ostride] = tmp194 - tmp195;
	  output[56 * ostride] = tmp194 + tmp195;
     }
     {
	  fftw_real tmp100;
	  fftw_real tmp140;
	  fftw_real tmp146;
	  fftw_real tmp150;
	  fftw_real tmp139;
	  fftw_real tmp147;
	  fftw_real tmp143;
	  fftw_real tmp148;
	  fftw_real tmp151;
	  fftw_real tmp149;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp80;
	       fftw_real tmp99;
	       fftw_real tmp144;
	       fftw_real tmp145;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp80 = tmp70 + tmp79;
	       tmp99 = (K1_847759065 * tmp89) - (K765366864 * tmp98);
	       tmp100 = tmp80 + tmp99;
	       tmp140 = tmp80 - tmp99;
	       tmp144 = tmp70 - tmp79;
	       tmp145 = (K765366864 * tmp89) + (K1_847759065 * tmp98);
	       tmp146 = tmp144 - tmp145;
	       tmp150 = tmp144 + tmp145;
	  }
	  {
	       fftw_real tmp119;
	       fftw_real tmp138;