Subversion Repositories shark

Rev

Rev 3 | Details | Compare with Previous | 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:57 EDT 1999 */
22
 
107 pj 23
#include <fftw-int.h>
24
#include <fftw.h>
2 pj 25
 
26
/* Generated by: ./genfft -magic-alignment-check -magic-twiddle-load-all -magic-variables 4 -magic-loopi -notwiddleinv 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 K195090322 = FFTW_KONST(+0.195090322016128267848284868477022240927691618);
34
static const fftw_real K980785280 = FFTW_KONST(+0.980785280403230449126182236134239036973933731);
35
static const fftw_real K831469612 = FFTW_KONST(+0.831469612302545237078788377617905756738560812);
36
static const fftw_real K555570233 = FFTW_KONST(+0.555570233019602224742830813948532874374937191);
37
static const fftw_real K382683432 = FFTW_KONST(+0.382683432365089771728459984030398866761344562);
38
static const fftw_real K923879532 = FFTW_KONST(+0.923879532511286756128183189396788286822416626);
39
static const fftw_real K707106781 = FFTW_KONST(+0.707106781186547524400844362104849039284835938);
40
 
41
/*
42
 * Generator Id's :
107 pj 43
 * $Id: fni_32.c,v 1.2 2003-03-24 11:14:51 pj Exp $
44
 * $Id: fni_32.c,v 1.2 2003-03-24 11:14:51 pj Exp $
45
 * $Id: fni_32.c,v 1.2 2003-03-24 11:14:51 pj Exp $
2 pj 46
 */
47
 
48
void fftwi_no_twiddle_32(const fftw_complex *input, fftw_complex *output, int istride, int ostride)
49
{
50
     fftw_real tmp7;
51
     fftw_real tmp339;
52
     fftw_real tmp70;
53
     fftw_real tmp313;
54
     fftw_real tmp97;
55
     fftw_real tmp215;
56
     fftw_real tmp179;
57
     fftw_real tmp241;
58
     fftw_real tmp14;
59
     fftw_real tmp314;
60
     fftw_real tmp77;
61
     fftw_real tmp340;
62
     fftw_real tmp182;
63
     fftw_real tmp216;
64
     fftw_real tmp104;
65
     fftw_real tmp242;
66
     fftw_real tmp153;
67
     fftw_real tmp236;
68
     fftw_real tmp53;
69
     fftw_real tmp60;
70
     fftw_real tmp287;
71
     fftw_real tmp336;
72
     fftw_real tmp360;
73
     fftw_real tmp290;
74
     fftw_real tmp293;
75
     fftw_real tmp294;
76
     fftw_real tmp170;
77
     fftw_real tmp233;
78
     fftw_real tmp333;
79
     fftw_real tmp359;
80
     fftw_real tmp164;
81
     fftw_real tmp234;
82
     fftw_real tmp173;
83
     fftw_real tmp237;
84
     fftw_real tmp22;
85
     fftw_real tmp318;
86
     fftw_real tmp343;
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 tmp321;
94
     fftw_real tmp342;
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 tmp278;
105
     fftw_real tmp329;
106
     fftw_real tmp357;
107
     fftw_real tmp281;
108
     fftw_real tmp284;
109
     fftw_real tmp285;
110
     fftw_real tmp143;
111
     fftw_real tmp226;
112
     fftw_real tmp326;
113
     fftw_real tmp356;
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 tmp177;
122
          fftw_real tmp66;
123
          fftw_real tmp96;
124
          fftw_real tmp6;
125
          fftw_real tmp95;
126
          fftw_real tmp69;
127
          fftw_real tmp178;
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
               tmp177 = tmp1 - tmp2;
139
               tmp64 = c_im(input[0]);
140
               tmp65 = c_im(input[16 * istride]);
141
               tmp66 = tmp64 + tmp65;
142
               tmp96 = 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
               tmp95 = tmp4 - tmp5;
154
               tmp67 = c_im(input[8 * istride]);
155
               tmp68 = c_im(input[24 * istride]);
156
               tmp69 = tmp67 + tmp68;
157
               tmp178 = tmp67 - tmp68;
158
          }
159
          tmp7 = tmp3 + tmp6;
160
          tmp339 = tmp3 - tmp6;
161
          tmp70 = tmp66 + tmp69;
162
          tmp313 = tmp66 - tmp69;
