fn_64.c

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

	  {
	       fftw_real tmp317;
	       fftw_real tmp328;
	       fftw_real tmp517;
	       fftw_real tmp518;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp317 = tmp313 - tmp316;
	       tmp328 = K707106781 * (tmp322 - tmp327);
	       tmp329 = tmp317 - tmp328;
	       tmp429 = tmp317 + tmp328;
	       tmp517 = tmp313 + tmp316;
	       tmp518 = K707106781 * (tmp360 + tmp359);
	       tmp519 = tmp517 - tmp518;
	       tmp593 = tmp517 + tmp518;
	  }
	  {
	       fftw_real tmp358;
	       fftw_real tmp361;
	       fftw_real tmp528;
	       fftw_real tmp529;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp358 = tmp354 + tmp357;
	       tmp361 = K707106781 * (tmp359 - tmp360);
	       tmp362 = tmp358 - tmp361;
	       tmp432 = tmp358 + tmp361;
	       tmp528 = tmp357 - tmp354;
	       tmp529 = K707106781 * (tmp322 + tmp327);
	       tmp530 = tmp528 - tmp529;
	       tmp596 = tmp528 + tmp529;
	  }
     }
     {
	  fftw_real tmp18;
	  fftw_real tmp205;
	  fftw_real tmp145;
	  fftw_real tmp203;
	  fftw_real tmp21;
	  fftw_real tmp202;
	  fftw_real tmp148;
	  fftw_real tmp206;
	  fftw_real tmp25;
	  fftw_real tmp212;
	  fftw_real tmp152;
	  fftw_real tmp210;
	  fftw_real tmp28;
	  fftw_real tmp209;
	  fftw_real tmp155;
	  fftw_real tmp213;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp16;
	       fftw_real tmp17;
	       fftw_real tmp143;
	       fftw_real tmp144;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp16 = c_re(input[4 * istride]);
	       tmp17 = c_re(input[36 * istride]);
	       tmp18 = tmp16 + tmp17;
	       tmp205 = tmp16 - tmp17;
	       tmp143 = c_im(input[4 * istride]);
	       tmp144 = c_im(input[36 * istride]);
	       tmp145 = tmp143 + tmp144;
	       tmp203 = tmp143 - tmp144;
	  }
	  {
	       fftw_real tmp19;
	       fftw_real tmp20;
	       fftw_real tmp146;
	       fftw_real tmp147;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp19 = c_re(input[20 * istride]);
	       tmp20 = c_re(input[52 * istride]);
	       tmp21 = tmp19 + tmp20;
	       tmp202 = tmp19 - tmp20;
	       tmp146 = c_im(input[20 * istride]);
	       tmp147 = c_im(input[52 * istride]);
	       tmp148 = tmp146 + tmp147;
	       tmp206 = tmp146 - tmp147;
	  }
	  {
	       fftw_real tmp23;
	       fftw_real tmp24;
	       fftw_real tmp150;
	       fftw_real tmp151;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp23 = c_re(input[60 * istride]);
	       tmp24 = c_re(input[28 * istride]);
	       tmp25 = tmp23 + tmp24;
	       tmp212 = tmp23 - tmp24;
	       tmp150 = c_im(input[60 * istride]);
	       tmp151 = c_im(input[28 * istride]);
	       tmp152 = tmp150 + tmp151;
	       tmp210 = tmp150 - tmp151;
	  }
	  {
	       fftw_real tmp26;
	       fftw_real tmp27;
	       fftw_real tmp153;
	       fftw_real tmp154;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp26 = c_re(input[12 * istride]);
	       tmp27 = c_re(input[44 * istride]);
	       tmp28 = tmp26 + tmp27;
	       tmp209 = tmp26 - tmp27;
	       tmp153 = c_im(input[12 * istride]);
	       tmp154 = c_im(input[44 * istride]);
	       tmp155 = tmp153 + tmp154;
	       tmp213 = tmp153 - tmp154;
	  }
	  {
	       fftw_real tmp22;
	       fftw_real tmp29;
	       fftw_real tmp638;
	       fftw_real tmp639;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp22 = tmp18 + tmp21;
	       tmp29 = tmp25 + tmp28;
	       tmp30 = tmp22 + tmp29;
	       tmp850 = tmp29 - tmp22;
	       tmp638 = tmp145 - tmp148;
	       tmp639 = tmp18 - tmp21;
	       tmp640 = tmp638 - tmp639;
	       tmp721 = tmp639 + tmp638;
	  }
	  {
	       fftw_real tmp149;
	       fftw_real tmp156;
	       fftw_real tmp641;
	       fftw_real tmp642;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp149 = tmp145 + tmp148;
	       tmp156 = tmp152 + tmp155;
	       tmp157 = tmp149 + tmp156;
	       tmp816 = tmp149 - tmp156;
	       tmp641 = tmp25 - tmp28;
	       tmp642 = tmp152 - tmp155;
	       tmp643 = tmp641 + tmp642;
	       tmp720 = tmp641 - tmp642;
	  }
	  {
	       fftw_real tmp204;
	       fftw_real tmp207;
	       fftw_real tmp474;
	       fftw_real tmp475;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp204 = tmp202 + tmp203;
	       tmp207 = tmp205 - tmp206;
	       tmp208 = (K382683432 * tmp204) - (K923879532 * tmp207);
	       tmp377 = (K923879532 * tmp204) + (K382683432 * tmp207);
	       tmp474 = tmp203 - tmp202;
	       tmp475 = tmp205 + tmp206;
	       tmp476 = (K923879532 * tmp474) - (K382683432 * tmp475);
	       tmp541 = (K382683432 * tmp474) + (K923879532 * tmp475);
	  }
	  {
	       fftw_real tmp211;
	       fftw_real tmp214;
	       fftw_real tmp477;
	       fftw_real tmp478;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp211 = tmp209 + tmp210;
	       tmp214 = tmp212 - tmp213;
	       tmp215 = (K382683432 * tmp211) + (K923879532 * tmp214);
	       tmp376 = (K382683432 * tmp214) - (K923879532 * tmp211);
	       tmp477 = tmp210 - tmp209;
	       tmp478 = tmp212 + tmp213;
	       tmp479 = (K923879532 * tmp477) + (K382683432 * tmp478);
	       tmp540 = (K923879532 * tmp478) - (K382683432 * tmp477);
	  }
     }
     {
	  fftw_real tmp112;
	  fftw_real tmp694;
	  fftw_real tmp115;
	  fftw_real tmp695;
	  fftw_real tmp334;
	  fftw_real tmp520;
	  fftw_real tmp339;
	  fftw_real tmp521;
	  fftw_real tmp697;
	  fftw_real tmp696;
	  fftw_real tmp119;
	  fftw_real tmp700;
	  fftw_real tmp122;
	  fftw_real tmp701;
	  fftw_real tmp345;
	  fftw_real tmp523;
	  fftw_real tmp350;
	  fftw_real tmp524;
	  fftw_real tmp702;
	  fftw_real tmp699;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp335;
	       fftw_real tmp333;
	       fftw_real tmp330;
	       fftw_real tmp338;
	       ASSERT_ALIGNED_DOUBLE;
	       {
		    fftw_real tmp110;
		    fftw_real tmp111;
		    fftw_real tmp331;
		    fftw_real tmp332;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp110 = c_re(input[3 * istride]);
		    tmp111 = c_re(input[35 * istride]);
		    tmp112 = tmp110 + tmp111;
		    tmp335 = tmp110 - tmp111;
		    tmp331 = c_im(input[3 * istride]);
		    tmp332 = c_im(input[35 * istride]);
		    tmp333 = tmp331 - tmp332;
		    tmp694 = tmp331 + tmp332;
	       }
	       {
		    fftw_real tmp113;
		    fftw_real tmp114;
		    fftw_real tmp336;
		    fftw_real tmp337;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp113 = c_re(input[19 * istride]);
		    tmp114 = c_re(input[51 * istride]);
		    tmp115 = tmp113 + tmp114;
		    tmp330 = tmp113 - tmp114;
		    tmp336 = c_im(input[19 * istride]);
		    tmp337 = c_im(input[51 * istride]);
		    tmp338 = tmp336 - tmp337;
		    tmp695 = tmp336 + tmp337;
	       }
	       tmp334 = tmp330 + tmp333;
	       tmp520 = tmp333 - tmp330;
	       tmp339 = tmp335 - tmp338;
	       tmp521 = tmp335 + tmp338;
	       tmp697 = tmp112 - tmp115;
	       tmp696 = tmp694 - tmp695;
	  }
	  {
	       fftw_real tmp346;
	       fftw_real tmp344;
	       fftw_real tmp341;
	       fftw_real tmp349;
	       ASSERT_ALIGNED_DOUBLE;
	       {
		    fftw_real tmp117;
		    fftw_real tmp118;
		    fftw_real tmp342;
		    fftw_real tmp343;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp117 = c_re(input[59 * istride]);
		    tmp118 = c_re(input[27 * istride]);
		    tmp119 = tmp117 + tmp118;
		    tmp346 = tmp117 - tmp118;
		    tmp342 = c_im(input[59 * istride]);
		    tmp343 = c_im(input[27 * istride]);
		    tmp344 = tmp342 - tmp343;
		    tmp700 = tmp342 + tmp343;
	       }
	       {
		    fftw_real tmp120;
		    fftw_real tmp121;
		    fftw_real tmp347;
		    fftw_real tmp348;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp120 = c_re(input[11 * istride]);
		    tmp121 = c_re(input[43 * istride]);
		    tmp122 = tmp120 + tmp121;
		    tmp341 = tmp120 - tmp121;
		    tmp347 = c_im(input[11 * istride]);
		    tmp348 = c_im(input[43 * istride]);
		    tmp349 = tmp347 - tmp348;
		    tmp701 = tmp347 + tmp348;
	       }
	       tmp345 = tmp341 + tmp344;
	       tmp523 = tmp344 - tmp341;
	       tmp350 = tmp346 - tmp349;
	       tmp524 = tmp346 + tmp349;
	       tmp702 = tmp700 - tmp701;
	       tmp699 = tmp119 - tmp122;
	  }
	  {
	       fftw_real tmp116;
	       fftw_real tmp123;
	       fftw_real tmp363;
	       fftw_real tmp364;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp116 = tmp112 + tmp115;
	       tmp123 = tmp119 + tmp122;
	       tmp124 = tmp116 + tmp123;
	       tmp845 = tmp123 - tmp116;
	       tmp363 = (K382683432 * tmp350) - (K923879532 * tmp345);
	       tmp364 = (K923879532 * tmp334) + (K382683432 * tmp339);
	       tmp365 = tmp363 - tmp364;
	       tmp430 = tmp364 + tmp363;
	  }
	  {
	       fftw_real tmp340;
	       fftw_real tmp351;
	       fftw_real tmp838;
	       fftw_real tmp839;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp340 = (K382683432 * tmp334) - (K923879532 * tmp339);
	       tmp351 = (K382683432 * tmp345) + (K923879532 * tmp350);
	       tmp352 = tmp340 - tmp351;
	       tmp433 = tmp340 + tmp351;
	       tmp838 = tmp694 + tmp695;
	       tmp839 = tmp700 + tmp701;
	       tmp840 = tmp838 - tmp839;
	       tmp893 = tmp838 + tmp839;
	  }
	  {
	       fftw_real tmp522;
	       fftw_real tmp525;
	       fftw_real tmp531;
	       fftw_real tmp532;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp522 = (K923879532 * tmp520) - (K382683432 * tmp521);
	       tmp525 = (K923879532 * tmp523) + (K382683432 * tmp524);
	       tmp526 = tmp522 - tmp525;
	       tmp597 = tmp522 + tmp525;
	       tmp531 = (K923879532 * tmp524) - (K382683432 * tmp523);
	       tmp532 = (K382683432 * tmp520) + (K923879532 * tmp521);
	       tmp533 = tmp531 - tmp532;
	       tmp594 = tmp532 + tmp531;
	  }
	  {
	       fftw_real tmp698;
	       fftw_real tmp703;
	       fftw_real tmp711;
	       fftw_real tmp712;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp698 = tmp696 - tmp697;
	       tmp703 = tmp699 + tmp702;
	       tmp704 = K707106781 * (tmp698 - tmp703);
	       tmp777 = K707106781 * (tmp698 + tmp703);
	       tmp711 = tmp699 - tmp702;
	       tmp712 = tmp697 + tmp696;
	       tmp713 = K707106781 * (tmp711 - tmp712);
	       tmp774 = K707106781 * (tmp712 + tmp711);
	  }
     }
     {
	  fftw_real tmp34;
	  fftw_real tmp229;
	  fftw_real tmp161;
	  fftw_real tmp219;
	  fftw_real tmp37;
	  fftw_real tmp218;
	  fftw_real tmp164;
	  fftw_real tmp230;
	  fftw_real tmp44;
	  fftw_real tmp233;
	  fftw_real tmp223;
	  fftw_real tmp171;
	  fftw_real tmp41;
	  fftw_real tmp232;
	  fftw_real tmp226;
	  fftw_real tmp168;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp32;
	       fftw_real tmp33;
	       fftw_real tmp162;
	       fftw_real tmp163;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp32 = c_re(input[2 * istride]);
	       tmp33 = c_re(input[34 * istride]);
	       tmp34 = tmp32 + tmp33;
	       tmp229 = tmp32 - tmp33;
	       {
		    fftw_real tmp159;
		    fftw_real tmp160;
		    fftw_real tmp35;
		    fftw_real tmp36;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp159 = c_im(input[2 * istride]);
		    tmp160 = c_im(input[34 * istride]);
		    tmp161 = tmp159 + tmp160;
		    tmp219 = tmp159 - tmp160;
		    tmp35 = c_re(input[18 * istride]);
		    tmp36 = c_re(input[50 * istride]);
		    tmp37 = tmp35 + tmp36;
		    tmp218 = tmp35 - tmp36;
	       }
	       tmp162 = c_im(input[18 * istride]);
	       tmp163 = c_im(input[50 * istride]);
	       tmp164 = tmp162 + tmp163;
	       tmp230 = tmp162 - tmp163;
	       {
		    fftw_real tmp42;
		    fftw_real tmp43;
		    fftw_real tmp221;
		    fftw_real tmp169;
		    fftw_real tmp170;
		    fftw_real tmp222;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp42 = c_re(input[58 * istride]);
		    tmp43 = c_re(input[26 * istride]);
		    tmp221 = tmp42 - tmp43;
		    tmp169 = c_im(input[58 * istride]);
		    tmp170 = c_im(input[26 * istride]);
		    tmp222 = tmp169 - tmp170;
		    tmp44 = tmp42 + tmp43;
		    tmp233 = tmp221 + tmp222;
		    tmp223 = tmp221 - tmp222;
		    tmp171 = tmp169 + tmp170;
	       }
	       {
		    fftw_real tmp39;
		    fftw_real tmp40;
		    fftw_real tmp224;
		    fftw_real tmp166;
		    fftw_real tmp167;
		    fftw_real tmp225;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp39 = c_re(input[10 * istride]);
		    tmp40 = c_re(input[42 * istride]);
		    tmp224 = tmp39 - tmp40;
		    tmp166 = c_im(input[10 * istride]);
		    tmp167 = c_im(input[42 * istride]);
		    tmp225 = tmp166 - tmp167;
		    tmp41 = tmp39 + tmp40;
		    tmp232 = tmp225 - tmp224;
		    tmp226 = tmp224 + tmp225;
		    tmp168 = tmp166 + tmp167;
	       }
	  }
	  {
	       fftw_real tmp38;
	       fftw_real tmp45;
	       fftw_real tmp646;
	       fftw_real tmp647;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp38 = tmp34 + tmp37;
	       tmp45 = tmp41 + tmp44;
	       tmp46 = tmp38 + tmp45;
	       tmp819 = tmp38 - tmp45;
	       tmp646 = tmp161 - tmp164;
	       tmp647 = tmp44 - tmp41;
	       tmp648 = tmp646 - tmp647;
	       tmp758 = tmp647 + tmp646;
	  }
	  {
	       fftw_real tmp165;
	       fftw_real tmp172;
	       fftw_real tmp649;
	       fftw_real tmp650;
	       ASSERT_ALIGNED_DOUBLE;