| Bsoft 1.9 Bernard's software package

- Main Page
- Namespaces
- Classes
- Files

- File List
- File Members

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