163
          tmp97 = tmp95 + tmp96;
164
          tmp215 = tmp96 - tmp95;
165
          tmp179 = tmp177 - tmp178;
166
          tmp241 = tmp177 + tmp178;
167
     }
168
     {
169
          fftw_real tmp10;
170
          fftw_real tmp98;
171
          fftw_real tmp73;
172
          fftw_real tmp99;
173
          fftw_real tmp13;
174
          fftw_real tmp102;
175
          fftw_real tmp76;
176
          fftw_real tmp101;
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
               tmp98 = tmp8 - tmp9;
188
               tmp71 = c_im(input[4 * istride]);
189
               tmp72 = c_im(input[20 * istride]);
190
               tmp73 = tmp71 + tmp72;
191
               tmp99 = 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
               tmp102 = tmp11 - tmp12;
203
               tmp74 = c_im(input[28 * istride]);
204
               tmp75 = c_im(input[12 * istride]);
205
               tmp76 = tmp74 + tmp75;
206
               tmp101 = tmp74 - tmp75;
207
          }
208
          tmp14 = tmp10 + tmp13;
209
          tmp314 = tmp10 - tmp13;
210
          tmp77 = tmp73 + tmp76;
211
          tmp340 = tmp76 - tmp73;
212
          {
213
               fftw_real tmp180;
214
               fftw_real tmp181;
215
               fftw_real tmp100;
216
               fftw_real tmp103;
217
               ASSERT_ALIGNED_DOUBLE();
218
               tmp180 = tmp98 - tmp99;
219
               tmp181 = tmp102 + tmp101;
220
               tmp182 = K707106781 * (tmp180 + tmp181);
221
               tmp216 = K707106781 * (tmp180 - tmp181);
222
               tmp100 = tmp98 + tmp99;
223
               tmp103 = tmp101 - tmp102;
224
               tmp104 = K707106781 * (tmp100 + tmp103);
225
               tmp242 = K707106781 * (tmp103 - tmp100);
226
          }
227
     }
228
     {
229
          fftw_real tmp49;
230
          fftw_real tmp149;
231
          fftw_real tmp169;
232
          fftw_real tmp288;
233
          fftw_real tmp52;
234
          fftw_real tmp166;
235
          fftw_real tmp152;
236
          fftw_real tmp289;
237
          fftw_real tmp56;
238
          fftw_real tmp154;
239
          fftw_real tmp157;
240
          fftw_real tmp291;
241
          fftw_real tmp59;
242
          fftw_real tmp159;
243
          fftw_real tmp162;
244
          fftw_real tmp292;
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
               tmp288 = 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
               tmp289 = tmp150 + tmp151;
275
          }
276
          {
277
               fftw_real tmp54;
278
               fftw_real tmp55;
279
               fftw_real tmp155;
280
               fftw_real tmp156;
281
               ASSERT_ALIGNED_DOUBLE();
282
               tmp54 = c_re(input[3 * istride]);
283
               tmp55 = c_re(input[19 * istride]);
284
               tmp56 = tmp54 + tmp55;
285
               tmp154 = tmp54 - tmp55;
286
               tmp155 = c_im(input[3 * istride]);
287
               tmp156 = c_im(input[19 * istride]);
288
               tmp157 = tmp155 - tmp156;
289
               tmp291 = tmp155 + tmp156;
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
               tmp292 = tmp160 + tmp161;
305
          }
