| Bsoft 1.9 Bernard's software package

- Main Page
- Namespaces
- Classes
- Files

- File List
- File Members

fn_14.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:05:41 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 -notwiddle 14 */
27
28 /*
29 * This function contains 148 FP additions, 72 FP multiplications,
30 * (or, 148 additions, 72 multiplications, 0 fused multiply/add),
31 * 36 stack variables, and 56 memory accesses
32 */
33 static const fftw_real K900968867 =
34 FFTW_KONST(+0.900968867902419126236102319507445051165919162);
35 static const fftw_real K222520933 =
36 FFTW_KONST(+0.222520933956314404288902564496794759466355569);
37 static const fftw_real K623489801 =
38 FFTW_KONST(+0.623489801858733530525004884004239810632274731);
39 static const fftw_real K781831482 =
40 FFTW_KONST(+0.781831482468029808708444526674057750232334519);
41 static const fftw_real K974927912 =
42 FFTW_KONST(+0.974927912181823607018131682993931217232785801);
43 static const fftw_real K433883739 =
44 FFTW_KONST(+0.433883739117558120475768332848358754609990728);
45
46 /*
47 * Generator Id's :
48 * $Id: exprdag.ml,v 1.43 2003/03/16 23:43:46 stevenj Exp $
49 * $Id: fft.ml,v 1.44 2003/03/16 23:43:46 stevenj Exp $
50 * $Id: to_c.ml,v 1.26 2003/03/16 23:43:46 stevenj Exp $
51 */
52
53 void fftw_no_twiddle_14(const fftw_complex *input, fftw_complex *output,
54 int istride, int ostride)
55 {
56 fftw_real tmp3;
57 fftw_real tmp25;
58 fftw_real tmp84;
59 fftw_real tmp93;
60 fftw_real tmp10;
61 fftw_real tmp77;
62 fftw_real tmp28;
63 fftw_real tmp97;
64 fftw_real tmp42;
65 fftw_real tmp86;
66 fftw_real tmp65;
67 fftw_real tmp92;
68 fftw_real tmp17;
69 fftw_real tmp79;
70 fftw_real tmp31;
71 fftw_real tmp99;
72 fftw_real tmp56;
73 fftw_real tmp81;
74 fftw_real tmp68;
75 fftw_real tmp94;
76 fftw_real tmp24;
77 fftw_real tmp78;
78 fftw_real tmp34;
79 fftw_real tmp98;
80 fftw_real tmp49;
81 fftw_real tmp85;
82 fftw_real tmp71;
83 fftw_real tmp95;
84 ASSERT_ALIGNED_DOUBLE;
85 {
86 fftw_real tmp1;
87 fftw_real tmp2;
88 fftw_real tmp82;
89 fftw_real tmp83;
90 ASSERT_ALIGNED_DOUBLE;
91 tmp1 = c_re(input[0]);
92 tmp2 = c_re(input[7 * istride]);
93 tmp3 = tmp1 - tmp2;
94 tmp25 = tmp1 + tmp2;
95 tmp82 = c_im(input[0]);
96 tmp83 = c_im(input[7 * istride]);
97 tmp84 = tmp82 - tmp83;
98 tmp93 = tmp82 + tmp83;
99 }
100 {
101 fftw_real tmp6;
102 fftw_real tmp26;
103 fftw_real tmp9;
104 fftw_real tmp27;
105 ASSERT_ALIGNED_DOUBLE;
106 {
107 fftw_real tmp4;
108 fftw_real tmp5;
109 fftw_real tmp7;
110 fftw_real tmp8;
111 ASSERT_ALIGNED_DOUBLE;
112 tmp4 = c_re(input[2 * istride]);
113 tmp5 = c_re(input[9 * istride]);
114 tmp6 = tmp4 - tmp5;
115 tmp26 = tmp4 + tmp5;
116 tmp7 = c_re(input[12 * istride]);
117 tmp8 = c_re(input[5 * istride]);
118 tmp9 = tmp7 - tmp8;
119 tmp27 = tmp7 + tmp8;
120 }
121 tmp10 = tmp6 + tmp9;
122 tmp77 = tmp9 - tmp6;
123 tmp28 = tmp26 + tmp27;
124 tmp97 = tmp27 - tmp26;
125 }
126 {
127 fftw_real tmp38;
128 fftw_real tmp63;
129 fftw_real tmp41;
130 fftw_real tmp64;
131 ASSERT_ALIGNED_DOUBLE;
132 {
133 fftw_real tmp36;
134 fftw_real tmp37;
135 fftw_real tmp39;
136 fftw_real tmp40;
137 ASSERT_ALIGNED_DOUBLE;
138 tmp36 = c_im(input[2 * istride]);
139 tmp37 = c_im(input[9 * istride]);
140 tmp38 = tmp36 - tmp37;
141 tmp63 = tmp36 + tmp37;
142 tmp39 = c_im(input[12 * istride]);
143 tmp40 = c_im(input[5 * istride]);
144 tmp41 = tmp39 - tmp40;
145 tmp64 = tmp39 + tmp40;
146 }
147 tmp42 = tmp38 - tmp41;
148 tmp86 = tmp38 + tmp41;
149 tmp65 = tmp63 - tmp64;
150 tmp92 = tmp63 + tmp64;
151 }
152 {
153 fftw_real tmp13;
154 fftw_real tmp29;
155 fftw_real tmp16;
156 fftw_real tmp30;
157 ASSERT_ALIGNED_DOUBLE;
158 {
159 fftw_real tmp11;
160 fftw_real tmp12;
161 fftw_real tmp14;
162 fftw_real tmp15;
163 ASSERT_ALIGNED_DOUBLE;
164 tmp11 = c_re(input[4 * istride]);
165 tmp12 = c_re(input[11 * istride]);
166 tmp13 = tmp11 - tmp12;
167 tmp29 = tmp11 + tmp12;
168 tmp14 = c_re(input[10 * istride]);
169 tmp15 = c_re(input[3 * istride]);
170 tmp16 = tmp14 - tmp15;
171 tmp30 = tmp14 + tmp15;
172 }
173 tmp17 = tmp13 + tmp16;
174 tmp79 = tmp16 - tmp13;
175 tmp31 = tmp29 + tmp30;
176 tmp99 = tmp29 - tmp30;
177 }
178 {
179 fftw_real tmp52;
180 fftw_real tmp67;
181 fftw_real tmp55;
182 fftw_real tmp66;
183 ASSERT_ALIGNED_DOUBLE;
184 {
185 fftw_real tmp50;
186 fftw_real tmp51;
187 fftw_real tmp53;
188 fftw_real tmp54;
189 ASSERT_ALIGNED_DOUBLE;
190 tmp50 = c_im(input[4 * istride]);
191 tmp51 = c_im(input[11 * istride]);
192 tmp52 = tmp50 - tmp51;
193 tmp67 = tmp50 + tmp51;
194 tmp53 = c_im(input[10 * istride]);
195 tmp54 = c_im(input[3 * istride]);
196 tmp55 = tmp53 - tmp54;
197 tmp66 = tmp53 + tmp54;
198 }
199 tmp56 = tmp52 - tmp55;
200 tmp81 = tmp52 + tmp55;
201 tmp68 = tmp66 - tmp67;
202 tmp94 = tmp67 + tmp66;
203 }
204 {
205 fftw_real tmp20;
206 fftw_real tmp32;
207 fftw_real tmp23;
208 fftw_real tmp33;
209 ASSERT_ALIGNED_DOUBLE;
210 {
211 fftw_real tmp18;
212 fftw_real tmp19;
213 fftw_real tmp21;
214 fftw_real tmp22;
215 ASSERT_ALIGNED_DOUBLE;
216 tmp18 = c_re(input[6 * istride]);
217 tmp19 = c_re(input[13 * istride]);
218 tmp20 = tmp18 - tmp19;
219 tmp32 = tmp18 + tmp19;
220 tmp21 = c_re(input[8 * istride]);
221 tmp22 = c_re(input[istride]);
222 tmp23 = tmp21 - tmp22;
223 tmp33 = tmp21 + tmp22;
224 }
225 tmp24 = tmp20 + tmp23;
226 tmp78 = tmp23 - tmp20;
227 tmp34 = tmp32 + tmp33;
228 tmp98 = tmp32 - tmp33;
229 }
230 {
231 fftw_real tmp45;
232 fftw_real tmp70;
233 fftw_real tmp48;
234 fftw_real tmp69;
235 ASSERT_ALIGNED_DOUBLE;
236 {
237 fftw_real tmp43;
238 fftw_real tmp44;
239 fftw_real tmp46;
240 fftw_real tmp47;
241 ASSERT_ALIGNED_DOUBLE;
242 tmp43 = c_im(input[6 * istride]);
243 tmp44 = c_im(input[13 * istride]);
244 tmp45 = tmp43 - tmp44;
245 tmp70 = tmp43 + tmp44;
246 tmp46 = c_im(input[8 * istride]);
247 tmp47 = c_im(input[istride]);
248 tmp48 = tmp46 - tmp47;
249 tmp69 = tmp46 + tmp47;
250 }
251 tmp49 = tmp45 - tmp48;
252 tmp85 = tmp45 + tmp48;
253 tmp71 = tmp69 - tmp70;
254 tmp95 = tmp70 + tmp69;
255 }
256 {
257 fftw_real tmp57;
258 fftw_real tmp35;
259 fftw_real tmp72;
260 fftw_real tmp62;
261 ASSERT_ALIGNED_DOUBLE;
262 c_re(output[7 * ostride]) = tmp3 + tmp10 + tmp17 + tmp24;
263 tmp57 =
264 (K433883739 * tmp42) + (K974927912 * tmp49) -
265 (K781831482 * tmp56);
266 tmp35 =
267 tmp3 + (K623489801 * tmp17) - (K222520933 * tmp24) -
268 (K900968867 * tmp10);
269 c_re(output[11 * ostride]) = tmp35 - tmp57;
270 c_re(output[3 * ostride]) = tmp35 + tmp57;
271 {
272 fftw_real tmp59;
273 fftw_real tmp58;
274 fftw_real tmp61;
275 fftw_real tmp60;
276 ASSERT_ALIGNED_DOUBLE;
277 tmp59 =
278 (K974927912 * tmp42) - (K433883739 * tmp56) -
279 (K781831482 * tmp49);
280 tmp58 =
281 tmp3 + (K623489801 * tmp24) - (K900968867 * tmp17) -
282 (K222520933 * tmp10);
283 c_re(output[5 * ostride]) = tmp58 - tmp59;
284 c_re(output[9 * ostride]) = tmp58 + tmp59;
285 tmp61 =
286 (K781831482 * tmp42) + (K433883739 * tmp49) +
287 (K974927912 * tmp56);
288 tmp60 =
289 tmp3 + (K623489801 * tmp10) - (K900968867 * tmp24) -
290 (K222520933 * tmp17);
291 c_re(output[13 * ostride]) = tmp60 - tmp61;
292 c_re(output[ostride]) = tmp60 + tmp61;
293 }
294 c_re(output[0]) = tmp25 + tmp28 + tmp31 + tmp34;
295 tmp72 =
296 (K781831482 * tmp65) - (K974927912 * tmp68) -
297 (K433883739 * tmp71);
298 tmp62 =
299 tmp25 + (K623489801 * tmp28) - (K900968867 * tmp34) -
300 (K222520933 * tmp31);
301 c_re(output[6 * ostride]) = tmp62 - tmp72;
302 c_re(output[8 * ostride]) = tmp62 + tmp72;
303 {
304 fftw_real tmp74;
305 fftw_real tmp73;
306 fftw_real tmp76;
307 fftw_real tmp75;
308 ASSERT_ALIGNED_DOUBLE;
309 tmp74 =
310 (K433883739 * tmp65) + (K781831482 * tmp68) -
311 (K974927912 * tmp71);
312 tmp73 =
313 tmp25 + (K623489801 * tmp31) - (K222520933 * tmp34) -
314 (K900968867 * tmp28);
315 c_re(output[4 * ostride]) = tmp73 - tmp74;
316 c_re(output[10 * ostride]) = tmp73 + tmp74;
317 tmp76 =
318 (K974927912 * tmp65) + (K781831482 * tmp71) +
319 (K433883739 * tmp68);
320 tmp75 =
321 tmp25 + (K623489801 * tmp34) - (K900968867 * tmp31) -
322 (K222520933 * tmp28);
323 c_re(output[12 * ostride]) = tmp75 - tmp76;
324 c_re(output[2 * ostride]) = tmp75 + tmp76;
325 }
326 }
327 {
328 fftw_real tmp91;
329 fftw_real tmp90;
330 fftw_real tmp103;
331 fftw_real tmp104;
332 ASSERT_ALIGNED_DOUBLE;
333 c_im(output[7 * ostride]) = tmp86 + tmp85 + tmp81 + tmp84;
334 tmp91 =
335 (K974927912 * tmp77) - (K781831482 * tmp78) -
336 (K433883739 * tmp79);
337 tmp90 =
338 (K623489801 * tmp85) + tmp84 - (K900968867 * tmp81) -
339 (K222520933 * tmp86);
340 c_im(output[5 * ostride]) = tmp90 - tmp91;
341 c_im(output[9 * ostride]) = tmp91 + tmp90;
342 {
343 fftw_real tmp88;
344 fftw_real tmp89;
345 fftw_real tmp80;
346 fftw_real tmp87;
347 ASSERT_ALIGNED_DOUBLE;
348 tmp88 =
349 (K781831482 * tmp77) + (K974927912 * tmp79) +
350 (K433883739 * tmp78);
351 tmp89 =
352 (K623489801 * tmp86) + tmp84 - (K222520933 * tmp81) -
353 (K900968867 * tmp85);
354 c_im(output[ostride]) = tmp88 + tmp89;
355 c_im(output[13 * ostride]) = tmp89 - tmp88;
356 tmp80 =
357 (K433883739 * tmp77) + (K974927912 * tmp78) -
358 (K781831482 * tmp79);
359 tmp87 =
360 (K623489801 * tmp81) + tmp84 - (K222520933 * tmp85) -
361 (K900968867 * tmp86);
362 c_im(output[3 * ostride]) = tmp80 + tmp87;
363 c_im(output[11 * ostride]) = tmp87 - tmp80;
364 }
365 c_im(output[0]) = tmp92 + tmp95 + tmp94 + tmp93;
366 tmp103 =
367 (K974927912 * tmp97) + (K433883739 * tmp99) +
368 (K781831482 * tmp98);
369 tmp104 =
370 (K623489801 * tmp95) + tmp93 - (K900968867 * tmp94) -
371 (K222520933 * tmp92);
372 c_im(output[2 * ostride]) = tmp103 + tmp104;
373 c_im(output[12 * ostride]) = tmp104 - tmp103;
374 {
375 fftw_real tmp100;
376 fftw_real tmp96;
377 fftw_real tmp102;
378 fftw_real tmp101;
379 ASSERT_ALIGNED_DOUBLE;
380 tmp100 =
381 (K781831482 * tmp97) - (K433883739 * tmp98) -
382 (K974927912 * tmp99);
383 tmp96 =
384 (K623489801 * tmp92) + tmp93 - (K222520933 * tmp94) -
385 (K900968867 * tmp95);
386 c_im(output[6 * ostride]) = tmp96 - tmp100;
387 c_im(output[8 * ostride]) = tmp100 + tmp96;
388 tmp102 =
389 (K433883739 * tmp97) + (K781831482 * tmp99) -
390 (K974927912 * tmp98);
391 tmp101 =
392 (K623489801 * tmp94) + tmp93 - (K222520933 * tmp95) -
393 (K900968867 * tmp92);
394 c_im(output[4 * ostride]) = tmp101 - tmp102;
395 c_im(output[10 * ostride]) = tmp102 + tmp101;
396 }
397 }
398 }
399
400 fftw_codelet_desc fftw_no_twiddle_14_desc = {
401 "fftw_no_twiddle_14",
402 (void (*)()) fftw_no_twiddle_14,
403 14,
404 FFTW_FORWARD,
405 FFTW_NOTW,
406 309,
407 0,
408 (const int *) 0,
409 };
fftw_codelet_desc
Definition: fftw.h:126
fftw-int.h
fftw_complex
Definition: fftw.h:57
FFTW_KONST
#define FFTW_KONST(x)
Definition: fftw-int.h:178
fftw.h
fftw_no_twiddle_14_desc
fftw_codelet_desc fftw_no_twiddle_14_desc
Definition: fn_14.c:400
c_re
#define c_re(c)
Definition: fftw.h:60
fftw_no_twiddle_14
void fftw_no_twiddle_14(const fftw_complex *input, fftw_complex *output, int istride, int ostride)
Definition: fn_14.c:53
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
FFTW_FORWARD
Definition: fftw.h:64