fn_32.c

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

	  fftw_real tmp345;
	  fftw_real tmp361;
	  fftw_real tmp364;
	  fftw_real tmp366;
	  fftw_real tmp356;
	  fftw_real tmp360;
	  fftw_real tmp359;
	  fftw_real tmp365;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp343;
	       fftw_real tmp344;
	       fftw_real tmp362;
	       fftw_real tmp363;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp343 = tmp7 - tmp14;
	       tmp344 = tmp85 - tmp92;
	       tmp345 = tmp343 + tmp344;
	       tmp361 = tmp343 - tmp344;
	       tmp362 = tmp349 - tmp346;
	       tmp363 = tmp351 + tmp354;
	       tmp364 = K707106781 * (tmp362 - tmp363);
	       tmp366 = K707106781 * (tmp362 + tmp363);
	  }
	  {
	       fftw_real tmp350;
	       fftw_real tmp355;
	       fftw_real tmp357;
	       fftw_real tmp358;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp350 = tmp346 + tmp349;
	       tmp355 = tmp351 - tmp354;
	       tmp356 = K707106781 * (tmp350 + tmp355);
	       tmp360 = K707106781 * (tmp355 - tmp350);
	       tmp357 = tmp70 - tmp77;
	       tmp358 = tmp29 - tmp22;
	       tmp359 = tmp357 - tmp358;
	       tmp365 = tmp358 + tmp357;
	  }
	  c_re(output[20 * ostride]) = tmp345 - tmp356;
	  c_re(output[4 * ostride]) = tmp345 + tmp356;
	  c_im(output[28 * ostride]) = tmp359 - tmp360;
	  c_im(output[12 * ostride]) = tmp359 + tmp360;
	  c_re(output[28 * ostride]) = tmp361 - tmp364;
	  c_re(output[12 * ostride]) = tmp361 + tmp364;
	  c_im(output[20 * ostride]) = tmp365 - tmp366;
	  c_im(output[4 * ostride]) = tmp365 + tmp366;
     }
     {
	  fftw_real tmp31;
	  fftw_real tmp367;
	  fftw_real tmp370;
	  fftw_real tmp372;
	  fftw_real tmp62;
	  fftw_real tmp63;
	  fftw_real tmp94;
	  fftw_real tmp371;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp15;
	       fftw_real tmp30;
	       fftw_real tmp368;
	       fftw_real tmp369;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp15 = tmp7 + tmp14;
	       tmp30 = tmp22 + tmp29;
	       tmp31 = tmp15 + tmp30;
	       tmp367 = tmp15 - tmp30;
	       tmp368 = tmp347 + tmp348;
	       tmp369 = tmp352 + tmp353;
	       tmp370 = tmp368 - tmp369;
	       tmp372 = tmp368 + tmp369;
	  }
	  {
	       fftw_real tmp46;
	       fftw_real tmp61;
	       fftw_real tmp78;
	       fftw_real tmp93;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp46 = tmp38 + tmp45;
	       tmp61 = tmp53 + tmp60;
	       tmp62 = tmp46 + tmp61;
	       tmp63 = tmp61 - tmp46;
	       tmp78 = tmp70 + tmp77;
	       tmp93 = tmp85 + tmp92;
	       tmp94 = tmp78 - tmp93;
	       tmp371 = tmp78 + tmp93;
	  }
	  c_re(output[16 * ostride]) = tmp31 - tmp62;
	  c_re(output[0]) = tmp31 + tmp62;
	  c_im(output[8 * ostride]) = tmp63 + tmp94;
	  c_im(output[24 * ostride]) = tmp94 - tmp63;
	  c_re(output[24 * ostride]) = tmp367 - tmp370;
	  c_re(output[8 * ostride]) = tmp367 + tmp370;
	  c_im(output[16 * ostride]) = tmp371 - tmp372;
	  c_im(output[0]) = tmp371 + tmp372;
     }
     {
	  fftw_real tmp121;
	  fftw_real tmp189;
	  fftw_real tmp187;
	  fftw_real tmp193;
	  fftw_real tmp148;
	  fftw_real tmp190;
	  fftw_real tmp175;
	  fftw_real tmp191;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp105;
	       fftw_real tmp120;
	       fftw_real tmp183;
	       fftw_real tmp186;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp105 = tmp97 - tmp104;
	       tmp120 = tmp112 - tmp119;
	       tmp121 = tmp105 + tmp120;
	       tmp189 = tmp105 - tmp120;
	       tmp183 = tmp179 - tmp182;
	       tmp186 = tmp184 - tmp185;
	       tmp187 = tmp183 - tmp186;
	       tmp193 = tmp183 + tmp186;
	  }
	  {
	       fftw_real tmp138;
	       fftw_real tmp147;
	       fftw_real tmp165;
	       fftw_real tmp174;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp138 = tmp126 - tmp137;
	       tmp147 = tmp143 - tmp146;
	       tmp148 = (K980785280 * tmp138) + (K195090322 * tmp147);
	       tmp190 = (K195090322 * tmp138) - (K980785280 * tmp147);
	       tmp165 = tmp153 - tmp164;
	       tmp174 = tmp170 - tmp173;
	       tmp175 = (K195090322 * tmp165) - (K980785280 * tmp174);
	       tmp191 = (K195090322 * tmp174) + (K980785280 * tmp165);
	  }
	  {
	       fftw_real tmp176;
	       fftw_real tmp188;
	       fftw_real tmp192;
	       fftw_real tmp194;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp176 = tmp148 + tmp175;
	       c_re(output[23 * ostride]) = tmp121 - tmp176;
	       c_re(output[7 * ostride]) = tmp121 + tmp176;
	       tmp188 = tmp175 - tmp148;
	       c_im(output[31 * ostride]) = tmp187 - tmp188;
	       c_im(output[15 * ostride]) = tmp187 + tmp188;
	       tmp192 = tmp190 - tmp191;
	       c_re(output[31 * ostride]) = tmp189 - tmp192;
	       c_re(output[15 * ostride]) = tmp189 + tmp192;
	       tmp194 = tmp190 + tmp191;
	       c_im(output[23 * ostride]) = tmp193 - tmp194;
	       c_im(output[7 * ostride]) = tmp193 + tmp194;
	  }
     }
     {
	  fftw_real tmp197;
	  fftw_real tmp209;
	  fftw_real tmp207;
	  fftw_real tmp213;
	  fftw_real tmp200;
	  fftw_real tmp210;
	  fftw_real tmp203;
	  fftw_real tmp211;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp195;
	       fftw_real tmp196;
	       fftw_real tmp205;
	       fftw_real tmp206;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp195 = tmp97 + tmp104;
	       tmp196 = tmp185 + tmp184;
	       tmp197 = tmp195 + tmp196;
	       tmp209 = tmp195 - tmp196;
	       tmp205 = tmp179 + tmp182;
	       tmp206 = tmp112 + tmp119;
	       tmp207 = tmp205 - tmp206;
	       tmp213 = tmp205 + tmp206;
	  }
	  {
	       fftw_real tmp198;
	       fftw_real tmp199;
	       fftw_real tmp201;
	       fftw_real tmp202;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp198 = tmp126 + tmp137;
	       tmp199 = tmp143 + tmp146;
	       tmp200 = (K555570233 * tmp198) + (K831469612 * tmp199);
	       tmp210 = (K831469612 * tmp198) - (K555570233 * tmp199);
	       tmp201 = tmp153 + tmp164;
	       tmp202 = tmp170 + tmp173;
	       tmp203 = (K831469612 * tmp201) - (K555570233 * tmp202);
	       tmp211 = (K831469612 * tmp202) + (K555570233 * tmp201);
	  }
	  {
	       fftw_real tmp204;
	       fftw_real tmp208;
	       fftw_real tmp212;
	       fftw_real tmp214;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp204 = tmp200 + tmp203;
	       c_re(output[19 * ostride]) = tmp197 - tmp204;
	       c_re(output[3 * ostride]) = tmp197 + tmp204;
	       tmp208 = tmp203 - tmp200;
	       c_im(output[27 * ostride]) = tmp207 - tmp208;
	       c_im(output[11 * ostride]) = tmp207 + tmp208;
	       tmp212 = tmp210 - tmp211;
	       c_re(output[27 * ostride]) = tmp209 - tmp212;
	       c_re(output[11 * ostride]) = tmp209 + tmp212;
	       tmp214 = tmp210 + tmp211;
	       c_im(output[19 * ostride]) = tmp213 - tmp214;
	       c_im(output[3 * ostride]) = tmp213 + tmp214;
	  }
     }
     {
	  fftw_real tmp225;
	  fftw_real tmp249;
	  fftw_real tmp247;
	  fftw_real tmp253;
	  fftw_real tmp232;
	  fftw_real tmp250;
	  fftw_real tmp239;
	  fftw_real tmp251;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp217;
	       fftw_real tmp224;
	       fftw_real tmp243;
	       fftw_real tmp246;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp217 = tmp215 - tmp216;
	       tmp224 = tmp220 - tmp223;
	       tmp225 = tmp217 + tmp224;
	       tmp249 = tmp217 - tmp224;
	       tmp243 = tmp241 - tmp242;
	       tmp246 = tmp244 - tmp245;
	       tmp247 = tmp243 - tmp246;
	       tmp253 = tmp243 + tmp246;
	  }
	  {
	       fftw_real tmp228;
	       fftw_real tmp231;
	       fftw_real tmp235;
	       fftw_real tmp238;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp228 = tmp226 - tmp227;
	       tmp231 = tmp229 - tmp230;
	       tmp232 = (K555570233 * tmp228) + (K831469612 * tmp231);
	       tmp250 = (K555570233 * tmp231) - (K831469612 * tmp228);
	       tmp235 = tmp233 - tmp234;
	       tmp238 = tmp236 - tmp237;
	       tmp239 = (K555570233 * tmp235) - (K831469612 * tmp238);
	       tmp251 = (K831469612 * tmp235) + (K555570233 * tmp238);
	  }
	  {
	       fftw_real tmp240;
	       fftw_real tmp248;
	       fftw_real tmp252;
	       fftw_real tmp254;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp240 = tmp232 + tmp239;
	       c_re(output[21 * ostride]) = tmp225 - tmp240;
	       c_re(output[5 * ostride]) = tmp225 + tmp240;
	       tmp248 = tmp239 - tmp232;
	       c_im(output[29 * ostride]) = tmp247 - tmp248;
	       c_im(output[13 * ostride]) = tmp247 + tmp248;
	       tmp252 = tmp250 - tmp251;
	       c_re(output[29 * ostride]) = tmp249 - tmp252;
	       c_re(output[13 * ostride]) = tmp249 + tmp252;
	       tmp254 = tmp250 + tmp251;
	       c_im(output[21 * ostride]) = tmp253 - tmp254;
	       c_im(output[5 * ostride]) = tmp253 + tmp254;
	  }
     }
     {
	  fftw_real tmp257;
	  fftw_real tmp269;
	  fftw_real tmp267;
	  fftw_real tmp273;
	  fftw_real tmp260;
	  fftw_real tmp270;
	  fftw_real tmp263;
	  fftw_real tmp271;
	  ASSERT_ALIGNED_DOUBLE;
	  {
	       fftw_real tmp255;
	       fftw_real tmp256;
	       fftw_real tmp265;
	       fftw_real tmp266;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp255 = tmp215 + tmp216;
	       tmp256 = tmp245 + tmp244;
	       tmp257 = tmp255 + tmp256;
	       tmp269 = tmp255 - tmp256;
	       tmp265 = tmp241 + tmp242;
	       tmp266 = tmp220 + tmp223;
	       tmp267 = tmp265 - tmp266;
	       tmp273 = tmp265 + tmp266;
	  }
	  {
	       fftw_real tmp258;
	       fftw_real tmp259;
	       fftw_real tmp261;
	       fftw_real tmp262;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp258 = tmp226 + tmp227;
	       tmp259 = tmp229 + tmp230;
	       tmp260 = (K980785280 * tmp258) + (K195090322 * tmp259);
	       tmp270 = (K980785280 * tmp259) - (K195090322 * tmp258);
	       tmp261 = tmp233 + tmp234;
	       tmp262 = tmp236 + tmp237;
	       tmp263 = (K980785280 * tmp261) - (K195090322 * tmp262);
	       tmp271 = (K195090322 * tmp261) + (K980785280 * tmp262);
	  }
	  {
	       fftw_real tmp264;
	       fftw_real tmp268;
	       fftw_real tmp272;
	       fftw_real tmp274;
	       ASSERT_ALIGNED_DOUBLE;
	       tmp264 = tmp260 + tmp263;
	       c_re(output[17 * ostride]) = tmp257 - tmp264;
	       c_re(output[ostride]) = tmp257 + tmp264;
	       tmp268 = tmp263 - tmp260;
	       c_im(output[25 * ostride]) = tmp267 - tmp268;
	       c_im(output[9 * ostride]) = tmp267 + tmp268;
	       tmp272 = tmp270 - tmp271;
	       c_re(output[25 * ostride]) = tmp269 - tmp272;
	       c_re(output[9 * ostride]) = tmp269 + tmp272;
	       tmp274 = tmp270 + tmp271;
	       c_im(output[17 * ostride]) = tmp273 - tmp274;
	       c_im(output[ostride]) = tmp273 + tmp274;
	  }
     }
}

fftw_codelet_desc fftw_no_twiddle_32_desc = {
     "fftw_no_twiddle_32",
     (void (*)()) fftw_no_twiddle_32,
     32,
     FFTW_FORWARD,
     FFTW_NOTW,
     705,
     0,
     (const int *) 0,
};