306
          {
307
               fftw_real tmp334;
308
               fftw_real tmp335;
309
               fftw_real tmp331;
310
               fftw_real tmp332;
311
               ASSERT_ALIGNED_DOUBLE();
312
               tmp153 = tmp149 - tmp152;
313
               tmp236 = tmp149 + tmp152;
314
               tmp53 = tmp49 + tmp52;
315
               tmp60 = tmp56 + tmp59;
316
               tmp287 = tmp53 - tmp60;
317
               tmp334 = tmp49 - tmp52;
318
               tmp335 = tmp292 - tmp291;
319
               tmp336 = tmp334 - tmp335;
320
               tmp360 = tmp334 + tmp335;
321
               tmp290 = tmp288 + tmp289;
322
               tmp293 = tmp291 + tmp292;
323
               tmp294 = tmp290 - tmp293;
324
               tmp170 = tmp166 + tmp169;
325
               tmp233 = tmp169 - tmp166;
326
               tmp331 = tmp288 - tmp289;
327
               tmp332 = tmp56 - tmp59;
328
               tmp333 = tmp331 - tmp332;
329
               tmp359 = tmp332 + tmp331;
330
               {
331
                    fftw_real tmp158;
332
                    fftw_real tmp163;
333
                    fftw_real tmp171;
334
                    fftw_real tmp172;
335
                    ASSERT_ALIGNED_DOUBLE();
336
                    tmp158 = tmp154 - tmp157;
337
                    tmp163 = tmp159 + tmp162;
338
                    tmp164 = K707106781 * (tmp158 + tmp163);
339
                    tmp234 = K707106781 * (tmp158 - tmp163);
340
                    tmp171 = tmp154 + tmp157;
341
                    tmp172 = tmp162 - tmp159;
342
                    tmp173 = K707106781 * (tmp171 + tmp172);
343
                    tmp237 = K707106781 * (tmp172 - tmp171);
344
               }
345
          }
346
     }
347
     {
348
          fftw_real tmp18;
349
          fftw_real tmp106;
350
          fftw_real tmp81;
351
          fftw_real tmp110;
352
          fftw_real tmp21;
353
          fftw_real tmp109;
354
          fftw_real tmp84;
355
          fftw_real tmp107;
356
          fftw_real tmp316;
357
          fftw_real tmp317;
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
               tmp106 = tmp16 - tmp17;
369
               tmp79 = c_im(input[2 * istride]);
370
               tmp80 = c_im(input[18 * istride]);
371
               tmp81 = tmp79 + tmp80;
372
               tmp110 = 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
               tmp109 = tmp19 - tmp20;
384
               tmp82 = c_im(input[10 * istride]);
385
               tmp83 = c_im(input[26 * istride]);
386
               tmp84 = tmp82 + tmp83;
387
               tmp107 = tmp82 - tmp83;
388
          }
389
          tmp22 = tmp18 + tmp21;
390
          tmp316 = tmp18 - tmp21;
391
          tmp317 = tmp81 - tmp84;
392
          tmp318 = tmp316 - tmp317;
393
          tmp343 = tmp316 + tmp317;
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 = (K923879532 * tmp108) - (K382683432 * tmp111);
404
               tmp185 = (K923879532 * tmp111) + (K382683432 * tmp108);
405
               tmp218 = tmp106 + tmp107;
406
               tmp219 = tmp110 - tmp109;
407
               tmp220 = (K382683432 * tmp218) - (K923879532 * tmp219);
408
               tmp245 = (K382683432 * tmp219) + (K923879532 * tmp218);
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 tmp319;
421
          fftw_real tmp320;
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
          tmp319 = tmp25 - tmp28;
455
          tmp320 = tmp88 - tmp91;
456
          tmp321 = tmp319 + tmp320;
457
          tmp342 = tmp320 - tmp319;
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 = (K923879532 * tmp115) - (K382683432 * tmp118);
469
               tmp221 = tmp114 - tmp113;
470
               tmp222 = tmp116 + tmp117;
471
               tmp223 = (K923879532 * tmp221) + (K382683432 * tmp222);
472
               tmp244 = (K382683432 * tmp221) - (K923879532 * tmp222);
473
          }
474
     }
