?  search:
summary | log | graphiclog1 | graphiclog2 | commit | commitdiff | tree | refs | edit | fork
blame | history | raw | HEAD
changed reading hint
[gromacs/adressmacs.git] / src / fftw / fhf_32.c
blobd6fc9d641baeb960328413ea8dfb24431062b5d4
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:55:39 EDT 1999 */
22
23 #include <fftw-int.h>
24 #include <fftw.h>
25
26 /* Generated by: ./genfft -magic-alignment-check -magic-twiddle-load-all -magic-variables 4 -magic-loopi -hc2hc-forward 32 */
27
28 /*
29 * This function contains 764 FP additions, 340 FP multiplications,
30 * (or, 618 additions, 194 multiplications, 146 fused multiply/add),
31 * 91 stack variables, and 256 memory accesses
32 */
33 static const fftw_real K145142338 = FFTW_KONST(+0.145142338627231183818096187908697637345738139);
34 static const fftw_real K478470167 = FFTW_KONST(+0.478470167866104432467898943490134984741424603);
35 static const fftw_real K235698368 = FFTW_KONST(+0.235698368412998824278193812952627188828730159);
36 static const fftw_real K440960632 = FFTW_KONST(+0.440960632174177514856378431830194174754221310);
37 static const fftw_real K317196642 = FFTW_KONST(+0.317196642081822749107585806612746685337843547);
38 static const fftw_real K386505226 = FFTW_KONST(+0.386505226681368480405453304879234900485520646);
39 static const fftw_real K497592363 = FFTW_KONST(+0.497592363336098443122418476554739960787737434);
40 static const fftw_real K049008570 = FFTW_KONST(+0.049008570164780300997097781944320922930568337);
41 static const fftw_real K277785116 = FFTW_KONST(+0.277785116509801112371415406974266437187468595);
42 static const fftw_real K415734806 = FFTW_KONST(+0.415734806151272618539394188808952878369280406);
43 static const fftw_real K097545161 = FFTW_KONST(+0.097545161008064133924142434238511120463845809);
44 static const fftw_real K490392640 = FFTW_KONST(+0.490392640201615224563091118067119518486966865);
45 static const fftw_real K765366864 = FFTW_KONST(+0.765366864730179543456919968060797733522689125);
46 static const fftw_real K1_847759065 = FFTW_KONST(+1.847759065022573512256366378793576573644833252);
47 static const fftw_real K1_414213562 = FFTW_KONST(+1.414213562373095048801688724209698078569671875);
48 static const fftw_real K2_000000000 = FFTW_KONST(+2.000000000000000000000000000000000000000000000);
49 static const fftw_real K195090322 = FFTW_KONST(+0.195090322016128267848284868477022240927691618);
50 static const fftw_real K980785280 = FFTW_KONST(+0.980785280403230449126182236134239036973933731);
51 static const fftw_real K555570233 = FFTW_KONST(+0.555570233019602224742830813948532874374937191);
52 static const fftw_real K831469612 = FFTW_KONST(+0.831469612302545237078788377617905756738560812);
53 static const fftw_real K382683432 = FFTW_KONST(+0.382683432365089771728459984030398866761344562);
54 static const fftw_real K923879532 = FFTW_KONST(+0.923879532511286756128183189396788286822416626);
55 static const fftw_real K707106781 = FFTW_KONST(+0.707106781186547524400844362104849039284835938);
56
57 /*
58 * Generator Id's :
59 * $Id$
60 * $Id$
61 * $Id$
62 */
63
64 void fftw_hc2hc_forward_32(fftw_real *A, const fftw_complex *W, int iostride, int m, int dist)
65 {
66 int i;
67 fftw_real *X;
68 fftw_real *Y;
69 X = A;
70 Y = A + (32 * iostride);
71 {
72 fftw_real tmp685;
73 fftw_real tmp813;
74 fftw_real tmp709;
75 fftw_real tmp761;
76 fftw_real tmp692;
77 fftw_real tmp826;
78 fftw_real tmp712;
79 fftw_real tmp760;
80 fftw_real tmp801;
81 fftw_real tmp821;
82 fftw_real tmp749;
83 fftw_real tmp777;
84 fftw_real tmp804;
85 fftw_real tmp822;
86 fftw_real tmp754;
87 fftw_real tmp778;
88 fftw_real tmp700;
89 fftw_real tmp814;
90 fftw_real tmp716;
91 fftw_real tmp758;
92 fftw_real tmp707;
93 fftw_real tmp815;
94 fftw_real tmp719;
95 fftw_real tmp757;
96 fftw_real tmp794;
97 fftw_real tmp818;
98 fftw_real tmp732;
99 fftw_real tmp774;
100 fftw_real tmp797;
101 fftw_real tmp819;
102 fftw_real tmp737;
103 fftw_real tmp775;
104 ASSERT_ALIGNED_DOUBLE();
105 {
106 fftw_real tmp679;
107 fftw_real tmp680;
108 fftw_real tmp681;
109 fftw_real tmp682;
110 fftw_real tmp683;
111 fftw_real tmp684;
112 ASSERT_ALIGNED_DOUBLE();
113 tmp679 = X[0];
114 tmp680 = X[16 * iostride];
115 tmp681 = tmp679 + tmp680;
116 tmp682 = X[8 * iostride];
117 tmp683 = X[24 * iostride];
118 tmp684 = tmp682 + tmp683;
119 tmp685 = tmp681 + tmp684;
120 tmp813 = tmp681 - tmp684;
121 tmp709 = tmp679 - tmp680;
122 tmp761 = tmp682 - tmp683;
123 }
124 {
125 fftw_real tmp688;
126 fftw_real tmp710;
127 fftw_real tmp691;
128 fftw_real tmp711;
129 ASSERT_ALIGNED_DOUBLE();
130 {
131 fftw_real tmp686;
132 fftw_real tmp687;
133 fftw_real tmp689;
134 fftw_real tmp690;
135 ASSERT_ALIGNED_DOUBLE();
136 tmp686 = X[4 * iostride];
137 tmp687 = X[20 * iostride];
138 tmp688 = tmp686 + tmp687;
139 tmp710 = tmp686 - tmp687;
140 tmp689 = X[28 * iostride];
141 tmp690 = X[12 * iostride];
142 tmp691 = tmp689 + tmp690;
143 tmp711 = tmp689 - tmp690;
144 }
145 tmp692 = tmp688 + tmp691;
146 tmp826 = tmp691 - tmp688;
147 tmp712 = K707106781 * (tmp710 + tmp711);
148 tmp760 = K707106781 * (tmp711 - tmp710);
149 }
150 {
151 fftw_real tmp741;
152 fftw_real tmp799;
153 fftw_real tmp753;
154 fftw_real tmp800;
155 fftw_real tmp744;
156 fftw_real tmp802;
157 fftw_real tmp747;
158 fftw_real tmp803;
159 fftw_real tmp748;
160 fftw_real tmp750;
161 ASSERT_ALIGNED_DOUBLE();
162 {
163 fftw_real tmp739;
164 fftw_real tmp740;
165 fftw_real tmp751;
166 fftw_real tmp752;
167 ASSERT_ALIGNED_DOUBLE();
168 tmp739 = X[31 * iostride];
169 tmp740 = X[15 * iostride];
170 tmp741 = tmp739 - tmp740;
171 tmp799 = tmp739 + tmp740;
172 tmp751 = X[7 * iostride];
173 tmp752 = X[23 * iostride];
174 tmp753 = tmp751 - tmp752;
175 tmp800 = tmp751 + tmp752;
176 }
177 {
178 fftw_real tmp742;
179 fftw_real tmp743;
180 fftw_real tmp745;
181 fftw_real tmp746;
182 ASSERT_ALIGNED_DOUBLE();
183 tmp742 = X[3 * iostride];
184 tmp743 = X[19 * iostride];
185 tmp744 = tmp742 - tmp743;
186 tmp802 = tmp742 + tmp743;
187 tmp745 = X[27 * iostride];
188 tmp746 = X[11 * iostride];
189 tmp747 = tmp745 - tmp746;
190 tmp803 = tmp745 + tmp746;
191 }
192 tmp801 = tmp799 + tmp800;
193 tmp821 = tmp799 - tmp800;
194 tmp748 = K707106781 * (tmp744 + tmp747);
195 tmp749 = tmp741 + tmp748;
196 tmp777 = tmp741 - tmp748;
197 tmp804 = tmp802 + tmp803;
198 tmp822 = tmp803 - tmp802;
199 tmp750 = K707106781 * (tmp747 - tmp744);
200 tmp754 = tmp750 - tmp753;
201 tmp778 = tmp753 + tmp750;
202 }
203 {
204 fftw_real tmp696;
205 fftw_real tmp714;
206 fftw_real tmp699;
207 fftw_real tmp715;
208 ASSERT_ALIGNED_DOUBLE();
209 {
210 fftw_real tmp694;
211 fftw_real tmp695;
212 fftw_real tmp697;
213 fftw_real tmp698;
214 ASSERT_ALIGNED_DOUBLE();
215 tmp694 = X[2 * iostride];
216 tmp695 = X[18 * iostride];
217 tmp696 = tmp694 + tmp695;
218 tmp714 = tmp694 - tmp695;
219 tmp697 = X[10 * iostride];
220 tmp698 = X[26 * iostride];
221 tmp699 = tmp697 + tmp698;
222 tmp715 = tmp697 - tmp698;
223 }
224 tmp700 = tmp696 + tmp699;
225 tmp814 = tmp696 - tmp699;
226 tmp716 = (K923879532 * tmp714) - (K382683432 * tmp715);
227 tmp758 = (K382683432 * tmp714) + (K923879532 * tmp715);
228 }
229 {
230 fftw_real tmp703;
231 fftw_real tmp717;
232 fftw_real tmp706;
233 fftw_real tmp718;
234 ASSERT_ALIGNED_DOUBLE();
235 {
236 fftw_real tmp701;
237 fftw_real tmp702;
238 fftw_real tmp704;
239 fftw_real tmp705;
240 ASSERT_ALIGNED_DOUBLE();
241 tmp701 = X[30 * iostride];
242 tmp702 = X[14 * iostride];
243 tmp703 = tmp701 + tmp702;
244 tmp717 = tmp701 - tmp702;
245 tmp704 = X[6 * iostride];
246 tmp705 = X[22 * iostride];
247 tmp706 = tmp704 + tmp705;
248 tmp718 = tmp704 - tmp705;
249 }
250 tmp707 = tmp703 + tmp706;
251 tmp815 = tmp703 - tmp706;
252 tmp719 = (K923879532 * tmp717) + (K382683432 * tmp718);
253 tmp757 = (K382683432 * tmp717) - (K923879532 * tmp718);
254 }
255 {
256 fftw_real tmp724;
257 fftw_real tmp792;
258 fftw_real tmp736;
259 fftw_real tmp793;
260 fftw_real tmp727;
261 fftw_real tmp795;
262 fftw_real tmp730;
263 fftw_real tmp796;
264 fftw_real tmp731;
265 fftw_real tmp733;
266 ASSERT_ALIGNED_DOUBLE();
267 {
268 fftw_real tmp722;
269 fftw_real tmp723;
270 fftw_real tmp734;
271 fftw_real tmp735;
272 ASSERT_ALIGNED_DOUBLE();
273 tmp722 = X[iostride];
274 tmp723 = X[17 * iostride];
275 tmp724 = tmp722 - tmp723;
276 tmp792 = tmp722 + tmp723;
277 tmp734 = X[9 * iostride];
278 tmp735 = X[25 * iostride];
279 tmp736 = tmp734 - tmp735;
280 tmp793 = tmp734 + tmp735;
281 }
282 {
283 fftw_real tmp725;
284 fftw_real tmp726;
285 fftw_real tmp728;
286 fftw_real tmp729;
287 ASSERT_ALIGNED_DOUBLE();
288 tmp725 = X[5 * iostride];
289 tmp726 = X[21 * iostride];
290 tmp727 = tmp725 - tmp726;
291 tmp795 = tmp725 + tmp726;
292 tmp728 = X[29 * iostride];
293 tmp729 = X[13 * iostride];
294 tmp730 = tmp728 - tmp729;
295 tmp796 = tmp728 + tmp729;
296 }
297 tmp794 = tmp792 + tmp793;
298 tmp818 = tmp792 - tmp793;
299 tmp731 = K707106781 * (tmp727 + tmp730);
300 tmp732 = tmp724 + tmp731;
301 tmp774 = tmp724 - tmp731;
302 tmp797 = tmp795 + tmp796;
303 tmp819 = tmp796 - tmp795;
304 tmp733 = K707106781 * (tmp730 - tmp727);
305 tmp737 = tmp733 - tmp736;
306 tmp775 = tmp736 + tmp733;
307 }
308 {
309 fftw_real tmp693;
310 fftw_real tmp708;
311 fftw_real tmp809;
312 fftw_real tmp810;
313 fftw_real tmp811;
314 fftw_real tmp812;
315 ASSERT_ALIGNED_DOUBLE();
316 tmp693 = tmp685 + tmp692;
317 tmp708 = tmp700 + tmp707;
318 tmp809 = tmp693 + tmp708;
319 tmp810 = tmp794 + tmp797;
320 tmp811 = tmp801 + tmp804;
321 tmp812 = tmp810 + tmp811;
322 X[8 * iostride] = tmp693 - tmp708;
323 Y[-8 * iostride] = tmp811 - tmp810;
324 X[16 * iostride] = tmp809 - tmp812;
325 X[0] = tmp809 + tmp812;
326 }
327 {
328 fftw_real tmp791;
329 fftw_real tmp807;
330 fftw_real tmp806;
331 fftw_real tmp808;
332 fftw_real tmp798;
333 fftw_real tmp805;
334 ASSERT_ALIGNED_DOUBLE();
335 tmp791 = tmp685 - tmp692;
336 tmp807 = tmp707 - tmp700;
337 tmp798 = tmp794 - tmp797;
338 tmp805 = tmp801 - tmp804;
339 tmp806 = K707106781 * (tmp798 + tmp805);
340 tmp808 = K707106781 * (tmp805 - tmp798);
341 X[12 * iostride] = tmp791 - tmp806;
342 X[4 * iostride] = tmp791 + tmp806;
343 Y[-4 * iostride] = tmp807 + tmp808;
344 Y[-12 * iostride] = tmp808 - tmp807;
345 }
346 {
347 fftw_real tmp817;
348 fftw_real tmp833;
349 fftw_real tmp827;
350 fftw_real tmp829;
351 fftw_real tmp824;
352 fftw_real tmp828;
353 fftw_real tmp832;
354 fftw_real tmp834;
355 fftw_real tmp816;
356 fftw_real tmp825;
357 ASSERT_ALIGNED_DOUBLE();
358 tmp816 = K707106781 * (tmp814 + tmp815);
359 tmp817 = tmp813 + tmp816;
360 tmp833 = tmp813 - tmp816;
361 tmp825 = K707106781 * (tmp815 - tmp814);
362 tmp827 = tmp825 - tmp826;
363 tmp829 = tmp826 + tmp825;
364 {
365 fftw_real tmp820;
366 fftw_real tmp823;
367 fftw_real tmp830;
368 fftw_real tmp831;
369 ASSERT_ALIGNED_DOUBLE();
370 tmp820 = (K923879532 * tmp818) + (K382683432 * tmp819);
371 tmp823 = (K923879532 * tmp821) - (K382683432 * tmp822);
372 tmp824 = tmp820 + tmp823;
373 tmp828 = tmp823 - tmp820;
374 tmp830 = (K923879532 * tmp819) - (K382683432 * tmp818);
375 tmp831 = (K382683432 * tmp821) + (K923879532 * tmp822);
376 tmp832 = tmp830 + tmp831;
377 tmp834 = tmp831 - tmp830;
378 }
379 X[14 * iostride] = tmp817 - tmp824;
380 X[2 * iostride] = tmp817 + tmp824;
381 Y[-6 * iostride] = tmp827 + tmp828;
382 Y[-10 * iostride] = tmp828 - tmp827;
383 Y[-2 * iostride] = tmp829 + tmp832;
384 Y[-14 * iostride] = tmp832 - tmp829;
385 X[10 * iostride] = tmp833 - tmp834;
386 X[6 * iostride] = tmp833 + tmp834;
387 }
388 {
389 fftw_real tmp773;
390 fftw_real tmp789;
391 fftw_real tmp788;
392 fftw_real tmp790;
393 fftw_real tmp780;
394 fftw_real tmp784;
395 fftw_real tmp783;
396 fftw_real tmp785;
397 ASSERT_ALIGNED_DOUBLE();
398 {
399 fftw_real tmp771;
400 fftw_real tmp772;
401 fftw_real tmp786;
402 fftw_real tmp787;
403 ASSERT_ALIGNED_DOUBLE();
404 tmp771 = tmp709 - tmp712;
405 tmp772 = tmp758 + tmp757;
406 tmp773 = tmp771 + tmp772;
407 tmp789 = tmp771 - tmp772;
408 tmp786 = (K831469612 * tmp775) - (K555570233 * tmp774);
409 tmp787 = (K555570233 * tmp777) + (K831469612 * tmp778);
410 tmp788 = tmp786 + tmp787;
411 tmp790 = tmp787 - tmp786;
412 }
413 {
414 fftw_real tmp776;
415 fftw_real tmp779;
416 fftw_real tmp781;
417 fftw_real tmp782;
418 ASSERT_ALIGNED_DOUBLE();
419 tmp776 = (K831469612 * tmp774) + (K555570233 * tmp775);
420 tmp779 = (K831469612 * tmp777) - (K555570233 * tmp778);
421 tmp780 = tmp776 + tmp779;
422 tmp784 = tmp779 - tmp776;
423 tmp781 = tmp719 - tmp716;
424 tmp782 = tmp761 + tmp760;
425 tmp783 = tmp781 - tmp782;
426 tmp785 = tmp782 + tmp781;
427 }
428 X[13 * iostride] = tmp773 - tmp780;
429 X[3 * iostride] = tmp773 + tmp780;
430 Y[-5 * iostride] = tmp783 + tmp784;
431 Y[-11 * iostride] = tmp784 - tmp783;
432 Y[-3 * iostride] = tmp785 + tmp788;
433 Y[-13 * iostride] = tmp788 - tmp785;
434 X[11 * iostride] = tmp789 - tmp790;
435 X[5 * iostride] = tmp789 + tmp790;
436 }
437 {
438 fftw_real tmp721;
439 fftw_real tmp769;
440 fftw_real tmp768;
441 fftw_real tmp770;
442 fftw_real tmp756;
443 fftw_real tmp764;
444 fftw_real tmp763;
445 fftw_real tmp765;
446 ASSERT_ALIGNED_DOUBLE();
447 {
448 fftw_real tmp713;
449 fftw_real tmp720;
450 fftw_real tmp766;
451 fftw_real tmp767;
452 ASSERT_ALIGNED_DOUBLE();
453 tmp713 = tmp709 + tmp712;
454 tmp720 = tmp716 + tmp719;
455 tmp721 = tmp713 + tmp720;
456 tmp769 = tmp713 - tmp720;
457 tmp766 = (K980785280 * tmp737) - (K195090322 * tmp732);
458 tmp767 = (K195090322 * tmp749) + (K980785280 * tmp754);
459 tmp768 = tmp766 + tmp767;
460 tmp770 = tmp767 - tmp766;
461 }
462 {
463 fftw_real tmp738;
464 fftw_real tmp755;
465 fftw_real tmp759;
466 fftw_real tmp762;
467 ASSERT_ALIGNED_DOUBLE();
468 tmp738 = (K980785280 * tmp732) + (K195090322 * tmp737);
469 tmp755 = (K980785280 * tmp749) - (K195090322 * tmp754);
470 tmp756 = tmp738 + tmp755;
471 tmp764 = tmp755 - tmp738;
472 tmp759 = tmp757 - tmp758;
473 tmp762 = tmp760 - tmp761;
474 tmp763 = tmp759 - tmp762;
475 tmp765 = tmp762 + tmp759;
476 }
477 X[15 * iostride] = tmp721 - tmp756;
478 X[iostride] = tmp721 + tmp756;
479 Y[-7 * iostride] = tmp763 + tmp764;
480 Y[-9 * iostride] = tmp764 - tmp763;
481 Y[-iostride] = tmp765 + tmp768;
482 Y[-15 * iostride] = tmp768 - tmp765;
483 X[9 * iostride] = tmp769 - tmp770;
484 X[7 * iostride] = tmp769 + tmp770;
485 }
486 }
487 X = X + dist;
488 Y = Y - dist;
489 for (i = 2; i < m; i = i + 2, X = X + dist, Y = Y - dist, W = W + 31) {
490 fftw_real tmp201;
491 fftw_real tmp533;
492 fftw_real tmp653;
493 fftw_real tmp667;
494 fftw_real tmp623;
495 fftw_real tmp637;
496 fftw_real tmp373;
497 fftw_real tmp485;
498 fftw_real tmp343;
499 fftw_real tmp561;
500 fftw_real tmp458;
501 fftw_real tmp508;
502 fftw_real tmp568;
503 fftw_real tmp604;
504 fftw_real tmp441;
505 fftw_real tmp505;
506 fftw_real tmp224;
507 fftw_real tmp636;
508 fftw_real tmp383;
509 fftw_real tmp487;
510 fftw_real tmp536;
511 fftw_real tmp618;
512 fftw_real tmp378;
513 fftw_real tmp486;
514 fftw_real tmp366;
515 fftw_real tmp569;
516 fftw_real tmp564;
517 fftw_real tmp605;
518 fftw_real tmp452;
519 fftw_real tmp509;
520 fftw_real tmp461;
521 fftw_real tmp506;
522 fftw_real tmp248;
523 fftw_real tmp541;
524 fftw_real tmp395;
525 fftw_real tmp491;
526 fftw_real tmp540;
527 fftw_real tmp594;
528 fftw_real tmp390;
529 fftw_real tmp490;
530 fftw_real tmp296;
531 fftw_real tmp555;
532 fftw_real tmp431;
533 fftw_real tmp498;
534 fftw_real tmp552;
535 fftw_real tmp599;
536 fftw_real tmp414;
537 fftw_real tmp501;
538 fftw_real tmp271;
539 fftw_real tmp543;
540 fftw_real tmp406;
541 fftw_real tmp494;
542 fftw_real tmp546;
543 fftw_real tmp595;
544 fftw_real tmp401;
545 fftw_real tmp493;
546 fftw_real tmp319;
547 fftw_real tmp553;
548 fftw_real tmp558;
549 fftw_real tmp600;
550 fftw_real tmp425;
551 fftw_real tmp499;
552 fftw_real tmp434;
553 fftw_real tmp502;
554 ASSERT_ALIGNED_DOUBLE();
555 {
556 fftw_real tmp183;
557 fftw_real tmp621;
558 fftw_real tmp188;
559 fftw_real tmp620;
560 fftw_real tmp194;
561 fftw_real tmp370;
562 fftw_real tmp199;
563 fftw_real tmp371;
564 ASSERT_ALIGNED_DOUBLE();
565 tmp183 = X[0];
566 tmp621 = Y[-31 * iostride];
567 {
568 fftw_real tmp185;
569 fftw_real tmp187;
570 fftw_real tmp184;
571 fftw_real tmp186;
572 ASSERT_ALIGNED_DOUBLE();
573 tmp185 = X[16 * iostride];
574 tmp187 = Y[-15 * iostride];
575 tmp184 = c_re(W[15]);
576 tmp186 = c_im(W[15]);
577 tmp188 = (tmp184 * tmp185) - (tmp186 * tmp187);
578 tmp620 = (tmp186 * tmp185) + (tmp184 * tmp187);
579 }
580 {
581 fftw_real tmp191;
582 fftw_real tmp193;
583 fftw_real tmp190;
584 fftw_real tmp192;
585 ASSERT_ALIGNED_DOUBLE();
586 tmp191 = X[8 * iostride];
587 tmp193 = Y[-23 * iostride];
588 tmp190 = c_re(W[7]);
589 tmp192 = c_im(W[7]);
590 tmp194 = (tmp190 * tmp191) - (tmp192 * tmp193);
591 tmp370 = (tmp192 * tmp191) + (tmp190 * tmp193);
592 }
593 {
594 fftw_real tmp196;
595 fftw_real tmp198;
596 fftw_real tmp195;
597 fftw_real tmp197;
598 ASSERT_ALIGNED_DOUBLE();
599 tmp196 = X[24 * iostride];
600 tmp198 = Y[-7 * iostride];
601 tmp195 = c_re(W[23]);
602 tmp197 = c_im(W[23]);
603 tmp199 = (tmp195 * tmp196) - (tmp197 * tmp198);
604 tmp371 = (tmp197 * tmp196) + (tmp195 * tmp198);
605 }
606 {
607 fftw_real tmp189;
608 fftw_real tmp200;
609 fftw_real tmp651;
610 fftw_real tmp652;
611 ASSERT_ALIGNED_DOUBLE();
612 tmp189 = tmp183 + tmp188;
613 tmp200 = tmp194 + tmp199;
614 tmp201 = tmp189 + tmp200;
615 tmp533 = tmp189 - tmp200;
616 tmp651 = tmp621 - tmp620;
617 tmp652 = tmp194 - tmp199;
618 tmp653 = tmp651 - tmp652;
619 tmp667 = tmp652 + tmp651;
620 }
621 {
622 fftw_real tmp619;
623 fftw_real tmp622;
624 fftw_real tmp369;
625 fftw_real tmp372;
626 ASSERT_ALIGNED_DOUBLE();
627 tmp619 = tmp370 + tmp371;
628 tmp622 = tmp620 + tmp621;
629 tmp623 = tmp619 + tmp622;
630 tmp637 = tmp622 - tmp619;
631 tmp369 = tmp183 - tmp188;
632 tmp372 = tmp370 - tmp371;
633 tmp373 = tmp369 - tmp372;
634 tmp485 = tmp369 + tmp372;
635 }
636 }
637 {
638 fftw_real tmp325;
639 fftw_real tmp454;
640 fftw_real tmp341;
641 fftw_real tmp439;
642 fftw_real tmp330;
643 fftw_real tmp455;
644 fftw_real tmp336;
645 fftw_real tmp438;
646 ASSERT_ALIGNED_DOUBLE();
647 {
648 fftw_real tmp322;
649 fftw_real tmp324;
650 fftw_real tmp321;
651 fftw_real tmp323;
652 ASSERT_ALIGNED_DOUBLE();
653 tmp322 = X[31 * iostride];
654 tmp324 = Y[0];
655 tmp321 = c_re(W[30]);
656 tmp323 = c_im(W[30]);
657 tmp325 = (tmp321 * tmp322) - (tmp323 * tmp324);
658 tmp454 = (tmp323 * tmp322) + (tmp321 * tmp324);
659 }
660 {
661 fftw_real tmp338;
662 fftw_real tmp340;
663 fftw_real tmp337;
664 fftw_real tmp339;
665 ASSERT_ALIGNED_DOUBLE();
666 tmp338 = X[23 * iostride];
667 tmp340 = Y[-8 * iostride];
668 tmp337 = c_re(W[22]);
669 tmp339 = c_im(W[22]);
670 tmp341 = (tmp337 * tmp338) - (tmp339 * tmp340);
671 tmp439 = (tmp339 * tmp338) + (tmp337 * tmp340);
672 }
673 {
674 fftw_real tmp327;
675 fftw_real tmp329;
676 fftw_real tmp326;
677 fftw_real tmp328;
678 ASSERT_ALIGNED_DOUBLE();
679 tmp327 = X[15 * iostride];
680 tmp329 = Y[-16 * iostride];
681 tmp326 = c_re(W[14]);
682 tmp328 = c_im(W[14]);
683 tmp330 = (tmp326 * tmp327) - (tmp328 * tmp329);
684 tmp455 = (tmp328 * tmp327) + (tmp326 * tmp329);
685 }
686 {
687 fftw_real tmp333;
688 fftw_real tmp335;
689 fftw_real tmp332;
690 fftw_real tmp334;
691 ASSERT_ALIGNED_DOUBLE();
692 tmp333 = X[7 * iostride];
693 tmp335 = Y[-24 * iostride];
694 tmp332 = c_re(W[6]);
695 tmp334 = c_im(W[6]);
696 tmp336 = (tmp332 * tmp333) - (tmp334 * tmp335);
697 tmp438 = (tmp334 * tmp333) + (tmp332 * tmp335);
698 }
699 {
700 fftw_real tmp331;
701 fftw_real tmp342;
702 fftw_real tmp456;
703 fftw_real tmp457;
704 ASSERT_ALIGNED_DOUBLE();
705 tmp331 = tmp325 + tmp330;
706 tmp342 = tmp336 + tmp341;
707 tmp343 = tmp331 + tmp342;
708 tmp561 = tmp331 - tmp342;
709 tmp456 = tmp454 - tmp455;
710 tmp457 = tmp336 - tmp341;
711 tmp458 = tmp456 + tmp457;
712 tmp508 = tmp456 - tmp457;
713 }
714 {
715 fftw_real tmp566;
716 fftw_real tmp567;
717 fftw_real tmp437;
718 fftw_real tmp440;
719 ASSERT_ALIGNED_DOUBLE();
720 tmp566 = tmp454 + tmp455;
721 tmp567 = tmp438 + tmp439;
722 tmp568 = tmp566 - tmp567;
723 tmp604 = tmp566 + tmp567;
724 tmp437 = tmp325 - tmp330;
725 tmp440 = tmp438 - tmp439;
726 tmp441 = tmp437 - tmp440;
727 tmp505 = tmp437 + tmp440;
728 }
729 }
730 {
731 fftw_real tmp206;
732 fftw_real tmp374;
733 fftw_real tmp222;
734 fftw_real tmp381;
735 fftw_real tmp211;
736 fftw_real tmp375;
737 fftw_real tmp217;
738 fftw_real tmp380;
739 ASSERT_ALIGNED_DOUBLE();
740 {
741 fftw_real tmp203;
742 fftw_real tmp205;
743 fftw_real tmp202;
744 fftw_real tmp204;
745 ASSERT_ALIGNED_DOUBLE();
746 tmp203 = X[4 * iostride];
747 tmp205 = Y[-27 * iostride];
748 tmp202 = c_re(W[3]);
749 tmp204 = c_im(W[3]);
750 tmp206 = (tmp202 * tmp203) - (tmp204 * tmp205);
751 tmp374 = (tmp204 * tmp203) + (tmp202 * tmp205);
752 }
753 {
754 fftw_real tmp219;
755 fftw_real tmp221;
756 fftw_real tmp218;
757 fftw_real tmp220;
758 ASSERT_ALIGNED_DOUBLE();
759 tmp219 = X[12 * iostride];
760 tmp221 = Y[-19 * iostride];
761 tmp218 = c_re(W[11]);
762 tmp220 = c_im(W[11]);
763 tmp222 = (tmp218 * tmp219) - (tmp220 * tmp221);
764 tmp381 = (tmp220 * tmp219) + (tmp218 * tmp221);
765 }
766 {
767 fftw_real tmp208;
768 fftw_real tmp210;
769 fftw_real tmp207;
770 fftw_real tmp209;
771 ASSERT_ALIGNED_DOUBLE();
772 tmp208 = X[20 * iostride];
773 tmp210 = Y[-11 * iostride];
774 tmp207 = c_re(W[19]);
775 tmp209 = c_im(W[19]);
776 tmp211 = (tmp207 * tmp208) - (tmp209 * tmp210);
777 tmp375 = (tmp209 * tmp208) + (tmp207 * tmp210);
778 }
779 {
780 fftw_real tmp214;
781 fftw_real tmp216;
782 fftw_real tmp213;
783 fftw_real tmp215;
784 ASSERT_ALIGNED_DOUBLE();
785 tmp214 = X[28 * iostride];
786 tmp216 = Y[-3 * iostride];
787 tmp213 = c_re(W[27]);
788 tmp215 = c_im(W[27]);
789 tmp217 = (tmp213 * tmp214) - (tmp215 * tmp216);
790 tmp380 = (tmp215 * tmp214) + (tmp213 * tmp216);
791 }
792 {
793 fftw_real tmp212;
794 fftw_real tmp223;
795 fftw_real tmp379;
796 fftw_real tmp382;
797 ASSERT_ALIGNED_DOUBLE();
798 tmp212 = tmp206 + tmp211;
799 tmp223 = tmp217 + tmp222;
800 tmp224 = tmp212 + tmp223;
801 tmp636 = tmp223 - tmp212;
802 tmp379 = tmp217 - tmp222;
803 tmp382 = tmp380 - tmp381;
804 tmp383 = tmp379 + tmp382;
805 tmp487 = tmp379 - tmp382;
806 }
807 {
808 fftw_real tmp534;
809 fftw_real tmp535;
810 fftw_real tmp376;
811 fftw_real tmp377;
812 ASSERT_ALIGNED_DOUBLE();
813 tmp534 = tmp374 + tmp375;
814 tmp535 = tmp380 + tmp381;
815 tmp536 = tmp534 - tmp535;
816 tmp618 = tmp534 + tmp535;
817 tmp376 = tmp374 - tmp375;
818 tmp377 = tmp206 - tmp211;
819 tmp378 = tmp376 - tmp377;
820 tmp486 = tmp377 + tmp376;
821 }
822 }
823 {
824 fftw_real tmp348;
825 fftw_real tmp442;
826 fftw_real tmp353;
827 fftw_real tmp443;
828 fftw_real tmp444;
829 fftw_real tmp445;
830 fftw_real tmp359;
831 fftw_real tmp448;
832 fftw_real tmp364;
833 fftw_real tmp449;
834 fftw_real tmp447;
835 fftw_real tmp450;
836 ASSERT_ALIGNED_DOUBLE();
837 {
838 fftw_real tmp345;
839 fftw_real tmp347;
840 fftw_real tmp344;
841 fftw_real tmp346;
842 ASSERT_ALIGNED_DOUBLE();
843 tmp345 = X[3 * iostride];
844 tmp347 = Y[-28 * iostride];
845 tmp344 = c_re(W[2]);
846 tmp346 = c_im(W[2]);
847 tmp348 = (tmp344 * tmp345) - (tmp346 * tmp347);
848 tmp442 = (tmp346 * tmp345) + (tmp344 * tmp347);
849 }
850 {
851 fftw_real tmp350;
852 fftw_real tmp352;
853 fftw_real tmp349;
854 fftw_real tmp351;
855 ASSERT_ALIGNED_DOUBLE();
856 tmp350 = X[19 * iostride];
857 tmp352 = Y[-12 * iostride];
858 tmp349 = c_re(W[18]);
859 tmp351 = c_im(W[18]);
860 tmp353 = (tmp349 * tmp350) - (tmp351 * tmp352);
861 tmp443 = (tmp351 * tmp350) + (tmp349 * tmp352);
862 }
863 tmp444 = tmp442 - tmp443;
864 tmp445 = tmp348 - tmp353;
865 {
866 fftw_real tmp356;
867 fftw_real tmp358;
868 fftw_real tmp355;
869 fftw_real tmp357;
870 ASSERT_ALIGNED_DOUBLE();
871 tmp356 = X[27 * iostride];
872 tmp358 = Y[-4 * iostride];
873 tmp355 = c_re(W[26]);
874 tmp357 = c_im(W[26]);
875 tmp359 = (tmp355 * tmp356) - (tmp357 * tmp358);
876 tmp448 = (tmp357 * tmp356) + (tmp355 * tmp358);
877 }
878 {
879 fftw_real tmp361;
880 fftw_real tmp363;
881 fftw_real tmp360;
882 fftw_real tmp362;
883 ASSERT_ALIGNED_DOUBLE();
884 tmp361 = X[11 * iostride];
885 tmp363 = Y[-20 * iostride];
886 tmp360 = c_re(W[10]);
887 tmp362 = c_im(W[10]);
888 tmp364 = (tmp360 * tmp361) - (tmp362 * tmp363);
889 tmp449 = (tmp362 * tmp361) + (tmp360 * tmp363);
890 }
891 tmp447 = tmp359 - tmp364;
892 tmp450 = tmp448 - tmp449;
893 {
894 fftw_real tmp354;
895 fftw_real tmp365;
896 fftw_real tmp562;
897 fftw_real tmp563;
898 ASSERT_ALIGNED_DOUBLE();
899 tmp354 = tmp348 + tmp353;
900 tmp365 = tmp359 + tmp364;
901 tmp366 = tmp354 + tmp365;
902 tmp569 = tmp365 - tmp354;
903 tmp562 = tmp442 + tmp443;
904 tmp563 = tmp448 + tmp449;
905 tmp564 = tmp562 - tmp563;
906 tmp605 = tmp562 + tmp563;
907 }
908 {
909 fftw_real tmp446;
910 fftw_real tmp451;
911 fftw_real tmp459;
912 fftw_real tmp460;
913 ASSERT_ALIGNED_DOUBLE();
914 tmp446 = tmp444 - tmp445;
915 tmp451 = tmp447 + tmp450;
916 tmp452 = K707106781 * (tmp446 - tmp451);
917 tmp509 = K707106781 * (tmp446 + tmp451);
918 tmp459 = tmp447 - tmp450;
919 tmp460 = tmp445 + tmp444;
920 tmp461 = K707106781 * (tmp459 - tmp460);
921 tmp506 = K707106781 * (tmp460 + tmp459);
922 }
923 }
924 {
925 fftw_real tmp230;
926 fftw_real tmp386;
927 fftw_real tmp246;
928 fftw_real tmp393;
929 fftw_real tmp235;
930 fftw_real tmp387;
931 fftw_real tmp241;
932 fftw_real tmp392;
933 ASSERT_ALIGNED_DOUBLE();
934 {
935 fftw_real tmp227;
936 fftw_real tmp229;
937 fftw_real tmp226;
938 fftw_real tmp228;
939 ASSERT_ALIGNED_DOUBLE();
940 tmp227 = X[2 * iostride];
941 tmp229 = Y[-29 * iostride];
942 tmp226 = c_re(W[1]);
943 tmp228 = c_im(W[1]);
944 tmp230 = (tmp226 * tmp227) - (tmp228 * tmp229);
945 tmp386 = (tmp228 * tmp227) + (tmp226 * tmp229);
946 }
947 {
948 fftw_real tmp243;
949 fftw_real tmp245;
950 fftw_real tmp242;
951 fftw_real tmp244;
952 ASSERT_ALIGNED_DOUBLE();
953 tmp243 = X[26 * iostride];
954 tmp245 = Y[-5 * iostride];
955 tmp242 = c_re(W[25]);
956 tmp244 = c_im(W[25]);
957 tmp246 = (tmp242 * tmp243) - (tmp244 * tmp245);
958 tmp393 = (tmp244 * tmp243) + (tmp242 * tmp245);
959 }
960 {
961 fftw_real tmp232;
962 fftw_real tmp234;
963 fftw_real tmp231;
964 fftw_real tmp233;
965 ASSERT_ALIGNED_DOUBLE();
966 tmp232 = X[18 * iostride];
967 tmp234 = Y[-13 * iostride];
968 tmp231 = c_re(W[17]);
969 tmp233 = c_im(W[17]);
970 tmp235 = (tmp231 * tmp232) - (tmp233 * tmp234);
971 tmp387 = (tmp233 * tmp232) + (tmp231 * tmp234);
972 }
973 {
974 fftw_real tmp238;
975 fftw_real tmp240;
976 fftw_real tmp237;
977 fftw_real tmp239;
978 ASSERT_ALIGNED_DOUBLE();
979 tmp238 = X[10 * iostride];
980 tmp240 = Y[-21 * iostride];
981 tmp237 = c_re(W[9]);
982 tmp239 = c_im(W[9]);
983 tmp241 = (tmp237 * tmp238) - (tmp239 * tmp240);
984 tmp392 = (tmp239 * tmp238) + (tmp237 * tmp240);
985 }
986 {
987 fftw_real tmp236;
988 fftw_real tmp247;
989 fftw_real tmp391;
990 fftw_real tmp394;
991 ASSERT_ALIGNED_DOUBLE();
992 tmp236 = tmp230 + tmp235;
993 tmp247 = tmp241 + tmp246;
994 tmp248 = tmp236 + tmp247;
995 tmp541 = tmp236 - tmp247;
996 tmp391 = tmp230 - tmp235;
997 tmp394 = tmp392 - tmp393;
998 tmp395 = tmp391 - tmp394;
999 tmp491 = tmp391 + tmp394;
1000 }
1001 {
1002 fftw_real tmp538;
1003 fftw_real tmp539;
1004 fftw_real tmp388;
1005 fftw_real tmp389;
1006 ASSERT_ALIGNED_DOUBLE();
1007 tmp538 = tmp386 + tmp387;
1008 tmp539 = tmp392 + tmp393;
1009 tmp540 = tmp538 - tmp539;
1010 tmp594 = tmp538 + tmp539;
1011 tmp388 = tmp386 - tmp387;
1012 tmp389 = tmp241 - tmp246;
1013 tmp390 = tmp388 + tmp389;
1014 tmp490 = tmp388 - tmp389;
1015 }
1016 }
1017 {
1018 fftw_real tmp278;
1019 fftw_real tmp410;
1020 fftw_real tmp294;
1021 fftw_real tmp429;
1022 fftw_real tmp283;
1023 fftw_real tmp411;
1024 fftw_real tmp289;
1025 fftw_real tmp428;
1026 ASSERT_ALIGNED_DOUBLE();
1027 {
1028 fftw_real tmp275;
1029 fftw_real tmp277;
1030 fftw_real tmp274;
1031 fftw_real tmp276;
1032 ASSERT_ALIGNED_DOUBLE();
1033 tmp275 = X[iostride];
1034 tmp277 = Y[-30 * iostride];
1035 tmp274 = c_re(W[0]);
1036 tmp276 = c_im(W[0]);
1037 tmp278 = (tmp274 * tmp275) - (tmp276 * tmp277);
1038 tmp410 = (tmp276 * tmp275) + (tmp274 * tmp277);
1039 }
1040 {
1041 fftw_real tmp291;
1042 fftw_real tmp293;
1043 fftw_real tmp290;
1044 fftw_real tmp292;
1045 ASSERT_ALIGNED_DOUBLE();
1046 tmp291 = X[25 * iostride];
1047 tmp293 = Y[-6 * iostride];
1048 tmp290 = c_re(W[24]);
1049 tmp292 = c_im(W[24]);
1050 tmp294 = (tmp290 * tmp291) - (tmp292 * tmp293);
1051 tmp429 = (tmp292 * tmp291) + (tmp290 * tmp293);
1052 }
1053 {
1054 fftw_real tmp280;
1055 fftw_real tmp282;
1056 fftw_real tmp279;
1057 fftw_real tmp281;
1058 ASSERT_ALIGNED_DOUBLE();
1059 tmp280 = X[17 * iostride];
1060 tmp282 = Y[-14 * iostride];
1061 tmp279 = c_re(W[16]);
1062 tmp281 = c_im(W[16]);
1063 tmp283 = (tmp279 * tmp280) - (tmp281 * tmp282);
1064 tmp411 = (tmp281 * tmp280) + (tmp279 * tmp282);
1065 }
1066 {
1067 fftw_real tmp286;
1068 fftw_real tmp288;
1069 fftw_real tmp285;
1070 fftw_real tmp287;
1071 ASSERT_ALIGNED_DOUBLE();
1072 tmp286 = X[9 * iostride];
1073 tmp288 = Y[-22 * iostride];
1074 tmp285 = c_re(W[8]);
1075 tmp287 = c_im(W[8]);
1076 tmp289 = (tmp285 * tmp286) - (tmp287 * tmp288);
1077 tmp428 = (tmp287 * tmp286) + (tmp285 * tmp288);
1078 }
1079 {
1080 fftw_real tmp284;
1081 fftw_real tmp295;
1082 fftw_real tmp427;
1083 fftw_real tmp430;
1084 ASSERT_ALIGNED_DOUBLE();
1085 tmp284 = tmp278 + tmp283;
1086 tmp295 = tmp289 + tmp294;
1087 tmp296 = tmp284 + tmp295;
1088 tmp555 = tmp284 - tmp295;
1089 tmp427 = tmp278 - tmp283;
1090 tmp430 = tmp428 - tmp429;
1091 tmp431 = tmp427 - tmp430;
1092 tmp498 = tmp427 + tmp430;
1093 }
1094 {
1095 fftw_real tmp550;
1096 fftw_real tmp551;
1097 fftw_real tmp412;
1098 fftw_real tmp413;
1099 ASSERT_ALIGNED_DOUBLE();
1100 tmp550 = tmp410 + tmp411;
1101 tmp551 = tmp428 + tmp429;
1102 tmp552 = tmp550 - tmp551;
1103 tmp599 = tmp550 + tmp551;
1104 tmp412 = tmp410 - tmp411;
1105 tmp413 = tmp289 - tmp294;
1106 tmp414 = tmp412 + tmp413;
1107 tmp501 = tmp412 - tmp413;
1108 }
1109 }
1110 {
1111 fftw_real tmp253;
1112 fftw_real tmp397;
1113 fftw_real tmp269;
1114 fftw_real tmp404;
1115 fftw_real tmp258;
1116 fftw_real tmp398;
1117 fftw_real tmp264;
1118 fftw_real tmp403;
1119 ASSERT_ALIGNED_DOUBLE();
1120 {
1121 fftw_real tmp250;
1122 fftw_real tmp252;
1123 fftw_real tmp249;
1124 fftw_real tmp251;
1125 ASSERT_ALIGNED_DOUBLE();
1126 tmp250 = X[30 * iostride];
1127 tmp252 = Y[-iostride];
1128 tmp249 = c_re(W[29]);
1129 tmp251 = c_im(W[29]);
1130 tmp253 = (tmp249 * tmp250) - (tmp251 * tmp252);
1131 tmp397 = (tmp251 * tmp250) + (tmp249 * tmp252);
1132 }
1133 {
1134 fftw_real tmp266;
1135 fftw_real tmp268;
1136 fftw_real tmp265;
1137 fftw_real tmp267;
1138 ASSERT_ALIGNED_DOUBLE();
1139 tmp266 = X[22 * iostride];
1140 tmp268 = Y[-9 * iostride];
1141 tmp265 = c_re(W[21]);
1142 tmp267 = c_im(W[21]);
1143 tmp269 = (tmp265 * tmp266) - (tmp267 * tmp268);
1144 tmp404 = (tmp267 * tmp266) + (tmp265 * tmp268);
1145 }
1146 {
1147 fftw_real tmp255;
1148 fftw_real tmp257;
1149 fftw_real tmp254;
1150 fftw_real tmp256;
1151 ASSERT_ALIGNED_DOUBLE();
1152 tmp255 = X[14 * iostride];
1153 tmp257 = Y[-17 * iostride];
1154 tmp254 = c_re(W[13]);
1155 tmp256 = c_im(W[13]);
1156 tmp258 = (tmp254 * tmp255) - (tmp256 * tmp257);
1157 tmp398 = (tmp256 * tmp255) + (tmp254 * tmp257);
1158 }
1159 {
1160 fftw_real tmp261;
1161 fftw_real tmp263;
1162 fftw_real tmp260;
1163 fftw_real tmp262;
1164 ASSERT_ALIGNED_DOUBLE();
1165 tmp261 = X[6 * iostride];
1166 tmp263 = Y[-25 * iostride];
1167 tmp260 = c_re(W[5]);
1168 tmp262 = c_im(W[5]);
1169 tmp264 = (tmp260 * tmp261) - (tmp262 * tmp263);
1170 tmp403 = (tmp262 * tmp261) + (tmp260 * tmp263);
1171 }
1172 {
1173 fftw_real tmp259;
1174 fftw_real tmp270;
1175 fftw_real tmp402;
1176 fftw_real tmp405;
1177 ASSERT_ALIGNED_DOUBLE();
1178 tmp259 = tmp253 + tmp258;
1179 tmp270 = tmp264 + tmp269;
1180 tmp271 = tmp259 + tmp270;
1181 tmp543 = tmp259 - tmp270;
1182 tmp402 = tmp253 - tmp258;
1183 tmp405 = tmp403 - tmp404;
1184 tmp406 = tmp402 - tmp405;
1185 tmp494 = tmp402 + tmp405;
1186 }
1187 {
1188 fftw_real tmp544;
1189 fftw_real tmp545;
1190 fftw_real tmp399;
1191 fftw_real tmp400;
1192 ASSERT_ALIGNED_DOUBLE();
1193 tmp544 = tmp397 + tmp398;
1194 tmp545 = tmp403 + tmp404;
1195 tmp546 = tmp544 - tmp545;
1196 tmp595 = tmp544 + tmp545;
1197 tmp399 = tmp397 - tmp398;
1198 tmp400 = tmp264 - tmp269;
1199 tmp401 = tmp399 + tmp400;
1200 tmp493 = tmp399 - tmp400;
1201 }
1202 }
1203 {
1204 fftw_real tmp301;
1205 fftw_real tmp421;
1206 fftw_real tmp306;
1207 fftw_real tmp422;
1208 fftw_real tmp420;
1209 fftw_real tmp423;
1210 fftw_real tmp312;
1211 fftw_real tmp416;
1212 fftw_real tmp317;
1213 fftw_real tmp417;
1214 fftw_real tmp415;
1215 fftw_real tmp418;
1216 ASSERT_ALIGNED_DOUBLE();
1217 {
1218 fftw_real tmp298;
1219 fftw_real tmp300;
1220 fftw_real tmp297;
1221 fftw_real tmp299;
1222 ASSERT_ALIGNED_DOUBLE();
1223 tmp298 = X[5 * iostride];
1224 tmp300 = Y[-26 * iostride];
1225 tmp297 = c_re(W[4]);
1226 tmp299 = c_im(W[4]);
1227 tmp301 = (tmp297 * tmp298) - (tmp299 * tmp300);
1228 tmp421 = (tmp299 * tmp298) + (tmp297 * tmp300);
1229 }
1230 {
1231 fftw_real tmp303;
1232 fftw_real tmp305;
1233 fftw_real tmp302;
1234 fftw_real tmp304;
1235 ASSERT_ALIGNED_DOUBLE();
1236 tmp303 = X[21 * iostride];
1237 tmp305 = Y[-10 * iostride];
1238 tmp302 = c_re(W[20]);
1239 tmp304 = c_im(W[20]);
1240 tmp306 = (tmp302 * tmp303) - (tmp304 * tmp305);
1241 tmp422 = (tmp304 * tmp303) + (tmp302 * tmp305);
1242 }
1243 tmp420 = tmp301 - tmp306;
1244 tmp423 = tmp421 - tmp422;
1245 {
1246 fftw_real tmp309;
1247 fftw_real tmp311;
1248 fftw_real tmp308;
1249 fftw_real tmp310;
1250 ASSERT_ALIGNED_DOUBLE();
1251 tmp309 = X[29 * iostride];
1252 tmp311 = Y[-2 * iostride];
1253 tmp308 = c_re(W[28]);
1254 tmp310 = c_im(W[28]);
1255 tmp312 = (tmp308 * tmp309) - (tmp310 * tmp311);
1256 tmp416 = (tmp310 * tmp309) + (tmp308 * tmp311);
1257 }
1258 {
1259 fftw_real tmp314;
1260 fftw_real tmp316;
1261 fftw_real tmp313;
1262 fftw_real tmp315;
1263 ASSERT_ALIGNED_DOUBLE();
1264 tmp314 = X[13 * iostride];
1265 tmp316 = Y[-18 * iostride];
1266 tmp313 = c_re(W[12]);
1267 tmp315 = c_im(W[12]);
1268 tmp317 = (tmp313 * tmp314) - (tmp315 * tmp316);
1269 tmp417 = (tmp315 * tmp314) + (tmp313 * tmp316);
1270 }
1271 tmp415 = tmp312 - tmp317;
1272 tmp418 = tmp416 - tmp417;
1273 {
1274 fftw_real tmp307;
1275 fftw_real tmp318;
1276 fftw_real tmp556;
1277 fftw_real tmp557;
1278 ASSERT_ALIGNED_DOUBLE();
1279 tmp307 = tmp301 + tmp306;
1280 tmp318 = tmp312 + tmp317;
1281 tmp319 = tmp307 + tmp318;
1282 tmp553 = tmp318 - tmp307;
1283 tmp556 = tmp421 + tmp422;
1284 tmp557 = tmp416 + tmp417;
1285 tmp558 = tmp556 - tmp557;
1286 tmp600 = tmp556 + tmp557;
1287 }
1288 {
1289 fftw_real tmp419;
1290 fftw_real tmp424;
1291 fftw_real tmp432;
1292 fftw_real tmp433;
1293 ASSERT_ALIGNED_DOUBLE();
1294 tmp419 = tmp415 - tmp418;
1295 tmp424 = tmp420 + tmp423;
1296 tmp425 = K707106781 * (tmp419 - tmp424);
1297 tmp499 = K707106781 * (tmp424 + tmp419);
1298 tmp432 = tmp423 - tmp420;
1299 tmp433 = tmp415 + tmp418;
1300 tmp434 = K707106781 * (tmp432 - tmp433);
1301 tmp502 = K707106781 * (tmp432 + tmp433);
1302 }
1303 }
1304 {
1305 fftw_real tmp273;
1306 fftw_real tmp613;
1307 fftw_real tmp625;
1308 fftw_real tmp627;
1309 fftw_real tmp368;
1310 fftw_real tmp628;
1311 fftw_real tmp616;
1312 fftw_real tmp626;
1313 ASSERT_ALIGNED_DOUBLE();
1314 {
1315 fftw_real tmp225;
1316 fftw_real tmp272;
1317 fftw_real tmp617;
1318 fftw_real tmp624;
1319 ASSERT_ALIGNED_DOUBLE();
1320 tmp225 = tmp201 + tmp224;
1321 tmp272 = tmp248 + tmp271;
1322 tmp273 = tmp225 + tmp272;
1323 tmp613 = tmp225 - tmp272;
1324 tmp617 = tmp594 + tmp595;
1325 tmp624 = tmp618 + tmp623;
1326 tmp625 = tmp617 + tmp624;
1327 tmp627 = tmp624 - tmp617;
1328 }
1329 {
1330 fftw_real tmp320;
1331 fftw_real tmp367;
1332 fftw_real tmp614;
1333 fftw_real tmp615;
1334 ASSERT_ALIGNED_DOUBLE();
1335 tmp320 = tmp296 + tmp319;
1336 tmp367 = tmp343 + tmp366;
1337 tmp368 = tmp320 + tmp367;
1338 tmp628 = tmp367 - tmp320;
1339 tmp614 = tmp599 + tmp600;
1340 tmp615 = tmp604 + tmp605;
1341 tmp616 = tmp614 - tmp615;
1342 tmp626 = tmp614 + tmp615;
1343 }
1344 Y[-16 * iostride] = tmp273 - tmp368;
1345 X[0] = tmp273 + tmp368;
1346 Y[-24 * iostride] = tmp613 - tmp616;
1347 X[8 * iostride] = tmp613 + tmp616;
1348 X[16 * iostride] = -(tmp625 - tmp626);
1349 Y[0] = tmp626 + tmp625;
1350 X[24 * iostride] = -(tmp627 - tmp628);
1351 Y[-8 * iostride] = tmp628 + tmp627;
1352 }
1353 {
1354 fftw_real tmp597;
1355 fftw_real tmp609;
1356 fftw_real tmp631;
1357 fftw_real tmp633;
1358 fftw_real tmp602;
1359 fftw_real tmp610;
1360 fftw_real tmp607;
1361 fftw_real tmp611;
1362 ASSERT_ALIGNED_DOUBLE();
1363 {
1364 fftw_real tmp593;
1365 fftw_real tmp596;
1366 fftw_real tmp629;
1367 fftw_real tmp630;
1368 ASSERT_ALIGNED_DOUBLE();
1369 tmp593 = tmp201 - tmp224;
1370 tmp596 = tmp594 - tmp595;
1371 tmp597 = tmp593 + tmp596;
1372 tmp609 = tmp593 - tmp596;
1373 tmp629 = tmp271 - tmp248;
1374 tmp630 = tmp623 - tmp618;
1375 tmp631 = tmp629 + tmp630;
1376 tmp633 = tmp630 - tmp629;
1377 }
1378 {
1379 fftw_real tmp598;
1380 fftw_real tmp601;
1381 fftw_real tmp603;
1382 fftw_real tmp606;
1383 ASSERT_ALIGNED_DOUBLE();
1384 tmp598 = tmp296 - tmp319;
1385 tmp601 = tmp599 - tmp600;
1386 tmp602 = tmp598 + tmp601;
1387 tmp610 = tmp601 - tmp598;
1388 tmp603 = tmp343 - tmp366;
1389 tmp606 = tmp604 - tmp605;
1390 tmp607 = tmp603 - tmp606;
1391 tmp611 = tmp603 + tmp606;
1392 }
1393 {
1394 fftw_real tmp608;
1395 fftw_real tmp634;
1396 fftw_real tmp612;
1397 fftw_real tmp632;
1398 ASSERT_ALIGNED_DOUBLE();
1399 tmp608 = K707106781 * (tmp602 + tmp607);
1400 Y[-20 * iostride] = tmp597 - tmp608;
1401 X[4 * iostride] = tmp597 + tmp608;
1402 tmp634 = K707106781 * (tmp607 - tmp602);
1403 X[28 * iostride] = -(tmp633 - tmp634);
1404 Y[-12 * iostride] = tmp634 + tmp633;
1405 tmp612 = K707106781 * (tmp610 - tmp611);
1406 Y[-28 * iostride] = tmp609 - tmp612;
1407 X[12 * iostride] = tmp609 + tmp612;
1408 tmp632 = K707106781 * (tmp610 + tmp611);
1409 X[20 * iostride] = -(tmp631 - tmp632);
1410 Y[-4 * iostride] = tmp632 + tmp631;
1411 }
1412 }
1413 {
1414 fftw_real tmp537;
1415 fftw_real tmp577;
1416 fftw_real tmp548;
1417 fftw_real tmp635;
1418 fftw_real tmp580;
1419 fftw_real tmp643;
1420 fftw_real tmp560;
1421 fftw_real tmp574;
1422 fftw_real tmp638;
1423 fftw_real tmp644;
1424 fftw_real tmp584;
1425 fftw_real tmp590;
1426 fftw_real tmp571;
1427 fftw_real tmp575;
1428 fftw_real tmp587;
1429 fftw_real tmp591;
1430 ASSERT_ALIGNED_DOUBLE();
1431 {
1432 fftw_real tmp542;
1433 fftw_real tmp547;
1434 fftw_real tmp582;
1435 fftw_real tmp583;
1436 ASSERT_ALIGNED_DOUBLE();
1437 tmp537 = tmp533 - tmp536;
1438 tmp577 = tmp533 + tmp536;
1439 tmp542 = tmp540 - tmp541;
1440 tmp547 = tmp543 + tmp546;
1441 tmp548 = K707106781 * (tmp542 - tmp547);
1442 tmp635 = K707106781 * (tmp542 + tmp547);
1443 {
1444 fftw_real tmp578;
1445 fftw_real tmp579;
1446 fftw_real tmp554;
1447 fftw_real tmp559;
1448 ASSERT_ALIGNED_DOUBLE();
1449 tmp578 = tmp541 + tmp540;
1450 tmp579 = tmp543 - tmp546;
1451 tmp580 = K707106781 * (tmp578 + tmp579);
1452 tmp643 = K707106781 * (tmp579 - tmp578);
1453 tmp554 = tmp552 - tmp553;
1454 tmp559 = tmp555 - tmp558;
1455 tmp560 = (K923879532 * tmp554) + (K382683432 * tmp559);
1456 tmp574 = (K382683432 * tmp554) - (K923879532 * tmp559);
1457 }
1458 tmp638 = tmp636 + tmp637;
1459 tmp644 = tmp637 - tmp636;
1460 tmp582 = tmp552 + tmp553;
1461 tmp583 = tmp555 + tmp558;
1462 tmp584 = (K382683432 * tmp582) + (K923879532 * tmp583);
1463 tmp590 = (K923879532 * tmp582) - (K382683432 * tmp583);
1464 {
1465 fftw_real tmp565;
1466 fftw_real tmp570;
1467 fftw_real tmp585;
1468 fftw_real tmp586;
1469 ASSERT_ALIGNED_DOUBLE();
1470 tmp565 = tmp561 - tmp564;
1471 tmp570 = tmp568 - tmp569;
1472 tmp571 = (K382683432 * tmp565) - (K923879532 * tmp570);
1473 tmp575 = (K382683432 * tmp570) + (K923879532 * tmp565);
1474 tmp585 = tmp561 + tmp564;
1475 tmp586 = tmp568 + tmp569;
1476 tmp587 = (K923879532 * tmp585) - (K382683432 * tmp586);
1477 tmp591 = (K923879532 * tmp586) + (K382683432 * tmp585);
1478 }
1479 }
1480 {
1481 fftw_real tmp549;
1482 fftw_real tmp572;
1483 fftw_real tmp573;
1484 fftw_real tmp576;
1485 ASSERT_ALIGNED_DOUBLE();
1486 tmp549 = tmp537 + tmp548;
1487 tmp572 = tmp560 + tmp571;
1488 Y[-22 * iostride] = tmp549 - tmp572;
1489 X[6 * iostride] = tmp549 + tmp572;
1490 tmp573 = tmp537 - tmp548;
1491 tmp576 = tmp574 - tmp575;
1492 Y[-30 * iostride] = tmp573 - tmp576;
1493 X[14 * iostride] = tmp573 + tmp576;
1494 }
1495 {
1496 fftw_real tmp645;
1497 fftw_real tmp646;
1498 fftw_real tmp647;
1499 fftw_real tmp648;
1500 ASSERT_ALIGNED_DOUBLE();
1501 tmp645 = tmp643 + tmp644;
1502 tmp646 = tmp574 + tmp575;
1503 X[22 * iostride] = -(tmp645 - tmp646);
1504 Y[-6 * iostride] = tmp646 + tmp645;
1505 tmp647 = tmp644 - tmp643;
1506 tmp648 = tmp571 - tmp560;
1507 X[30 * iostride] = -(tmp647 - tmp648);
1508 Y[-14 * iostride] = tmp648 + tmp647;
1509 }
1510 {
1511 fftw_real tmp581;
1512 fftw_real tmp588;
1513 fftw_real tmp589;
1514 fftw_real tmp592;
1515 ASSERT_ALIGNED_DOUBLE();
1516 tmp581 = tmp577 + tmp580;
1517 tmp588 = tmp584 + tmp587;
1518 Y[-18 * iostride] = tmp581 - tmp588;
1519 X[2 * iostride] = tmp581 + tmp588;
1520 tmp589 = tmp577 - tmp580;
1521 tmp592 = tmp590 - tmp591;
1522 Y[-26 * iostride] = tmp589 - tmp592;
1523 X[10 * iostride] = tmp589 + tmp592;
1524 }
1525 {
1526 fftw_real tmp639;
1527 fftw_real tmp640;
1528 fftw_real tmp641;
1529 fftw_real tmp642;
1530 ASSERT_ALIGNED_DOUBLE();
1531 tmp639 = tmp635 + tmp638;
1532 tmp640 = tmp590 + tmp591;
1533 X[18 * iostride] = -(tmp639 - tmp640);
1534 Y[-2 * iostride] = tmp640 + tmp639;
1535 tmp641 = tmp638 - tmp635;
1536 tmp642 = tmp587 - tmp584;
1537 X[26 * iostride] = -(tmp641 - tmp642);
1538 Y[-10 * iostride] = tmp642 + tmp641;
1539 }
1540 }
1541 {
1542 fftw_real tmp489;
1543 fftw_real tmp517;
1544 fftw_real tmp520;
1545 fftw_real tmp659;
1546 fftw_real tmp654;
1547 fftw_real tmp660;
1548 fftw_real tmp496;
1549 fftw_real tmp649;
1550 fftw_real tmp504;
1551 fftw_real tmp514;
1552 fftw_real tmp524;
1553 fftw_real tmp530;
1554 fftw_real tmp511;
1555 fftw_real tmp515;
1556 fftw_real tmp527;
1557 fftw_real tmp531;
1558 ASSERT_ALIGNED_DOUBLE();
1559 {
1560 fftw_real tmp488;
1561 fftw_real tmp518;
1562 fftw_real tmp519;
1563 fftw_real tmp650;
1564 fftw_real tmp492;
1565 fftw_real tmp495;
1566 ASSERT_ALIGNED_DOUBLE();
1567 tmp488 = K707106781 * (tmp486 + tmp487);
1568 tmp489 = tmp485 - tmp488;
1569 tmp517 = tmp485 + tmp488;
1570 tmp518 = (K382683432 * tmp490) + (K923879532 * tmp491);
1571 tmp519 = (K923879532 * tmp494) - (K382683432 * tmp493);
1572 tmp520 = tmp518 + tmp519;
1573 tmp659 = tmp519 - tmp518;
1574 tmp650 = K707106781 * (tmp378 + tmp383);
1575 tmp654 = tmp650 + tmp653;
1576 tmp660 = tmp653 - tmp650;
1577 tmp492 = (K923879532 * tmp490) - (K382683432 * tmp491);
1578 tmp495 = (K923879532 * tmp493) + (K382683432 * tmp494);
1579 tmp496 = tmp492 - tmp495;
1580 tmp649 = tmp492 + tmp495;
1581 }
1582 {
1583 fftw_real tmp500;
1584 fftw_real tmp503;
1585 fftw_real tmp522;
1586 fftw_real tmp523;
1587 ASSERT_ALIGNED_DOUBLE();
1588 tmp500 = tmp498 - tmp499;
1589 tmp503 = tmp501 - tmp502;
1590 tmp504 = (K555570233 * tmp500) + (K831469612 * tmp503);
1591 tmp514 = (K555570233 * tmp503) - (K831469612 * tmp500);
1592 tmp522 = tmp498 + tmp499;
1593 tmp523 = tmp501 + tmp502;
1594 tmp524 = (K980785280 * tmp522) + (K195090322 * tmp523);
1595 tmp530 = (K980785280 * tmp523) - (K195090322 * tmp522);
1596 }
1597 {
1598 fftw_real tmp507;
1599 fftw_real tmp510;
1600 fftw_real tmp525;
1601 fftw_real tmp526;
1602 ASSERT_ALIGNED_DOUBLE();
1603 tmp507 = tmp505 - tmp506;
1604 tmp510 = tmp508 - tmp509;
1605 tmp511 = (K555570233 * tmp507) - (K831469612 * tmp510);
1606 tmp515 = (K831469612 * tmp507) + (K555570233 * tmp510);
1607 tmp525 = tmp505 + tmp506;
1608 tmp526 = tmp508 + tmp509;
1609 tmp527 = (K980785280 * tmp525) - (K195090322 * tmp526);
1610 tmp531 = (K195090322 * tmp525) + (K980785280 * tmp526);
1611 }
1612 {
1613 fftw_real tmp497;
1614 fftw_real tmp512;
1615 fftw_real tmp513;
1616 fftw_real tmp516;
1617 ASSERT_ALIGNED_DOUBLE();
1618 tmp497 = tmp489 + tmp496;
1619 tmp512 = tmp504 + tmp511;
1620 Y[-21 * iostride] = tmp497 - tmp512;
1621 X[5 * iostride] = tmp497 + tmp512;
1622 tmp513 = tmp489 - tmp496;
1623 tmp516 = tmp514 - tmp515;
1624 Y[-29 * iostride] = tmp513 - tmp516;
1625 X[13 * iostride] = tmp513 + tmp516;
1626 }
1627 {
1628 fftw_real tmp661;
1629 fftw_real tmp662;
1630 fftw_real tmp663;
1631 fftw_real tmp664;
1632 ASSERT_ALIGNED_DOUBLE();
1633 tmp661 = tmp659 + tmp660;
1634 tmp662 = tmp514 + tmp515;
1635 X[21 * iostride] = -(tmp661 - tmp662);
1636 Y[-5 * iostride] = tmp662 + tmp661;
1637 tmp663 = tmp660 - tmp659;
1638 tmp664 = tmp511 - tmp504;
1639 X[29 * iostride] = -(tmp663 - tmp664);
1640 Y[-13 * iostride] = tmp664 + tmp663;
1641 }
1642 {
1643 fftw_real tmp521;
1644 fftw_real tmp528;
1645 fftw_real tmp529;
1646 fftw_real tmp532;
1647 ASSERT_ALIGNED_DOUBLE();
1648 tmp521 = tmp517 + tmp520;
1649 tmp528 = tmp524 + tmp527;
1650 Y[-17 * iostride] = tmp521 - tmp528;
1651 X[iostride] = tmp521 + tmp528;
1652 tmp529 = tmp517 - tmp520;
1653 tmp532 = tmp530 - tmp531;
1654 Y[-25 * iostride] = tmp529 - tmp532;
1655 X[9 * iostride] = tmp529 + tmp532;
1656 }
1657 {
1658 fftw_real tmp655;
1659 fftw_real tmp656;
1660 fftw_real tmp657;
1661 fftw_real tmp658;
1662 ASSERT_ALIGNED_DOUBLE();
1663 tmp655 = tmp649 + tmp654;
1664 tmp656 = tmp530 + tmp531;
1665 X[17 * iostride] = -(tmp655 - tmp656);
1666 Y[-iostride] = tmp656 + tmp655;
1667 tmp657 = tmp654 - tmp649;
1668 tmp658 = tmp527 - tmp524;
1669 X[25 * iostride] = -(tmp657 - tmp658);
1670 Y[-9 * iostride] = tmp658 + tmp657;
1671 }
1672 }
1673 {
1674 fftw_real tmp385;
1675 fftw_real tmp469;
1676 fftw_real tmp472;
1677 fftw_real tmp673;
1678 fftw_real tmp668;
1679 fftw_real tmp674;
1680 fftw_real tmp408;
1681 fftw_real tmp665;
1682 fftw_real tmp436;
1683 fftw_real tmp466;
1684 fftw_real tmp476;
1685 fftw_real tmp482;
1686 fftw_real tmp463;
1687 fftw_real tmp467;
1688 fftw_real tmp479;
1689 fftw_real tmp483;
1690 ASSERT_ALIGNED_DOUBLE();
1691 {
1692 fftw_real tmp384;
1693 fftw_real tmp470;
1694 fftw_real tmp471;
1695 fftw_real tmp666;
1696 fftw_real tmp396;
1697 fftw_real tmp407;
1698 ASSERT_ALIGNED_DOUBLE();
1699 tmp384 = K707106781 * (tmp378 - tmp383);
1700 tmp385 = tmp373 - tmp384;
1701 tmp469 = tmp373 + tmp384;
1702 tmp470 = (K923879532 * tmp390) + (K382683432 * tmp395);
1703 tmp471 = (K382683432 * tmp406) - (K923879532 * tmp401);
1704 tmp472 = tmp470 + tmp471;
1705 tmp673 = tmp471 - tmp470;
1706 tmp666 = K707106781 * (tmp487 - tmp486);
1707 tmp668 = tmp666 + tmp667;
1708 tmp674 = tmp667 - tmp666;
1709 tmp396 = (K382683432 * tmp390) - (K923879532 * tmp395);
1710 tmp407 = (K382683432 * tmp401) + (K923879532 * tmp406);
1711 tmp408 = tmp396 - tmp407;
1712 tmp665 = tmp396 + tmp407;
1713 }
1714 {
1715 fftw_real tmp426;
1716 fftw_real tmp435;
1717 fftw_real tmp474;
1718 fftw_real tmp475;
1719 ASSERT_ALIGNED_DOUBLE();
1720 tmp426 = tmp414 - tmp425;
1721 tmp435 = tmp431 - tmp434;
1722 tmp436 = (K980785280 * tmp426) + (K195090322 * tmp435);
1723 tmp466 = (K195090322 * tmp426) - (K980785280 * tmp435);
1724 tmp474 = tmp414 + tmp425;
1725 tmp475 = tmp431 + tmp434;
1726 tmp476 = (K555570233 * tmp474) + (K831469612 * tmp475);
1727 tmp482 = (K831469612 * tmp474) - (K555570233 * tmp475);
1728 }
1729 {
1730 fftw_real tmp453;
1731 fftw_real tmp462;
1732 fftw_real tmp477;
1733 fftw_real tmp478;
1734 ASSERT_ALIGNED_DOUBLE();
1735 tmp453 = tmp441 - tmp452;
1736 tmp462 = tmp458 - tmp461;
1737 tmp463 = (K195090322 * tmp453) - (K980785280 * tmp462);
1738 tmp467 = (K195090322 * tmp462) + (K980785280 * tmp453);
1739 tmp477 = tmp441 + tmp452;
1740 tmp478 = tmp458 + tmp461;
1741 tmp479 = (K831469612 * tmp477) - (K555570233 * tmp478);
1742 tmp483 = (K831469612 * tmp478) + (K555570233 * tmp477);
1743 }
1744 {
1745 fftw_real tmp409;
1746 fftw_real tmp464;
1747 fftw_real tmp465;
1748 fftw_real tmp468;
1749 ASSERT_ALIGNED_DOUBLE();
1750 tmp409 = tmp385 + tmp408;
1751 tmp464 = tmp436 + tmp463;
1752 Y[-23 * iostride] = tmp409 - tmp464;
1753 X[7 * iostride] = tmp409 + tmp464;
1754 tmp465 = tmp385 - tmp408;
1755 tmp468 = tmp466 - tmp467;
1756 Y[-31 * iostride] = tmp465 - tmp468;
1757 X[15 * iostride] = tmp465 + tmp468;
1758 }
1759 {
1760 fftw_real tmp675;
1761 fftw_real tmp676;
1762 fftw_real tmp677;
1763 fftw_real tmp678;
1764 ASSERT_ALIGNED_DOUBLE();
1765 tmp675 = tmp673 + tmp674;
1766 tmp676 = tmp466 + tmp467;
1767 X[23 * iostride] = -(tmp675 - tmp676);
1768 Y[-7 * iostride] = tmp676 + tmp675;
1769 tmp677 = tmp674 - tmp673;
1770 tmp678 = tmp463 - tmp436;
1771 X[31 * iostride] = -(tmp677 - tmp678);
1772 Y[-15 * iostride] = tmp678 + tmp677;
1773 }
1774 {
1775 fftw_real tmp473;
1776 fftw_real tmp480;
1777 fftw_real tmp481;
1778 fftw_real tmp484;
1779 ASSERT_ALIGNED_DOUBLE();
1780 tmp473 = tmp469 + tmp472;
1781 tmp480 = tmp476 + tmp479;
1782 Y[-19 * iostride] = tmp473 - tmp480;
1783 X[3 * iostride] = tmp473 + tmp480;
1784 tmp481 = tmp469 - tmp472;
1785 tmp484 = tmp482 - tmp483;
1786 Y[-27 * iostride] = tmp481 - tmp484;
1787 X[11 * iostride] = tmp481 + tmp484;
1788 }
1789 {
1790 fftw_real tmp669;
1791 fftw_real tmp670;
1792 fftw_real tmp671;
1793 fftw_real tmp672;
1794 ASSERT_ALIGNED_DOUBLE();
1795 tmp669 = tmp665 + tmp668;
1796 tmp670 = tmp482 + tmp483;
1797 X[19 * iostride] = -(tmp669 - tmp670);
1798 Y[-3 * iostride] = tmp670 + tmp669;
1799 tmp671 = tmp668 - tmp665;
1800 tmp672 = tmp479 - tmp476;
1801 X[27 * iostride] = -(tmp671 - tmp672);
1802 Y[-11 * iostride] = tmp672 + tmp671;
1803 }
1804 }
1805 }
1806 if (i == m) {
1807 fftw_real tmp5;
1808 fftw_real tmp105;
1809 fftw_real tmp158;
1810 fftw_real tmp171;
1811 fftw_real tmp12;
1812 fftw_real tmp170;
1813 fftw_real tmp108;
1814 fftw_real tmp155;
1815 fftw_real tmp74;
1816 fftw_real tmp97;
1817 fftw_real tmp130;
1818 fftw_real tmp146;
1819 fftw_real tmp82;
1820 fftw_real tmp98;
1821 fftw_real tmp127;
1822 fftw_real tmp145;
1823 fftw_real tmp24;
1824 fftw_real tmp90;
1825 fftw_real tmp115;
1826 fftw_real tmp138;
1827 fftw_real tmp35;
1828 fftw_real tmp91;
1829 fftw_real tmp112;
1830 fftw_real tmp139;
1831 fftw_real tmp51;
1832 fftw_real tmp94;
1833 fftw_real tmp123;
1834 fftw_real tmp143;
1835 fftw_real tmp59;
1836 fftw_real tmp95;
1837 fftw_real tmp120;
1838 fftw_real tmp142;
1839 ASSERT_ALIGNED_DOUBLE();
1840 {
1841 fftw_real tmp1;
1842 fftw_real tmp157;
1843 fftw_real tmp4;
1844 fftw_real tmp156;
1845 fftw_real tmp2;
1846 fftw_real tmp3;
1847 ASSERT_ALIGNED_DOUBLE();
1848 tmp1 = X[0];
1849 tmp157 = X[16 * iostride];
1850 tmp2 = X[8 * iostride];
1851 tmp3 = X[24 * iostride];
1852 tmp4 = K707106781 * (tmp2 - tmp3);
1853 tmp156 = K707106781 * (tmp2 + tmp3);
1854 tmp5 = tmp1 + tmp4;
1855 tmp105 = tmp1 - tmp4;
1856 tmp158 = tmp156 + tmp157;
1857 tmp171 = tmp157 - tmp156;
1858 }
1859 {
1860 fftw_real tmp8;
1861 fftw_real tmp106;
1862 fftw_real tmp11;
1863 fftw_real tmp107;
1864 ASSERT_ALIGNED_DOUBLE();
1865 {
1866 fftw_real tmp6;
1867 fftw_real tmp7;
1868 fftw_real tmp9;
1869 fftw_real tmp10;
1870 ASSERT_ALIGNED_DOUBLE();
1871 tmp6 = X[4 * iostride];
1872 tmp7 = X[20 * iostride];
1873 tmp8 = (K923879532 * tmp6) - (K382683432 * tmp7);
1874 tmp106 = (K382683432 * tmp6) + (K923879532 * tmp7);
1875 tmp9 = X[12 * iostride];
1876 tmp10 = X[28 * iostride];
1877 tmp11 = (K382683432 * tmp9) - (K923879532 * tmp10);
1878 tmp107 = (K923879532 * tmp9) + (K382683432 * tmp10);
1879 }
1880 tmp12 = tmp8 + tmp11;
1881 tmp170 = tmp11 - tmp8;
1882 tmp108 = tmp106 - tmp107;
1883 tmp155 = tmp106 + tmp107;
1884 }
1885 {
1886 fftw_real tmp65;
1887 fftw_real tmp80;
1888 fftw_real tmp63;
1889 fftw_real tmp78;
1890 fftw_real tmp69;
1891 fftw_real tmp75;
1892 fftw_real tmp72;
1893 fftw_real tmp76;
1894 ASSERT_ALIGNED_DOUBLE();
1895 {
1896 fftw_real tmp64;
1897 fftw_real tmp79;
1898 fftw_real tmp61;
1899 fftw_real tmp62;
1900 ASSERT_ALIGNED_DOUBLE();
1901 tmp64 = X[31 * iostride];
1902 tmp65 = K2_000000000 * tmp64;
1903 tmp79 = X[15 * iostride];
1904 tmp80 = K2_000000000 * tmp79;
1905 tmp61 = X[7 * iostride];
1906 tmp62 = X[23 * iostride];
1907 tmp63 = K1_414213562 * (tmp61 - tmp62);
1908 tmp78 = K1_414213562 * (tmp61 + tmp62);
1909 {
1910 fftw_real tmp67;
1911 fftw_real tmp68;
1912 fftw_real tmp70;
1913 fftw_real tmp71;
1914 ASSERT_ALIGNED_DOUBLE();
1915 tmp67 = X[3 * iostride];
1916 tmp68 = X[19 * iostride];
1917 tmp69 = (K1_847759065 * tmp67) - (K765366864 * tmp68);
1918 tmp75 = (K765366864 * tmp67) + (K1_847759065 * tmp68);
1919 tmp70 = X[11 * iostride];
1920 tmp71 = X[27 * iostride];
1921 tmp72 = (K765366864 * tmp70) - (K1_847759065 * tmp71);
1922 tmp76 = (K1_847759065 * tmp70) + (K765366864 * tmp71);
1923 }
1924 }
1925 {
1926 fftw_real tmp66;
1927 fftw_real tmp73;
1928 fftw_real tmp128;
1929 fftw_real tmp129;
1930 ASSERT_ALIGNED_DOUBLE();
1931 tmp66 = tmp63 - tmp65;
1932 tmp73 = tmp69 + tmp72;
1933 tmp74 = tmp66 + tmp73;
1934 tmp97 = tmp66 - tmp73;
1935 tmp128 = tmp72 - tmp69;
1936 tmp129 = tmp80 - tmp78;
1937 tmp130 = tmp128 - tmp129;
1938 tmp146 = tmp128 + tmp129;
1939 }
1940 {
1941 fftw_real tmp77;
1942 fftw_real tmp81;
1943 fftw_real tmp125;
1944 fftw_real tmp126;
1945 ASSERT_ALIGNED_DOUBLE();
1946 tmp77 = tmp75 + tmp76;
1947 tmp81 = tmp78 + tmp80;
1948 tmp82 = tmp77 + tmp81;
1949 tmp98 = tmp81 - tmp77;
1950 tmp125 = tmp63 + tmp65;
1951 tmp126 = tmp75 - tmp76;
1952 tmp127 = tmp125 + tmp126;
1953 tmp145 = tmp126 - tmp125;
1954 }
1955 }
1956 {
1957 fftw_real tmp15;
1958 fftw_real tmp22;
1959 fftw_real tmp18;
1960 fftw_real tmp20;
1961 ASSERT_ALIGNED_DOUBLE();
1962 {
1963 fftw_real tmp14;
1964 fftw_real tmp21;
1965 fftw_real tmp16;
1966 fftw_real tmp17;
1967 ASSERT_ALIGNED_DOUBLE();
1968 tmp14 = X[2 * iostride];
1969 tmp15 = K2_000000000 * tmp14;
1970 tmp21 = X[18 * iostride];
1971 tmp22 = K2_000000000 * tmp21;
1972 tmp16 = X[10 * iostride];
1973 tmp17 = X[26 * iostride];
1974 tmp18 = K1_414213562 * (tmp16 - tmp17);
1975 tmp20 = K1_414213562 * (tmp16 + tmp17);
1976 }
1977 {
1978 fftw_real tmp19;
1979 fftw_real tmp23;
1980 fftw_real tmp113;
1981 fftw_real tmp114;
1982 ASSERT_ALIGNED_DOUBLE();
1983 tmp19 = tmp15 + tmp18;
1984 tmp23 = tmp20 + tmp22;
1985 tmp24 = (K490392640 * tmp19) - (K097545161 * tmp23);
1986 tmp90 = (K097545161 * tmp19) + (K490392640 * tmp23);
1987 tmp113 = tmp22 - tmp20;
1988 tmp114 = tmp15 - tmp18;
1989 tmp115 = (K415734806 * tmp113) - (K277785116 * tmp114);
1990 tmp138 = (K415734806 * tmp114) + (K277785116 * tmp113);
1991 }
1992 }
1993 {
1994 fftw_real tmp29;
1995 fftw_real tmp33;
1996 fftw_real tmp27;
1997 fftw_real tmp31;
1998 ASSERT_ALIGNED_DOUBLE();
1999 {
2000 fftw_real tmp28;
2001 fftw_real tmp32;
2002 fftw_real tmp25;
2003 fftw_real tmp26;
2004 ASSERT_ALIGNED_DOUBLE();
2005 tmp28 = X[30 * iostride];
2006 tmp29 = K2_000000000 * tmp28;
2007 tmp32 = X[14 * iostride];
2008 tmp33 = K2_000000000 * tmp32;
2009 tmp25 = X[6 * iostride];
2010 tmp26 = X[22 * iostride];
2011 tmp27 = K1_414213562 * (tmp25 - tmp26);
2012 tmp31 = K1_414213562 * (tmp25 + tmp26);
2013 }
2014 {
2015 fftw_real tmp30;
2016 fftw_real tmp34;
2017 fftw_real tmp110;
2018 fftw_real tmp111;
2019 ASSERT_ALIGNED_DOUBLE();
2020 tmp30 = tmp27 - tmp29;
2021 tmp34 = tmp31 + tmp33;
2022 tmp35 = (K490392640 * tmp30) + (K097545161 * tmp34);
2023 tmp91 = (K097545161 * tmp30) - (K490392640 * tmp34);
2024 tmp110 = tmp33 - tmp31;
2025 tmp111 = tmp27 + tmp29;
2026 tmp112 = (K415734806 * tmp110) - (K277785116 * tmp111);
2027 tmp139 = (K415734806 * tmp111) + (K277785116 * tmp110);
2028 }
2029 }
2030 {
2031 fftw_real tmp39;
2032 fftw_real tmp57;
2033 fftw_real tmp42;
2034 fftw_real tmp55;
2035 fftw_real tmp46;
2036 fftw_real tmp52;
2037 fftw_real tmp49;
2038 fftw_real tmp53;
2039 ASSERT_ALIGNED_DOUBLE();
2040 {
2041 fftw_real tmp38;
2042 fftw_real tmp56;
2043 fftw_real tmp40;
2044 fftw_real tmp41;
2045 ASSERT_ALIGNED_DOUBLE();
2046 tmp38 = X[iostride];
2047 tmp39 = K2_000000000 * tmp38;
2048 tmp56 = X[17 * iostride];
2049 tmp57 = K2_000000000 * tmp56;
2050 tmp40 = X[9 * iostride];
2051 tmp41 = X[25 * iostride];
2052 tmp42 = K1_414213562 * (tmp40 - tmp41);
2053 tmp55 = K1_414213562 * (tmp40 + tmp41);
2054 {
2055 fftw_real tmp44;
2056 fftw_real tmp45;
2057 fftw_real tmp47;
2058 fftw_real tmp48;
2059 ASSERT_ALIGNED_DOUBLE();
2060 tmp44 = X[5 * iostride];
2061 tmp45 = X[21 * iostride];
2062 tmp46 = (K1_847759065 * tmp44) - (K765366864 * tmp45);
2063 tmp52 = (K765366864 * tmp44) + (K1_847759065 * tmp45);
2064 tmp47 = X[13 * iostride];
2065 tmp48 = X[29 * iostride];
2066 tmp49 = (K765366864 * tmp47) - (K1_847759065 * tmp48);
2067 tmp53 = (K1_847759065 * tmp47) + (K765366864 * tmp48);
2068 }
2069 }
2070 {
2071 fftw_real tmp43;
2072 fftw_real tmp50;
2073 fftw_real tmp121;
2074 fftw_real tmp122;
2075 ASSERT_ALIGNED_DOUBLE();
2076 tmp43 = tmp39 + tmp42;
2077 tmp50 = tmp46 + tmp49;
2078 tmp51 = tmp43 + tmp50;
2079 tmp94 = tmp43 - tmp50;
2080 tmp121 = tmp49 - tmp46;
2081 tmp122 = tmp57 - tmp55;
2082 tmp123 = tmp121 - tmp122;
2083 tmp143 = tmp121 + tmp122;
2084 }
2085 {
2086 fftw_real tmp54;
2087 fftw_real tmp58;
2088 fftw_real tmp118;
2089 fftw_real tmp119;
2090 ASSERT_ALIGNED_DOUBLE();
2091 tmp54 = tmp52 + tmp53;
2092 tmp58 = tmp55 + tmp57;
2093 tmp59 = tmp54 + tmp58;
2094 tmp95 = tmp58 - tmp54;
2095 tmp118 = tmp39 - tmp42;
2096 tmp119 = tmp52 - tmp53;
2097 tmp120 = tmp118 - tmp119;
2098 tmp142 = tmp118 + tmp119;
2099 }
2100 }
2101 {
2102 fftw_real tmp37;
2103 fftw_real tmp85;
2104 fftw_real tmp160;
2105 fftw_real tmp162;
2106 fftw_real tmp84;
2107 fftw_real tmp153;
2108 fftw_real tmp88;
2109 fftw_real tmp161;
2110 ASSERT_ALIGNED_DOUBLE();
2111 {
2112 fftw_real tmp13;
2113 fftw_real tmp36;
2114 fftw_real tmp154;
2115 fftw_real tmp159;
2116 ASSERT_ALIGNED_DOUBLE();
2117 tmp13 = tmp5 + tmp12;
2118 tmp36 = tmp24 + tmp35;
2119 tmp37 = tmp13 - tmp36;
2120 tmp85 = tmp13 + tmp36;
2121 tmp154 = tmp91 - tmp90;
2122 tmp159 = tmp155 + tmp158;
2123 tmp160 = tmp154 - tmp159;
2124 tmp162 = tmp154 + tmp159;
2125 }
2126 {
2127 fftw_real tmp60;
2128 fftw_real tmp83;
2129 fftw_real tmp86;
2130 fftw_real tmp87;
2131 ASSERT_ALIGNED_DOUBLE();
2132 tmp60 = (K049008570 * tmp51) + (K497592363 * tmp59);
2133 tmp83 = (K049008570 * tmp74) - (K497592363 * tmp82);
2134 tmp84 = tmp60 + tmp83;
2135 tmp153 = tmp83 - tmp60;
2136 tmp86 = (K497592363 * tmp51) - (K049008570 * tmp59);
2137 tmp87 = (K497592363 * tmp74) + (K049008570 * tmp82);
2138 tmp88 = tmp86 + tmp87;
2139 tmp161 = tmp87 - tmp86;
2140 }
2141 X[8 * iostride] = tmp37 - tmp84;
2142 X[7 * iostride] = tmp37 + tmp84;
2143 X[15 * iostride] = tmp85 - tmp88;
2144 X[0] = tmp85 + tmp88;
2145 Y[-15 * iostride] = tmp153 - tmp160;
2146 Y[0] = tmp153 + tmp160;
2147 Y[-8 * iostride] = tmp161 - tmp162;
2148 Y[-7 * iostride] = tmp161 + tmp162;
2149 }
2150 {
2151 fftw_real tmp93;
2152 fftw_real tmp101;
2153 fftw_real tmp166;
2154 fftw_real tmp168;
2155 fftw_real tmp100;
2156 fftw_real tmp163;
2157 fftw_real tmp104;
2158 fftw_real tmp167;
2159 ASSERT_ALIGNED_DOUBLE();
2160 {
2161 fftw_real tmp89;
2162 fftw_real tmp92;
2163 fftw_real tmp164;
2164 fftw_real tmp165;
2165 ASSERT_ALIGNED_DOUBLE();
2166 tmp89 = tmp5 - tmp12;
2167 tmp92 = tmp90 + tmp91;
2168 tmp93 = tmp89 + tmp92;
2169 tmp101 = tmp89 - tmp92;
2170 tmp164 = tmp35 - tmp24;
2171 tmp165 = tmp158 - tmp155;
2172 tmp166 = tmp164 - tmp165;
2173 tmp168 = tmp164 + tmp165;
2174 }
2175 {
2176 fftw_real tmp96;
2177 fftw_real tmp99;
2178 fftw_real tmp102;
2179 fftw_real tmp103;
2180 ASSERT_ALIGNED_DOUBLE();
2181 tmp96 = (K386505226 * tmp94) + (K317196642 * tmp95);
2182 tmp99 = (K386505226 * tmp97) - (K317196642 * tmp98);
2183 tmp100 = tmp96 + tmp99;
2184 tmp163 = tmp99 - tmp96;
2185 tmp102 = (K317196642 * tmp97) + (K386505226 * tmp98);
2186 tmp103 = (K386505226 * tmp95) - (K317196642 * tmp94);
2187 tmp104 = tmp102 - tmp103;
2188 tmp167 = tmp103 + tmp102;
2189 }
2190 X[12 * iostride] = tmp93 - tmp100;
2191 X[3 * iostride] = tmp93 + tmp100;
2192 X[11 * iostride] = tmp101 - tmp104;
2193 X[4 * iostride] = tmp101 + tmp104;
2194 Y[-11 * iostride] = tmp163 - tmp166;
2195 Y[-4 * iostride] = tmp163 + tmp166;
2196 Y[-12 * iostride] = tmp167 - tmp168;
2197 Y[-3 * iostride] = tmp167 + tmp168;
2198 }
2199 {
2200 fftw_real tmp117;
2201 fftw_real tmp133;
2202 fftw_real tmp174;
2203 fftw_real tmp176;
2204 fftw_real tmp132;
2205 fftw_real tmp175;
2206 fftw_real tmp136;
2207 fftw_real tmp169;
2208 ASSERT_ALIGNED_DOUBLE();
2209 {
2210 fftw_real tmp109;
2211 fftw_real tmp116;
2212 fftw_real tmp172;
2213 fftw_real tmp173;
2214 ASSERT_ALIGNED_DOUBLE();
2215 tmp109 = tmp105 - tmp108;
2216 tmp116 = tmp112 - tmp115;
2217 tmp117 = tmp109 + tmp116;
2218 tmp133 = tmp109 - tmp116;
2219 tmp172 = tmp170 - tmp171;
2220 tmp173 = tmp138 + tmp139;
2221 tmp174 = tmp172 - tmp173;
2222 tmp176 = tmp173 + tmp172;
2223 }
2224 {
2225 fftw_real tmp124;
2226 fftw_real tmp131;
2227 fftw_real tmp134;
2228 fftw_real tmp135;
2229 ASSERT_ALIGNED_DOUBLE();
2230 tmp124 = (K440960632 * tmp120) + (K235698368 * tmp123);
2231 tmp131 = (K440960632 * tmp127) + (K235698368 * tmp130);
2232 tmp132 = tmp124 - tmp131;
2233 tmp175 = tmp124 + tmp131;
2234 tmp134 = (K440960632 * tmp130) - (K235698368 * tmp127);
2235 tmp135 = (K440960632 * tmp123) - (K235698368 * tmp120);
2236 tmp136 = tmp134 - tmp135;
2237 tmp169 = tmp135 + tmp134;
2238 }
2239 X[13 * iostride] = tmp117 - tmp132;
2240 X[2 * iostride] = tmp117 + tmp132;
2241 X[10 * iostride] = tmp133 - tmp136;
2242 X[5 * iostride] = tmp133 + tmp136;
2243 Y[-13 * iostride] = tmp169 - tmp174;
2244 Y[-2 * iostride] = tmp169 + tmp174;
2245 Y[-5 * iostride] = -(tmp175 + tmp176);
2246 Y[-10 * iostride] = tmp176 - tmp175;
2247 }
2248 {
2249 fftw_real tmp141;
2250 fftw_real tmp149;
2251 fftw_real tmp180;
2252 fftw_real tmp182;
2253 fftw_real tmp148;
2254 fftw_real tmp177;
2255 fftw_real tmp152;
2256 fftw_real tmp181;
2257 ASSERT_ALIGNED_DOUBLE();
2258 {
2259 fftw_real tmp137;
2260 fftw_real tmp140;
2261 fftw_real tmp178;
2262 fftw_real tmp179;
2263 ASSERT_ALIGNED_DOUBLE();
2264 tmp137 = tmp105 + tmp108;
2265 tmp140 = tmp138 - tmp139;
2266 tmp141 = tmp137 + tmp140;
2267 tmp149 = tmp137 - tmp140;
2268 tmp178 = tmp115 + tmp112;
2269 tmp179 = tmp170 + tmp171;
2270 tmp180 = tmp178 - tmp179;
2271 tmp182 = tmp178 + tmp179;
2272 }
2273 {
2274 fftw_real tmp144;
2275 fftw_real tmp147;
2276 fftw_real tmp150;
2277 fftw_real tmp151;
2278 ASSERT_ALIGNED_DOUBLE();
2279 tmp144 = (K478470167 * tmp142) + (K145142338 * tmp143);
2280 tmp147 = (K478470167 * tmp145) - (K145142338 * tmp146);
2281 tmp148 = tmp144 + tmp147;
2282 tmp177 = tmp147 - tmp144;
2283 tmp150 = (K145142338 * tmp145) + (K478470167 * tmp146);
2284 tmp151 = (K478470167 * tmp143) - (K145142338 * tmp142);
2285 tmp152 = tmp150 - tmp151;
2286 tmp181 = tmp151 + tmp150;
2287 }
2288 X[14 * iostride] = tmp141 - tmp148;
2289 X[iostride] = tmp141 + tmp148;
2290 X[9 * iostride] = tmp149 - tmp152;
2291 X[6 * iostride] = tmp149 + tmp152;
2292 Y[-9 * iostride] = tmp177 - tmp180;
2293 Y[-6 * iostride] = tmp177 + tmp180;
2294 Y[-14 * iostride] = tmp181 - tmp182;
2295 Y[-iostride] = tmp181 + tmp182;
2296 }
2297 }
2298 }
2299
2300 static const int twiddle_order[] =
2301 {1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31};
2302 fftw_codelet_desc fftw_hc2hc_forward_32_desc =
2303 {
2304 "fftw_hc2hc_forward_32",
2305 (void (*)()) fftw_hc2hc_forward_32,
2306 32,
2307 FFTW_FORWARD,
2308 FFTW_HC2HC,
2309 707,
2310 31,
2311 twiddle_order,
2312 };
got merge a long time ago