Mercurial > hg > js-dsp-test

comparison fft/fftw/fftw-3.3.4/rdft/scalar/r2cb/r2cbIII_15.c @ 19:26056e866c29

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add FFTW to comparison table
author Chris Cannam
date Tue, 06 Oct 2015 13:08:39 +0100
parents
children
comparison
equal deleted inserted replaced
18:8db794ca3e0b 19:26056e866c29
1 /*
2 * Copyright (c) 2003, 2007-14 Matteo Frigo
3 * Copyright (c) 2003, 2007-14 Massachusetts Institute of Technology
4 *
5 * This program is free software; you can redistribute it and/or modify
6 * it under the terms of the GNU General Public License as published by
7 * the Free Software Foundation; either version 2 of the License, or
8 * (at your option) any later version.
9 *
10 * This program is distributed in the hope that it will be useful,
11 * but WITHOUT ANY WARRANTY; without even the implied warranty of
12 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
13 * GNU General Public License for more details.
14 *
15 * You should have received a copy of the GNU General Public License
16 * along with this program; if not, write to the Free Software
17 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
18 *
19 */
20
21 /* This file was automatically generated --- DO NOT EDIT */
22 /* Generated on Tue Mar 4 13:50:34 EST 2014 */
23
24 #include "codelet-rdft.h"
25
26 #ifdef HAVE_FMA
27
28 /* Generated by: ../../../genfft/gen_r2cb.native -fma -reorder-insns -schedule-for-pipeline -compact -variables 4 -pipeline-latency 4 -sign 1 -n 15 -name r2cbIII_15 -dft-III -include r2cbIII.h */
29
30 /*
31 * This function contains 64 FP additions, 43 FP multiplications,
32 * (or, 21 additions, 0 multiplications, 43 fused multiply/add),
33 * 48 stack variables, 9 constants, and 30 memory accesses
34 */
35 #include "r2cbIII.h"
36
37 static void r2cbIII_15(R *R0, R *R1, R *Cr, R *Ci, stride rs, stride csr, stride csi, INT v, INT ivs, INT ovs)
38 {
39 DK(KP951056516, +0.951056516295153572116439333379382143405698634);
40 DK(KP1_732050807, +1.732050807568877293527446341505872366942805254);
41 DK(KP559016994, +0.559016994374947424102293417182819058860154590);
42 DK(KP250000000, +0.250000000000000000000000000000000000000000000);
43 DK(KP1_902113032, +1.902113032590307144232878666758764286811397268);
44 DK(KP1_118033988, +1.118033988749894848204586834365638117720309180);
45 DK(KP2_000000000, +2.000000000000000000000000000000000000000000000);
46 DK(KP500000000, +0.500000000000000000000000000000000000000000000);
47 DK(KP618033988, +0.618033988749894848204586834365638117720309180);
48 {
49 INT i;
50 for (i = v; i > 0; i = i - 1, R0 = R0 + ovs, R1 = R1 + ovs, Cr = Cr + ivs, Ci = Ci + ivs, MAKE_VOLATILE_STRIDE(60, rs), MAKE_VOLATILE_STRIDE(60, csr), MAKE_VOLATILE_STRIDE(60, csi)) {
51 E TX, Tv, To, TW, Tl, Tx, Ty, Tw;
52 {
53 E TA, Tk, T6, T5, Tz, Th, TI, Tp, Tu, TK, TR, Tn, Td, Tq;
54 {
55 E T1, T2, T3, Ti, Tj;
56 Ti = Ci[WS(csi, 4)];
57 Tj = Ci[WS(csi, 1)];
58 T1 = Cr[WS(csr, 7)];
59 T2 = Cr[WS(csr, 4)];
60 T3 = Cr[WS(csr, 1)];
61 TA = FNMS(KP618033988, Ti, Tj);
62 Tk = FMA(KP618033988, Tj, Ti);
63 {
64 E T7, TP, Tc, T8;
65 T6 = Cr[WS(csr, 2)];
66 {
67 E T4, Tg, Ta, Tb, Tf;
68 T4 = T2 + T3;
69 Tg = T2 - T3;
70 Ta = Cr[WS(csr, 3)];
71 Tb = Cr[WS(csr, 6)];
72 T7 = Cr[0];
73 Tf = FNMS(KP500000000, T4, T1);
74 T5 = FMA(KP2_000000000, T4, T1);
75 TP = Ta - Tb;
76 Tc = Ta + Tb;
77 Tz = FNMS(KP1_118033988, Tg, Tf);
78 Th = FMA(KP1_118033988, Tg, Tf);
79 T8 = Cr[WS(csr, 5)];
80 }
81 TI = Ci[WS(csi, 2)];
82 {
83 E Ts, Tt, TQ, T9;
84 Ts = Ci[WS(csi, 3)];
85 Tt = Ci[WS(csi, 6)];
86 TQ = T7 - T8;
87 T9 = T7 + T8;
88 Tp = Ci[0];
89 Tu = Ts - Tt;
90 TK = Ts + Tt;
91 TX = FMA(KP618033988, TP, TQ);
92 TR = FNMS(KP618033988, TQ, TP);
93 Tn = T9 - Tc;
94 Td = T9 + Tc;
95 Tq = Ci[WS(csi, 5)];
96 }
97 }
98 }
99 {
100 E TB, TF, TO, TG, TE;
101 {
102 E Tm, T11, TN, TD, TM, T12, TC;
103 TB = FNMS(KP1_902113032, TA, Tz);
104 TF = FMA(KP1_902113032, TA, Tz);
105 {
106 E Te, Tr, TJ, TL;
107 Tm = FNMS(KP250000000, Td, T6);
108 Te = T6 + Td;
109 Tr = Tp + Tq;
110 TJ = Tq - Tp;
111 R0[0] = FMA(KP2_000000000, Te, T5);
112 T11 = Te - T5;
113 TN = TJ + TK;
114 TL = TJ - TK;
115 Tv = FMA(KP618033988, Tu, Tr);
116 TD = FNMS(KP618033988, Tr, Tu);
117 TM = FNMS(KP250000000, TL, TI);
118 T12 = TL + TI;
119 }
120 TC = FNMS(KP559016994, Tn, Tm);
121 To = FMA(KP559016994, Tn, Tm);
122 R1[WS(rs, 2)] = FMA(KP1_732050807, T12, T11);
123 R0[WS(rs, 5)] = FMS(KP1_732050807, T12, T11);
124 TW = FMA(KP559016994, TN, TM);
125 TO = FNMS(KP559016994, TN, TM);
126 TG = FNMS(KP951056516, TD, TC);
127 TE = FMA(KP951056516, TD, TC);
128 }
129 Tl = FNMS(KP1_902113032, Tk, Th);
130 Tx = FMA(KP1_902113032, Tk, Th);
131 {
132 E TS, TU, TT, TH;
133 TS = FMA(KP951056516, TR, TO);
134 TU = FNMS(KP951056516, TR, TO);
135 TT = TF - TG;
136 R1[WS(rs, 1)] = -(FMA(KP2_000000000, TG, TF));
137 TH = TB - TE;
138 R0[WS(rs, 6)] = FMA(KP2_000000000, TE, TB);
139 R1[WS(rs, 6)] = -(FMA(KP1_732050807, TU, TT));
140 R0[WS(rs, 4)] = FNMS(KP1_732050807, TU, TT);
141 R1[WS(rs, 3)] = -(FMA(KP1_732050807, TS, TH));
142 R0[WS(rs, 1)] = FNMS(KP1_732050807, TS, TH);
143 }
144 }
145 }
146 Ty = FNMS(KP951056516, Tv, To);
147 Tw = FMA(KP951056516, Tv, To);
148 {
149 E T10, TY, TV, TZ;
150 T10 = FMA(KP951056516, TX, TW);
151 TY = FNMS(KP951056516, TX, TW);
152 TV = Ty - Tx;
153 R0[WS(rs, 3)] = FMA(KP2_000000000, Ty, Tx);
154 TZ = Tl - Tw;
155 R1[WS(rs, 4)] = -(FMA(KP2_000000000, Tw, Tl));
156 R1[WS(rs, 5)] = FMA(KP1_732050807, TY, TV);
157 R1[0] = FNMS(KP1_732050807, TY, TV);
158 R0[WS(rs, 2)] = FMA(KP1_732050807, T10, TZ);
159 R0[WS(rs, 7)] = FNMS(KP1_732050807, T10, TZ);
160 }
161 }
162 }
163 }
164
165 static const kr2c_desc desc = { 15, "r2cbIII_15", {21, 0, 43, 0}, &GENUS };
166
167 void X(codelet_r2cbIII_15) (planner *p) {
168 X(kr2c_register) (p, r2cbIII_15, &desc);
169 }
170
171 #else /* HAVE_FMA */
172
173 /* Generated by: ../../../genfft/gen_r2cb.native -compact -variables 4 -pipeline-latency 4 -sign 1 -n 15 -name r2cbIII_15 -dft-III -include r2cbIII.h */
174
175 /*
176 * This function contains 64 FP additions, 26 FP multiplications,
177 * (or, 49 additions, 11 multiplications, 15 fused multiply/add),
178 * 47 stack variables, 14 constants, and 30 memory accesses
179 */
180 #include "r2cbIII.h"
181
182 static void r2cbIII_15(R *R0, R *R1, R *Cr, R *Ci, stride rs, stride csr, stride csi, INT v, INT ivs, INT ovs)
183 {
184 DK(KP1_732050807, +1.732050807568877293527446341505872366942805254);
185 DK(KP433012701, +0.433012701892219323381861585376468091735701313);
186 DK(KP968245836, +0.968245836551854221294816349945599902708230426);
187 DK(KP587785252, +0.587785252292473129168705954639072768597652438);
188 DK(KP951056516, +0.951056516295153572116439333379382143405698634);
189 DK(KP250000000, +0.250000000000000000000000000000000000000000000);
190 DK(KP1_647278207, +1.647278207092663851754840078556380006059321028);
191 DK(KP1_018073920, +1.018073920910254366901961726787815297021466329);
192 DK(KP559016994, +0.559016994374947424102293417182819058860154590);
193 DK(KP500000000, +0.500000000000000000000000000000000000000000000);
194 DK(KP2_000000000, +2.000000000000000000000000000000000000000000000);
195 DK(KP1_118033988, +1.118033988749894848204586834365638117720309180);
196 DK(KP1_175570504, +1.175570504584946258337411909278145537195304875);
197 DK(KP1_902113032, +1.902113032590307144232878666758764286811397268);
198 {
199 INT i;
200 for (i = v; i > 0; i = i - 1, R0 = R0 + ovs, R1 = R1 + ovs, Cr = Cr + ivs, Ci = Ci + ivs, MAKE_VOLATILE_STRIDE(60, rs), MAKE_VOLATILE_STRIDE(60, csr), MAKE_VOLATILE_STRIDE(60, csi)) {
201 E Tv, TD, T5, Ts, TC, T6, Tf, TW, TK, Td, Tg, TP, To, TN, TA;
202 E TO, TQ, Tt, Tu, T12, Te, T11;
203 Tt = Ci[WS(csi, 4)];
204 Tu = Ci[WS(csi, 1)];
205 Tv = FMA(KP1_902113032, Tt, KP1_175570504 * Tu);
206 TD = FNMS(KP1_175570504, Tt, KP1_902113032 * Tu);
207 {
208 E T1, T4, Tq, T2, T3, Tr;
209 T1 = Cr[WS(csr, 7)];
210 T2 = Cr[WS(csr, 4)];
211 T3 = Cr[WS(csr, 1)];
212 T4 = T2 + T3;
213 Tq = KP1_118033988 * (T2 - T3);
214 T5 = FMA(KP2_000000000, T4, T1);
215 Tr = FNMS(KP500000000, T4, T1);
216 Ts = Tq + Tr;
217 TC = Tr - Tq;
218 }
219 {
220 E Tc, TJ, T9, TI;
221 T6 = Cr[WS(csr, 2)];
222 {
223 E Ta, Tb, T7, T8;
224 Ta = Cr[WS(csr, 3)];
225 Tb = Cr[WS(csr, 6)];
226 Tc = Ta + Tb;
227 TJ = Ta - Tb;
228 T7 = Cr[0];
229 T8 = Cr[WS(csr, 5)];
230 T9 = T7 + T8;
231 TI = T7 - T8;
232 }
233 Tf = KP559016994 * (T9 - Tc);
234 TW = FNMS(KP1_647278207, TJ, KP1_018073920 * TI);
235 TK = FMA(KP1_647278207, TI, KP1_018073920 * TJ);
236 Td = T9 + Tc;
237 Tg = FNMS(KP250000000, Td, T6);
238 }
239 {
240 E Tn, TM, Tk, TL;
241 TP = Ci[WS(csi, 2)];
242 {
243 E Tl, Tm, Ti, Tj;
244 Tl = Ci[WS(csi, 3)];
245 Tm = Ci[WS(csi, 6)];
246 Tn = Tl - Tm;
247 TM = Tl + Tm;
248 Ti = Ci[0];
249 Tj = Ci[WS(csi, 5)];
250 Tk = Ti + Tj;
251 TL = Ti - Tj;
252 }
253 To = FMA(KP951056516, Tk, KP587785252 * Tn);
254 TN = KP968245836 * (TL - TM);
255 TA = FNMS(KP587785252, Tk, KP951056516 * Tn);
256 TO = TL + TM;
257 TQ = FMA(KP433012701, TO, KP1_732050807 * TP);
258 }
259 T12 = KP1_732050807 * (TP - TO);
260 Te = T6 + Td;
261 T11 = Te - T5;
262 R0[0] = FMA(KP2_000000000, Te, T5);
263 R0[WS(rs, 5)] = T12 - T11;
264 R1[WS(rs, 2)] = T11 + T12;
265 {
266 E TE, TG, TB, TF, TY, T10, Tz, TX, TV, TZ;
267 TE = TC - TD;
268 TG = TC + TD;
269 Tz = Tg - Tf;
270 TB = Tz + TA;
271 TF = TA - Tz;
272 TX = TN + TQ;
273 TY = TW - TX;
274 T10 = TW + TX;
275 R0[WS(rs, 6)] = FMA(KP2_000000000, TB, TE);
276 R1[WS(rs, 1)] = FMS(KP2_000000000, TF, TG);
277 TV = TE - TB;
278 R0[WS(rs, 1)] = TV + TY;
279 R1[WS(rs, 3)] = TY - TV;
280 TZ = TF + TG;
281 R0[WS(rs, 4)] = TZ - T10;
282 R1[WS(rs, 6)] = -(TZ + T10);
283 }
284 {
285 E Tw, Ty, Tp, Tx, TS, TU, Th, TR, TH, TT;
286 Tw = Ts - Tv;
287 Ty = Ts + Tv;
288 Th = Tf + Tg;
289 Tp = Th + To;
290 Tx = Th - To;
291 TR = TN - TQ;
292 TS = TK + TR;
293 TU = TR - TK;
294 R1[WS(rs, 4)] = -(FMA(KP2_000000000, Tp, Tw));
295 R0[WS(rs, 3)] = FMA(KP2_000000000, Tx, Ty);
296 TH = Tx - Ty;
297 R1[WS(rs, 5)] = TH - TS;
298 R1[0] = TH + TS;
299 TT = Tw - Tp;
300 R0[WS(rs, 2)] = TT - TU;
301 R0[WS(rs, 7)] = TT + TU;
302 }
303 }
304 }
305 }
306
307 static const kr2c_desc desc = { 15, "r2cbIII_15", {49, 11, 15, 0}, &GENUS };
308
309 void X(codelet_r2cbIII_15) (planner *p) {
310 X(kr2c_register) (p, r2cbIII_15, &desc);
311 }
312
313 #endif /* HAVE_FMA */