475
     {
476
          fftw_real tmp34;
477
          fftw_real tmp122;
478
          fftw_real tmp142;
479
          fftw_real tmp279;
480
          fftw_real tmp37;
481
          fftw_real tmp139;
482
          fftw_real tmp125;
483
          fftw_real tmp280;
484
          fftw_real tmp41;
485
          fftw_real tmp127;
486
          fftw_real tmp130;
487
          fftw_real tmp282;
488
          fftw_real tmp44;
489
          fftw_real tmp132;
490
          fftw_real tmp135;
491
          fftw_real tmp283;
492
          ASSERT_ALIGNED_DOUBLE();
493
          {
494
               fftw_real tmp32;
495
               fftw_real tmp33;
496
               fftw_real tmp140;
497
               fftw_real tmp141;
498
               ASSERT_ALIGNED_DOUBLE();
499
               tmp32 = c_re(input[istride]);
500
               tmp33 = c_re(input[17 * istride]);
501
               tmp34 = tmp32 + tmp33;
502
               tmp122 = tmp32 - tmp33;
503
               tmp140 = c_im(input[istride]);
504
               tmp141 = c_im(input[17 * istride]);
505
               tmp142 = tmp140 - tmp141;
506
               tmp279 = tmp140 + tmp141;
507
          }
508
          {
509
               fftw_real tmp35;
510
               fftw_real tmp36;
511
               fftw_real tmp123;
512
               fftw_real tmp124;
513
               ASSERT_ALIGNED_DOUBLE();
514
               tmp35 = c_re(input[9 * istride]);
515
               tmp36 = c_re(input[25 * istride]);
516
               tmp37 = tmp35 + tmp36;
517
               tmp139 = tmp35 - tmp36;
518
               tmp123 = c_im(input[9 * istride]);
519
               tmp124 = c_im(input[25 * istride]);
520
               tmp125 = tmp123 - tmp124;
521
               tmp280 = tmp123 + tmp124;
522
          }
523
          {
524
               fftw_real tmp39;
525
               fftw_real tmp40;
526
               fftw_real tmp128;
527
               fftw_real tmp129;
528
               ASSERT_ALIGNED_DOUBLE();
529
               tmp39 = c_re(input[5 * istride]);
530
               tmp40 = c_re(input[21 * istride]);
531
               tmp41 = tmp39 + tmp40;
532
               tmp127 = tmp39 - tmp40;
533
               tmp128 = c_im(input[5 * istride]);
534
               tmp129 = c_im(input[21 * istride]);
535
               tmp130 = tmp128 - tmp129;
536
               tmp282 = tmp128 + tmp129;
537
          }
538
          {
539
               fftw_real tmp42;
540
               fftw_real tmp43;
541
               fftw_real tmp133;
542
               fftw_real tmp134;
543
               ASSERT_ALIGNED_DOUBLE();
544
               tmp42 = c_re(input[29 * istride]);
545
               tmp43 = c_re(input[13 * istride]);
546
               tmp44 = tmp42 + tmp43;
547
               tmp132 = tmp42 - tmp43;
548
               tmp133 = c_im(input[29 * istride]);
549
               tmp134 = c_im(input[13 * istride]);
550
               tmp135 = tmp133 - tmp134;
551
               tmp283 = tmp133 + tmp134;
552
          }
553
          {
554
               fftw_real tmp327;
555
               fftw_real tmp328;
556
               fftw_real tmp324;
557
               fftw_real tmp325;
558
               ASSERT_ALIGNED_DOUBLE();
559
               tmp126 = tmp122 - tmp125;
560
               tmp229 = tmp122 + tmp125;
561
               tmp38 = tmp34 + tmp37;
562
               tmp45 = tmp41 + tmp44;
563
               tmp278 = tmp38 - tmp45;
564
               tmp327 = tmp34 - tmp37;
565
               tmp328 = tmp283 - tmp282;
566
               tmp329 = tmp327 - tmp328;
567
               tmp357 = tmp327 + tmp328;
568
               tmp281 = tmp279 + tmp280;
569
               tmp284 = tmp282 + tmp283;
570
               tmp285 = tmp281 - tmp284;
571
               tmp143 = tmp139 + tmp142;
572
               tmp226 = tmp142 - tmp139;
573
               tmp324 = tmp279 - tmp280;
574
               tmp325 = tmp41 - tmp44;
575
               tmp326 = tmp324 - tmp325;
576
               tmp356 = tmp325 + tmp324;
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 * (tmp131 - tmp136);
587
                    tmp144 = tmp127 + tmp130;
588
                    tmp145 = tmp135 - tmp132;
589
                    tmp146 = K707106781 * (tmp144 + tmp145);
590
                    tmp230 = K707106781 * (tmp145 - tmp144);
591
               }
592
          }
593
     }
