src/fftw/frc_64.c

   1 /*
   2  * Copyright (c) 1997-1999 Massachusetts Institute of Technology
   3  *
   4  * This program is free software; you can redistribute it and/or modify
   5  * it under the terms of the GNU General Public License as published by
   6  * the Free Software Foundation; either version 2 of the License, or
   7  * (at your option) any later version.
   8  *
   9  * This program is distributed in the hope that it will be useful,
  10  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  11  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  12  * GNU General Public License for more details.
  13  *
  14  * You should have received a copy of the GNU General Public License
  15  * along with this program; if not, write to the Free Software
  16  * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
  17  *
  18  */
  19
  20 /* This file was automatically generated --- DO NOT EDIT */
  21 /* Generated on Tue May 18 13:54:40 EDT 1999 */
  22
  23 #include <fftw-int.h>
  24 #include <fftw.h>
  25
  26 /* Generated by: ./genfft -magic-alignment-check -magic-twiddle-load-all -magic-variables 4 -magic-loopi -real2hc 64 */
  27
  28 /*
  29  * This function contains 394 FP additions, 124 FP multiplications,
  30  * (or, 342 additions, 72 multiplications, 52 fused multiply/add),
  31  * 90 stack variables, and 128 memory accesses
  32  */
  33 static const fftw_real K471396736 = FFTW_KONST(+0.471396736825997648556387625905254377657460319);
  34 static const fftw_real K881921264 = FFTW_KONST(+0.881921264348355029712756863660388349508442621);
  35 static const fftw_real K290284677 = FFTW_KONST(+0.290284677254462367636192375817395274691476278);
  36 static const fftw_real K956940335 = FFTW_KONST(+0.956940335732208864935797886980269969482849206);
  37 static const fftw_real K555570233 = FFTW_KONST(+0.555570233019602224742830813948532874374937191);
  38 static const fftw_real K831469612 = FFTW_KONST(+0.831469612302545237078788377617905756738560812);
  39 static const fftw_real K995184726 = FFTW_KONST(+0.995184726672196886244836953109479921575474869);
  40 static const fftw_real K098017140 = FFTW_KONST(+0.098017140329560601994195563888641845861136673);
  41 static const fftw_real K773010453 = FFTW_KONST(+0.773010453362736960810906609758469800971041293);
  42 static const fftw_real K634393284 = FFTW_KONST(+0.634393284163645498215171613225493370675687095);
  43 static const fftw_real K980785280 = FFTW_KONST(+0.980785280403230449126182236134239036973933731);
  44 static const fftw_real K195090322 = FFTW_KONST(+0.195090322016128267848284868477022240927691618);
  45 static const fftw_real K382683432 = FFTW_KONST(+0.382683432365089771728459984030398866761344562);
  46 static const fftw_real K923879532 = FFTW_KONST(+0.923879532511286756128183189396788286822416626);
  47 static const fftw_real K707106781 = FFTW_KONST(+0.707106781186547524400844362104849039284835938);
  48
  49 /*
  50  * Generator Id's :
  51  * $Id$
  52  * $Id$
  53  * $Id$
  54  */
  55
  56 void fftw_real2hc_64(const fftw_real *input, fftw_real *real_output, fftw_real *imag_output, int istride, int real_ostride, int imag_ostride)
  57 {
  58      fftw_real tmp269;
  59      fftw_real tmp7;
  60      fftw_real tmp14;
  61      fftw_real tmp15;
  62      fftw_real tmp320;
  63      fftw_real tmp67;
  64      fftw_real tmp209;
  65      fftw_real tmp179;
  66      fftw_real tmp238;
  67      fftw_real tmp301;
  68      fftw_real tmp361;
  69      fftw_real tmp381;
  70      fftw_real tmp312;
  71      fftw_real tmp143;
  72      fftw_real tmp227;
  73      fftw_real tmp167;
  74      fftw_real tmp230;
  75      fftw_real tmp158;
  76      fftw_real tmp231;
  77      fftw_real tmp162;
  78      fftw_real tmp228;
  79      fftw_real tmp308;
  80      fftw_real tmp313;
  81      fftw_real tmp364;
  82      fftw_real tmp382;
  83      fftw_real tmp22;
  84      fftw_real tmp29;
  85      fftw_real tmp30;
  86      fftw_real tmp272;
  87      fftw_real tmp321;
  88      fftw_real tmp74;
  89      fftw_real tmp239;
  90      fftw_real tmp176;
  91      fftw_real tmp210;
  92      fftw_real tmp284;
  93      fftw_real tmp354;
  94      fftw_real tmp378;
  95      fftw_real tmp295;
  96      fftw_real tmp121;
  97      fftw_real tmp223;
  98      fftw_real tmp127;
  99      fftw_real tmp220;
 100      fftw_real tmp110;
 101      fftw_real tmp221;
 102      fftw_real tmp130;
 103      fftw_real tmp224;
 104      fftw_real tmp291;
 105      fftw_real tmp296;
 106      fftw_real tmp357;
 107      fftw_real tmp379;
 108      fftw_real tmp89;
 109      fftw_real tmp213;
 110      fftw_real tmp92;
 111      fftw_real tmp212;
 112      fftw_real tmp46;
 113      fftw_real tmp374;
 114      fftw_real tmp276;
 115      fftw_real tmp317;
 116      fftw_real tmp80;
 117      fftw_real tmp215;
 118      fftw_real tmp83;
 119      fftw_real tmp216;
 120      fftw_real tmp61;
 121      fftw_real tmp375;
 122      fftw_real tmp279;
 123      fftw_real tmp318;
 124      ASSERT_ALIGNED_DOUBLE();
 125      {
 126           fftw_real tmp3;
 127           fftw_real tmp63;
 128           fftw_real tmp13;
 129           fftw_real tmp65;
 130           fftw_real tmp6;
 131           fftw_real tmp178;
 132           fftw_real tmp10;
 133           fftw_real tmp64;
 134           fftw_real tmp66;
 135           fftw_real tmp177;
 136           ASSERT_ALIGNED_DOUBLE();
 137           {
 138                fftw_real tmp1;
 139                fftw_real tmp2;
 140                fftw_real tmp11;
 141                fftw_real tmp12;
 142                ASSERT_ALIGNED_DOUBLE();
 143                tmp1 = input[0];
 144                tmp2 = input[32 * istride];
 145                tmp3 = tmp1 + tmp2;
 146                tmp63 = tmp1 - tmp2;
 147                tmp11 = input[56 * istride];
 148                tmp12 = input[24 * istride];
 149                tmp13 = tmp11 + tmp12;
 150                tmp65 = tmp11 - tmp12;
 151           }
 152           {
 153                fftw_real tmp4;
 154                fftw_real tmp5;
 155                fftw_real tmp8;
 156                fftw_real tmp9;
 157                ASSERT_ALIGNED_DOUBLE();
 158                tmp4 = input[16 * istride];
 159                tmp5 = input[48 * istride];
 160                tmp6 = tmp4 + tmp5;
 161                tmp178 = tmp4 - tmp5;
 162                tmp8 = input[8 * istride];
 163                tmp9 = input[40 * istride];
 164                tmp10 = tmp8 + tmp9;
 165                tmp64 = tmp8 - tmp9;
 166           }
 167           tmp269 = tmp3 - tmp6;
 168           tmp7 = tmp3 + tmp6;
 169           tmp14 = tmp10 + tmp13;
 170           tmp15 = tmp7 + tmp14;
 171           tmp320 = tmp13 - tmp10;
 172           tmp66 = K707106781 * (tmp64 + tmp65);
 173           tmp67 = tmp63 + tmp66;
 174           tmp209 = tmp63 - tmp66;
 175           tmp177 = K707106781 * (tmp65 - tmp64);
 176           tmp179 = tmp177 - tmp178;
 177           tmp238 = tmp178 + tmp177;
 178      }
 179      {
 180           fftw_real tmp135;
 181           fftw_real tmp299;
 182           fftw_real tmp166;
 183           fftw_real tmp300;
 184           fftw_real tmp138;
 185           fftw_real tmp311;
 186           fftw_real tmp141;
 187           fftw_real tmp310;
 188           ASSERT_ALIGNED_DOUBLE();
 189           {
 190                fftw_real tmp133;
 191                fftw_real tmp134;
 192                fftw_real tmp164;
 193                fftw_real tmp165;
 194                ASSERT_ALIGNED_DOUBLE();
 195                tmp133 = input[63 * istride];
 196                tmp134 = input[31 * istride];
 197                tmp135 = tmp133 - tmp134;
 198                tmp299 = tmp133 + tmp134;
 199                tmp164 = input[15 * istride];
 200                tmp165 = input[47 * istride];
 201                tmp166 = tmp164 - tmp165;
 202                tmp300 = tmp164 + tmp165;
 203           }
 204           {
 205                fftw_real tmp136;
 206                fftw_real tmp137;
 207                fftw_real tmp139;
 208                fftw_real tmp140;
 209                ASSERT_ALIGNED_DOUBLE();
 210                tmp136 = input[7 * istride];
 211                tmp137 = input[39 * istride];
 212                tmp138 = tmp136 - tmp137;
 213                tmp311 = tmp136 + tmp137;
 214                tmp139 = input[55 * istride];
 215                tmp140 = input[23 * istride];
 216                tmp141 = tmp139 - tmp140;
 217                tmp310 = tmp139 + tmp140;
 218           }
 219           {
 220                fftw_real tmp359;
 221                fftw_real tmp360;
 222                fftw_real tmp142;
 223                fftw_real tmp163;
 224                ASSERT_ALIGNED_DOUBLE();
 225                tmp301 = tmp299 - tmp300;
 226                tmp359 = tmp299 + tmp300;
 227                tmp360 = tmp311 + tmp310;
 228                tmp361 = tmp359 + tmp360;
 229                tmp381 = tmp359 - tmp360;
 230                tmp312 = tmp310 - tmp311;
 231                tmp142 = K707106781 * (tmp138 + tmp141);
 232                tmp143 = tmp135 + tmp142;
 233                tmp227 = tmp135 - tmp142;
 234                tmp163 = K707106781 * (tmp141 - tmp138);
 235                tmp167 = tmp163 - tmp166;
 236                tmp230 = tmp166 + tmp163;
 237           }
 238      }
 239      {
 240           fftw_real tmp146;
 241           fftw_real tmp302;
 242           fftw_real tmp156;
 243           fftw_real tmp306;
 244           fftw_real tmp149;
 245           fftw_real tmp303;
 246           fftw_real tmp153;
 247           fftw_real tmp305;
 248           ASSERT_ALIGNED_DOUBLE();
 249           {
 250                fftw_real tmp144;
 251                fftw_real tmp145;
 252                fftw_real tmp154;
 253                fftw_real tmp155;
 254                ASSERT_ALIGNED_DOUBLE();
 255                tmp144 = input[3 * istride];
 256                tmp145 = input[35 * istride];
 257                tmp146 = tmp144 - tmp145;
 258                tmp302 = tmp144 + tmp145;
 259                tmp154 = input[11 * istride];
 260                tmp155 = input[43 * istride];
 261                tmp156 = tmp154 - tmp155;
 262                tmp306 = tmp154 + tmp155;
 263           }
 264           {
 265                fftw_real tmp147;
 266                fftw_real tmp148;
 267                fftw_real tmp151;
 268                fftw_real tmp152;
 269                ASSERT_ALIGNED_DOUBLE();
 270                tmp147 = input[19 * istride];
 271                tmp148 = input[51 * istride];
 272                tmp149 = tmp147 - tmp148;
 273                tmp303 = tmp147 + tmp148;
 274                tmp151 = input[59 * istride];
 275                tmp152 = input[27 * istride];
 276                tmp153 = tmp151 - tmp152;
 277                tmp305 = tmp151 + tmp152;
 278           }
 279           {
 280                fftw_real tmp150;
 281                fftw_real tmp157;
 282                fftw_real tmp160;
 283                fftw_real tmp161;
 284                ASSERT_ALIGNED_DOUBLE();
 285                tmp150 = (K923879532 * tmp146) - (K382683432 * tmp149);
 286                tmp157 = (K923879532 * tmp153) + (K382683432 * tmp156);
 287                tmp158 = tmp150 + tmp157;
 288                tmp231 = tmp157 - tmp150;
 289                tmp160 = (K382683432 * tmp153) - (K923879532 * tmp156);
 290                tmp161 = (K382683432 * tmp146) + (K923879532 * tmp149);
 291                tmp162 = tmp160 - tmp161;
 292                tmp228 = tmp161 + tmp160;
 293           }
 294           {
 295                fftw_real tmp304;
 296                fftw_real tmp307;
 297                fftw_real tmp362;
 298                fftw_real tmp363;
 299                ASSERT_ALIGNED_DOUBLE();
 300                tmp304 = tmp302 - tmp303;
 301                tmp307 = tmp305 - tmp306;
 302                tmp308 = K707106781 * (tmp304 + tmp307);
 303                tmp313 = K707106781 * (tmp307 - tmp304);
 304                tmp362 = tmp302 + tmp303;
 305                tmp363 = tmp305 + tmp306;
 306                tmp364 = tmp362 + tmp363;
 307                tmp382 = tmp363 - tmp362;
 308           }
 309      }
 310      {
 311           fftw_real tmp18;
 312           fftw_real tmp68;
 313           fftw_real tmp28;
 314           fftw_real tmp72;
 315           fftw_real tmp21;
 316           fftw_real tmp69;
 317           fftw_real tmp25;
 318           fftw_real tmp71;
 319           fftw_real tmp270;
 320           fftw_real tmp271;
 321           ASSERT_ALIGNED_DOUBLE();
 322           {
 323                fftw_real tmp16;
 324                fftw_real tmp17;
 325                fftw_real tmp26;
 326                fftw_real tmp27;
 327                ASSERT_ALIGNED_DOUBLE();
 328                tmp16 = input[4 * istride];
 329                tmp17 = input[36 * istride];
 330                tmp18 = tmp16 + tmp17;
 331                tmp68 = tmp16 - tmp17;
 332                tmp26 = input[12 * istride];
 333                tmp27 = input[44 * istride];
 334                tmp28 = tmp26 + tmp27;
 335                tmp72 = tmp26 - tmp27;
 336           }
 337           {
 338                fftw_real tmp19;
 339                fftw_real tmp20;
 340                fftw_real tmp23;
 341                fftw_real tmp24;
 342                ASSERT_ALIGNED_DOUBLE();
 343                tmp19 = input[20 * istride];
 344                tmp20 = input[52 * istride];
 345                tmp21 = tmp19 + tmp20;
 346                tmp69 = tmp19 - tmp20;
 347                tmp23 = input[60 * istride];
 348                tmp24 = input[28 * istride];
 349                tmp25 = tmp23 + tmp24;
 350                tmp71 = tmp23 - tmp24;
 351           }
 352           tmp22 = tmp18 + tmp21;
 353           tmp29 = tmp25 + tmp28;
 354           tmp30 = tmp22 + tmp29;
 355           tmp270 = tmp18 - tmp21;
 356           tmp271 = tmp25 - tmp28;
 357           tmp272 = K707106781 * (tmp270 + tmp271);
 358           tmp321 = K707106781 * (tmp271 - tmp270);
 359           {
 360                fftw_real tmp70;
 361                fftw_real tmp73;
 362                fftw_real tmp174;
 363                fftw_real tmp175;
 364                ASSERT_ALIGNED_DOUBLE();
 365                tmp70 = (K923879532 * tmp68) - (K382683432 * tmp69);
 366                tmp73 = (K923879532 * tmp71) + (K382683432 * tmp72);
 367                tmp74 = tmp70 + tmp73;
 368                tmp239 = tmp73 - tmp70;
 369                tmp174 = (K382683432 * tmp71) - (K923879532 * tmp72);
 370                tmp175 = (K382683432 * tmp68) + (K923879532 * tmp69);
 371                tmp176 = tmp174 - tmp175;
 372                tmp210 = tmp175 + tmp174;
 373           }
 374      }
 375      {
 376           fftw_real tmp113;
 377           fftw_real tmp293;
 378           fftw_real tmp125;
 379           fftw_real tmp282;
 380           fftw_real tmp116;
 381           fftw_real tmp294;
 382           fftw_real tmp120;
 383           fftw_real tmp283;
 384           ASSERT_ALIGNED_DOUBLE();
 385           {
 386                fftw_real tmp111;
 387                fftw_real tmp112;
 388                fftw_real tmp123;
 389                fftw_real tmp124;
 390                ASSERT_ALIGNED_DOUBLE();
 391                tmp111 = input[57 * istride];
 392                tmp112 = input[25 * istride];
 393                tmp113 = tmp111 - tmp112;
 394                tmp293 = tmp111 + tmp112;
 395                tmp123 = input[istride];
 396                tmp124 = input[33 * istride];
 397                tmp125 = tmp123 - tmp124;
 398                tmp282 = tmp123 + tmp124;
 399           }
 400           {
 401                fftw_real tmp114;
 402                fftw_real tmp115;
 403                fftw_real tmp118;
 404                fftw_real tmp119;
 405                ASSERT_ALIGNED_DOUBLE();
 406                tmp114 = input[9 * istride];
 407                tmp115 = input[41 * istride];
 408                tmp116 = tmp114 - tmp115;
 409                tmp294 = tmp114 + tmp115;
 410                tmp118 = input[17 * istride];
 411                tmp119 = input[49 * istride];
 412                tmp120 = tmp118 - tmp119;
 413                tmp283 = tmp118 + tmp119;
 414           }
 415           {
 416                fftw_real tmp352;
 417                fftw_real tmp353;
 418                fftw_real tmp117;
 419                fftw_real tmp126;
 420                ASSERT_ALIGNED_DOUBLE();
 421                tmp284 = tmp282 - tmp283;
 422                tmp352 = tmp282 + tmp283;
 423                tmp353 = tmp294 + tmp293;
 424                tmp354 = tmp352 + tmp353;
 425                tmp378 = tmp352 - tmp353;
 426                tmp295 = tmp293 - tmp294;
 427                tmp117 = K707106781 * (tmp113 - tmp116);
 428                tmp121 = tmp117 - tmp120;
 429                tmp223 = tmp120 + tmp117;
 430                tmp126 = K707106781 * (tmp116 + tmp113);
 431                tmp127 = tmp125 + tmp126;
 432                tmp220 = tmp125 - tmp126;
 433           }
 434      }
 435      {
 436           fftw_real tmp98;
 437           fftw_real tmp288;
 438           fftw_real tmp108;
 439           fftw_real tmp286;
 440           fftw_real tmp101;
 441           fftw_real tmp289;
 442           fftw_real tmp105;
 443           fftw_real tmp285;
 444           ASSERT_ALIGNED_DOUBLE();
 445           {
 446                fftw_real tmp96;
 447                fftw_real tmp97;
 448                fftw_real tmp106;
 449                fftw_real tmp107;
 450                ASSERT_ALIGNED_DOUBLE();
 451                tmp96 = input[61 * istride];
 452                tmp97 = input[29 * istride];
 453                tmp98 = tmp96 - tmp97;
 454                tmp288 = tmp96 + tmp97;
 455                tmp106 = input[21 * istride];
 456                tmp107 = input[53 * istride];
 457                tmp108 = tmp106 - tmp107;
 458                tmp286 = tmp106 + tmp107;
 459           }
 460           {
 461                fftw_real tmp99;
 462                fftw_real tmp100;
 463                fftw_real tmp103;
 464                fftw_real tmp104;
 465                ASSERT_ALIGNED_DOUBLE();
 466                tmp99 = input[13 * istride];
 467                tmp100 = input[45 * istride];
 468                tmp101 = tmp99 - tmp100;
 469                tmp289 = tmp99 + tmp100;
 470                tmp103 = input[5 * istride];
 471                tmp104 = input[37 * istride];
 472                tmp105 = tmp103 - tmp104;
 473                tmp285 = tmp103 + tmp104;
 474           }
 475           {
 476                fftw_real tmp102;
 477                fftw_real tmp109;
 478                fftw_real tmp128;
 479                fftw_real tmp129;
 480                ASSERT_ALIGNED_DOUBLE();
 481                tmp102 = (K382683432 * tmp98) - (K923879532 * tmp101);
 482                tmp109 = (K382683432 * tmp105) + (K923879532 * tmp108);
 483                tmp110 = tmp102 - tmp109;
 484                tmp221 = tmp109 + tmp102;
 485                tmp128 = (K923879532 * tmp105) - (K382683432 * tmp108);
 486                tmp129 = (K923879532 * tmp98) + (K382683432 * tmp101);
 487                tmp130 = tmp128 + tmp129;
 488                tmp224 = tmp129 - tmp128;
 489           }
 490           {
 491                fftw_real tmp287;
 492                fftw_real tmp290;
 493                fftw_real tmp355;
 494                fftw_real tmp356;
 495                ASSERT_ALIGNED_DOUBLE();
 496                tmp287 = tmp285 - tmp286;
 497                tmp290 = tmp288 - tmp289;
 498                tmp291 = K707106781 * (tmp287 + tmp290);
 499                tmp296 = K707106781 * (tmp290 - tmp287);
 500                tmp355 = tmp285 + tmp286;
 501                tmp356 = tmp288 + tmp289;
 502                tmp357 = tmp355 + tmp356;
 503                tmp379 = tmp356 - tmp355;
 504           }
 505      }
 506      {
 507           fftw_real tmp34;
 508           fftw_real tmp90;
 509           fftw_real tmp44;
 510           fftw_real tmp85;
 511           fftw_real tmp37;
 512           fftw_real tmp88;
 513           fftw_real tmp41;
 514           fftw_real tmp86;
 515           fftw_real tmp87;
 516           fftw_real tmp91;
 517           ASSERT_ALIGNED_DOUBLE();
 518           {
 519                fftw_real tmp32;
 520                fftw_real tmp33;
 521                fftw_real tmp42;
 522                fftw_real tmp43;
 523                ASSERT_ALIGNED_DOUBLE();
 524                tmp32 = input[2 * istride];
 525                tmp33 = input[34 * istride];
 526                tmp34 = tmp32 + tmp33;
 527                tmp90 = tmp32 - tmp33;
 528                tmp42 = input[58 * istride];
 529                tmp43 = input[26 * istride];
 530                tmp44 = tmp42 + tmp43;
 531                tmp85 = tmp42 - tmp43;
 532           }
 533           {
 534                fftw_real tmp35;
 535                fftw_real tmp36;
 536                fftw_real tmp39;
 537                fftw_real tmp40;
 538                ASSERT_ALIGNED_DOUBLE();
 539                tmp35 = input[18 * istride];
 540                tmp36 = input[50 * istride];
 541                tmp37 = tmp35 + tmp36;
 542                tmp88 = tmp35 - tmp36;
 543                tmp39 = input[10 * istride];
 544                tmp40 = input[42 * istride];
 545                tmp41 = tmp39 + tmp40;
 546                tmp86 = tmp39 - tmp40;
 547           }
 548           tmp87 = K707106781 * (tmp85 - tmp86);
 549           tmp89 = tmp87 - tmp88;
 550           tmp213 = tmp88 + tmp87;
 551           tmp91 = K707106781 * (tmp86 + tmp85);
 552           tmp92 = tmp90 + tmp91;
 553           tmp212 = tmp90 - tmp91;
 554           {
 555                fftw_real tmp38;
 556                fftw_real tmp45;
 557                fftw_real tmp274;
 558                fftw_real tmp275;
 559                ASSERT_ALIGNED_DOUBLE();
 560                tmp38 = tmp34 + tmp37;
 561                tmp45 = tmp41 + tmp44;
 562                tmp46 = tmp38 + tmp45;
 563                tmp374 = tmp38 - tmp45;
 564                tmp274 = tmp34 - tmp37;
 565                tmp275 = tmp44 - tmp41;
 566                tmp276 = (K923879532 * tmp274) + (K382683432 * tmp275);
 567                tmp317 = (K923879532 * tmp275) - (K382683432 * tmp274);
 568           }
 569      }
 570      {
 571           fftw_real tmp49;
 572           fftw_real tmp76;
 573           fftw_real tmp59;
 574           fftw_real tmp78;
 575           fftw_real tmp52;
 576           fftw_real tmp82;
 577           fftw_real tmp56;
 578           fftw_real tmp77;
 579           fftw_real tmp79;
 580           fftw_real tmp81;
 581           ASSERT_ALIGNED_DOUBLE();
 582           {
 583                fftw_real tmp47;
 584                fftw_real tmp48;
 585                fftw_real tmp57;
 586                fftw_real tmp58;
 587                ASSERT_ALIGNED_DOUBLE();
 588                tmp47 = input[62 * istride];
 589                tmp48 = input[30 * istride];
 590                tmp49 = tmp47 + tmp48;
 591                tmp76 = tmp47 - tmp48;
 592                tmp57 = input[54 * istride];
 593                tmp58 = input[22 * istride];
 594                tmp59 = tmp57 + tmp58;
 595                tmp78 = tmp57 - tmp58;
 596           }
 597           {
 598                fftw_real tmp50;
 599                fftw_real tmp51;
 600                fftw_real tmp54;
 601                fftw_real tmp55;
 602                ASSERT_ALIGNED_DOUBLE();
 603                tmp50 = input[14 * istride];
 604                tmp51 = input[46 * istride];
 605                tmp52 = tmp50 + tmp51;
 606                tmp82 = tmp50 - tmp51;
 607                tmp54 = input[6 * istride];
 608                tmp55 = input[38 * istride];
 609                tmp56 = tmp54 + tmp55;
 610                tmp77 = tmp54 - tmp55;
 611           }
 612           tmp79 = K707106781 * (tmp77 + tmp78);
 613           tmp80 = tmp76 + tmp79;
 614           tmp215 = tmp76 - tmp79;
 615           tmp81 = K707106781 * (tmp78 - tmp77);
 616           tmp83 = tmp81 - tmp82;
 617           tmp216 = tmp82 + tmp81;
 618           {
 619                fftw_real tmp53;
 620                fftw_real tmp60;
 621                fftw_real tmp277;
 622                fftw_real tmp278;
 623                ASSERT_ALIGNED_DOUBLE();
 624                tmp53 = tmp49 + tmp52;
 625                tmp60 = tmp56 + tmp59;
 626                tmp61 = tmp53 + tmp60;
 627                tmp375 = tmp53 - tmp60;
 628                tmp277 = tmp49 - tmp52;
 629                tmp278 = tmp59 - tmp56;
 630                tmp279 = (K923879532 * tmp277) - (K382683432 * tmp278);
 631                tmp318 = (K382683432 * tmp277) + (K923879532 * tmp278);
 632           }
 633      }
 634      {
 635           fftw_real tmp31;
 636           fftw_real tmp62;
 637           fftw_real tmp369;
 638           fftw_real tmp370;
 639           fftw_real tmp371;
 640           fftw_real tmp372;
 641           ASSERT_ALIGNED_DOUBLE();
 642           tmp31 = tmp15 + tmp30;
 643           tmp62 = tmp46 + tmp61;
 644           tmp369 = tmp31 + tmp62;
 645           tmp370 = tmp354 + tmp357;
 646           tmp371 = tmp361 + tmp364;
 647           tmp372 = tmp370 + tmp371;
 648           real_output[16 * real_ostride] = tmp31 - tmp62;
 649           imag_output[16 * imag_ostride] = tmp371 - tmp370;
 650           real_output[32 * real_ostride] = tmp369 - tmp372;
 651           real_output[0] = tmp369 + tmp372;
 652      }
 653      {
 654           fftw_real tmp351;
 655           fftw_real tmp367;
 656           fftw_real tmp366;
 657           fftw_real tmp368;
 658           fftw_real tmp358;
 659           fftw_real tmp365;
 660           ASSERT_ALIGNED_DOUBLE();
 661           tmp351 = tmp15 - tmp30;
 662           tmp367 = tmp61 - tmp46;
 663           tmp358 = tmp354 - tmp357;
 664           tmp365 = tmp361 - tmp364;
 665           tmp366 = K707106781 * (tmp358 + tmp365);
 666           tmp368 = K707106781 * (tmp365 - tmp358);
 667           real_output[24 * real_ostride] = tmp351 - tmp366;
 668           real_output[8 * real_ostride] = tmp351 + tmp366;
 669           imag_output[8 * imag_ostride] = tmp367 + tmp368;
 670           imag_output[24 * imag_ostride] = tmp368 - tmp367;
 671      }
 672      {
 673           fftw_real tmp75;
 674           fftw_real tmp189;
 675           fftw_real tmp94;
 676           fftw_real tmp199;
 677           fftw_real tmp173;
 678           fftw_real tmp190;
 679           fftw_real tmp132;
 680           fftw_real tmp184;
 681           fftw_real tmp180;
 682           fftw_real tmp200;
 683           fftw_real tmp194;
 684           fftw_real tmp204;
 685           fftw_real tmp169;
 686           fftw_real tmp185;
 687           fftw_real tmp197;
 688           fftw_real tmp205;
 689           ASSERT_ALIGNED_DOUBLE();
 690           {
 691                fftw_real tmp84;
 692                fftw_real tmp93;
 693                fftw_real tmp192;
 694                fftw_real tmp193;
 695                ASSERT_ALIGNED_DOUBLE();
 696                tmp75 = tmp67 - tmp74;
 697                tmp189 = tmp67 + tmp74;
 698                tmp84 = (K195090322 * tmp80) + (K980785280 * tmp83);
 699                tmp93 = (K980785280 * tmp89) - (K195090322 * tmp92);
 700                tmp94 = tmp84 - tmp93;
 701                tmp199 = tmp93 + tmp84;
 702                {
 703                     fftw_real tmp171;
 704                     fftw_real tmp172;
 705                     fftw_real tmp122;
 706                     fftw_real tmp131;
 707                     ASSERT_ALIGNED_DOUBLE();
 708                     tmp171 = (K980785280 * tmp80) - (K195090322 * tmp83);
 709                     tmp172 = (K980785280 * tmp92) + (K195090322 * tmp89);
 710                     tmp173 = tmp171 - tmp172;
 711                     tmp190 = tmp172 + tmp171;
 712                     tmp122 = tmp110 - tmp121;
 713                     tmp131 = tmp127 - tmp130;
 714                     tmp132 = (K634393284 * tmp122) + (K773010453 * tmp131);
 715                     tmp184 = (K773010453 * tmp122) - (K634393284 * tmp131);
 716                }
 717                tmp180 = tmp176 - tmp179;
 718                tmp200 = tmp179 + tmp176;
 719                tmp192 = tmp121 + tmp110;
 720                tmp193 = tmp127 + tmp130;
 721                tmp194 = (K098017140 * tmp192) + (K995184726 * tmp193);
 722                tmp204 = (K995184726 * tmp192) - (K098017140 * tmp193);
 723                {
 724                     fftw_real tmp159;
 725                     fftw_real tmp168;
 726                     fftw_real tmp195;
 727                     fftw_real tmp196;
 728                     ASSERT_ALIGNED_DOUBLE();
 729                     tmp159 = tmp143 - tmp158;
 730                     tmp168 = tmp162 - tmp167;
 731                     tmp169 = (K773010453 * tmp159) - (K634393284 * tmp168);
 732                     tmp185 = (K773010453 * tmp168) + (K634393284 * tmp159);
 733                     tmp195 = tmp143 + tmp158;
 734                     tmp196 = tmp167 + tmp162;
 735                     tmp197 = (K995184726 * tmp195) - (K098017140 * tmp196);
 736                     tmp205 = (K995184726 * tmp196) + (K098017140 * tmp195);
 737                }
 738           }
 739           {
 740                fftw_real tmp95;
 741                fftw_real tmp170;
 742                fftw_real tmp181;
 743                fftw_real tmp182;
 744                ASSERT_ALIGNED_DOUBLE();
 745                tmp95 = tmp75 + tmp94;
 746                tmp170 = tmp132 + tmp169;
 747                real_output[25 * real_ostride] = tmp95 - tmp170;
 748                real_output[7 * real_ostride] = tmp95 + tmp170;
 749                tmp181 = tmp173 - tmp180;
 750                tmp182 = tmp169 - tmp132;
 751                imag_output[9 * imag_ostride] = tmp181 + tmp182;
 752                imag_output[23 * imag_ostride] = tmp182 - tmp181;
 753           }
 754           {
 755                fftw_real tmp183;
 756                fftw_real tmp186;
 757                fftw_real tmp187;
 758                fftw_real tmp188;
 759                ASSERT_ALIGNED_DOUBLE();
 760                tmp183 = tmp180 + tmp173;
 761                tmp186 = tmp184 + tmp185;
 762                imag_output[7 * imag_ostride] = tmp183 + tmp186;
 763                imag_output[25 * imag_ostride] = tmp186 - tmp183;
 764                tmp187 = tmp75 - tmp94;
 765                tmp188 = tmp185 - tmp184;
 766                real_output[23 * real_ostride] = tmp187 - tmp188;
 767                real_output[9 * real_ostride] = tmp187 + tmp188;
 768           }
 769           {
 770                fftw_real tmp191;
 771                fftw_real tmp198;
 772                fftw_real tmp201;
 773                fftw_real tmp202;
 774                ASSERT_ALIGNED_DOUBLE();
 775                tmp191 = tmp189 + tmp190;
 776                tmp198 = tmp194 + tmp197;
 777                real_output[31 * real_ostride] = tmp191 - tmp198;
 778                real_output[real_ostride] = tmp191 + tmp198;
 779                tmp201 = tmp199 - tmp200;
 780                tmp202 = tmp197 - tmp194;
 781                imag_output[15 * imag_ostride] = tmp201 + tmp202;
 782                imag_output[17 * imag_ostride] = tmp202 - tmp201;
 783           }
 784           {
 785                fftw_real tmp203;
 786                fftw_real tmp206;
 787                fftw_real tmp207;
 788                fftw_real tmp208;
 789                ASSERT_ALIGNED_DOUBLE();
 790                tmp203 = tmp200 + tmp199;
 791                tmp206 = tmp204 + tmp205;
 792                imag_output[imag_ostride] = tmp203 + tmp206;
 793                imag_output[31 * imag_ostride] = tmp206 - tmp203;
 794                tmp207 = tmp189 - tmp190;
 795                tmp208 = tmp205 - tmp204;
 796                real_output[17 * real_ostride] = tmp207 - tmp208;
 797                real_output[15 * real_ostride] = tmp207 + tmp208;
 798           }
 799      }
 800      {
 801           fftw_real tmp377;
 802           fftw_real tmp393;
 803           fftw_real tmp387;
 804           fftw_real tmp389;
 805           fftw_real tmp384;
 806           fftw_real tmp388;
 807           fftw_real tmp392;
 808           fftw_real tmp394;
 809           ASSERT_ALIGNED_DOUBLE();
 810           {
 811                fftw_real tmp373;
 812                fftw_real tmp376;
 813                fftw_real tmp385;
 814                fftw_real tmp386;
 815                ASSERT_ALIGNED_DOUBLE();
 816                tmp373 = tmp7 - tmp14;
 817                tmp376 = K707106781 * (tmp374 + tmp375);
 818                tmp377 = tmp373 + tmp376;
 819                tmp393 = tmp373 - tmp376;
 820                tmp385 = K707106781 * (tmp375 - tmp374);
 821                tmp386 = tmp29 - tmp22;
 822                tmp387 = tmp385 - tmp386;
 823                tmp389 = tmp386 + tmp385;
 824           }
 825           {
 826                fftw_real tmp380;
 827                fftw_real tmp383;
 828                fftw_real tmp390;
 829                fftw_real tmp391;
 830                ASSERT_ALIGNED_DOUBLE();
 831                tmp380 = (K923879532 * tmp378) + (K382683432 * tmp379);
 832                tmp383 = (K923879532 * tmp381) - (K382683432 * tmp382);
 833                tmp384 = tmp380 + tmp383;
 834                tmp388 = tmp383 - tmp380;
 835                tmp390 = (K923879532 * tmp379) - (K382683432 * tmp378);
 836                tmp391 = (K382683432 * tmp381) + (K923879532 * tmp382);
 837                tmp392 = tmp390 + tmp391;
 838                tmp394 = tmp391 - tmp390;
 839           }
 840           real_output[28 * real_ostride] = tmp377 - tmp384;
 841           real_output[4 * real_ostride] = tmp377 + tmp384;
 842           imag_output[12 * imag_ostride] = tmp387 + tmp388;
 843           imag_output[20 * imag_ostride] = tmp388 - tmp387;
 844           imag_output[4 * imag_ostride] = tmp389 + tmp392;
 845           imag_output[28 * imag_ostride] = tmp392 - tmp389;
 846           real_output[20 * real_ostride] = tmp393 - tmp394;
 847           real_output[12 * real_ostride] = tmp393 + tmp394;
 848      }
 849      {
 850           fftw_real tmp281;
 851           fftw_real tmp329;
 852           fftw_real tmp323;
 853           fftw_real tmp325;
 854           fftw_real tmp298;
 855           fftw_real tmp326;
 856           fftw_real tmp315;
 857           fftw_real tmp327;
 858           ASSERT_ALIGNED_DOUBLE();
 859           {
 860                fftw_real tmp273;
 861                fftw_real tmp280;
 862                fftw_real tmp319;
 863                fftw_real tmp322;
 864                ASSERT_ALIGNED_DOUBLE();
 865                tmp273 = tmp269 + tmp272;
 866                tmp280 = tmp276 + tmp279;
 867                tmp281 = tmp273 + tmp280;
 868                tmp329 = tmp273 - tmp280;
 869                tmp319 = tmp317 + tmp318;
 870                tmp322 = tmp320 + tmp321;
 871                tmp323 = tmp319 - tmp322;
 872                tmp325 = tmp322 + tmp319;
 873           }
 874           {
 875                fftw_real tmp292;
 876                fftw_real tmp297;
 877                fftw_real tmp309;
 878                fftw_real tmp314;
 879                ASSERT_ALIGNED_DOUBLE();
 880                tmp292 = tmp284 + tmp291;
 881                tmp297 = tmp295 + tmp296;
 882                tmp298 = (K980785280 * tmp292) + (K195090322 * tmp297);
 883                tmp326 = (K980785280 * tmp297) - (K195090322 * tmp292);
 884                tmp309 = tmp301 + tmp308;
 885                tmp314 = tmp312 + tmp313;
 886                tmp315 = (K980785280 * tmp309) - (K195090322 * tmp314);
 887                tmp327 = (K195090322 * tmp309) + (K980785280 * tmp314);
 888           }
 889           {
 890                fftw_real tmp316;
 891                fftw_real tmp324;
 892                fftw_real tmp328;
 893                fftw_real tmp330;
 894                ASSERT_ALIGNED_DOUBLE();
 895                tmp316 = tmp298 + tmp315;
 896                real_output[30 * real_ostride] = tmp281 - tmp316;
 897                real_output[2 * real_ostride] = tmp281 + tmp316;
 898                tmp324 = tmp315 - tmp298;
 899                imag_output[14 * imag_ostride] = tmp323 + tmp324;
 900                imag_output[18 * imag_ostride] = tmp324 - tmp323;
 901                tmp328 = tmp326 + tmp327;
 902                imag_output[2 * imag_ostride] = tmp325 + tmp328;
 903                imag_output[30 * imag_ostride] = tmp328 - tmp325;
 904                tmp330 = tmp327 - tmp326;
 905                real_output[18 * real_ostride] = tmp329 - tmp330;
 906                real_output[14 * real_ostride] = tmp329 + tmp330;
 907           }
 908      }
 909      {
 910           fftw_real tmp333;
 911           fftw_real tmp349;
 912           fftw_real tmp343;
 913           fftw_real tmp345;
 914           fftw_real tmp336;
 915           fftw_real tmp346;
 916           fftw_real tmp339;
 917           fftw_real tmp347;
 918           ASSERT_ALIGNED_DOUBLE();
 919           {
 920                fftw_real tmp331;
 921                fftw_real tmp332;
 922                fftw_real tmp341;
 923                fftw_real tmp342;
 924                ASSERT_ALIGNED_DOUBLE();
 925                tmp331 = tmp269 - tmp272;
 926                tmp332 = tmp318 - tmp317;
 927                tmp333 = tmp331 + tmp332;
 928                tmp349 = tmp331 - tmp332;
 929                tmp341 = tmp279 - tmp276;
 930                tmp342 = tmp321 - tmp320;
 931                tmp343 = tmp341 - tmp342;
 932                tmp345 = tmp342 + tmp341;
 933           }
 934           {
 935                fftw_real tmp334;
 936                fftw_real tmp335;
 937                fftw_real tmp337;
 938                fftw_real tmp338;
 939                ASSERT_ALIGNED_DOUBLE();
 940                tmp334 = tmp284 - tmp291;
 941                tmp335 = tmp296 - tmp295;
 942                tmp336 = (K831469612 * tmp334) + (K555570233 * tmp335);
 943                tmp346 = (K831469612 * tmp335) - (K555570233 * tmp334);
 944                tmp337 = tmp301 - tmp308;
 945                tmp338 = tmp313 - tmp312;
 946                tmp339 = (K831469612 * tmp337) - (K555570233 * tmp338);
 947                tmp347 = (K555570233 * tmp337) + (K831469612 * tmp338);
 948           }
 949           {
 950                fftw_real tmp340;
 951                fftw_real tmp344;
 952                fftw_real tmp348;
 953                fftw_real tmp350;
 954                ASSERT_ALIGNED_DOUBLE();
 955                tmp340 = tmp336 + tmp339;
 956                real_output[26 * real_ostride] = tmp333 - tmp340;
 957                real_output[6 * real_ostride] = tmp333 + tmp340;
 958                tmp344 = tmp339 - tmp336;
 959                imag_output[10 * imag_ostride] = tmp343 + tmp344;
 960                imag_output[22 * imag_ostride] = tmp344 - tmp343;
 961                tmp348 = tmp346 + tmp347;
 962                imag_output[6 * imag_ostride] = tmp345 + tmp348;
 963                imag_output[26 * imag_ostride] = tmp348 - tmp345;
 964                tmp350 = tmp347 - tmp346;
 965                real_output[22 * real_ostride] = tmp349 - tmp350;
 966                real_output[10 * real_ostride] = tmp349 + tmp350;
 967           }
 968      }
 969      {
 970           fftw_real tmp211;
 971           fftw_real tmp249;
 972           fftw_real tmp218;
 973           fftw_real tmp259;
 974           fftw_real tmp237;
 975           fftw_real tmp250;
 976           fftw_real tmp226;
 977           fftw_real tmp244;
 978           fftw_real tmp240;
 979           fftw_real tmp260;
 980           fftw_real tmp254;
 981           fftw_real tmp264;
 982           fftw_real tmp233;
 983           fftw_real tmp245;
 984           fftw_real tmp257;
 985           fftw_real tmp265;
 986           ASSERT_ALIGNED_DOUBLE();
 987           {
 988                fftw_real tmp214;
 989                fftw_real tmp217;
 990                fftw_real tmp252;
 991                fftw_real tmp253;
 992                ASSERT_ALIGNED_DOUBLE();
 993                tmp211 = tmp209 + tmp210;
 994                tmp249 = tmp209 - tmp210;
 995                tmp214 = (K831469612 * tmp212) + (K555570233 * tmp213);
 996                tmp217 = (K831469612 * tmp215) - (K555570233 * tmp216);
 997                tmp218 = tmp214 + tmp217;
 998                tmp259 = tmp217 - tmp214;
 999                {
1000                     fftw_real tmp235;
1001                     fftw_real tmp236;
1002                     fftw_real tmp222;
1003                     fftw_real tmp225;
1004                     ASSERT_ALIGNED_DOUBLE();
1005                     tmp235 = (K831469612 * tmp213) - (K555570233 * tmp212);
1006                     tmp236 = (K555570233 * tmp215) + (K831469612 * tmp216);
1007                     tmp237 = tmp235 + tmp236;
1008                     tmp250 = tmp236 - tmp235;
1009                     tmp222 = tmp220 + tmp221;
1010                     tmp225 = tmp223 + tmp224;
1011                     tmp226 = (K956940335 * tmp222) + (K290284677 * tmp225);
1012                     tmp244 = (K956940335 * tmp225) - (K290284677 * tmp222);
1013                }
1014                tmp240 = tmp238 + tmp239;
1015                tmp260 = tmp239 - tmp238;
1016                tmp252 = tmp220 - tmp221;
1017                tmp253 = tmp224 - tmp223;
1018                tmp254 = (K881921264 * tmp252) + (K471396736 * tmp253);
1019                tmp264 = (K881921264 * tmp253) - (K471396736 * tmp252);
1020                {
1021                     fftw_real tmp229;
1022                     fftw_real tmp232;
1023                     fftw_real tmp255;
1024                     fftw_real tmp256;
1025                     ASSERT_ALIGNED_DOUBLE();
1026                     tmp229 = tmp227 + tmp228;
1027                     tmp232 = tmp230 + tmp231;
1028                     tmp233 = (K956940335 * tmp229) - (K290284677 * tmp232);
1029                     tmp245 = (K290284677 * tmp229) + (K956940335 * tmp232);
1030                     tmp255 = tmp227 - tmp228;
1031                     tmp256 = tmp231 - tmp230;
1032                     tmp257 = (K881921264 * tmp255) - (K471396736 * tmp256);
1033                     tmp265 = (K471396736 * tmp255) + (K881921264 * tmp256);
1034                }
1035           }
1036           {
1037                fftw_real tmp219;
1038                fftw_real tmp234;
1039                fftw_real tmp241;
1040                fftw_real tmp242;
1041                ASSERT_ALIGNED_DOUBLE();
1042                tmp219 = tmp211 + tmp218;
1043                tmp234 = tmp226 + tmp233;
1044                real_output[29 * real_ostride] = tmp219 - tmp234;
1045                real_output[3 * real_ostride] = tmp219 + tmp234;
1046                tmp241 = tmp237 - tmp240;
1047                tmp242 = tmp233 - tmp226;
1048                imag_output[13 * imag_ostride] = tmp241 + tmp242;
1049                imag_output[19 * imag_ostride] = tmp242 - tmp241;
1050           }
1051           {
1052                fftw_real tmp243;
1053                fftw_real tmp246;
1054                fftw_real tmp247;
1055                fftw_real tmp248;
1056                ASSERT_ALIGNED_DOUBLE();
1057                tmp243 = tmp240 + tmp237;
1058                tmp246 = tmp244 + tmp245;
1059                imag_output[3 * imag_ostride] = tmp243 + tmp246;
1060                imag_output[29 * imag_ostride] = tmp246 - tmp243;
1061                tmp247 = tmp211 - tmp218;
1062                tmp248 = tmp245 - tmp244;
1063                real_output[19 * real_ostride] = tmp247 - tmp248;
1064                real_output[13 * real_ostride] = tmp247 + tmp248;
1065           }
1066           {
1067                fftw_real tmp251;
1068                fftw_real tmp258;
1069                fftw_real tmp261;
1070                fftw_real tmp262;
1071                ASSERT_ALIGNED_DOUBLE();
1072                tmp251 = tmp249 + tmp250;
1073                tmp258 = tmp254 + tmp257;
1074                real_output[27 * real_ostride] = tmp251 - tmp258;
1075                real_output[5 * real_ostride] = tmp251 + tmp258;
1076                tmp261 = tmp259 - tmp260;
1077                tmp262 = tmp257 - tmp254;
1078                imag_output[11 * imag_ostride] = tmp261 + tmp262;
1079                imag_output[21 * imag_ostride] = tmp262 - tmp261;
1080           }
1081           {
1082                fftw_real tmp263;
1083                fftw_real tmp266;
1084                fftw_real tmp267;
1085                fftw_real tmp268;
1086                ASSERT_ALIGNED_DOUBLE();
1087                tmp263 = tmp260 + tmp259;
1088                tmp266 = tmp264 + tmp265;
1089                imag_output[5 * imag_ostride] = tmp263 + tmp266;
1090                imag_output[27 * imag_ostride] = tmp266 - tmp263;
1091                tmp267 = tmp249 - tmp250;
1092                tmp268 = tmp265 - tmp264;
1093                real_output[21 * real_ostride] = tmp267 - tmp268;
1094                real_output[11 * real_ostride] = tmp267 + tmp268;
1095           }
1096      }
1097 }
1098
1099 fftw_codelet_desc fftw_real2hc_64_desc =
1100 {
1101      "fftw_real2hc_64",
1102      (void (*)()) fftw_real2hc_64,
1103      64,
1104      FFTW_FORWARD,
1105      FFTW_REAL2HC,
1106      1410,
1107      0,
1108      (const int *) 0,
1109 };