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