594
     {
595
          fftw_real tmp277;
596
          fftw_real tmp301;
597
          fftw_real tmp304;
598
          fftw_real tmp306;
599
          fftw_real tmp296;
600
          fftw_real tmp300;
601
          fftw_real tmp299;
602
          fftw_real tmp305;
603
          ASSERT_ALIGNED_DOUBLE();
604
          {
605
               fftw_real tmp275;
606
               fftw_real tmp276;
607
               fftw_real tmp302;
608
               fftw_real tmp303;
609
               ASSERT_ALIGNED_DOUBLE();
610
               tmp275 = tmp70 - tmp77;
611
               tmp276 = tmp22 - tmp29;
612
               tmp277 = tmp275 - tmp276;
613
               tmp301 = tmp276 + tmp275;
614
               tmp302 = tmp278 + tmp285;
615
               tmp303 = tmp294 - tmp287;
616
               tmp304 = K707106781 * (tmp302 + tmp303);
617
               tmp306 = K707106781 * (tmp303 - tmp302);
618
          }
619
          {
620
               fftw_real tmp286;
621
               fftw_real tmp295;
622
               fftw_real tmp297;
623
               fftw_real tmp298;
624
               ASSERT_ALIGNED_DOUBLE();
625
               tmp286 = tmp278 - tmp285;
626
               tmp295 = tmp287 + tmp294;
627
               tmp296 = K707106781 * (tmp286 - tmp295);
628
               tmp300 = K707106781 * (tmp286 + tmp295);
629
               tmp297 = tmp7 - tmp14;
630
               tmp298 = tmp92 - tmp85;
631
               tmp299 = tmp297 + tmp298;
632
               tmp305 = tmp297 - tmp298;
633
          }
634
          c_im(output[28 * ostride]) = tmp277 - tmp296;
635
          c_im(output[12 * ostride]) = tmp277 + tmp296;
636
          c_re(output[20 * ostride]) = tmp299 - tmp300;
637
          c_re(output[4 * ostride]) = tmp299 + tmp300;
638
          c_im(output[20 * ostride]) = tmp301 - tmp304;
639
          c_im(output[4 * ostride]) = tmp301 + tmp304;
640
          c_re(output[28 * ostride]) = tmp305 - tmp306;
641
          c_re(output[12 * ostride]) = tmp305 + tmp306;
642
     }
643
     {
644
          fftw_real tmp31;
645
          fftw_real tmp311;
646
          fftw_real tmp310;
647
          fftw_real tmp312;
648
          fftw_real tmp62;
649
          fftw_real tmp63;
650
          fftw_real tmp94;
651
          fftw_real tmp307;
652
          ASSERT_ALIGNED_DOUBLE();
653
          {
654
               fftw_real tmp15;
655
               fftw_real tmp30;
656
               fftw_real tmp308;
657
               fftw_real tmp309;
658
               ASSERT_ALIGNED_DOUBLE();
659
               tmp15 = tmp7 + tmp14;
660
               tmp30 = tmp22 + tmp29;
661
               tmp31 = tmp15 + tmp30;
662
               tmp311 = tmp15 - tmp30;
663
               tmp308 = tmp281 + tmp284;
664
               tmp309 = tmp290 + tmp293;
665
               tmp310 = tmp308 + tmp309;
666
               tmp312 = tmp309 - tmp308;
667
          }
668
          {
669
               fftw_real tmp46;
670
               fftw_real tmp61;
671
               fftw_real tmp78;
672
               fftw_real tmp93;
673
               ASSERT_ALIGNED_DOUBLE();
674
               tmp46 = tmp38 + tmp45;
675
               tmp61 = tmp53 + tmp60;
676
               tmp62 = tmp46 + tmp61;
677
               tmp63 = tmp46 - tmp61;
678
               tmp78 = tmp70 + tmp77;
679
               tmp93 = tmp85 + tmp92;
680
               tmp94 = tmp78 - tmp93;
681
               tmp307 = tmp78 + tmp93;
682
          }
683
          c_re(output[16 * ostride]) = tmp31 - tmp62;
684
          c_re(output[0]) = tmp31 + tmp62;
685
          c_im(output[8 * ostride]) = tmp63 + tmp94;
686
          c_im(output[24 * ostride]) = tmp94 - tmp63;
687
          c_im(output[16 * ostride]) = tmp307 - tmp310;
688
          c_im(output[0]) = tmp307 + tmp310;
689
          c_re(output[24 * ostride]) = tmp311 - tmp312;
690
          c_re(output[8 * ostride]) = tmp311 + tmp312;
691
     }
