fn_32.c

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

		    tmp234 = K707106781 * (tmp172 + tmp171);
	       }
	  }
     }
     {
	  fftw_real tmp18;
	  fftw_real tmp109;
	  fftw_real tmp81;
	  fftw_real tmp107;
	  fftw_real tmp21;
	  fftw_real tmp106;
	  fftw_real tmp84;
	  fftw_real tmp110;
	  fftw_real tmp278;
	  fftw_real tmp279;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp16;
	       fftw_real tmp17;
	       fftw_real tmp79;
	       fftw_real tmp80;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp16 = c_re(input[2 * istride]);
	       tmp17 = c_re(input[18 * istride]);
	       tmp18 = tmp16 + tmp17;
	       tmp109 = tmp16 - tmp17;
	       tmp79 = c_im(input[2 * istride]);
	       tmp80 = c_im(input[18 * istride]);
	       tmp81 = tmp79 + tmp80;
	       tmp107 = tmp79 - tmp80;
	  }
	  {
	       fftw_real tmp19;
	       fftw_real tmp20;
	       fftw_real tmp82;
	       fftw_real tmp83;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp19 = c_re(input[10 * istride]);
	       tmp20 = c_re(input[26 * istride]);
	       tmp21 = tmp19 + tmp20;
	       tmp106 = tmp19 - tmp20;
	       tmp82 = c_im(input[10 * istride]);
	       tmp83 = c_im(input[26 * istride]);
	       tmp84 = tmp82 + tmp83;
	       tmp110 = tmp82 - tmp83;
	  }
	  tmp22 = tmp18 + tmp21;
	  tmp278 = tmp81 - tmp84;
	  tmp279 = tmp18 - tmp21;
	  tmp280 = tmp278 - tmp279;
	  tmp313 = tmp279 + tmp278;
	  tmp85 = tmp81 + tmp84;
	  {
	       fftw_real tmp108;
	       fftw_real tmp111;
	       fftw_real tmp218;
	       fftw_real tmp219;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp108 = tmp106 + tmp107;
	       tmp111 = tmp109 - tmp110;
	       tmp112 = (K382683432 * tmp108) - (K923879532 * tmp111);
	       tmp185 = (K923879532 * tmp108) + (K382683432 * tmp111);
	       tmp218 = tmp107 - tmp106;
	       tmp219 = tmp109 + tmp110;
	       tmp220 = (K923879532 * tmp218) - (K382683432 * tmp219);
	       tmp245 = (K382683432 * tmp218) + (K923879532 * tmp219);
	  }
     }
     {
	  fftw_real tmp25;
	  fftw_real tmp116;
	  fftw_real tmp88;
	  fftw_real tmp114;
	  fftw_real tmp28;
	  fftw_real tmp113;
	  fftw_real tmp91;
	  fftw_real tmp117;
	  fftw_real tmp281;
	  fftw_real tmp282;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp23;
	       fftw_real tmp24;
	       fftw_real tmp86;
	       fftw_real tmp87;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp23 = c_re(input[30 * istride]);
	       tmp24 = c_re(input[14 * istride]);
	       tmp25 = tmp23 + tmp24;
	       tmp116 = tmp23 - tmp24;
	       tmp86 = c_im(input[30 * istride]);
	       tmp87 = c_im(input[14 * istride]);
	       tmp88 = tmp86 + tmp87;
	       tmp114 = tmp86 - tmp87;
	  }
	  {
	       fftw_real tmp26;
	       fftw_real tmp27;
	       fftw_real tmp89;
	       fftw_real tmp90;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp26 = c_re(input[6 * istride]);
	       tmp27 = c_re(input[22 * istride]);
	       tmp28 = tmp26 + tmp27;
	       tmp113 = tmp26 - tmp27;
	       tmp89 = c_im(input[6 * istride]);
	       tmp90 = c_im(input[22 * istride]);
	       tmp91 = tmp89 + tmp90;
	       tmp117 = tmp89 - tmp90;
	  }
	  tmp29 = tmp25 + tmp28;
	  tmp281 = tmp25 - tmp28;
	  tmp282 = tmp88 - tmp91;
	  tmp283 = tmp281 + tmp282;
	  tmp312 = tmp281 - tmp282;
	  tmp92 = tmp88 + tmp91;
	  {
	       fftw_real tmp115;
	       fftw_real tmp118;
	       fftw_real tmp221;
	       fftw_real tmp222;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp115 = tmp113 + tmp114;
	       tmp118 = tmp116 - tmp117;
	       tmp119 = (K382683432 * tmp115) + (K923879532 * tmp118);
	       tmp184 = (K382683432 * tmp118) - (K923879532 * tmp115);
	       tmp221 = tmp114 - tmp113;
	       tmp222 = tmp116 + tmp117;
	       tmp223 = (K923879532 * tmp221) + (K382683432 * tmp222);
	       tmp244 = (K923879532 * tmp222) - (K382683432 * tmp221);
	  }
     }
     {
	  fftw_real tmp34;
	  fftw_real tmp139;
	  fftw_real tmp125;
	  fftw_real tmp286;
	  fftw_real tmp37;
	  fftw_real tmp122;
	  fftw_real tmp142;
	  fftw_real tmp287;
	  fftw_real tmp41;
	  fftw_real tmp132;
	  fftw_real tmp135;
	  fftw_real tmp292;
	  fftw_real tmp44;
	  fftw_real tmp127;
	  fftw_real tmp130;
	  fftw_real tmp293;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp32;
	       fftw_real tmp33;
	       fftw_real tmp123;
	       fftw_real tmp124;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp32 = c_re(input[istride]);
	       tmp33 = c_re(input[17 * istride]);
	       tmp34 = tmp32 + tmp33;
	       tmp139 = tmp32 - tmp33;
	       tmp123 = c_im(input[istride]);
	       tmp124 = c_im(input[17 * istride]);
	       tmp125 = tmp123 - tmp124;
	       tmp286 = tmp123 + tmp124;
	  }
	  {
	       fftw_real tmp35;
	       fftw_real tmp36;
	       fftw_real tmp140;
	       fftw_real tmp141;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp35 = c_re(input[9 * istride]);
	       tmp36 = c_re(input[25 * istride]);
	       tmp37 = tmp35 + tmp36;
	       tmp122 = tmp35 - tmp36;
	       tmp140 = c_im(input[9 * istride]);
	       tmp141 = c_im(input[25 * istride]);
	       tmp142 = tmp140 - tmp141;
	       tmp287 = tmp140 + tmp141;
	  }
	  {
	       fftw_real tmp39;
	       fftw_real tmp40;
	       fftw_real tmp133;
	       fftw_real tmp134;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp39 = c_re(input[5 * istride]);
	       tmp40 = c_re(input[21 * istride]);
	       tmp41 = tmp39 + tmp40;
	       tmp132 = tmp39 - tmp40;
	       tmp133 = c_im(input[5 * istride]);
	       tmp134 = c_im(input[21 * istride]);
	       tmp135 = tmp133 - tmp134;
	       tmp292 = tmp133 + tmp134;
	  }
	  {
	       fftw_real tmp42;
	       fftw_real tmp43;
	       fftw_real tmp128;
	       fftw_real tmp129;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp42 = c_re(input[29 * istride]);
	       tmp43 = c_re(input[13 * istride]);
	       tmp44 = tmp42 + tmp43;
	       tmp127 = tmp42 - tmp43;
	       tmp128 = c_im(input[29 * istride]);
	       tmp129 = c_im(input[13 * istride]);
	       tmp130 = tmp128 - tmp129;
	       tmp293 = tmp128 + tmp129;
	  }
	  {
	       fftw_real tmp291;
	       fftw_real tmp294;
	       fftw_real tmp288;
	       fftw_real tmp289;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp126 = tmp122 + tmp125;
	       tmp229 = tmp125 - tmp122;
	       tmp38 = tmp34 + tmp37;
	       tmp45 = tmp41 + tmp44;
	       tmp346 = tmp38 - tmp45;
	       tmp291 = tmp34 - tmp37;
	       tmp294 = tmp292 - tmp293;
	       tmp295 = tmp291 - tmp294;
	       tmp327 = tmp291 + tmp294;
	       tmp347 = tmp286 + tmp287;
	       tmp348 = tmp292 + tmp293;
	       tmp349 = tmp347 - tmp348;
	       tmp143 = tmp139 - tmp142;
	       tmp226 = tmp139 + tmp142;
	       tmp288 = tmp286 - tmp287;
	       tmp289 = tmp44 - tmp41;
	       tmp290 = tmp288 - tmp289;
	       tmp326 = tmp289 + tmp288;
	       {
		    fftw_real tmp131;
		    fftw_real tmp136;
		    fftw_real tmp144;
		    fftw_real tmp145;
		    ASSERT_ALIGNED_DOUBLE;
		    tmp131 = tmp127 - tmp130;
		    tmp136 = tmp132 + tmp135;
		    tmp137 = K707106781 * (tmp131 - tmp136);
		    tmp227 = K707106781 * (tmp136 + tmp131);
		    tmp144 = tmp135 - tmp132;
		    tmp145 = tmp127 + tmp130;
		    tmp146 = K707106781 * (tmp144 - tmp145);
		    tmp230 = K707106781 * (tmp144 + tmp145);
	       }
	  }
     }
     {
	  fftw_real tmp285;
	  fftw_real tmp317;
	  fftw_real tmp320;
	  fftw_real tmp322;
	  fftw_real tmp308;
	  fftw_real tmp316;
	  fftw_real tmp315;
	  fftw_real tmp321;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp277;
	       fftw_real tmp284;
	       fftw_real tmp318;
	       fftw_real tmp319;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp277 = tmp275 - tmp276;
	       tmp284 = K707106781 * (tmp280 - tmp283);
	       tmp285 = tmp277 + tmp284;
	       tmp317 = tmp277 - tmp284;
	       tmp318 = (K382683432 * tmp290) - (K923879532 * tmp295);
	       tmp319 = (K382683432 * tmp306) + (K923879532 * tmp301);
	       tmp320 = tmp318 - tmp319;
	       tmp322 = tmp318 + tmp319;
	  }
	  {
	       fftw_real tmp296;
	       fftw_real tmp307;
	       fftw_real tmp311;
	       fftw_real tmp314;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp296 = (K923879532 * tmp290) + (K382683432 * tmp295);
	       tmp307 = (K382683432 * tmp301) - (K923879532 * tmp306);
	       tmp308 = tmp296 + tmp307;
	       tmp316 = tmp307 - tmp296;
	       tmp311 = tmp309 - tmp310;
	       tmp314 = K707106781 * (tmp312 - tmp313);
	       tmp315 = tmp311 - tmp314;
	       tmp321 = tmp311 + tmp314;
	  }
	  c_re(output[22 * ostride]) = tmp285 - tmp308;
	  c_re(output[6 * ostride]) = tmp285 + tmp308;
	  c_im(output[30 * ostride]) = tmp315 - tmp316;
	  c_im(output[14 * ostride]) = tmp315 + tmp316;
	  c_re(output[30 * ostride]) = tmp317 - tmp320;
	  c_re(output[14 * ostride]) = tmp317 + tmp320;
	  c_im(output[22 * ostride]) = tmp321 - tmp322;
	  c_im(output[6 * ostride]) = tmp321 + tmp322;
     }
     {
	  fftw_real tmp325;
	  fftw_real tmp337;
	  fftw_real tmp340;
	  fftw_real tmp342;
	  fftw_real tmp332;
	  fftw_real tmp336;
	  fftw_real tmp335;
	  fftw_real tmp341;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp323;
	       fftw_real tmp324;
	       fftw_real tmp338;
	       fftw_real tmp339;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp323 = tmp275 + tmp276;
	       tmp324 = K707106781 * (tmp313 + tmp312);
	       tmp325 = tmp323 + tmp324;
	       tmp337 = tmp323 - tmp324;
	       tmp338 = (K923879532 * tmp326) - (K382683432 * tmp327);
	       tmp339 = (K923879532 * tmp330) + (K382683432 * tmp329);
	       tmp340 = tmp338 - tmp339;
	       tmp342 = tmp338 + tmp339;
	  }
	  {
	       fftw_real tmp328;
	       fftw_real tmp331;
	       fftw_real tmp333;
	       fftw_real tmp334;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp328 = (K382683432 * tmp326) + (K923879532 * tmp327);
	       tmp331 = (K923879532 * tmp329) - (K382683432 * tmp330);
	       tmp332 = tmp328 + tmp331;
	       tmp336 = tmp331 - tmp328;
	       tmp333 = tmp310 + tmp309;
	       tmp334 = K707106781 * (tmp280 + tmp283);
	       tmp335 = tmp333 - tmp334;
	       tmp341 = tmp333 + tmp334;
	  }
	  c_re(output[18 * ostride]) = tmp325 - tmp332;
	  c_re(output[2 * ostride]) = tmp325 + tmp332;
	  c_im(output[26 * ostride]) = tmp335 - tmp336;
	  c_im(output[10 * ostride]) = tmp335 + tmp336;
	  c_re(output[26 * ostride]) = tmp337 - tmp340;
	  c_re(output[10 * ostride]) = tmp337 + tmp340;
	  c_im(output[18 * ostride]) = tmp341 - tmp342;
	  c_im(output[2 * ostride]) = tmp341 + tmp342;
     }
     {