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