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