692
     {
693
          fftw_real tmp121;
694
          fftw_real tmp189;
695
          fftw_real tmp187;
696
          fftw_real tmp193;
697
          fftw_real tmp148;
698
          fftw_real tmp190;
699
          fftw_real tmp175;
700
          fftw_real tmp191;
701
          ASSERT_ALIGNED_DOUBLE();
702
          {
703
               fftw_real tmp105;
704
               fftw_real tmp120;
705
               fftw_real tmp183;
706
               fftw_real tmp186;
707
               ASSERT_ALIGNED_DOUBLE();
708
               tmp105 = tmp97 - tmp104;
709
               tmp120 = tmp112 - tmp119;
710
               tmp121 = tmp105 - tmp120;
711
               tmp189 = tmp105 + tmp120;
712
               tmp183 = tmp179 - tmp182;
713
               tmp186 = tmp184 - tmp185;
714
               tmp187 = tmp183 + tmp186;
715
               tmp193 = tmp183 - tmp186;
716
          }
717
          {
718
               fftw_real tmp138;
719
               fftw_real tmp147;
720
               fftw_real tmp165;
721
               fftw_real tmp174;
722
               ASSERT_ALIGNED_DOUBLE();
723
               tmp138 = tmp126 - tmp137;
724
               tmp147 = tmp143 - tmp146;
725
               tmp148 = (K555570233 * tmp138) - (K831469612 * tmp147);
726
               tmp190 = (K831469612 * tmp138) + (K555570233 * tmp147);
727
               tmp165 = tmp153 - tmp164;
728
               tmp174 = tmp170 - tmp173;
729
               tmp175 = (K555570233 * tmp165) + (K831469612 * tmp174);
730
               tmp191 = (K555570233 * tmp174) - (K831469612 * tmp165);
731
          }
732
          {
733
               fftw_real tmp176;
734
               fftw_real tmp188;
735
               fftw_real tmp192;
736
               fftw_real tmp194;
737
               ASSERT_ALIGNED_DOUBLE();
738
               tmp176 = tmp148 - tmp175;
739
               c_im(output[29 * ostride]) = tmp121 - tmp176;
740
               c_im(output[13 * ostride]) = tmp121 + tmp176;
741
               tmp188 = tmp148 + tmp175;
742
               c_re(output[21 * ostride]) = tmp187 - tmp188;
743
               c_re(output[5 * ostride]) = tmp187 + tmp188;
744
               tmp192 = tmp190 + tmp191;
745
               c_im(output[21 * ostride]) = tmp189 - tmp192;
746
               c_im(output[5 * ostride]) = tmp189 + tmp192;
747
               tmp194 = tmp191 - tmp190;
748
               c_re(output[29 * ostride]) = tmp193 - tmp194;
749
               c_re(output[13 * ostride]) = tmp193 + tmp194;
750
          }
751
     }
752
     {
753
          fftw_real tmp197;
754
          fftw_real tmp209;
755
          fftw_real tmp207;
756
          fftw_real tmp213;
757
          fftw_real tmp200;
758
          fftw_real tmp210;
759
          fftw_real tmp203;
760
          fftw_real tmp211;
761
          ASSERT_ALIGNED_DOUBLE();
762
          {
763
               fftw_real tmp195;
764
               fftw_real tmp196;
765
               fftw_real tmp205;
766
               fftw_real tmp206;
767
               ASSERT_ALIGNED_DOUBLE();
768
               tmp195 = tmp97 + tmp104;
769
               tmp196 = tmp185 + tmp184;
770
               tmp197 = tmp195 - tmp196;
771
               tmp209 = tmp195 + tmp196;
772
               tmp205 = tmp179 + tmp182;
773
               tmp206 = tmp112 + tmp119;
774
               tmp207 = tmp205 + tmp206;
775
               tmp213 = tmp205 - tmp206;
776
          }
777
          {
778
               fftw_real tmp198;
779
               fftw_real tmp199;
780
               fftw_real tmp201;
781
               fftw_real tmp202;
782
               ASSERT_ALIGNED_DOUBLE();
783
               tmp198 = tmp126 + tmp137;
784
               tmp199 = tmp143 + tmp146;
785
               tmp200 = (K980785280 * tmp198) - (K195090322 * tmp199);
786
               tmp210 = (K195090322 * tmp198) + (K980785280 * tmp199);
787
               tmp201 = tmp153 + tmp164;
788
               tmp202 = tmp170 + tmp173;
789
               tmp203 = (K980785280 * tmp201) + (K195090322 * tmp202);
790
               tmp211 = (K980785280 * tmp202) - (K195090322 * tmp201);
791
          }
792
          {
793
               fftw_real tmp204;
794
               fftw_real tmp208;
795
               fftw_real tmp212;
796
               fftw_real tmp214;
797
               ASSERT_ALIGNED_DOUBLE();
798
               tmp204 = tmp200 - tmp203;
799
               c_im(output[25 * ostride]) = tmp197 - tmp204;
800
               c_im(output[9 * ostride]) = tmp197 + tmp204;
801
               tmp208 = tmp200 + tmp203;
802
               c_re(output[17 * ostride]) = tmp207 - tmp208;
803
               c_re(output[ostride]) = tmp207 + tmp208;
804
               tmp212 = tmp210 + tmp211;
805
               c_im(output[17 * ostride]) = tmp209 - tmp212;
806
               c_im(output[ostride]) = tmp209 + tmp212;
807
               tmp214 = tmp211 - tmp210;
808
               c_re(output[25 * ostride]) = tmp213 - tmp214;
809
               c_re(output[9 * ostride]) = tmp213 + tmp214;
810
          }
811
     }
