frc_64.c

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

	       fftw_real tmp70;
	       fftw_real tmp73;
	       fftw_real tmp174;
	       fftw_real tmp175;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp70 = (K923879532 * tmp68) - (K382683432 * tmp69);
	       tmp73 = (K923879532 * tmp71) + (K382683432 * tmp72);
	       tmp74 = tmp70 + tmp73;
	       tmp239 = tmp73 - tmp70;
	       tmp174 = (K382683432 * tmp71) - (K923879532 * tmp72);
	       tmp175 = (K382683432 * tmp68) + (K923879532 * tmp69);
	       tmp176 = tmp174 - tmp175;
	       tmp210 = tmp175 + tmp174;
	  }
     }
     {
	  fftw_real tmp113;
	  fftw_real tmp293;
	  fftw_real tmp125;
	  fftw_real tmp282;
	  fftw_real tmp116;
	  fftw_real tmp294;
	  fftw_real tmp120;
	  fftw_real tmp283;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp111;
	       fftw_real tmp112;
	       fftw_real tmp123;
	       fftw_real tmp124;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp111 = input[57 * istride];
	       tmp112 = input[25 * istride];
	       tmp113 = tmp111 - tmp112;
	       tmp293 = tmp111 + tmp112;
	       tmp123 = input[istride];
	       tmp124 = input[33 * istride];
	       tmp125 = tmp123 - tmp124;
	       tmp282 = tmp123 + tmp124;
	  }
	  {
	       fftw_real tmp114;
	       fftw_real tmp115;
	       fftw_real tmp118;
	       fftw_real tmp119;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp114 = input[9 * istride];
	       tmp115 = input[41 * istride];
	       tmp116 = tmp114 - tmp115;
	       tmp294 = tmp114 + tmp115;
	       tmp118 = input[17 * istride];
	       tmp119 = input[49 * istride];
	       tmp120 = tmp118 - tmp119;
	       tmp283 = tmp118 + tmp119;
	  }
	  {
	       fftw_real tmp352;
	       fftw_real tmp353;
	       fftw_real tmp117;
	       fftw_real tmp126;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp284 = tmp282 - tmp283;
	       tmp352 = tmp282 + tmp283;
	       tmp353 = tmp294 + tmp293;
	       tmp354 = tmp352 + tmp353;
	       tmp378 = tmp352 - tmp353;
	       tmp295 = tmp293 - tmp294;
	       tmp117 = K707106781 * (tmp113 - tmp116);
	       tmp121 = tmp117 - tmp120;
	       tmp223 = tmp120 + tmp117;
	       tmp126 = K707106781 * (tmp116 + tmp113);
	       tmp127 = tmp125 + tmp126;
	       tmp220 = tmp125 - tmp126;
	  }
     }
     {
	  fftw_real tmp98;
	  fftw_real tmp288;
	  fftw_real tmp108;
	  fftw_real tmp286;
	  fftw_real tmp101;
	  fftw_real tmp289;
	  fftw_real tmp105;
	  fftw_real tmp285;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp96;
	       fftw_real tmp97;
	       fftw_real tmp106;
	       fftw_real tmp107;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp96 = input[61 * istride];
	       tmp97 = input[29 * istride];
	       tmp98 = tmp96 - tmp97;
	       tmp288 = tmp96 + tmp97;
	       tmp106 = input[21 * istride];
	       tmp107 = input[53 * istride];
	       tmp108 = tmp106 - tmp107;
	       tmp286 = tmp106 + tmp107;
	  }
	  {
	       fftw_real tmp99;
	       fftw_real tmp100;
	       fftw_real tmp103;
	       fftw_real tmp104;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp99 = input[13 * istride];
	       tmp100 = input[45 * istride];
	       tmp101 = tmp99 - tmp100;
	       tmp289 = tmp99 + tmp100;
	       tmp103 = input[5 * istride];
	       tmp104 = input[37 * istride];
	       tmp105 = tmp103 - tmp104;
	       tmp285 = tmp103 + tmp104;
	  }
	  {
	       fftw_real tmp102;
	       fftw_real tmp109;
	       fftw_real tmp128;
	       fftw_real tmp129;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp102 = (K382683432 * tmp98) - (K923879532 * tmp101);
	       tmp109 = (K382683432 * tmp105) + (K923879532 * tmp108);
	       tmp110 = tmp102 - tmp109;
	       tmp221 = tmp109 + tmp102;
	       tmp128 = (K923879532 * tmp105) - (K382683432 * tmp108);
	       tmp129 = (K923879532 * tmp98) + (K382683432 * tmp101);
	       tmp130 = tmp128 + tmp129;
	       tmp224 = tmp129 - tmp128;
	  }
	  {
	       fftw_real tmp287;
	       fftw_real tmp290;
	       fftw_real tmp355;
	       fftw_real tmp356;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp287 = tmp285 - tmp286;
	       tmp290 = tmp288 - tmp289;
	       tmp291 = K707106781 * (tmp287 + tmp290);
	       tmp296 = K707106781 * (tmp290 - tmp287);
	       tmp355 = tmp285 + tmp286;
	       tmp356 = tmp288 + tmp289;
	       tmp357 = tmp355 + tmp356;
	       tmp379 = tmp356 - tmp355;
	  }
     }
     {
	  fftw_real tmp34;
	  fftw_real tmp90;
	  fftw_real tmp44;
	  fftw_real tmp85;
	  fftw_real tmp37;
	  fftw_real tmp88;
	  fftw_real tmp41;
	  fftw_real tmp86;
	  fftw_real tmp87;
	  fftw_real tmp91;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp32;
	       fftw_real tmp33;
	       fftw_real tmp42;
	       fftw_real tmp43;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp32 = input[2 * istride];
	       tmp33 = input[34 * istride];
	       tmp34 = tmp32 + tmp33;
	       tmp90 = tmp32 - tmp33;
	       tmp42 = input[58 * istride];
	       tmp43 = input[26 * istride];
	       tmp44 = tmp42 + tmp43;
	       tmp85 = tmp42 - tmp43;
	  }
	  {
	       fftw_real tmp35;
	       fftw_real tmp36;
	       fftw_real tmp39;
	       fftw_real tmp40;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp35 = input[18 * istride];
	       tmp36 = input[50 * istride];
	       tmp37 = tmp35 + tmp36;
	       tmp88 = tmp35 - tmp36;
	       tmp39 = input[10 * istride];
	       tmp40 = input[42 * istride];
	       tmp41 = tmp39 + tmp40;
	       tmp86 = tmp39 - tmp40;
	  }
	  tmp87 = K707106781 * (tmp85 - tmp86);
	  tmp89 = tmp87 - tmp88;
	  tmp213 = tmp88 + tmp87;
	  tmp91 = K707106781 * (tmp86 + tmp85);
	  tmp92 = tmp90 + tmp91;
	  tmp212 = tmp90 - tmp91;
	  {
	       fftw_real tmp38;
	       fftw_real tmp45;
	       fftw_real tmp274;
	       fftw_real tmp275;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp38 = tmp34 + tmp37;
	       tmp45 = tmp41 + tmp44;
	       tmp46 = tmp38 + tmp45;
	       tmp374 = tmp38 - tmp45;
	       tmp274 = tmp34 - tmp37;
	       tmp275 = tmp44 - tmp41;
	       tmp276 = (K923879532 * tmp274) + (K382683432 * tmp275);
	       tmp317 = (K923879532 * tmp275) - (K382683432 * tmp274);
	  }
     }
     {
	  fftw_real tmp49;
	  fftw_real tmp76;
	  fftw_real tmp59;
	  fftw_real tmp78;
	  fftw_real tmp52;
	  fftw_real tmp82;
	  fftw_real tmp56;
	  fftw_real tmp77;
	  fftw_real tmp79;
	  fftw_real tmp81;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp47;
	       fftw_real tmp48;
	       fftw_real tmp57;
	       fftw_real tmp58;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp47 = input[62 * istride];
	       tmp48 = input[30 * istride];
	       tmp49 = tmp47 + tmp48;
	       tmp76 = tmp47 - tmp48;
	       tmp57 = input[54 * istride];
	       tmp58 = input[22 * istride];
	       tmp59 = tmp57 + tmp58;
	       tmp78 = tmp57 - tmp58;
	  }
	  {
	       fftw_real tmp50;
	       fftw_real tmp51;
	       fftw_real tmp54;
	       fftw_real tmp55;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp50 = input[14 * istride];
	       tmp51 = input[46 * istride];
	       tmp52 = tmp50 + tmp51;
	       tmp82 = tmp50 - tmp51;
	       tmp54 = input[6 * istride];
	       tmp55 = input[38 * istride];
	       tmp56 = tmp54 + tmp55;
	       tmp77 = tmp54 - tmp55;
	  }
	  tmp79 = K707106781 * (tmp77 + tmp78);
	  tmp80 = tmp76 + tmp79;
	  tmp215 = tmp76 - tmp79;
	  tmp81 = K707106781 * (tmp78 - tmp77);
	  tmp83 = tmp81 - tmp82;
	  tmp216 = tmp82 + tmp81;
	  {
	       fftw_real tmp53;
	       fftw_real tmp60;
	       fftw_real tmp277;
	       fftw_real tmp278;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp53 = tmp49 + tmp52;
	       tmp60 = tmp56 + tmp59;
	       tmp61 = tmp53 + tmp60;
	       tmp375 = tmp53 - tmp60;
	       tmp277 = tmp49 - tmp52;
	       tmp278 = tmp59 - tmp56;
	       tmp279 = (K923879532 * tmp277) - (K382683432 * tmp278);
	       tmp318 = (K382683432 * tmp277) + (K923879532 * tmp278);
	  }
     }
     {
	  fftw_real tmp31;
	  fftw_real tmp62;
	  fftw_real tmp369;
	  fftw_real tmp370;
	  fftw_real tmp371;
	  fftw_real tmp372;
	  ASSERT_ALIGNED_DOUBLE;
	  tmp31 = tmp15 + tmp30;
	  tmp62 = tmp46 + tmp61;
	  tmp369 = tmp31 + tmp62;
	  tmp370 = tmp354 + tmp357;
	  tmp371 = tmp361 + tmp364;
	  tmp372 = tmp370 + tmp371;
	  real_output[16 * real_ostride] = tmp31 - tmp62;
	  imag_output[16 * imag_ostride] = tmp371 - tmp370;
	  real_output[32 * real_ostride] = tmp369 - tmp372;
	  real_output[0] = tmp369 + tmp372;
     }
     {
	  fftw_real tmp351;
	  fftw_real tmp367;
	  fftw_real tmp366;
	  fftw_real tmp368;
	  fftw_real tmp358;
	  fftw_real tmp365;
	  ASSERT_ALIGNED_DOUBLE;
	  tmp351 = tmp15 - tmp30;
	  tmp367 = tmp61 - tmp46;
	  tmp358 = tmp354 - tmp357;
	  tmp365 = tmp361 - tmp364;
	  tmp366 = K707106781 * (tmp358 + tmp365);
	  tmp368 = K707106781 * (tmp365 - tmp358);
	  real_output[24 * real_ostride] = tmp351 - tmp366;
	  real_output[8 * real_ostride] = tmp351 + tmp366;
	  imag_output[8 * imag_ostride] = tmp367 + tmp368;
	  imag_output[24 * imag_ostride] = tmp368 - tmp367;
     }
     {
	  fftw_real tmp75;
	  fftw_real tmp189;
	  fftw_real tmp94;
	  fftw_real tmp199;
	  fftw_real tmp173;
	  fftw_real tmp190;
	  fftw_real tmp132;
	  fftw_real tmp184;
	  fftw_real tmp180;
	  fftw_real tmp200;
	  fftw_real tmp194;
	  fftw_real tmp204;
	  fftw_real tmp169;
	  fftw_real tmp185;
	  fftw_real tmp197;
	  fftw_real tmp205;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp84;
	       fftw_real tmp93;
	       fftw_real tmp192;
	       fftw_real tmp193;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp75 = tmp67 - tmp74;
	       tmp189 = tmp67 + tmp74;
	       tmp84 = (K195090322 * tmp80) + (K980785280 * tmp83);
	       tmp93 = (K980785280 * tmp89) - (K195090322 * tmp92);
	       tmp94 = tmp84 - tmp93;
	       tmp199 = tmp93 + tmp84;
	       {
		    fftw_real tmp171;
		    fftw_real tmp172;
		    fftw_real tmp122;
		    fftw_real tmp131;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp171 = (K980785280 * tmp80) - (K195090322 * tmp83);
		    tmp172 = (K980785280 * tmp92) + (K195090322 * tmp89);
		    tmp173 = tmp171 - tmp172;
		    tmp190 = tmp172 + tmp171;
		    tmp122 = tmp110 - tmp121;
		    tmp131 = tmp127 - tmp130;
		    tmp132 = (K634393284 * tmp122) + (K773010453 * tmp131);
		    tmp184 = (K773010453 * tmp122) - (K634393284 * tmp131);
	       }
	       tmp180 = tmp176 - tmp179;
	       tmp200 = tmp179 + tmp176;
	       tmp192 = tmp121 + tmp110;
	       tmp193 = tmp127 + tmp130;
	       tmp194 = (K098017140 * tmp192) + (K995184726 * tmp193);
	       tmp204 = (K995184726 * tmp192) - (K098017140 * tmp193);
	       {
		    fftw_real tmp159;
		    fftw_real tmp168;
		    fftw_real tmp195;
		    fftw_real tmp196;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp159 = tmp143 - tmp158;
		    tmp168 = tmp162 - tmp167;
		    tmp169 = (K773010453 * tmp159) - (K634393284 * tmp168);
		    tmp185 = (K773010453 * tmp168) + (K634393284 * tmp159);
		    tmp195 = tmp143 + tmp158;
		    tmp196 = tmp167 + tmp162;
		    tmp197 = (K995184726 * tmp195) - (K098017140 * tmp196);