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:20 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 7 */
27
28 /*
29 * This function contains 60 FP additions, 36 FP multiplications,
30 * (or, 60 additions, 36 multiplications, 0 fused multiply/add),
31 * 22 stack variables, and 28 memory accesses
32 */
33 static const fftw_real K222520933 =
34 FFTW_KONST(+0.222520933956314404288902564496794759466355569);
35 static const fftw_real K900968867 =
36 FFTW_KONST(+0.900968867902419126236102319507445051165919162);
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 K433883739 =
42 FFTW_KONST(+0.433883739117558120475768332848358754609990728);
43 static const fftw_real K974927912 =
44 FFTW_KONST(+0.974927912181823607018131682993931217232785801);
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
fftwi_no_twiddle_7(const fftw_complex * input,fftw_complex * output,int istride,int ostride)53 void fftwi_no_twiddle_7(const fftw_complex *input, fftw_complex *output,
54 int istride, int ostride)
55 {
56 fftw_real tmp1;
57 fftw_real tmp15;
58 fftw_real tmp4;
59 fftw_real tmp11;
60 fftw_real tmp21;
61 fftw_real tmp31;
62 fftw_real tmp7;
63 fftw_real tmp13;
64 fftw_real tmp24;
65 fftw_real tmp33;
66 fftw_real tmp10;
67 fftw_real tmp12;
68 fftw_real tmp18;
69 fftw_real tmp32;
70 ASSERT_ALIGNED_DOUBLE;
71 tmp1 = c_re(input[0]);
72 tmp15 = c_im(input[0]);
73 {
74 fftw_real tmp2;
75 fftw_real tmp3;
76 fftw_real tmp19;
77 fftw_real tmp20;
78 ASSERT_ALIGNED_DOUBLE;
79 tmp2 = c_re(input[istride]);
80 tmp3 = c_re(input[6 * istride]);
81 tmp4 = tmp2 + tmp3;
82 tmp11 = tmp2 - tmp3;
83 tmp19 = c_im(input[istride]);
84 tmp20 = c_im(input[6 * istride]);
85 tmp21 = tmp19 + tmp20;
86 tmp31 = tmp20 - tmp19;
87 }
88 {
89 fftw_real tmp5;
90 fftw_real tmp6;
91 fftw_real tmp22;
92 fftw_real tmp23;
93 ASSERT_ALIGNED_DOUBLE;
94 tmp5 = c_re(input[2 * istride]);
95 tmp6 = c_re(input[5 * istride]);
96 tmp7 = tmp5 + tmp6;
97 tmp13 = tmp5 - tmp6;
98 tmp22 = c_im(input[2 * istride]);
99 tmp23 = c_im(input[5 * istride]);
100 tmp24 = tmp22 + tmp23;
101 tmp33 = tmp23 - tmp22;
102 }
103 {
104 fftw_real tmp8;
105 fftw_real tmp9;
106 fftw_real tmp16;
107 fftw_real tmp17;
108 ASSERT_ALIGNED_DOUBLE;
109 tmp8 = c_re(input[3 * istride]);
110 tmp9 = c_re(input[4 * istride]);
111 tmp10 = tmp8 + tmp9;
112 tmp12 = tmp8 - tmp9;
113 tmp16 = c_im(input[3 * istride]);
114 tmp17 = c_im(input[4 * istride]);
115 tmp18 = tmp16 + tmp17;
116 tmp32 = tmp17 - tmp16;
117 }
118 {
119 fftw_real tmp36;
120 fftw_real tmp35;
121 fftw_real tmp26;
122 fftw_real tmp27;
123 ASSERT_ALIGNED_DOUBLE;
124 c_re(output[0]) = tmp1 + tmp4 + tmp7 + tmp10;
125 tmp36 =
126 (K974927912 * tmp33) + (K433883739 * tmp32) +
127 (K781831482 * tmp31);
128 tmp35 =
129 tmp1 + (K623489801 * tmp4) - (K900968867 * tmp10) -
130 (K222520933 * tmp7);
131 c_re(output[6 * ostride]) = tmp35 - tmp36;
132 c_re(output[ostride]) = tmp35 + tmp36;
133 {
134 fftw_real tmp38;
135 fftw_real tmp37;
136 fftw_real tmp34;
137 fftw_real tmp30;
138 ASSERT_ALIGNED_DOUBLE;
139 tmp38 =
140 (K974927912 * tmp32) + (K433883739 * tmp31) -
141 (K781831482 * tmp33);
142 tmp37 =
143 tmp1 + (K623489801 * tmp7) - (K222520933 * tmp10) -
144 (K900968867 * tmp4);
145 c_re(output[4 * ostride]) = tmp37 - tmp38;
146 c_re(output[3 * ostride]) = tmp37 + tmp38;
147 tmp34 =
148 (K974927912 * tmp31) - (K781831482 * tmp32) -
149 (K433883739 * tmp33);
150 tmp30 =
151 tmp1 + (K623489801 * tmp10) - (K900968867 * tmp7) -
152 (K222520933 * tmp4);
153 c_re(output[5 * ostride]) = tmp30 - tmp34;
154 c_re(output[2 * ostride]) = tmp30 + tmp34;
155 }
156 c_im(output[0]) = tmp15 + tmp24 + tmp18 + tmp21;
157 tmp26 =
158 (K433883739 * tmp11) + (K974927912 * tmp12) -
159 (K781831482 * tmp13);
160 tmp27 =
161 tmp15 + (K623489801 * tmp24) - (K900968867 * tmp21) -
162 (K222520933 * tmp18);
163 c_im(output[3 * ostride]) = tmp26 + tmp27;
164 c_im(output[4 * ostride]) = tmp27 - tmp26;
165 {
166 fftw_real tmp14;
167 fftw_real tmp25;
168 fftw_real tmp28;
169 fftw_real tmp29;
170 ASSERT_ALIGNED_DOUBLE;
171 tmp14 =
172 (K974927912 * tmp11) - (K781831482 * tmp12) -
173 (K433883739 * tmp13);
174 tmp25 =
175 tmp15 + (K623489801 * tmp18) - (K222520933 * tmp21) -
176 (K900968867 * tmp24);
177 c_im(output[2 * ostride]) = tmp14 + tmp25;
178 c_im(output[5 * ostride]) = tmp25 - tmp14;
179 tmp28 =
180 (K781831482 * tmp11) + (K974927912 * tmp13) +
181 (K433883739 * tmp12);
182 tmp29 =
183 tmp15 + (K623489801 * tmp21) - (K900968867 * tmp18) -
184 (K222520933 * tmp24);
185 c_im(output[ostride]) = tmp28 + tmp29;
186 c_im(output[6 * ostride]) = tmp29 - tmp28;
187 }
188 }
189 }
190
191 fftw_codelet_desc fftwi_no_twiddle_7_desc = {
192 "fftwi_no_twiddle_7",
193 (void (*)()) fftwi_no_twiddle_7,
194 7,
195 FFTW_BACKWARD,
196 FFTW_NOTW,
197 166,
198 0,
199 (const int *) 0,
200 };
201