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