Subversion Repositories shark

Rev

Details | Last modification | View Log | RSS feed

Rev Author Line No. Line
2 pj 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:21 EDT 1999 */
22
 
23
#include <ports/fftw-int.h>
24
#include <ports/fftw.h>
25
 
26
/* Generated by: ./genfft -magic-alignment-check -magic-twiddle-load-all -magic-variables 4 -magic-loopi -twiddle 8 */
27
 
28
/*
29
 * This function contains 66 FP additions, 32 FP multiplications,
30
 * (or, 52 additions, 18 multiplications, 14 fused multiply/add),
31
 * 28 stack variables, and 32 memory accesses
32
 */
33
static const fftw_real K707106781 = FFTW_KONST(+0.707106781186547524400844362104849039284835938);
34
 
35
/*
36
 * Generator Id's :
37
 * $Id: ftw_8.c,v 1.1.1.1 2002-03-29 14:12:57 pj Exp $
38
 * $Id: ftw_8.c,v 1.1.1.1 2002-03-29 14:12:57 pj Exp $
39
 * $Id: ftw_8.c,v 1.1.1.1 2002-03-29 14:12:57 pj Exp $
40
 */
41
 
42
void fftw_twiddle_8(fftw_complex *A, const fftw_complex *W, int iostride, int m, int dist)
43
{
44
     int i;
45
     fftw_complex *inout;
46
     inout = A;
47
     for (i = m; i > 0; i = i - 1, inout = inout + dist, W = W + 7) {
48
          fftw_real tmp7;
49
          fftw_real tmp43;
50
          fftw_real tmp71;
51
          fftw_real tmp76;
52
          fftw_real tmp41;
53
          fftw_real tmp53;
54
          fftw_real tmp56;
55
          fftw_real tmp65;
56
          fftw_real tmp18;
57
          fftw_real tmp77;
58
          fftw_real tmp46;
59
          fftw_real tmp68;
60
          fftw_real tmp30;
61
          fftw_real tmp48;
62
          fftw_real tmp51;
63
          fftw_real tmp64;
64
          ASSERT_ALIGNED_DOUBLE();
65
          {
66
               fftw_real tmp1;
67
               fftw_real tmp70;
68
               fftw_real tmp6;
69
               fftw_real tmp69;
70
               ASSERT_ALIGNED_DOUBLE();
71
               tmp1 = c_re(inout[0]);
72
               tmp70 = c_im(inout[0]);
73
               {
74
                    fftw_real tmp3;
75
                    fftw_real tmp5;
76
                    fftw_real tmp2;
77
                    fftw_real tmp4;
78
                    ASSERT_ALIGNED_DOUBLE();
79
                    tmp3 = c_re(inout[4 * iostride]);
80
                    tmp5 = c_im(inout[4 * iostride]);
81
                    tmp2 = c_re(W[3]);
82
                    tmp4 = c_im(W[3]);
83
                    tmp6 = (tmp2 * tmp3) - (tmp4 * tmp5);
84
                    tmp69 = (tmp4 * tmp3) + (tmp2 * tmp5);
85
               }
86
               tmp7 = tmp1 + tmp6;
87
               tmp43 = tmp1 - tmp6;
88
               tmp71 = tmp69 + tmp70;
89
               tmp76 = tmp70 - tmp69;
90
          }
91
          {
92
               fftw_real tmp35;
93
               fftw_real tmp54;
94
               fftw_real tmp40;
95
               fftw_real tmp55;
96
               ASSERT_ALIGNED_DOUBLE();
97
               {
98
                    fftw_real tmp32;
99
                    fftw_real tmp34;
100
                    fftw_real tmp31;
101
                    fftw_real tmp33;
102
                    ASSERT_ALIGNED_DOUBLE();
103
                    tmp32 = c_re(inout[7 * iostride]);
104
                    tmp34 = c_im(inout[7 * iostride]);
105
                    tmp31 = c_re(W[6]);
106
                    tmp33 = c_im(W[6]);
107
                    tmp35 = (tmp31 * tmp32) - (tmp33 * tmp34);
108
                    tmp54 = (tmp33 * tmp32) + (tmp31 * tmp34);
109
               }
110
               {
111
                    fftw_real tmp37;
112
                    fftw_real tmp39;
113
                    fftw_real tmp36;
114
                    fftw_real tmp38;
115
                    ASSERT_ALIGNED_DOUBLE();
116
                    tmp37 = c_re(inout[3 * iostride]);
117
                    tmp39 = c_im(inout[3 * iostride]);
118
                    tmp36 = c_re(W[2]);
119
                    tmp38 = c_im(W[2]);
120
                    tmp40 = (tmp36 * tmp37) - (tmp38 * tmp39);
121
                    tmp55 = (tmp38 * tmp37) + (tmp36 * tmp39);
122
               }
123
               tmp41 = tmp35 + tmp40;
124
               tmp53 = tmp35 - tmp40;
125
               tmp56 = tmp54 - tmp55;
126
               tmp65 = tmp54 + tmp55;
127
          }
128
          {
129
               fftw_real tmp12;
130
               fftw_real tmp44;
131
               fftw_real tmp17;
132
               fftw_real tmp45;
133
               ASSERT_ALIGNED_DOUBLE();
134
               {
135
                    fftw_real tmp9;
136
                    fftw_real tmp11;
137
                    fftw_real tmp8;
138
                    fftw_real tmp10;
139
                    ASSERT_ALIGNED_DOUBLE();
140
                    tmp9 = c_re(inout[2 * iostride]);
141
                    tmp11 = c_im(inout[2 * iostride]);
142
                    tmp8 = c_re(W[1]);
143
                    tmp10 = c_im(W[1]);
144
                    tmp12 = (tmp8 * tmp9) - (tmp10 * tmp11);
145
                    tmp44 = (tmp10 * tmp9) + (tmp8 * tmp11);
146
               }
147
               {
148
                    fftw_real tmp14;
149
                    fftw_real tmp16;
150
                    fftw_real tmp13;
151
                    fftw_real tmp15;
152
                    ASSERT_ALIGNED_DOUBLE();
153
                    tmp14 = c_re(inout[6 * iostride]);
154
                    tmp16 = c_im(inout[6 * iostride]);
155
                    tmp13 = c_re(W[5]);
156
                    tmp15 = c_im(W[5]);
157
                    tmp17 = (tmp13 * tmp14) - (tmp15 * tmp16);
158
                    tmp45 = (tmp15 * tmp14) + (tmp13 * tmp16);
159
               }
160
               tmp18 = tmp12 + tmp17;
161
               tmp77 = tmp12 - tmp17;
162
               tmp46 = tmp44 - tmp45;
163
               tmp68 = tmp44 + tmp45;
164
          }
165
          {
166
               fftw_real tmp24;
167
               fftw_real tmp49;
168
               fftw_real tmp29;
169
               fftw_real tmp50;
170
               ASSERT_ALIGNED_DOUBLE();
171
               {
172
                    fftw_real tmp21;
173
                    fftw_real tmp23;
174
                    fftw_real tmp20;
175
                    fftw_real tmp22;
176
                    ASSERT_ALIGNED_DOUBLE();
177
                    tmp21 = c_re(inout[iostride]);
178
                    tmp23 = c_im(inout[iostride]);
179
                    tmp20 = c_re(W[0]);
180
                    tmp22 = c_im(W[0]);
181
                    tmp24 = (tmp20 * tmp21) - (tmp22 * tmp23);
182
                    tmp49 = (tmp22 * tmp21) + (tmp20 * tmp23);
183
               }
184
               {
185
                    fftw_real tmp26;
186
                    fftw_real tmp28;
187
                    fftw_real tmp25;
188
                    fftw_real tmp27;
189
                    ASSERT_ALIGNED_DOUBLE();
190
                    tmp26 = c_re(inout[5 * iostride]);
191
                    tmp28 = c_im(inout[5 * iostride]);
192
                    tmp25 = c_re(W[4]);
193
                    tmp27 = c_im(W[4]);
194
                    tmp29 = (tmp25 * tmp26) - (tmp27 * tmp28);
195
                    tmp50 = (tmp27 * tmp26) + (tmp25 * tmp28);
196
               }
197
               tmp30 = tmp24 + tmp29;
198
               tmp48 = tmp24 - tmp29;
199
               tmp51 = tmp49 - tmp50;
200
               tmp64 = tmp49 + tmp50;
201
          }
202
          {
203
               fftw_real tmp19;
204
               fftw_real tmp42;
205
               fftw_real tmp63;
206
               fftw_real tmp66;
207
               ASSERT_ALIGNED_DOUBLE();
208
               tmp19 = tmp7 + tmp18;
209
               tmp42 = tmp30 + tmp41;
210
               c_re(inout[4 * iostride]) = tmp19 - tmp42;
211
               c_re(inout[0]) = tmp19 + tmp42;
212
               {
213
                    fftw_real tmp73;
214
                    fftw_real tmp74;
215
                    fftw_real tmp67;
216
                    fftw_real tmp72;
217
                    ASSERT_ALIGNED_DOUBLE();
218
                    tmp73 = tmp41 - tmp30;
219
                    tmp74 = tmp71 - tmp68;
220
                    c_im(inout[2 * iostride]) = tmp73 + tmp74;
221
                    c_im(inout[6 * iostride]) = tmp74 - tmp73;
222
                    tmp67 = tmp64 + tmp65;
223
                    tmp72 = tmp68 + tmp71;
224
                    c_im(inout[0]) = tmp67 + tmp72;
225
                    c_im(inout[4 * iostride]) = tmp72 - tmp67;
226
               }
227
               tmp63 = tmp7 - tmp18;
228
               tmp66 = tmp64 - tmp65;
229
               c_re(inout[6 * iostride]) = tmp63 - tmp66;
230
               c_re(inout[2 * iostride]) = tmp63 + tmp66;
231
               {
232
                    fftw_real tmp59;
233
                    fftw_real tmp78;
234
                    fftw_real tmp62;
235
                    fftw_real tmp75;
236
                    fftw_real tmp60;
237
                    fftw_real tmp61;
238
                    ASSERT_ALIGNED_DOUBLE();
239
                    tmp59 = tmp43 - tmp46;
240
                    tmp78 = tmp76 - tmp77;
241
                    tmp60 = tmp51 - tmp48;
242
                    tmp61 = tmp53 + tmp56;
243
                    tmp62 = K707106781 * (tmp60 - tmp61);
244
                    tmp75 = K707106781 * (tmp60 + tmp61);
245
                    c_re(inout[7 * iostride]) = tmp59 - tmp62;
246
                    c_re(inout[3 * iostride]) = tmp59 + tmp62;
247
                    c_im(inout[iostride]) = tmp75 + tmp78;
248
                    c_im(inout[5 * iostride]) = tmp78 - tmp75;
249
               }
250
               {
251
                    fftw_real tmp47;
252
                    fftw_real tmp80;
253
                    fftw_real tmp58;
254
                    fftw_real tmp79;
255
                    fftw_real tmp52;
256
                    fftw_real tmp57;
257
                    ASSERT_ALIGNED_DOUBLE();
258
                    tmp47 = tmp43 + tmp46;
259
                    tmp80 = tmp77 + tmp76;
260
                    tmp52 = tmp48 + tmp51;
261
                    tmp57 = tmp53 - tmp56;
262
                    tmp58 = K707106781 * (tmp52 + tmp57);
263
                    tmp79 = K707106781 * (tmp57 - tmp52);
264
                    c_re(inout[5 * iostride]) = tmp47 - tmp58;
265
                    c_re(inout[iostride]) = tmp47 + tmp58;
266
                    c_im(inout[3 * iostride]) = tmp79 + tmp80;
267
                    c_im(inout[7 * iostride]) = tmp80 - tmp79;
268
               }
269
          }
270
     }
271
}
272
 
273
static const int twiddle_order[] =
274
{1, 2, 3, 4, 5, 6, 7};
275
fftw_codelet_desc fftw_twiddle_8_desc =
276
{
277
     "fftw_twiddle_8",
278
     (void (*)()) fftw_twiddle_8,
279
     8,
280
     FFTW_FORWARD,
281
     FFTW_TWIDDLE,
282
     176,
283
     7,
284
     twiddle_order,
285
};