812
     {
813
          fftw_real tmp323;
814
          fftw_real tmp347;
815
          fftw_real tmp350;
816
          fftw_real tmp352;
817
          fftw_real tmp338;
818
          fftw_real tmp346;
819
          fftw_real tmp345;
820
          fftw_real tmp351;
821
          ASSERT_ALIGNED_DOUBLE();
822
          {
823
               fftw_real tmp315;
824
               fftw_real tmp322;
825
               fftw_real tmp348;
826
               fftw_real tmp349;
827
               ASSERT_ALIGNED_DOUBLE();
828
               tmp315 = tmp313 - tmp314;
829
               tmp322 = K707106781 * (tmp318 - tmp321);
830
               tmp323 = tmp315 + tmp322;
831
               tmp347 = tmp315 - tmp322;
832
               tmp348 = (K382683432 * tmp329) - (K923879532 * tmp326);
833
               tmp349 = (K923879532 * tmp333) + (K382683432 * tmp336);
834
               tmp350 = tmp348 - tmp349;
835
               tmp352 = tmp348 + tmp349;
836
          }
837
          {
838
               fftw_real tmp330;
839
               fftw_real tmp337;
840
               fftw_real tmp341;
841
               fftw_real tmp344;
842
               ASSERT_ALIGNED_DOUBLE();
843
               tmp330 = (K382683432 * tmp326) + (K923879532 * tmp329);
844
               tmp337 = (K382683432 * tmp333) - (K923879532 * tmp336);
845
               tmp338 = tmp330 + tmp337;
846
               tmp346 = tmp337 - tmp330;
847
               tmp341 = tmp339 - tmp340;
848
               tmp344 = K707106781 * (tmp342 - tmp343);
849
               tmp345 = tmp341 - tmp344;
850
               tmp351 = tmp341 + tmp344;
851
          }
852
          c_im(output[22 * ostride]) = tmp323 - tmp338;
853
          c_im(output[6 * ostride]) = tmp323 + tmp338;
854
          c_re(output[30 * ostride]) = tmp345 - tmp346;
855
          c_re(output[14 * ostride]) = tmp345 + tmp346;
856
          c_im(output[30 * ostride]) = tmp347 - tmp350;
857
          c_im(output[14 * ostride]) = tmp347 + tmp350;
858
          c_re(output[22 * ostride]) = tmp351 - tmp352;
859
          c_re(output[6 * ostride]) = tmp351 + tmp352;
860
     }
