src/fftw/ftwi_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:50 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 -twiddleinv 9 */
  27
  28 /*
  29  * This function contains 96 FP additions, 72 FP multiplications,
  30  * (or, 60 additions, 36 multiplications, 36 fused multiply/add),
  31  * 34 stack variables, and 36 memory accesses
  32  */
  33 static const fftw_real K642787609 = FFTW_KONST(+0.642787609686539326322643409907263432907559884);
  34 static const fftw_real K766044443 = FFTW_KONST(+0.766044443118978035202392650555416673935832457);
  35 static const fftw_real K939692620 = FFTW_KONST(+0.939692620785908384054109277324731469936208134);
  36 static const fftw_real K342020143 = FFTW_KONST(+0.342020143325668733044099614682259580763083368);
  37 static const fftw_real K984807753 = FFTW_KONST(+0.984807753012208059366743024589523013670643252);
  38 static const fftw_real K173648177 = FFTW_KONST(+0.173648177666930348851716626769314796000375677);
  39 static const fftw_real K500000000 = FFTW_KONST(+0.500000000000000000000000000000000000000000000);
  40 static const fftw_real K866025403 = FFTW_KONST(+0.866025403784438646763723170752936183471402627);
  41
  42 /*
  43  * Generator Id's :
  44  * $Id$
  45  * $Id$
  46  * $Id$
  47  */
  48
  49 void fftwi_twiddle_9(fftw_complex *A, const fftw_complex *W, int iostride, int m, int dist)
  50 {
  51      int i;
  52      fftw_complex *inout;
  53      inout = A;
  54      for (i = m; i > 0; i = i - 1, inout = inout + dist, W = W + 8) {
  55           fftw_real tmp1;
  56           fftw_real tmp99;
  57           fftw_real tmp64;
  58           fftw_real tmp98;
  59           fftw_real tmp105;
  60           fftw_real tmp104;
  61           fftw_real tmp12;
  62           fftw_real tmp61;
  63           fftw_real tmp47;
  64           fftw_real tmp78;
  65           fftw_real tmp89;
  66           fftw_real tmp54;
  67           fftw_real tmp75;
  68           fftw_real tmp90;
  69           fftw_real tmp30;
  70           fftw_real tmp68;
  71           fftw_real tmp86;
  72           fftw_real tmp59;
  73           fftw_real tmp71;
  74           fftw_real tmp87;
  75           ASSERT_ALIGNED_DOUBLE();
  76           {
  77                fftw_real tmp6;
  78                fftw_real tmp63;
  79                fftw_real tmp11;
  80                fftw_real tmp62;
  81                ASSERT_ALIGNED_DOUBLE();
  82                tmp1 = c_re(inout[0]);
  83                tmp99 = c_im(inout[0]);
  84                {
  85                     fftw_real tmp3;
  86                     fftw_real tmp5;
  87                     fftw_real tmp2;
  88                     fftw_real tmp4;
  89                     ASSERT_ALIGNED_DOUBLE();
  90                     tmp3 = c_re(inout[3 * iostride]);
  91                     tmp5 = c_im(inout[3 * iostride]);
  92                     tmp2 = c_re(W[2]);
  93                     tmp4 = c_im(W[2]);
  94                     tmp6 = (tmp2 * tmp3) + (tmp4 * tmp5);
  95                     tmp63 = (tmp2 * tmp5) - (tmp4 * tmp3);
  96                }
  97                {
  98                     fftw_real tmp8;
  99                     fftw_real tmp10;
 100                     fftw_real tmp7;
 101                     fftw_real tmp9;
 102                     ASSERT_ALIGNED_DOUBLE();
 103                     tmp8 = c_re(inout[6 * iostride]);
 104                     tmp10 = c_im(inout[6 * iostride]);
 105                     tmp7 = c_re(W[5]);
 106                     tmp9 = c_im(W[5]);
 107                     tmp11 = (tmp7 * tmp8) + (tmp9 * tmp10);
 108                     tmp62 = (tmp7 * tmp10) - (tmp9 * tmp8);
 109                }
 110                tmp64 = K866025403 * (tmp62 - tmp63);
 111                tmp98 = tmp63 + tmp62;
 112                tmp105 = tmp99 - (K500000000 * tmp98);
 113                tmp104 = K866025403 * (tmp6 - tmp11);
 114                tmp12 = tmp6 + tmp11;
 115                tmp61 = tmp1 - (K500000000 * tmp12);
 116           }
 117           {
 118                fftw_real tmp35;
 119                fftw_real tmp50;
 120                fftw_real tmp40;
 121                fftw_real tmp51;
 122                fftw_real tmp45;
 123                fftw_real tmp52;
 124                fftw_real tmp46;
 125                fftw_real tmp53;
 126                ASSERT_ALIGNED_DOUBLE();
 127                {
 128                     fftw_real tmp32;
 129                     fftw_real tmp34;
 130                     fftw_real tmp31;
 131                     fftw_real tmp33;
 132                     ASSERT_ALIGNED_DOUBLE();
 133                     tmp32 = c_re(inout[2 * iostride]);
 134                     tmp34 = c_im(inout[2 * iostride]);
 135                     tmp31 = c_re(W[1]);
 136                     tmp33 = c_im(W[1]);
 137                     tmp35 = (tmp31 * tmp32) + (tmp33 * tmp34);
 138                     tmp50 = (tmp31 * tmp34) - (tmp33 * tmp32);
 139                }
 140                {
 141                     fftw_real tmp37;
 142                     fftw_real tmp39;
 143                     fftw_real tmp36;
 144                     fftw_real tmp38;
 145                     ASSERT_ALIGNED_DOUBLE();
 146                     tmp37 = c_re(inout[5 * iostride]);
 147                     tmp39 = c_im(inout[5 * iostride]);
 148                     tmp36 = c_re(W[4]);
 149                     tmp38 = c_im(W[4]);
 150                     tmp40 = (tmp36 * tmp37) + (tmp38 * tmp39);
 151                     tmp51 = (tmp36 * tmp39) - (tmp38 * tmp37);
 152                }
 153                {
 154                     fftw_real tmp42;
 155                     fftw_real tmp44;
 156                     fftw_real tmp41;
 157                     fftw_real tmp43;
 158                     ASSERT_ALIGNED_DOUBLE();
 159                     tmp42 = c_re(inout[8 * iostride]);
 160                     tmp44 = c_im(inout[8 * iostride]);
 161                     tmp41 = c_re(W[7]);
 162                     tmp43 = c_im(W[7]);
 163                     tmp45 = (tmp41 * tmp42) + (tmp43 * tmp44);
 164                     tmp52 = (tmp41 * tmp44) - (tmp43 * tmp42);
 165                }
 166                tmp46 = tmp40 + tmp45;
 167                tmp53 = tmp51 + tmp52;
 168                {
 169                     fftw_real tmp76;
 170                     fftw_real tmp77;
 171                     fftw_real tmp73;
 172                     fftw_real tmp74;
 173                     ASSERT_ALIGNED_DOUBLE();
 174                     tmp47 = tmp35 + tmp46;
 175                     tmp76 = tmp35 - (K500000000 * tmp46);
 176                     tmp77 = K866025403 * (tmp52 - tmp51);
 177                     tmp78 = tmp76 - tmp77;
 178                     tmp89 = tmp76 + tmp77;
 179                     tmp54 = tmp50 + tmp53;
 180                     tmp73 = tmp50 - (K500000000 * tmp53);
 181                     tmp74 = K866025403 * (tmp40 - tmp45);
 182                     tmp75 = tmp73 - tmp74;
 183                     tmp90 = tmp74 + tmp73;
 184                }
 185           }
 186           {
 187                fftw_real tmp18;
 188                fftw_real tmp55;
 189                fftw_real tmp23;
 190                fftw_real tmp56;
 191                fftw_real tmp28;
 192                fftw_real tmp57;
 193                fftw_real tmp29;
 194                fftw_real tmp58;
 195                ASSERT_ALIGNED_DOUBLE();
 196                {
 197                     fftw_real tmp15;
 198                     fftw_real tmp17;
 199                     fftw_real tmp14;
 200                     fftw_real tmp16;
 201                     ASSERT_ALIGNED_DOUBLE();
 202                     tmp15 = c_re(inout[iostride]);
 203                     tmp17 = c_im(inout[iostride]);
 204                     tmp14 = c_re(W[0]);
 205                     tmp16 = c_im(W[0]);
 206                     tmp18 = (tmp14 * tmp15) + (tmp16 * tmp17);
 207                     tmp55 = (tmp14 * tmp17) - (tmp16 * tmp15);
 208                }
 209                {
 210                     fftw_real tmp20;
 211                     fftw_real tmp22;
 212                     fftw_real tmp19;
 213                     fftw_real tmp21;
 214                     ASSERT_ALIGNED_DOUBLE();
 215                     tmp20 = c_re(inout[4 * iostride]);
 216                     tmp22 = c_im(inout[4 * iostride]);
 217                     tmp19 = c_re(W[3]);
 218                     tmp21 = c_im(W[3]);
 219                     tmp23 = (tmp19 * tmp20) + (tmp21 * tmp22);
 220                     tmp56 = (tmp19 * tmp22) - (tmp21 * tmp20);
 221                }
 222                {
 223                     fftw_real tmp25;
 224                     fftw_real tmp27;
 225                     fftw_real tmp24;
 226                     fftw_real tmp26;
 227                     ASSERT_ALIGNED_DOUBLE();
 228                     tmp25 = c_re(inout[7 * iostride]);
 229                     tmp27 = c_im(inout[7 * iostride]);
 230                     tmp24 = c_re(W[6]);
 231                     tmp26 = c_im(W[6]);
 232                     tmp28 = (tmp24 * tmp25) + (tmp26 * tmp27);
 233                     tmp57 = (tmp24 * tmp27) - (tmp26 * tmp25);
 234                }
 235                tmp29 = tmp23 + tmp28;
 236                tmp58 = tmp56 + tmp57;
 237                {
 238                     fftw_real tmp66;
 239                     fftw_real tmp67;
 240                     fftw_real tmp69;
 241                     fftw_real tmp70;
 242                     ASSERT_ALIGNED_DOUBLE();
 243                     tmp30 = tmp18 + tmp29;
 244                     tmp66 = tmp18 - (K500000000 * tmp29);
 245                     tmp67 = K866025403 * (tmp57 - tmp56);
 246                     tmp68 = tmp66 - tmp67;
 247                     tmp86 = tmp66 + tmp67;
 248                     tmp59 = tmp55 + tmp58;
 249                     tmp69 = tmp55 - (K500000000 * tmp58);
 250                     tmp70 = K866025403 * (tmp23 - tmp28);
 251                     tmp71 = tmp69 - tmp70;
 252                     tmp87 = tmp70 + tmp69;
 253                }
 254           }
 255           {
 256                fftw_real tmp60;
 257                fftw_real tmp13;
 258                fftw_real tmp48;
 259                fftw_real tmp49;
 260                ASSERT_ALIGNED_DOUBLE();
 261                tmp60 = K866025403 * (tmp54 - tmp59);
 262                tmp13 = tmp1 + tmp12;
 263                tmp48 = tmp30 + tmp47;
 264                tmp49 = tmp13 - (K500000000 * tmp48);
 265                c_re(inout[0]) = tmp13 + tmp48;
 266                c_re(inout[3 * iostride]) = tmp49 + tmp60;
 267                c_re(inout[6 * iostride]) = tmp49 - tmp60;
 268           }
 269           {
 270                fftw_real tmp101;
 271                fftw_real tmp97;
 272                fftw_real tmp100;
 273                fftw_real tmp102;
 274                ASSERT_ALIGNED_DOUBLE();
 275                tmp101 = K866025403 * (tmp30 - tmp47);
 276                tmp97 = tmp59 + tmp54;
 277                tmp100 = tmp98 + tmp99;
 278                tmp102 = tmp100 - (K500000000 * tmp97);
 279                c_im(inout[0]) = tmp97 + tmp100;
 280                c_im(inout[6 * iostride]) = tmp102 - tmp101;
 281                c_im(inout[3 * iostride]) = tmp101 + tmp102;
 282           }
 283           {
 284                fftw_real tmp65;
 285                fftw_real tmp110;
 286                fftw_real tmp80;
 287                fftw_real tmp111;
 288                fftw_real tmp84;
 289                fftw_real tmp109;
 290                fftw_real tmp81;
 291                fftw_real tmp112;
 292                ASSERT_ALIGNED_DOUBLE();
 293                tmp65 = tmp61 - tmp64;
 294                tmp110 = tmp105 - tmp104;
 295                {
 296                     fftw_real tmp72;
 297                     fftw_real tmp79;
 298                     fftw_real tmp82;
 299                     fftw_real tmp83;
 300                     ASSERT_ALIGNED_DOUBLE();
 301                     tmp72 = (K173648177 * tmp68) - (K984807753 * tmp71);
 302                     tmp79 = (K342020143 * tmp75) + (K939692620 * tmp78);
 303                     tmp80 = tmp72 - tmp79;
 304                     tmp111 = K866025403 * (tmp72 + tmp79);
 305                     tmp82 = (K342020143 * tmp78) - (K939692620 * tmp75);
 306                     tmp83 = (K173648177 * tmp71) + (K984807753 * tmp68);
 307                     tmp84 = K866025403 * (tmp82 - tmp83);
 308                     tmp109 = tmp83 + tmp82;
 309                }
 310                c_re(inout[2 * iostride]) = tmp65 + tmp80;
 311                tmp81 = tmp65 - (K500000000 * tmp80);
 312                c_re(inout[8 * iostride]) = tmp81 - tmp84;
 313                c_re(inout[5 * iostride]) = tmp81 + tmp84;
 314                c_im(inout[2 * iostride]) = tmp109 + tmp110;
 315                tmp112 = tmp110 - (K500000000 * tmp109);
 316                c_im(inout[5 * iostride]) = tmp111 + tmp112;
 317                c_im(inout[8 * iostride]) = tmp112 - tmp111;
 318           }
 319           {
 320                fftw_real tmp85;
 321                fftw_real tmp106;
 322                fftw_real tmp92;
 323                fftw_real tmp107;
 324                fftw_real tmp96;
 325                fftw_real tmp103;
 326                fftw_real tmp93;
 327                fftw_real tmp108;
 328                ASSERT_ALIGNED_DOUBLE();
 329                tmp85 = tmp61 + tmp64;
 330                tmp106 = tmp104 + tmp105;
 331                {
 332                     fftw_real tmp88;
 333                     fftw_real tmp91;
 334                     fftw_real tmp94;
 335                     fftw_real tmp95;
 336                     ASSERT_ALIGNED_DOUBLE();
 337                     tmp88 = (K766044443 * tmp86) - (K642787609 * tmp87);
 338                     tmp91 = (K173648177 * tmp89) - (K984807753 * tmp90);
 339                     tmp92 = tmp88 + tmp91;
 340                     tmp107 = K866025403 * (tmp88 - tmp91);
 341                     tmp94 = (K173648177 * tmp90) + (K984807753 * tmp89);
 342                     tmp95 = (K766044443 * tmp87) + (K642787609 * tmp86);
 343                     tmp96 = K866025403 * (tmp94 - tmp95);
 344                     tmp103 = tmp95 + tmp94;
 345                }
 346                c_re(inout[iostride]) = tmp85 + tmp92;
 347                tmp93 = tmp85 - (K500000000 * tmp92);
 348                c_re(inout[7 * iostride]) = tmp93 - tmp96;
 349                c_re(inout[4 * iostride]) = tmp93 + tmp96;
 350                c_im(inout[iostride]) = tmp103 + tmp106;
 351                tmp108 = tmp106 - (K500000000 * tmp103);
 352                c_im(inout[4 * iostride]) = tmp107 + tmp108;
 353                c_im(inout[7 * iostride]) = tmp108 - tmp107;
 354           }
 355      }
 356 }
 357
 358 static const int twiddle_order[] =
 359 {1, 2, 3, 4, 5, 6, 7, 8};
 360 fftw_codelet_desc fftwi_twiddle_9_desc =
 361 {
 362      "fftwi_twiddle_9",
 363      (void (*)()) fftwi_twiddle_9,
 364      9,
 365      FFTW_BACKWARD,
 366      FFTW_TWIDDLE,
 367      209,
 368      8,
 369      twiddle_order,
 370 };