src/fftw/fhf_9.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:55:35 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 -hc2hc-forward 9 */
  27
  28 /*
  29  * This function contains 180 FP additions, 120 FP multiplications,
  30  * (or, 133 additions, 73 multiplications, 47 fused multiply/add),
  31  * 35 stack variables, and 72 memory accesses
  32  */
  33 static const fftw_real K342020143 = FFTW_KONST(+0.342020143325668733044099614682259580763083368);
  34 static const fftw_real K813797681 = FFTW_KONST(+0.813797681349373692844693217248393223289101568);
  35 static const fftw_real K939692620 = FFTW_KONST(+0.939692620785908384054109277324731469936208134);
  36 static const fftw_real K296198132 = FFTW_KONST(+0.296198132726023843175338011893050938967728390);
  37 static const fftw_real K852868531 = FFTW_KONST(+0.852868531952443209628250963940074071936020296);
  38 static const fftw_real K173648177 = FFTW_KONST(+0.173648177666930348851716626769314796000375677);
  39 static const fftw_real K556670399 = FFTW_KONST(+0.556670399226419366452912952047023132968291906);
  40 static const fftw_real K766044443 = FFTW_KONST(+0.766044443118978035202392650555416673935832457);
  41 static const fftw_real K984807753 = FFTW_KONST(+0.984807753012208059366743024589523013670643252);
  42 static const fftw_real K150383733 = FFTW_KONST(+0.150383733180435296639271897612501926072238258);
  43 static const fftw_real K642787609 = FFTW_KONST(+0.642787609686539326322643409907263432907559884);
  44 static const fftw_real K663413948 = FFTW_KONST(+0.663413948168938396205421319635891297216863310);
  45 static const fftw_real K866025403 = FFTW_KONST(+0.866025403784438646763723170752936183471402627);
  46 static const fftw_real K500000000 = FFTW_KONST(+0.500000000000000000000000000000000000000000000);
  47
  48 /*
  49  * Generator Id's :
  50  * $Id$
  51  * $Id$
  52  * $Id$
  53  */
  54
  55 void fftw_hc2hc_forward_9(fftw_real *A, const fftw_complex *W, int iostride, int m, int dist)
  56 {
  57      int i;
  58      fftw_real *X;
  59      fftw_real *Y;
  60      X = A;
  61      Y = A + (9 * iostride);
  62      {
  63           fftw_real tmp131;
  64           fftw_real tmp145;
  65           fftw_real tmp150;
  66           fftw_real tmp149;
  67           fftw_real tmp134;
  68           fftw_real tmp157;
  69           fftw_real tmp140;
  70           fftw_real tmp148;
  71           fftw_real tmp151;
  72           fftw_real tmp132;
  73           fftw_real tmp133;
  74           fftw_real tmp135;
  75           fftw_real tmp146;
  76           ASSERT_ALIGNED_DOUBLE();
  77           tmp131 = X[0];
  78           {
  79                fftw_real tmp141;
  80                fftw_real tmp142;
  81                fftw_real tmp143;
  82                fftw_real tmp144;
  83                ASSERT_ALIGNED_DOUBLE();
  84                tmp141 = X[2 * iostride];
  85                tmp142 = X[5 * iostride];
  86                tmp143 = X[8 * iostride];
  87                tmp144 = tmp142 + tmp143;
  88                tmp145 = tmp141 + tmp144;
  89                tmp150 = tmp141 - (K500000000 * tmp144);
  90                tmp149 = tmp143 - tmp142;
  91           }
  92           tmp132 = X[3 * iostride];
  93           tmp133 = X[6 * iostride];
  94           tmp134 = tmp132 + tmp133;
  95           tmp157 = tmp133 - tmp132;
  96           {
  97                fftw_real tmp136;
  98                fftw_real tmp137;
  99                fftw_real tmp138;
 100                fftw_real tmp139;
 101                ASSERT_ALIGNED_DOUBLE();
 102                tmp136 = X[iostride];
 103                tmp137 = X[4 * iostride];
 104                tmp138 = X[7 * iostride];
 105                tmp139 = tmp137 + tmp138;
 106                tmp140 = tmp136 + tmp139;
 107                tmp148 = tmp136 - (K500000000 * tmp139);
 108                tmp151 = tmp138 - tmp137;
 109           }
 110           Y[-3 * iostride] = K866025403 * (tmp145 - tmp140);
 111           tmp135 = tmp131 + tmp134;
 112           tmp146 = tmp140 + tmp145;
 113           X[3 * iostride] = tmp135 - (K500000000 * tmp146);
 114           X[0] = tmp135 + tmp146;
 115           {
 116                fftw_real tmp159;
 117                fftw_real tmp155;
 118                fftw_real tmp156;
 119                fftw_real tmp158;
 120                fftw_real tmp147;
 121                fftw_real tmp152;
 122                fftw_real tmp153;
 123                fftw_real tmp154;
 124                ASSERT_ALIGNED_DOUBLE();
 125                tmp159 = K866025403 * tmp157;
 126                tmp155 = (K663413948 * tmp151) - (K642787609 * tmp148);
 127                tmp156 = (K150383733 * tmp149) - (K984807753 * tmp150);
 128                tmp158 = tmp155 + tmp156;
 129                tmp147 = tmp131 - (K500000000 * tmp134);
 130                tmp152 = (K766044443 * tmp148) + (K556670399 * tmp151);
 131                tmp153 = (K173648177 * tmp150) + (K852868531 * tmp149);
 132                tmp154 = tmp152 + tmp153;
 133                X[iostride] = tmp147 + tmp154;
 134                X[4 * iostride] = tmp147 + (K866025403 * (tmp155 - tmp156)) - (K500000000 * tmp154);
 135                X[2 * iostride] = tmp147 + (K173648177 * tmp148) - (K296198132 * tmp149) - (K939692620 * tmp150) - (K852868531 * tmp151);
 136                Y[-iostride] = tmp159 + tmp158;
 137                Y[-4 * iostride] = (K866025403 * (tmp157 + (tmp153 - tmp152))) - (K500000000 * tmp158);
 138                Y[-2 * iostride] = (K813797681 * tmp149) - (K342020143 * tmp150) - (K150383733 * tmp151) - (K984807753 * tmp148) - tmp159;
 139           }
 140      }
 141      X = X + dist;
 142      Y = Y - dist;
 143      for (i = 2; i < m; i = i + 2, X = X + dist, Y = Y - dist, W = W + 8) {
 144           fftw_real tmp19;
 145           fftw_real tmp117;
 146           fftw_real tmp70;
 147           fftw_real tmp116;
 148           fftw_real tmp123;
 149           fftw_real tmp122;
 150           fftw_real tmp30;
 151           fftw_real tmp67;
 152           fftw_real tmp65;
 153           fftw_real tmp87;
 154           fftw_real tmp104;
 155           fftw_real tmp113;
 156           fftw_real tmp92;
 157           fftw_real tmp103;
 158           fftw_real tmp48;
 159           fftw_real tmp76;
 160           fftw_real tmp100;
 161           fftw_real tmp112;
 162           fftw_real tmp81;
 163           fftw_real tmp101;
 164           ASSERT_ALIGNED_DOUBLE();
 165           {
 166                fftw_real tmp24;
 167                fftw_real tmp68;
 168                fftw_real tmp29;
 169                fftw_real tmp69;
 170                ASSERT_ALIGNED_DOUBLE();
 171                tmp19 = X[0];
 172                tmp117 = Y[-8 * iostride];
 173                {
 174                     fftw_real tmp21;
 175                     fftw_real tmp23;
 176                     fftw_real tmp20;
 177                     fftw_real tmp22;
 178                     ASSERT_ALIGNED_DOUBLE();
 179                     tmp21 = X[3 * iostride];
 180                     tmp23 = Y[-5 * iostride];
 181                     tmp20 = c_re(W[2]);
 182                     tmp22 = c_im(W[2]);
 183                     tmp24 = (tmp20 * tmp21) - (tmp22 * tmp23);
 184                     tmp68 = (tmp22 * tmp21) + (tmp20 * tmp23);
 185                }
 186                {
 187                     fftw_real tmp26;
 188                     fftw_real tmp28;
 189                     fftw_real tmp25;
 190                     fftw_real tmp27;
 191                     ASSERT_ALIGNED_DOUBLE();
 192                     tmp26 = X[6 * iostride];
 193                     tmp28 = Y[-2 * iostride];
 194                     tmp25 = c_re(W[5]);
 195                     tmp27 = c_im(W[5]);
 196                     tmp29 = (tmp25 * tmp26) - (tmp27 * tmp28);
 197                     tmp69 = (tmp27 * tmp26) + (tmp25 * tmp28);
 198                }
 199                tmp70 = K866025403 * (tmp68 - tmp69);
 200                tmp116 = tmp68 + tmp69;
 201                tmp123 = tmp117 - (K500000000 * tmp116);
 202                tmp122 = K866025403 * (tmp29 - tmp24);
 203                tmp30 = tmp24 + tmp29;
 204                tmp67 = tmp19 - (K500000000 * tmp30);
 205           }
 206           {
 207                fftw_real tmp53;
 208                fftw_real tmp89;
 209                fftw_real tmp58;
 210                fftw_real tmp84;
 211                fftw_real tmp63;
 212                fftw_real tmp85;
 213                fftw_real tmp64;
 214                fftw_real tmp90;
 215                ASSERT_ALIGNED_DOUBLE();
 216                {
 217                     fftw_real tmp50;
 218                     fftw_real tmp52;
 219                     fftw_real tmp49;
 220                     fftw_real tmp51;
 221                     ASSERT_ALIGNED_DOUBLE();
 222                     tmp50 = X[2 * iostride];
 223                     tmp52 = Y[-6 * iostride];
 224                     tmp49 = c_re(W[1]);
 225                     tmp51 = c_im(W[1]);
 226                     tmp53 = (tmp49 * tmp50) - (tmp51 * tmp52);
 227                     tmp89 = (tmp51 * tmp50) + (tmp49 * tmp52);
 228                }
 229                {
 230                     fftw_real tmp55;
 231                     fftw_real tmp57;
 232                     fftw_real tmp54;
 233                     fftw_real tmp56;
 234                     ASSERT_ALIGNED_DOUBLE();
 235                     tmp55 = X[5 * iostride];
 236                     tmp57 = Y[-3 * iostride];
 237                     tmp54 = c_re(W[4]);
 238                     tmp56 = c_im(W[4]);
 239                     tmp58 = (tmp54 * tmp55) - (tmp56 * tmp57);
 240                     tmp84 = (tmp56 * tmp55) + (tmp54 * tmp57);
 241                }
 242                {
 243                     fftw_real tmp60;
 244                     fftw_real tmp62;
 245                     fftw_real tmp59;
 246                     fftw_real tmp61;
 247                     ASSERT_ALIGNED_DOUBLE();
 248                     tmp60 = X[8 * iostride];
 249                     tmp62 = Y[0];
 250                     tmp59 = c_re(W[7]);
 251                     tmp61 = c_im(W[7]);
 252                     tmp63 = (tmp59 * tmp60) - (tmp61 * tmp62);
 253                     tmp85 = (tmp61 * tmp60) + (tmp59 * tmp62);
 254                }
 255                tmp64 = tmp58 + tmp63;
 256                tmp90 = tmp84 + tmp85;
 257                {
 258                     fftw_real tmp83;
 259                     fftw_real tmp86;
 260                     fftw_real tmp88;
 261                     fftw_real tmp91;
 262                     ASSERT_ALIGNED_DOUBLE();
 263                     tmp65 = tmp53 + tmp64;
 264                     tmp83 = tmp53 - (K500000000 * tmp64);
 265                     tmp86 = K866025403 * (tmp84 - tmp85);
 266                     tmp87 = tmp83 + tmp86;
 267                     tmp104 = tmp83 - tmp86;
 268                     tmp113 = tmp89 + tmp90;
 269                     tmp88 = K866025403 * (tmp63 - tmp58);
 270                     tmp91 = tmp89 - (K500000000 * tmp90);
 271                     tmp92 = tmp88 + tmp91;
 272                     tmp103 = tmp91 - tmp88;
 273                }
 274           }
 275           {
 276                fftw_real tmp36;
 277                fftw_real tmp78;
 278                fftw_real tmp41;
 279                fftw_real tmp73;
 280                fftw_real tmp46;
 281                fftw_real tmp74;
 282                fftw_real tmp47;
 283                fftw_real tmp79;
 284                ASSERT_ALIGNED_DOUBLE();
 285                {
 286                     fftw_real tmp33;
 287                     fftw_real tmp35;
 288                     fftw_real tmp32;
 289                     fftw_real tmp34;
 290                     ASSERT_ALIGNED_DOUBLE();
 291                     tmp33 = X[iostride];
 292                     tmp35 = Y[-7 * iostride];
 293                     tmp32 = c_re(W[0]);
 294                     tmp34 = c_im(W[0]);
 295                     tmp36 = (tmp32 * tmp33) - (tmp34 * tmp35);
 296                     tmp78 = (tmp34 * tmp33) + (tmp32 * tmp35);
 297                }
 298                {
 299                     fftw_real tmp38;
 300                     fftw_real tmp40;
 301                     fftw_real tmp37;
 302                     fftw_real tmp39;
 303                     ASSERT_ALIGNED_DOUBLE();
 304                     tmp38 = X[4 * iostride];
 305                     tmp40 = Y[-4 * iostride];
 306                     tmp37 = c_re(W[3]);
 307                     tmp39 = c_im(W[3]);
 308                     tmp41 = (tmp37 * tmp38) - (tmp39 * tmp40);
 309                     tmp73 = (tmp39 * tmp38) + (tmp37 * tmp40);
 310                }
 311                {
 312                     fftw_real tmp43;
 313                     fftw_real tmp45;
 314                     fftw_real tmp42;
 315                     fftw_real tmp44;
 316                     ASSERT_ALIGNED_DOUBLE();
 317                     tmp43 = X[7 * iostride];
 318                     tmp45 = Y[-iostride];
 319                     tmp42 = c_re(W[6]);
 320                     tmp44 = c_im(W[6]);
 321                     tmp46 = (tmp42 * tmp43) - (tmp44 * tmp45);
 322                     tmp74 = (tmp44 * tmp43) + (tmp42 * tmp45);
 323                }
 324                tmp47 = tmp41 + tmp46;
 325                tmp79 = tmp73 + tmp74;
 326                {
 327                     fftw_real tmp72;
 328                     fftw_real tmp75;
 329                     fftw_real tmp77;
 330                     fftw_real tmp80;
 331                     ASSERT_ALIGNED_DOUBLE();
 332                     tmp48 = tmp36 + tmp47;
 333                     tmp72 = tmp36 - (K500000000 * tmp47);
 334                     tmp75 = K866025403 * (tmp73 - tmp74);
 335                     tmp76 = tmp72 + tmp75;
 336                     tmp100 = tmp72 - tmp75;
 337                     tmp112 = tmp78 + tmp79;
 338                     tmp77 = K866025403 * (tmp46 - tmp41);
 339                     tmp80 = tmp78 - (K500000000 * tmp79);
 340                     tmp81 = tmp77 + tmp80;
 341                     tmp101 = tmp80 - tmp77;
 342                }
 343           }
 344           {
 345                fftw_real tmp114;
 346                fftw_real tmp31;
 347                fftw_real tmp66;
 348                fftw_real tmp111;
 349                ASSERT_ALIGNED_DOUBLE();
 350                tmp114 = K866025403 * (tmp112 - tmp113);
 351                tmp31 = tmp19 + tmp30;
 352                tmp66 = tmp48 + tmp65;
 353                tmp111 = tmp31 - (K500000000 * tmp66);
 354                X[0] = tmp31 + tmp66;
 355                X[3 * iostride] = tmp111 + tmp114;
 356                Y[-6 * iostride] = tmp111 - tmp114;
 357           }
 358           {
 359                fftw_real tmp120;
 360                fftw_real tmp115;
 361                fftw_real tmp118;
 362                fftw_real tmp119;
 363                ASSERT_ALIGNED_DOUBLE();
 364                tmp120 = K866025403 * (tmp65 - tmp48);
 365                tmp115 = tmp112 + tmp113;
 366                tmp118 = tmp116 + tmp117;
 367                tmp119 = tmp118 - (K500000000 * tmp115);
 368                Y[0] = tmp115 + tmp118;
 369                Y[-3 * iostride] = tmp120 + tmp119;
 370                X[6 * iostride] = -(tmp119 - tmp120);
 371           }
 372           {
 373                fftw_real tmp71;
 374                fftw_real tmp124;
 375                fftw_real tmp94;
 376                fftw_real tmp126;
 377                fftw_real tmp98;
 378                fftw_real tmp121;
 379                fftw_real tmp95;
 380                fftw_real tmp125;
 381                ASSERT_ALIGNED_DOUBLE();
 382                tmp71 = tmp67 + tmp70;
 383                tmp124 = tmp122 + tmp123;
 384                {
 385                     fftw_real tmp82;
 386                     fftw_real tmp93;
 387                     fftw_real tmp96;
 388                     fftw_real tmp97;
 389                     ASSERT_ALIGNED_DOUBLE();
 390                     tmp82 = (K766044443 * tmp76) + (K642787609 * tmp81);
 391                     tmp93 = (K173648177 * tmp87) + (K984807753 * tmp92);
 392                     tmp94 = tmp82 + tmp93;
 393                     tmp126 = K866025403 * (tmp93 - tmp82);
 394                     tmp96 = (K766044443 * tmp81) - (K642787609 * tmp76);
 395                     tmp97 = (K173648177 * tmp92) - (K984807753 * tmp87);
 396                     tmp98 = K866025403 * (tmp96 - tmp97);
 397                     tmp121 = tmp96 + tmp97;
 398                }
 399                X[iostride] = tmp71 + tmp94;
 400                tmp95 = tmp71 - (K500000000 * tmp94);
 401                Y[-7 * iostride] = tmp95 - tmp98;
 402                X[4 * iostride] = tmp95 + tmp98;
 403                Y[-iostride] = tmp121 + tmp124;
 404                tmp125 = tmp124 - (K500000000 * tmp121);
 405                X[7 * iostride] = -(tmp125 - tmp126);
 406                Y[-4 * iostride] = tmp126 + tmp125;
 407           }
 408           {
 409                fftw_real tmp99;
 410                fftw_real tmp128;
 411                fftw_real tmp106;
 412                fftw_real tmp127;
 413                fftw_real tmp110;
 414                fftw_real tmp129;
 415                fftw_real tmp107;
 416                fftw_real tmp130;
 417                ASSERT_ALIGNED_DOUBLE();
 418                tmp99 = tmp67 - tmp70;
 419                tmp128 = tmp123 - tmp122;
 420                {
 421                     fftw_real tmp102;
 422                     fftw_real tmp105;
 423                     fftw_real tmp108;
 424                     fftw_real tmp109;
 425                     ASSERT_ALIGNED_DOUBLE();
 426                     tmp102 = (K173648177 * tmp100) + (K984807753 * tmp101);
 427                     tmp105 = (K342020143 * tmp103) - (K939692620 * tmp104);
 428                     tmp106 = tmp102 + tmp105;
 429                     tmp127 = K866025403 * (tmp105 - tmp102);
 430                     tmp108 = (K173648177 * tmp101) - (K984807753 * tmp100);
 431                     tmp109 = (K342020143 * tmp104) + (K939692620 * tmp103);
 432                     tmp110 = K866025403 * (tmp108 + tmp109);
 433                     tmp129 = tmp108 - tmp109;
 434                }
 435                X[2 * iostride] = tmp99 + tmp106;
 436                tmp107 = tmp99 - (K500000000 * tmp106);
 437                Y[-8 * iostride] = tmp107 - tmp110;
 438                Y[-5 * iostride] = tmp107 + tmp110;
 439                Y[-2 * iostride] = tmp129 + tmp128;
 440                tmp130 = tmp128 - (K500000000 * tmp129);
 441                X[5 * iostride] = -(tmp127 + tmp130);
 442                X[8 * iostride] = -(tmp130 - tmp127);
 443           }
 444      }
 445      if (i == m) {
 446           fftw_real tmp5;
 447           fftw_real tmp10;
 448           fftw_real tmp14;
 449           fftw_real tmp18;
 450           fftw_real tmp7;
 451           fftw_real tmp8;
 452           fftw_real tmp12;
 453           fftw_real tmp16;
 454           fftw_real tmp6;
 455           fftw_real tmp9;
 456           fftw_real tmp13;
 457           fftw_real tmp17;
 458           fftw_real tmp11;
 459           fftw_real tmp15;
 460           fftw_real tmp4;
 461           fftw_real tmp1;
 462           fftw_real tmp3;
 463           fftw_real tmp2;
 464           ASSERT_ALIGNED_DOUBLE();
 465           tmp5 = X[2 * iostride];
 466           tmp10 = X[7 * iostride];
 467           tmp14 = tmp5 + tmp10;
 468           tmp18 = tmp5 - tmp10;
 469           tmp7 = X[4 * iostride];
 470           tmp8 = X[5 * iostride];
 471           tmp12 = tmp8 + tmp7;
 472           tmp16 = tmp8 - tmp7;
 473           tmp6 = X[8 * iostride];
 474           tmp9 = X[iostride];
 475           tmp13 = tmp6 + tmp9;
 476           tmp17 = tmp6 - tmp9;
 477           tmp1 = X[0];
 478           tmp3 = X[3 * iostride];
 479           tmp2 = X[6 * iostride];
 480           tmp11 = K866025403 * (tmp3 + tmp2);
 481           tmp15 = tmp1 - (K500000000 * (tmp2 - tmp3));
 482           tmp4 = tmp1 + tmp2 - tmp3;
 483           Y[0] = -(tmp11 + (K984807753 * tmp12) + (K342020143 * tmp13) + (K642787609 * tmp14));
 484           Y[-3 * iostride] = (K342020143 * tmp12) + (K984807753 * tmp14) - (K642787609 * tmp13) - tmp11;
 485           Y[-2 * iostride] = tmp11 + (K342020143 * tmp14) - (K984807753 * tmp13) - (K642787609 * tmp12);
 486           X[2 * iostride] = tmp15 + (K173648177 * tmp17) - (K939692620 * tmp18) - (K766044443 * tmp16);
 487           X[3 * iostride] = tmp15 + (K939692620 * tmp16) + (K766044443 * tmp17) + (K173648177 * tmp18);
 488           X[0] = tmp15 + (K766044443 * tmp18) - (K939692620 * tmp17) - (K173648177 * tmp16);
 489           Y[-iostride] = K866025403 * (tmp8 + tmp7 - (tmp5 + tmp6 + tmp9 + tmp10));
 490           X[iostride] = tmp4 + (K500000000 * (tmp8 + tmp9 + tmp10 - (tmp5 + tmp6 + tmp7)));
 491           X[4 * iostride] = tmp4 + tmp5 + tmp6 + tmp7 - (tmp8 + tmp9 + tmp10);
 492      }
 493 }
 494
 495 static const int twiddle_order[] =
 496 {1, 2, 3, 4, 5, 6, 7, 8};
 497 fftw_codelet_desc fftw_hc2hc_forward_9_desc =
 498 {
 499      "fftw_hc2hc_forward_9",
 500      (void (*)()) fftw_hc2hc_forward_9,
 501      9,
 502      FFTW_FORWARD,
 503      FFTW_HC2HC,
 504      201,
 505      8,
 506      twiddle_order,
 507 };