861
     {
862
          fftw_real tmp355;
863
          fftw_real tmp367;
864
          fftw_real tmp370;
865
          fftw_real tmp372;
866
          fftw_real tmp362;
867
          fftw_real tmp366;
868
          fftw_real tmp365;
869
          fftw_real tmp371;
870
          ASSERT_ALIGNED_DOUBLE();
871
          {
872
               fftw_real tmp353;
873
               fftw_real tmp354;
874
               fftw_real tmp368;
875
               fftw_real tmp369;
876
               ASSERT_ALIGNED_DOUBLE();
877
               tmp353 = tmp314 + tmp313;
878
               tmp354 = K707106781 * (tmp343 + tmp342);
879
               tmp355 = tmp353 + tmp354;
880
               tmp367 = tmp353 - tmp354;
881
               tmp368 = (K923879532 * tmp357) - (K382683432 * tmp356);
882
               tmp369 = (K382683432 * tmp359) + (K923879532 * tmp360);
883
               tmp370 = tmp368 - tmp369;
884
               tmp372 = tmp368 + tmp369;
885
          }
886
          {
887
               fftw_real tmp358;
888
               fftw_real tmp361;
889
               fftw_real tmp363;
890
               fftw_real tmp364;
891
               ASSERT_ALIGNED_DOUBLE();
892
               tmp358 = (K923879532 * tmp356) + (K382683432 * tmp357);
893
               tmp361 = (K923879532 * tmp359) - (K382683432 * tmp360);
894
               tmp362 = tmp358 + tmp361;
895
               tmp366 = tmp361 - tmp358;
896
               tmp363 = tmp339 + tmp340;
897
               tmp364 = K707106781 * (tmp318 + tmp321);
898
               tmp365 = tmp363 - tmp364;
899
               tmp371 = tmp363 + tmp364;
900
          }
901
          c_im(output[18 * ostride]) = tmp355 - tmp362;
902
          c_im(output[2 * ostride]) = tmp355 + tmp362;
903
          c_re(output[26 * ostride]) = tmp365 - tmp366;
904
          c_re(output[10 * ostride]) = tmp365 + tmp366;
905
          c_im(output[26 * ostride]) = tmp367 - tmp370;
906
          c_im(output[10 * ostride]) = tmp367 + tmp370;
907
          c_re(output[18 * ostride]) = tmp371 - tmp372;
908
          c_re(output[2 * ostride]) = tmp371 + tmp372;
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 = (K195090322 * tmp228) + (K980785280 * tmp231);
944
               tmp250 = (K195090322 * tmp231) - (K980785280 * tmp228);
945
               tmp235 = tmp233 - tmp234;
946
               tmp238 = tmp236 - tmp237;
947
               tmp239 = (K195090322 * tmp235) - (K980785280 * tmp238);
948
               tmp251 = (K980785280 * tmp235) + (K195090322 * 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_im(output[23 * ostride]) = tmp225 - tmp240;
958
               c_im(output[7 * ostride]) = tmp225 + tmp240;
959
               tmp248 = tmp239 - tmp232;
960
               c_re(output[31 * ostride]) = tmp247 - tmp248;
961
               c_re(output[15 * ostride]) = tmp247 + tmp248;
962
               tmp252 = tmp250 - tmp251;
963
               c_im(output[31 * ostride]) = tmp249 - tmp252;
964
               c_im(output[15 * ostride]) = tmp249 + tmp252;
965
               tmp254 = tmp250 + tmp251;
966
               c_re(output[23 * ostride]) = tmp253 - tmp254;
967
               c_re(output[7 * 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 = (K831469612 * tmp258) + (K555570233 * tmp259);
1004
               tmp270 = (K831469612 * tmp259) - (K555570233 * tmp258);
1005
               tmp261 = tmp233 + tmp234;
1006
               tmp262 = tmp236 + tmp237;
1007
               tmp263 = (K831469612 * tmp261) - (K555570233 * tmp262);
1008
               tmp271 = (K555570233 * tmp261) + (K831469612 * 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_im(output[19 * ostride]) = tmp257 - tmp264;
1018
               c_im(output[3 * ostride]) = tmp257 + tmp264;
1019
               tmp268 = tmp263 - tmp260;
1020
               c_re(output[27 * ostride]) = tmp267 - tmp268;
1021
               c_re(output[11 * ostride]) = tmp267 + tmp268;
1022
               tmp272 = tmp270 - tmp271;
1023
               c_im(output[27 * ostride]) = tmp269 - tmp272;
1024
               c_im(output[11 * ostride]) = tmp269 + tmp272;
1025
               tmp274 = tmp270 + tmp271;
1026
               c_re(output[19 * ostride]) = tmp273 - tmp274;
1027
               c_re(output[3 * ostride]) = tmp273 + tmp274;
1028
          }
1029
     }
1030
}
1031
 
1032
fftw_codelet_desc fftwi_no_twiddle_32_desc =
1033
{
1034
     "fftwi_no_twiddle_32",
1035
     (void (*)()) fftwi_no_twiddle_32,
1036
     32,
1037
     FFTW_BACKWARD,
1038
     FFTW_NOTW,
1039
     716,
1040
     0,
1041
     (const int *) 0,
1042
};