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:24 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 16 */
27
 
28
/*
29
 * This function contains 174 FP additions, 84 FP multiplications,
30
 * (or, 136 additions, 46 multiplications, 38 fused multiply/add),
31
 * 50 stack variables, and 64 memory accesses
32
 */
33
static const fftw_real K382683432 = FFTW_KONST(+0.382683432365089771728459984030398866761344562);
34
static const fftw_real K923879532 = FFTW_KONST(+0.923879532511286756128183189396788286822416626);
35
static const fftw_real K707106781 = FFTW_KONST(+0.707106781186547524400844362104849039284835938);
36
 
37
/*
38
 * Generator Id's :
39
 * $Id: ftw_16.c,v 1.1.1.1 2002-03-29 14:12:56 pj Exp $
40
 * $Id: ftw_16.c,v 1.1.1.1 2002-03-29 14:12:56 pj Exp $
41
 * $Id: ftw_16.c,v 1.1.1.1 2002-03-29 14:12:56 pj Exp $
42
 */
43
 
44
void fftw_twiddle_16(fftw_complex *A, const fftw_complex *W, int iostride, int m, int dist)
45
{
46
     int i;
47
     fftw_complex *inout;
48
     inout = A;
49
     for (i = m; i > 0; i = i - 1, inout = inout + dist, W = W + 15) {
50
          fftw_real tmp7;
51
          fftw_real tmp91;
52
          fftw_real tmp180;
53
          fftw_real tmp193;
54
          fftw_real tmp18;
55
          fftw_real tmp194;
56
          fftw_real tmp94;
57
          fftw_real tmp177;
58
          fftw_real tmp77;
59
          fftw_real tmp88;
60
          fftw_real tmp161;
61
          fftw_real tmp128;
62
          fftw_real tmp144;
63
          fftw_real tmp162;
64
          fftw_real tmp163;
65
          fftw_real tmp164;
66
          fftw_real tmp123;
67
          fftw_real tmp143;
68
          fftw_real tmp30;
69
          fftw_real tmp152;
70
          fftw_real tmp100;
71
          fftw_real tmp136;
72
          fftw_real tmp41;
73
          fftw_real tmp153;
74
          fftw_real tmp105;
75
          fftw_real tmp137;
76
          fftw_real tmp54;
77
          fftw_real tmp65;
78
          fftw_real tmp156;
79
          fftw_real tmp117;
80
          fftw_real tmp141;
81
          fftw_real tmp157;
82
          fftw_real tmp158;
83
          fftw_real tmp159;
84
          fftw_real tmp112;
85
          fftw_real tmp140;
86
          ASSERT_ALIGNED_DOUBLE();
87
          {
88
               fftw_real tmp1;
89
               fftw_real tmp179;
90
               fftw_real tmp6;
91
               fftw_real tmp178;
92
               ASSERT_ALIGNED_DOUBLE();
93
               tmp1 = c_re(inout[0]);
94
               tmp179 = c_im(inout[0]);
95
               {
96
                    fftw_real tmp3;
97
                    fftw_real tmp5;
98
                    fftw_real tmp2;
99
                    fftw_real tmp4;
100
                    ASSERT_ALIGNED_DOUBLE();
101
                    tmp3 = c_re(inout[8 * iostride]);
102
                    tmp5 = c_im(inout[8 * iostride]);
103
                    tmp2 = c_re(W[7]);
104
                    tmp4 = c_im(W[7]);
105
                    tmp6 = (tmp2 * tmp3) - (tmp4 * tmp5);
106
                    tmp178 = (tmp4 * tmp3) + (tmp2 * tmp5);
107
               }
108
               tmp7 = tmp1 + tmp6;
109
               tmp91 = tmp1 - tmp6;
110
               tmp180 = tmp178 + tmp179;
111
               tmp193 = tmp179 - tmp178;
112
          }
113
          {
114
               fftw_real tmp12;
115
               fftw_real tmp92;
116
               fftw_real tmp17;
117
               fftw_real tmp93;
118
               ASSERT_ALIGNED_DOUBLE();
119
               {
120
                    fftw_real tmp9;
121
                    fftw_real tmp11;
122
                    fftw_real tmp8;
123
                    fftw_real tmp10;
124
                    ASSERT_ALIGNED_DOUBLE();
125
                    tmp9 = c_re(inout[4 * iostride]);
126
                    tmp11 = c_im(inout[4 * iostride]);
127
                    tmp8 = c_re(W[3]);
128
                    tmp10 = c_im(W[3]);
129
                    tmp12 = (tmp8 * tmp9) - (tmp10 * tmp11);
130
                    tmp92 = (tmp10 * tmp9) + (tmp8 * tmp11);
131
               }
132
               {
133
                    fftw_real tmp14;
134
                    fftw_real tmp16;
135
                    fftw_real tmp13;
136
                    fftw_real tmp15;
137
                    ASSERT_ALIGNED_DOUBLE();
138
                    tmp14 = c_re(inout[12 * iostride]);
139
                    tmp16 = c_im(inout[12 * iostride]);
140
                    tmp13 = c_re(W[11]);
141
                    tmp15 = c_im(W[11]);
142
                    tmp17 = (tmp13 * tmp14) - (tmp15 * tmp16);
143
                    tmp93 = (tmp15 * tmp14) + (tmp13 * tmp16);
144
               }
145
               tmp18 = tmp12 + tmp17;
146
               tmp194 = tmp12 - tmp17;
147
               tmp94 = tmp92 - tmp93;
148
               tmp177 = tmp92 + tmp93;
149
          }
150
          {
151
               fftw_real tmp71;
152
               fftw_real tmp124;
153
               fftw_real tmp87;
154
               fftw_real tmp121;
155
               fftw_real tmp76;
156
               fftw_real tmp125;
157
               fftw_real tmp82;
158
               fftw_real tmp120;
159
               ASSERT_ALIGNED_DOUBLE();
160
               {
161
                    fftw_real tmp68;
162
                    fftw_real tmp70;
163
                    fftw_real tmp67;
164
                    fftw_real tmp69;
165
                    ASSERT_ALIGNED_DOUBLE();
166
                    tmp68 = c_re(inout[15 * iostride]);
167
                    tmp70 = c_im(inout[15 * iostride]);
168
                    tmp67 = c_re(W[14]);
169
                    tmp69 = c_im(W[14]);
170
                    tmp71 = (tmp67 * tmp68) - (tmp69 * tmp70);
171
                    tmp124 = (tmp69 * tmp68) + (tmp67 * tmp70);
172
               }
173
               {
174
                    fftw_real tmp84;
175
                    fftw_real tmp86;
176
                    fftw_real tmp83;
177
                    fftw_real tmp85;
178
                    ASSERT_ALIGNED_DOUBLE();
179
                    tmp84 = c_re(inout[11 * iostride]);
180
                    tmp86 = c_im(inout[11 * iostride]);
181
                    tmp83 = c_re(W[10]);
182
                    tmp85 = c_im(W[10]);
183
                    tmp87 = (tmp83 * tmp84) - (tmp85 * tmp86);
184
                    tmp121 = (tmp85 * tmp84) + (tmp83 * tmp86);
185
               }
186
               {
187
                    fftw_real tmp73;
188
                    fftw_real tmp75;
189
                    fftw_real tmp72;
190
                    fftw_real tmp74;
191
                    ASSERT_ALIGNED_DOUBLE();
192
                    tmp73 = c_re(inout[7 * iostride]);
193
                    tmp75 = c_im(inout[7 * iostride]);
194
                    tmp72 = c_re(W[6]);
195
                    tmp74 = c_im(W[6]);
196
                    tmp76 = (tmp72 * tmp73) - (tmp74 * tmp75);
197
                    tmp125 = (tmp74 * tmp73) + (tmp72 * tmp75);
198
               }
199
               {
200
                    fftw_real tmp79;
201
                    fftw_real tmp81;
202
                    fftw_real tmp78;
203
                    fftw_real tmp80;
204
                    ASSERT_ALIGNED_DOUBLE();
205
                    tmp79 = c_re(inout[3 * iostride]);
206
                    tmp81 = c_im(inout[3 * iostride]);
207
                    tmp78 = c_re(W[2]);
208
                    tmp80 = c_im(W[2]);
209
                    tmp82 = (tmp78 * tmp79) - (tmp80 * tmp81);
210
                    tmp120 = (tmp80 * tmp79) + (tmp78 * tmp81);
211
               }
212
               {
213
                    fftw_real tmp126;
214
                    fftw_real tmp127;
215
                    fftw_real tmp119;
216
                    fftw_real tmp122;
217
                    ASSERT_ALIGNED_DOUBLE();
218
                    tmp77 = tmp71 + tmp76;
219
                    tmp88 = tmp82 + tmp87;
220
                    tmp161 = tmp77 - tmp88;
221
                    tmp126 = tmp124 - tmp125;
222
                    tmp127 = tmp82 - tmp87;
223
                    tmp128 = tmp126 + tmp127;
224
                    tmp144 = tmp126 - tmp127;
225
                    tmp162 = tmp124 + tmp125;
226
                    tmp163 = tmp120 + tmp121;
227
                    tmp164 = tmp162 - tmp163;
228
                    tmp119 = tmp71 - tmp76;
229
                    tmp122 = tmp120 - tmp121;
230
                    tmp123 = tmp119 - tmp122;
231
                    tmp143 = tmp119 + tmp122;
232
               }
233
          }
234
          {
235
               fftw_real tmp24;
236
               fftw_real tmp96;
237
               fftw_real tmp29;
238
               fftw_real tmp97;
239
               fftw_real tmp98;
240
               fftw_real tmp99;
241
               ASSERT_ALIGNED_DOUBLE();
242
               {
243
                    fftw_real tmp21;
244
                    fftw_real tmp23;
245
                    fftw_real tmp20;
246
                    fftw_real tmp22;
247
                    ASSERT_ALIGNED_DOUBLE();
248
                    tmp21 = c_re(inout[2 * iostride]);
249
                    tmp23 = c_im(inout[2 * iostride]);
250
                    tmp20 = c_re(W[1]);
251
                    tmp22 = c_im(W[1]);
252
                    tmp24 = (tmp20 * tmp21) - (tmp22 * tmp23);
253
                    tmp96 = (tmp22 * tmp21) + (tmp20 * tmp23);
254
               }
255
               {
256
                    fftw_real tmp26;
257
                    fftw_real tmp28;
258
                    fftw_real tmp25;
259
                    fftw_real tmp27;
260
                    ASSERT_ALIGNED_DOUBLE();
261
                    tmp26 = c_re(inout[10 * iostride]);
262
                    tmp28 = c_im(inout[10 * iostride]);
263
                    tmp25 = c_re(W[9]);
264
                    tmp27 = c_im(W[9]);
265
                    tmp29 = (tmp25 * tmp26) - (tmp27 * tmp28);
266
                    tmp97 = (tmp27 * tmp26) + (tmp25 * tmp28);
267
               }
268
               tmp30 = tmp24 + tmp29;
269
               tmp152 = tmp96 + tmp97;
270
               tmp98 = tmp96 - tmp97;
271
               tmp99 = tmp24 - tmp29;
272
               tmp100 = tmp98 - tmp99;
273
               tmp136 = tmp99 + tmp98;
274
          }
275
          {
276
               fftw_real tmp35;
277
               fftw_real tmp102;
278
               fftw_real tmp40;
279
               fftw_real tmp103;
280
               fftw_real tmp101;
281
               fftw_real tmp104;
282
               ASSERT_ALIGNED_DOUBLE();
283
               {
284
                    fftw_real tmp32;
285
                    fftw_real tmp34;
286
                    fftw_real tmp31;
287
                    fftw_real tmp33;
288
                    ASSERT_ALIGNED_DOUBLE();
289
                    tmp32 = c_re(inout[14 * iostride]);
290
                    tmp34 = c_im(inout[14 * iostride]);
291
                    tmp31 = c_re(W[13]);
292
                    tmp33 = c_im(W[13]);
293
                    tmp35 = (tmp31 * tmp32) - (tmp33 * tmp34);
294
                    tmp102 = (tmp33 * tmp32) + (tmp31 * tmp34);
295
               }
296
               {
297
                    fftw_real tmp37;
298
                    fftw_real tmp39;
299
                    fftw_real tmp36;
300
                    fftw_real tmp38;
301
                    ASSERT_ALIGNED_DOUBLE();
302
                    tmp37 = c_re(inout[6 * iostride]);
303
                    tmp39 = c_im(inout[6 * iostride]);
304
                    tmp36 = c_re(W[5]);
305
                    tmp38 = c_im(W[5]);
306
                    tmp40 = (tmp36 * tmp37) - (tmp38 * tmp39);
307
                    tmp103 = (tmp38 * tmp37) + (tmp36 * tmp39);
308
               }
309
               tmp41 = tmp35 + tmp40;
310
               tmp153 = tmp102 + tmp103;
311
               tmp101 = tmp35 - tmp40;
312
               tmp104 = tmp102 - tmp103;
313
               tmp105 = tmp101 + tmp104;
314
               tmp137 = tmp101 - tmp104;
315
          }
316
          {
317
               fftw_real tmp48;
318
               fftw_real tmp108;
319
               fftw_real tmp64;
320
               fftw_real tmp115;
321
               fftw_real tmp53;
322
               fftw_real tmp109;
323
               fftw_real tmp59;
324
               fftw_real tmp114;
325
               ASSERT_ALIGNED_DOUBLE();
326
               {
327
                    fftw_real tmp45;
328
                    fftw_real tmp47;
329
                    fftw_real tmp44;
330
                    fftw_real tmp46;
331
                    ASSERT_ALIGNED_DOUBLE();
332
                    tmp45 = c_re(inout[iostride]);
333
                    tmp47 = c_im(inout[iostride]);
334
                    tmp44 = c_re(W[0]);
335
                    tmp46 = c_im(W[0]);
336
                    tmp48 = (tmp44 * tmp45) - (tmp46 * tmp47);
337
                    tmp108 = (tmp46 * tmp45) + (tmp44 * tmp47);
338
               }
339
               {
340
                    fftw_real tmp61;
341
                    fftw_real tmp63;
342
                    fftw_real tmp60;
343
                    fftw_real tmp62;
344
                    ASSERT_ALIGNED_DOUBLE();
345
                    tmp61 = c_re(inout[13 * iostride]);
346
                    tmp63 = c_im(inout[13 * iostride]);
347
                    tmp60 = c_re(W[12]);
348
                    tmp62 = c_im(W[12]);
349
                    tmp64 = (tmp60 * tmp61) - (tmp62 * tmp63);
350
                    tmp115 = (tmp62 * tmp61) + (tmp60 * tmp63);
351
               }
352
               {
353
                    fftw_real tmp50;
354
                    fftw_real tmp52;
355
                    fftw_real tmp49;
356
                    fftw_real tmp51;
357
                    ASSERT_ALIGNED_DOUBLE();
358
                    tmp50 = c_re(inout[9 * iostride]);
359
                    tmp52 = c_im(inout[9 * iostride]);
360
                    tmp49 = c_re(W[8]);
361
                    tmp51 = c_im(W[8]);
362
                    tmp53 = (tmp49 * tmp50) - (tmp51 * tmp52);
363
                    tmp109 = (tmp51 * tmp50) + (tmp49 * tmp52);
364
               }
365
               {
366
                    fftw_real tmp56;
367
                    fftw_real tmp58;
368
                    fftw_real tmp55;
369
                    fftw_real tmp57;
370
                    ASSERT_ALIGNED_DOUBLE();
371
                    tmp56 = c_re(inout[5 * iostride]);
372
                    tmp58 = c_im(inout[5 * iostride]);
373
                    tmp55 = c_re(W[4]);
374
                    tmp57 = c_im(W[4]);
375
                    tmp59 = (tmp55 * tmp56) - (tmp57 * tmp58);
376
                    tmp114 = (tmp57 * tmp56) + (tmp55 * tmp58);
377
               }
378
               {
379
                    fftw_real tmp113;
380
                    fftw_real tmp116;
381
                    fftw_real tmp110;
382
                    fftw_real tmp111;
383
                    ASSERT_ALIGNED_DOUBLE();
384
                    tmp54 = tmp48 + tmp53;
385
                    tmp65 = tmp59 + tmp64;
386
                    tmp156 = tmp54 - tmp65;
387
                    tmp113 = tmp48 - tmp53;
388
                    tmp116 = tmp114 - tmp115;
389
                    tmp117 = tmp113 - tmp116;
390
                    tmp141 = tmp113 + tmp116;
391
                    tmp157 = tmp108 + tmp109;
392
                    tmp158 = tmp114 + tmp115;
393
                    tmp159 = tmp157 - tmp158;
394
                    tmp110 = tmp108 - tmp109;
395
                    tmp111 = tmp59 - tmp64;
396
                    tmp112 = tmp110 + tmp111;
397
                    tmp140 = tmp110 - tmp111;
398
               }
399
          }
400
          {
401
               fftw_real tmp107;
402
               fftw_real tmp131;
403
               fftw_real tmp202;
404
               fftw_real tmp204;
405
               fftw_real tmp130;
406
               fftw_real tmp203;
407
               fftw_real tmp134;
408
               fftw_real tmp199;
409
               ASSERT_ALIGNED_DOUBLE();
410
               {
411
                    fftw_real tmp95;
412
                    fftw_real tmp106;
413
                    fftw_real tmp200;
414
                    fftw_real tmp201;
415
                    ASSERT_ALIGNED_DOUBLE();
416
                    tmp95 = tmp91 - tmp94;
417
                    tmp106 = K707106781 * (tmp100 - tmp105);
418
                    tmp107 = tmp95 + tmp106;
419
                    tmp131 = tmp95 - tmp106;
420
                    tmp200 = K707106781 * (tmp137 - tmp136);
421
                    tmp201 = tmp194 + tmp193;
422
                    tmp202 = tmp200 + tmp201;
423
                    tmp204 = tmp201 - tmp200;
424
               }
425
               {
426
                    fftw_real tmp118;
427
                    fftw_real tmp129;
428
                    fftw_real tmp132;
429
                    fftw_real tmp133;
430
                    ASSERT_ALIGNED_DOUBLE();
431
                    tmp118 = (K923879532 * tmp112) + (K382683432 * tmp117);
432
                    tmp129 = (K382683432 * tmp123) - (K923879532 * tmp128);
433
                    tmp130 = tmp118 + tmp129;
434
                    tmp203 = tmp129 - tmp118;
435
                    tmp132 = (K382683432 * tmp112) - (K923879532 * tmp117);
436
                    tmp133 = (K382683432 * tmp128) + (K923879532 * tmp123);
437
                    tmp134 = tmp132 - tmp133;
438
                    tmp199 = tmp132 + tmp133;
439
               }
440
               c_re(inout[11 * iostride]) = tmp107 - tmp130;
441
               c_re(inout[3 * iostride]) = tmp107 + tmp130;
442
               c_re(inout[15 * iostride]) = tmp131 - tmp134;
443
               c_re(inout[7 * iostride]) = tmp131 + tmp134;
444
               c_im(inout[3 * iostride]) = tmp199 + tmp202;
445
               c_im(inout[11 * iostride]) = tmp202 - tmp199;
446
               c_im(inout[7 * iostride]) = tmp203 + tmp204;
447
               c_im(inout[15 * iostride]) = tmp204 - tmp203;
448
          }
449
          {
450
               fftw_real tmp139;
451
               fftw_real tmp147;
452
               fftw_real tmp196;
453
               fftw_real tmp198;
454
               fftw_real tmp146;
455
               fftw_real tmp197;
456
               fftw_real tmp150;
457
               fftw_real tmp191;
458
               ASSERT_ALIGNED_DOUBLE();
459
               {
460
                    fftw_real tmp135;
461
                    fftw_real tmp138;
462
                    fftw_real tmp192;
463
                    fftw_real tmp195;
464
                    ASSERT_ALIGNED_DOUBLE();
465
                    tmp135 = tmp91 + tmp94;
466
                    tmp138 = K707106781 * (tmp136 + tmp137);
467
                    tmp139 = tmp135 + tmp138;
468
                    tmp147 = tmp135 - tmp138;
469
                    tmp192 = K707106781 * (tmp100 + tmp105);
470
                    tmp195 = tmp193 - tmp194;
471
                    tmp196 = tmp192 + tmp195;
472
                    tmp198 = tmp195 - tmp192;
473
               }
474
               {
475
                    fftw_real tmp142;
476
                    fftw_real tmp145;
477
                    fftw_real tmp148;
478
                    fftw_real tmp149;
479
                    ASSERT_ALIGNED_DOUBLE();
480
                    tmp142 = (K382683432 * tmp140) + (K923879532 * tmp141);
481
                    tmp145 = (K923879532 * tmp143) - (K382683432 * tmp144);
482
                    tmp146 = tmp142 + tmp145;
483
                    tmp197 = tmp145 - tmp142;
484
                    tmp148 = (K923879532 * tmp140) - (K382683432 * tmp141);
485
                    tmp149 = (K923879532 * tmp144) + (K382683432 * tmp143);
486
                    tmp150 = tmp148 - tmp149;
487
                    tmp191 = tmp148 + tmp149;
488
               }
489
               c_re(inout[9 * iostride]) = tmp139 - tmp146;
490
               c_re(inout[iostride]) = tmp139 + tmp146;
491
               c_re(inout[13 * iostride]) = tmp147 - tmp150;
492
               c_re(inout[5 * iostride]) = tmp147 + tmp150;
493
               c_im(inout[iostride]) = tmp191 + tmp196;
494
               c_im(inout[9 * iostride]) = tmp196 - tmp191;
495
               c_im(inout[5 * iostride]) = tmp197 + tmp198;
496
               c_im(inout[13 * iostride]) = tmp198 - tmp197;
497
          }
498
          {
499
               fftw_real tmp155;
500
               fftw_real tmp167;
501
               fftw_real tmp188;
502
               fftw_real tmp190;
503
               fftw_real tmp166;
504
               fftw_real tmp189;
505
               fftw_real tmp170;
506
               fftw_real tmp185;
507
               ASSERT_ALIGNED_DOUBLE();
508
               {
509
                    fftw_real tmp151;
510
                    fftw_real tmp154;
511
                    fftw_real tmp186;
512
                    fftw_real tmp187;
513
                    ASSERT_ALIGNED_DOUBLE();
514
                    tmp151 = tmp7 - tmp18;
515
                    tmp154 = tmp152 - tmp153;
516
                    tmp155 = tmp151 + tmp154;
517
                    tmp167 = tmp151 - tmp154;
518
                    tmp186 = tmp41 - tmp30;
519
                    tmp187 = tmp180 - tmp177;
520
                    tmp188 = tmp186 + tmp187;
521
                    tmp190 = tmp187 - tmp186;
522
               }
523
               {
524
                    fftw_real tmp160;
525
                    fftw_real tmp165;
526
                    fftw_real tmp168;
527
                    fftw_real tmp169;
528
                    ASSERT_ALIGNED_DOUBLE();
529
                    tmp160 = tmp156 + tmp159;
530
                    tmp165 = tmp161 - tmp164;
531
                    tmp166 = K707106781 * (tmp160 + tmp165);
532
                    tmp189 = K707106781 * (tmp165 - tmp160);
533
                    tmp168 = tmp159 - tmp156;
534
                    tmp169 = tmp161 + tmp164;
535
                    tmp170 = K707106781 * (tmp168 - tmp169);
536
                    tmp185 = K707106781 * (tmp168 + tmp169);
537
               }
538
               c_re(inout[10 * iostride]) = tmp155 - tmp166;
539
               c_re(inout[2 * iostride]) = tmp155 + tmp166;
540
               c_re(inout[14 * iostride]) = tmp167 - tmp170;
541
               c_re(inout[6 * iostride]) = tmp167 + tmp170;
542
               c_im(inout[2 * iostride]) = tmp185 + tmp188;
543
               c_im(inout[10 * iostride]) = tmp188 - tmp185;
544
               c_im(inout[6 * iostride]) = tmp189 + tmp190;
545
               c_im(inout[14 * iostride]) = tmp190 - tmp189;
546
          }
547
          {
548
               fftw_real tmp43;
549
               fftw_real tmp171;
550
               fftw_real tmp182;
551
               fftw_real tmp184;
552
               fftw_real tmp90;
553
               fftw_real tmp183;
554
               fftw_real tmp174;
555
               fftw_real tmp175;
556
               ASSERT_ALIGNED_DOUBLE();
557
               {
558
                    fftw_real tmp19;
559
                    fftw_real tmp42;
560
                    fftw_real tmp176;
561
                    fftw_real tmp181;
562
                    ASSERT_ALIGNED_DOUBLE();
563
                    tmp19 = tmp7 + tmp18;
564
                    tmp42 = tmp30 + tmp41;
565
                    tmp43 = tmp19 + tmp42;
566
                    tmp171 = tmp19 - tmp42;
567
                    tmp176 = tmp152 + tmp153;
568
                    tmp181 = tmp177 + tmp180;
569
                    tmp182 = tmp176 + tmp181;
570
                    tmp184 = tmp181 - tmp176;
571
               }
572
               {
573
                    fftw_real tmp66;
574
                    fftw_real tmp89;
575
                    fftw_real tmp172;
576
                    fftw_real tmp173;
577
                    ASSERT_ALIGNED_DOUBLE();
578
                    tmp66 = tmp54 + tmp65;
579
                    tmp89 = tmp77 + tmp88;
580
                    tmp90 = tmp66 + tmp89;
581
                    tmp183 = tmp89 - tmp66;
582
                    tmp172 = tmp157 + tmp158;
583
                    tmp173 = tmp162 + tmp163;
584
                    tmp174 = tmp172 - tmp173;
585
                    tmp175 = tmp172 + tmp173;
586
               }
587
               c_re(inout[8 * iostride]) = tmp43 - tmp90;
588
               c_re(inout[0]) = tmp43 + tmp90;
589
               c_re(inout[12 * iostride]) = tmp171 - tmp174;
590
               c_re(inout[4 * iostride]) = tmp171 + tmp174;
591
               c_im(inout[0]) = tmp175 + tmp182;
592
               c_im(inout[8 * iostride]) = tmp182 - tmp175;
593
               c_im(inout[4 * iostride]) = tmp183 + tmp184;
594
               c_im(inout[12 * iostride]) = tmp184 - tmp183;
595
          }
596
     }
597
}
598
 
599
static const int twiddle_order[] =
600
{1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15};
601
fftw_codelet_desc fftw_twiddle_16_desc =
602
{
603
     "fftw_twiddle_16",
604
     (void (*)()) fftw_twiddle_16,
605
     16,
606
     FFTW_FORWARD,
607
     FFTW_TWIDDLE,
608
     352,
609
     15,
610
     twiddle_order,
611
};