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