/src/fftw3/rdft/scalar/r2cb/hb2_25.c
Line | Count | Source (jump to first uncovered line) |
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 Fri Jul 11 06:54:10 UTC 2025 */ |
23 | | |
24 | | #include "rdft/codelet-rdft.h" |
25 | | |
26 | | #if defined(ARCH_PREFERS_FMA) || defined(ISA_EXTENSION_PREFERS_FMA) |
27 | | |
28 | | /* Generated by: ../../../genfft/gen_hc2hc.native -fma -compact -variables 4 -pipeline-latency 4 -sign 1 -twiddle-log3 -precompute-twiddles -n 25 -dif -name hb2_25 -include rdft/scalar/hb.h */ |
29 | | |
30 | | /* |
31 | | * This function contains 440 FP additions, 434 FP multiplications, |
32 | | * (or, 84 additions, 78 multiplications, 356 fused multiply/add), |
33 | | * 206 stack variables, 47 constants, and 100 memory accesses |
34 | | */ |
35 | | #include "rdft/scalar/hb.h" |
36 | | |
37 | | static void hb2_25(R *cr, R *ci, const R *W, stride rs, INT mb, INT me, INT ms) |
38 | | { |
39 | | DK(KP921177326, +0.921177326965143320250447435415066029359282231); |
40 | | DK(KP833417178, +0.833417178328688677408962550243238843138996060); |
41 | | DK(KP541454447, +0.541454447536312777046285590082819509052033189); |
42 | | DK(KP968583161, +0.968583161128631119490168375464735813836012403); |
43 | | DK(KP242145790, +0.242145790282157779872542093866183953459003101); |
44 | | DK(KP904730450, +0.904730450839922351881287709692877908104763647); |
45 | | DK(KP683113946, +0.683113946453479238701949862233725244439656928); |
46 | | DK(KP559154169, +0.559154169276087864842202529084232643714075927); |
47 | | DK(KP831864738, +0.831864738706457140726048799369896829771167132); |
48 | | DK(KP871714437, +0.871714437527667770979999223229522602943903653); |
49 | | DK(KP803003575, +0.803003575438660414833440593570376004635464850); |
50 | | DK(KP554608978, +0.554608978404018097464974850792216217022558774); |
51 | | DK(KP992114701, +0.992114701314477831049793042785778521453036709); |
52 | | DK(KP248028675, +0.248028675328619457762448260696444630363259177); |
53 | | DK(KP851038619, +0.851038619207379630836264138867114231259902550); |
54 | | DK(KP525970792, +0.525970792408939708442463226536226366643874659); |
55 | | DK(KP726211448, +0.726211448929902658173535992263577167607493062); |
56 | | DK(KP912018591, +0.912018591466481957908415381764119056233607330); |
57 | | DK(KP912575812, +0.912575812670962425556968549836277086778922727); |
58 | | DK(KP943557151, +0.943557151597354104399655195398983005179443399); |
59 | | DK(KP994076283, +0.994076283785401014123185814696322018529298887); |
60 | | DK(KP614372930, +0.614372930789563808870829930444362096004872855); |
61 | | DK(KP621716863, +0.621716863012209892444754556304102309693593202); |
62 | | DK(KP772036680, +0.772036680810363904029489473607579825330539880); |
63 | | DK(KP734762448, +0.734762448793050413546343770063151342619912334); |
64 | | DK(KP949179823, +0.949179823508441261575555465843363271711583843); |
65 | | DK(KP860541664, +0.860541664367944677098261680920518816412804187); |
66 | | DK(KP557913902, +0.557913902031834264187699648465567037992437152); |
67 | | DK(KP249506682, +0.249506682107067890488084201715862638334226305); |
68 | | DK(KP998026728, +0.998026728428271561952336806863450553336905220); |
69 | | DK(KP906616052, +0.906616052148196230441134447086066874408359177); |
70 | | DK(KP560319534, +0.560319534973832390111614715371676131169633784); |
71 | | DK(KP681693190, +0.681693190061530575150324149145440022633095390); |
72 | | DK(KP845997307, +0.845997307939530944175097360758058292389769300); |
73 | | DK(KP968479752, +0.968479752739016373193524836781420152702090879); |
74 | | DK(KP062914667, +0.062914667253649757225485955897349402364686947); |
75 | | DK(KP827271945, +0.827271945972475634034355757144307982555673741); |
76 | | DK(KP126329378, +0.126329378446108174786050455341811215027378105); |
77 | | DK(KP470564281, +0.470564281212251493087595091036643380879947982); |
78 | | DK(KP634619297, +0.634619297544148100711287640319130485732531031); |
79 | | DK(KP256756360, +0.256756360367726783319498520922669048172391148); |
80 | | DK(KP939062505, +0.939062505817492352556001843133229685779824606); |
81 | | DK(KP549754652, +0.549754652192770074288023275540779861653779767); |
82 | | DK(KP951056516, +0.951056516295153572116439333379382143405698634); |
83 | | DK(KP559016994, +0.559016994374947424102293417182819058860154590); |
84 | | DK(KP250000000, +0.250000000000000000000000000000000000000000000); |
85 | | DK(KP618033988, +0.618033988749894848204586834365638117720309180); |
86 | | { |
87 | | INT m; |
88 | | for (m = mb, W = W + ((mb - 1) * 8); m < me; m = m + 1, cr = cr + ms, ci = ci - ms, W = W + 8, MAKE_VOLATILE_STRIDE(50, rs)) { |
89 | | E TN, TT, TO, TR, T23, T25, TQ, TS, T4l, TW, T4n, TX, T2e, T2y, T4z; |
90 | | E T2q, T76, T4o, T8d, T2u, T4e, T4i, T8a, T86, T71, T6Y, T6U, T26, T2a, T3U; |
91 | | E T8o, T8s, T4B, T4C, T4G, T2k, T5w, T5C, T6E, T5T, T4u, T7g, T7c, T1I, TY; |
92 | | E T5I, T8i, T5M; |
93 | | { |
94 | | E T2x, T2p, T85, T4d, T2t, T89, T4h, TU, T4m, T2j, T3T, TP, TV, T2d, T5v; |
95 | | E T5B; |
96 | | TN = W[0]; |
97 | | TT = W[4]; |
98 | | TO = W[2]; |
99 | | TR = W[3]; |
100 | | TP = TN * TO; |
101 | | T2x = TO * TT; |
102 | | TV = TN * TR; |
103 | | T2d = TN * TT; |
104 | | T23 = W[6]; |
105 | | T2p = TT * T23; |
106 | | T85 = TN * T23; |
107 | | T4d = TO * T23; |
108 | | T25 = W[7]; |
109 | | T2t = TT * T25; |
110 | | T89 = TN * T25; |
111 | | T4h = TO * T25; |
112 | | TQ = W[1]; |
113 | | TS = FNMS(TQ, TR, TP); |
114 | | T4l = FMA(TQ, TR, TP); |
115 | | TU = TS * TT; |
116 | | T4m = T4l * TT; |
117 | | TW = FMA(TQ, TO, TV); |
118 | | T4n = FNMS(TQ, TO, TV); |
119 | | TX = W[5]; |
120 | | T2j = TN * TX; |
121 | | T3T = TO * TX; |
122 | | T2e = FNMS(TQ, TX, T2d); |
123 | | T2y = FMA(TR, TX, T2x); |
124 | | T4z = FMA(TQ, TX, T2d); |
125 | | T2q = FMA(TX, T25, T2p); |
126 | | T76 = FMA(TR, TT, T3T); |
127 | | T4o = FNMS(T4n, TX, T4m); |
128 | | T8d = FNMS(TW, TX, TU); |
129 | | T2u = FNMS(TX, T23, T2t); |
130 | | T4e = FMA(TR, T25, T4d); |
131 | | T4i = FNMS(TR, T23, T4h); |
132 | | T8a = FNMS(TQ, T23, T89); |
133 | | T86 = FMA(TQ, T25, T85); |
134 | | { |
135 | | E T6X, T6T, T24, T29; |
136 | | T71 = FNMS(TR, TX, T2x); |
137 | | T6X = T4l * T25; |
138 | | T6Y = FNMS(T4n, T23, T6X); |
139 | | T6T = T4l * T23; |
140 | | T6U = FMA(T4n, T25, T6T); |
141 | | T24 = TS * T23; |
142 | | T26 = FMA(TW, T25, T24); |
143 | | T29 = TS * T25; |
144 | | T2a = FNMS(TW, T23, T29); |
145 | | } |
146 | | { |
147 | | E T8n, T8r, T4A, T4F; |
148 | | T8n = T2y * T23; |
149 | | T8r = T2y * T25; |
150 | | T3U = FNMS(TR, TT, T3T); |
151 | | T8o = FMA(T3U, T25, T8n); |
152 | | T8s = FNMS(T3U, T23, T8r); |
153 | | T4A = T4z * T23; |
154 | | T4F = T4z * T25; |
155 | | T4B = FNMS(TQ, TT, T2j); |
156 | | T4C = FMA(T4B, T25, T4A); |
157 | | T4G = FNMS(T4B, T23, T4F); |
158 | | } |
159 | | T5v = T2e * T23; |
160 | | T5B = T2e * T25; |
161 | | T2k = FMA(TQ, TT, T2j); |
162 | | T5w = FMA(T2k, T25, T5v); |
163 | | T5C = FNMS(T2k, T23, T5B); |
164 | | { |
165 | | E T4t, T7b, T7f, T1H, T5H, T5L; |
166 | | T4t = T4l * TX; |
167 | | T6E = FNMS(T4n, TT, T4t); |
168 | | T5T = FMA(T4n, TX, T4m); |
169 | | T7b = T5T * T23; |
170 | | T7f = T5T * T25; |
171 | | T4u = FMA(T4n, TT, T4t); |
172 | | T7g = FNMS(T6E, T23, T7f); |
173 | | T7c = FMA(T6E, T25, T7b); |
174 | | T1H = TS * TX; |
175 | | T1I = FNMS(TW, TT, T1H); |
176 | | TY = FMA(TW, TX, TU); |
177 | | T5H = TY * T23; |
178 | | T5L = TY * T25; |
179 | | T5I = FMA(T1I, T25, T5H); |
180 | | T8i = FMA(TW, TT, T1H); |
181 | | T5M = FNMS(T1I, T23, T5L); |
182 | | } |
183 | | } |
184 | | { |
185 | | E T9, T40, T1R, T6G, T6F, T3X, T6H, T2F, T7n, T4N, T5W, T1k, T1S, T1D, T1T; |
186 | | E Ti, Tr, Ts, TB, TK, TL, TM, T6p, T7K, T6w, T7A, T2U, T56, T3K, T4X; |
187 | | E T6i, T7J, T6v, T7x, T39, T57, T3L, T50, T3E, T59, T3O, T4Q, T63, T7H, T6y; |
188 | | E T7t, T3p, T5a, T3N, T4T, T6a, T7G, T6z, T7q; |
189 | | { |
190 | | E T1, T1J, T8, T3Z, T2A, T3Y, T1Q, T3W, T2C, T2D, T3V; |
191 | | T1 = cr[0]; |
192 | | T1J = ci[WS(rs, 24)]; |
193 | | { |
194 | | E T2, T3, T4, T5, T6, T7; |
195 | | T2 = cr[WS(rs, 5)]; |
196 | | T3 = ci[WS(rs, 4)]; |
197 | | T4 = T2 + T3; |
198 | | T5 = cr[WS(rs, 10)]; |
199 | | T6 = ci[WS(rs, 9)]; |
200 | | T7 = T5 + T6; |
201 | | T8 = T4 + T7; |
202 | | T3Z = T5 - T6; |
203 | | T2A = T4 - T7; |
204 | | T3Y = T2 - T3; |
205 | | } |
206 | | { |
207 | | E T1K, T1L, T1M, T1N, T1O, T1P; |
208 | | T1K = ci[WS(rs, 19)]; |
209 | | T1L = cr[WS(rs, 20)]; |
210 | | T1M = T1K - T1L; |
211 | | T1N = ci[WS(rs, 14)]; |
212 | | T1O = cr[WS(rs, 15)]; |
213 | | T1P = T1N - T1O; |
214 | | T1Q = T1M + T1P; |
215 | | T3W = T1M - T1P; |
216 | | T2C = T1K + T1L; |
217 | | T2D = T1N + T1O; |
218 | | } |
219 | | T9 = T1 + T8; |
220 | | T40 = FMA(KP618033988, T3Z, T3Y); |
221 | | T1R = T1J + T1Q; |
222 | | T6G = FNMS(KP618033988, T3Y, T3Z); |
223 | | T3V = FNMS(KP250000000, T1Q, T1J); |
224 | | T6F = FNMS(KP559016994, T3W, T3V); |
225 | | T3X = FMA(KP559016994, T3W, T3V); |
226 | | T6H = FNMS(KP951056516, T6G, T6F); |
227 | | { |
228 | | E T2E, T5V, T2B, T5U, T2z; |
229 | | T2E = FMA(KP618033988, T2D, T2C); |
230 | | T5V = FNMS(KP618033988, T2C, T2D); |
231 | | T2z = FNMS(KP250000000, T8, T1); |
232 | | T2B = FMA(KP559016994, T2A, T2z); |
233 | | T5U = FNMS(KP559016994, T2A, T2z); |
234 | | T2F = FNMS(KP951056516, T2E, T2B); |
235 | | T7n = FNMS(KP951056516, T5V, T5U); |
236 | | T4N = FMA(KP951056516, T2E, T2B); |
237 | | T5W = FMA(KP951056516, T5V, T5U); |
238 | | } |
239 | | } |
240 | | { |
241 | | E Ta, T2H, T6n, T2S, Th, T2G, TC, T3r, T5Y, T3C, TJ, T3q, Tj, T30, T6d; |
242 | | E T33, Tq, T32, T1u, T3v, T61, T3y, T1B, T3x, T12, T2L, T6k, T2O, T19, T2N; |
243 | | E T1b, T2W, T6g, T37, T1i, T2V, T1l, T3g, T68, T3j, T1s, T3i, Tt, T3c, T65; |
244 | | E T3n, TA, T3b; |
245 | | { |
246 | | E Tg, T2R, Td, T2Q; |
247 | | Ta = cr[WS(rs, 1)]; |
248 | | { |
249 | | E Te, Tf, Tb, Tc; |
250 | | Te = cr[WS(rs, 11)]; |
251 | | Tf = ci[WS(rs, 8)]; |
252 | | Tg = Te + Tf; |
253 | | T2R = Tf - Te; |
254 | | Tb = cr[WS(rs, 6)]; |
255 | | Tc = ci[WS(rs, 3)]; |
256 | | Td = Tb + Tc; |
257 | | T2Q = Tb - Tc; |
258 | | } |
259 | | T2H = Td - Tg; |
260 | | T6n = FMA(KP618033988, T2Q, T2R); |
261 | | T2S = FNMS(KP618033988, T2R, T2Q); |
262 | | Th = Td + Tg; |
263 | | T2G = FNMS(KP250000000, Th, Ta); |
264 | | } |
265 | | { |
266 | | E TI, T3B, TF, T3A; |
267 | | TC = cr[WS(rs, 3)]; |
268 | | { |
269 | | E TG, TH, TD, TE; |
270 | | TG = ci[WS(rs, 11)]; |
271 | | TH = ci[WS(rs, 6)]; |
272 | | TI = TG + TH; |
273 | | T3B = TG - TH; |
274 | | TD = cr[WS(rs, 8)]; |
275 | | TE = ci[WS(rs, 1)]; |
276 | | TF = TD + TE; |
277 | | T3A = TD - TE; |
278 | | } |
279 | | T3r = TI - TF; |
280 | | T5Y = FNMS(KP618033988, T3A, T3B); |
281 | | T3C = FMA(KP618033988, T3B, T3A); |
282 | | TJ = TF + TI; |
283 | | T3q = FNMS(KP250000000, TJ, TC); |
284 | | } |
285 | | { |
286 | | E Tp, T2Z, Tm, T2Y; |
287 | | Tj = cr[WS(rs, 4)]; |
288 | | { |
289 | | E Tn, To, Tk, Tl; |
290 | | Tn = ci[WS(rs, 10)]; |
291 | | To = ci[WS(rs, 5)]; |
292 | | Tp = Tn + To; |
293 | | T2Z = To - Tn; |
294 | | Tk = cr[WS(rs, 9)]; |
295 | | Tl = ci[0]; |
296 | | Tm = Tk + Tl; |
297 | | T2Y = Tl - Tk; |
298 | | } |
299 | | T30 = FMA(KP618033988, T2Z, T2Y); |
300 | | T6d = FNMS(KP618033988, T2Y, T2Z); |
301 | | T33 = Tm - Tp; |
302 | | Tq = Tm + Tp; |
303 | | T32 = FMS(KP250000000, Tq, Tj); |
304 | | } |
305 | | { |
306 | | E T1A, T3u, T1x, T3t; |
307 | | T1u = ci[WS(rs, 21)]; |
308 | | { |
309 | | E T1y, T1z, T1v, T1w; |
310 | | T1y = cr[WS(rs, 13)]; |
311 | | T1z = cr[WS(rs, 18)]; |
312 | | T1A = T1y + T1z; |
313 | | T3u = T1z - T1y; |
314 | | T1v = ci[WS(rs, 16)]; |
315 | | T1w = cr[WS(rs, 23)]; |
316 | | T1x = T1v - T1w; |
317 | | T3t = T1v + T1w; |
318 | | } |
319 | | T3v = FMA(KP618033988, T3u, T3t); |
320 | | T61 = FNMS(KP618033988, T3t, T3u); |
321 | | T3y = T1x + T1A; |
322 | | T1B = T1x - T1A; |
323 | | T3x = FMS(KP250000000, T1B, T1u); |
324 | | } |
325 | | { |
326 | | E T18, T2K, T15, T2J; |
327 | | T12 = ci[WS(rs, 23)]; |
328 | | { |
329 | | E T16, T17, T13, T14; |
330 | | T16 = ci[WS(rs, 13)]; |
331 | | T17 = cr[WS(rs, 16)]; |
332 | | T18 = T16 - T17; |
333 | | T2K = T16 + T17; |
334 | | T13 = ci[WS(rs, 18)]; |
335 | | T14 = cr[WS(rs, 21)]; |
336 | | T15 = T13 - T14; |
337 | | T2J = T13 + T14; |
338 | | } |
339 | | T2L = FMA(KP618033988, T2K, T2J); |
340 | | T6k = FNMS(KP618033988, T2J, T2K); |
341 | | T2O = T15 - T18; |
342 | | T19 = T15 + T18; |
343 | | T2N = FNMS(KP250000000, T19, T12); |
344 | | } |
345 | | { |
346 | | E T1h, T36, T1e, T35; |
347 | | T1b = ci[WS(rs, 20)]; |
348 | | { |
349 | | E T1f, T1g, T1c, T1d; |
350 | | T1f = cr[WS(rs, 14)]; |
351 | | T1g = cr[WS(rs, 19)]; |
352 | | T1h = T1f + T1g; |
353 | | T36 = T1g - T1f; |
354 | | T1c = ci[WS(rs, 15)]; |
355 | | T1d = cr[WS(rs, 24)]; |
356 | | T1e = T1c - T1d; |
357 | | T35 = T1c + T1d; |
358 | | } |
359 | | T2W = T1e + T1h; |
360 | | T6g = FNMS(KP618033988, T35, T36); |
361 | | T37 = FMA(KP618033988, T36, T35); |
362 | | T1i = T1e - T1h; |
363 | | T2V = FMS(KP250000000, T1i, T1b); |
364 | | } |
365 | | { |
366 | | E T1o, T3e, T1r, T3f; |
367 | | T1l = ci[WS(rs, 22)]; |
368 | | { |
369 | | E T1m, T1n, T1p, T1q; |
370 | | T1m = ci[WS(rs, 17)]; |
371 | | T1n = cr[WS(rs, 22)]; |
372 | | T1o = T1m - T1n; |
373 | | T3e = T1m + T1n; |
374 | | T1p = ci[WS(rs, 12)]; |
375 | | T1q = cr[WS(rs, 17)]; |
376 | | T1r = T1p - T1q; |
377 | | T3f = T1p + T1q; |
378 | | } |
379 | | T3g = FMA(KP618033988, T3f, T3e); |
380 | | T68 = FNMS(KP618033988, T3e, T3f); |
381 | | T3j = T1o - T1r; |
382 | | T1s = T1o + T1r; |
383 | | T3i = FMS(KP250000000, T1s, T1l); |
384 | | } |
385 | | { |
386 | | E Tw, T3l, Tz, T3m; |
387 | | Tt = cr[WS(rs, 2)]; |
388 | | { |
389 | | E Tu, Tv, Tx, Ty; |
390 | | Tu = cr[WS(rs, 7)]; |
391 | | Tv = ci[WS(rs, 2)]; |
392 | | Tw = Tu + Tv; |
393 | | T3l = Tu - Tv; |
394 | | Tx = cr[WS(rs, 12)]; |
395 | | Ty = ci[WS(rs, 7)]; |
396 | | Tz = Tx + Ty; |
397 | | T3m = Ty - Tx; |
398 | | } |
399 | | T3c = Tz - Tw; |
400 | | T65 = FMA(KP618033988, T3l, T3m); |
401 | | T3n = FNMS(KP618033988, T3m, T3l); |
402 | | TA = Tw + Tz; |
403 | | T3b = FNMS(KP250000000, TA, Tt); |
404 | | } |
405 | | { |
406 | | E T1a, T1j, T1t, T1C; |
407 | | T1a = T12 + T19; |
408 | | T1j = T1b + T1i; |
409 | | T1k = T1a - T1j; |
410 | | T1S = T1a + T1j; |
411 | | T1t = T1l + T1s; |
412 | | T1C = T1u + T1B; |
413 | | T1D = T1t - T1C; |
414 | | T1T = T1t + T1C; |
415 | | } |
416 | | Ti = Ta + Th; |
417 | | Tr = Tj + Tq; |
418 | | Ts = Ti + Tr; |
419 | | TB = Tt + TA; |
420 | | TK = TC + TJ; |
421 | | TL = TB + TK; |
422 | | TM = Ts + TL; |
423 | | { |
424 | | E T6l, T7y, T6o, T7z, T6j, T6m; |
425 | | T6j = FNMS(KP559016994, T2H, T2G); |
426 | | T6l = FMA(KP951056516, T6k, T6j); |
427 | | T7y = FNMS(KP951056516, T6k, T6j); |
428 | | T6m = FNMS(KP559016994, T2O, T2N); |
429 | | T6o = FMA(KP951056516, T6n, T6m); |
430 | | T7z = FNMS(KP951056516, T6n, T6m); |
431 | | T6p = FNMS(KP549754652, T6o, T6l); |
432 | | T7K = FMA(KP939062505, T7y, T7z); |
433 | | T6w = FMA(KP549754652, T6l, T6o); |
434 | | T7A = FNMS(KP939062505, T7z, T7y); |
435 | | } |
436 | | { |
437 | | E T2M, T4W, T2T, T4V, T2I, T2P; |
438 | | T2I = FMA(KP559016994, T2H, T2G); |
439 | | T2M = FNMS(KP951056516, T2L, T2I); |
440 | | T4W = FMA(KP951056516, T2L, T2I); |
441 | | T2P = FMA(KP559016994, T2O, T2N); |
442 | | T2T = FMA(KP951056516, T2S, T2P); |
443 | | T4V = FNMS(KP951056516, T2S, T2P); |
444 | | T2U = FNMS(KP256756360, T2T, T2M); |
445 | | T56 = FMA(KP634619297, T4V, T4W); |
446 | | T3K = FMA(KP256756360, T2M, T2T); |
447 | | T4X = FNMS(KP634619297, T4W, T4V); |
448 | | } |
449 | | { |
450 | | E T6e, T7w, T6h, T7v, T6c, T6f; |
451 | | T6c = FMA(KP559016994, T2W, T2V); |
452 | | T6e = FNMS(KP951056516, T6d, T6c); |
453 | | T7w = FMA(KP951056516, T6d, T6c); |
454 | | T6f = FMA(KP559016994, T33, T32); |
455 | | T6h = FNMS(KP951056516, T6g, T6f); |
456 | | T7v = FMA(KP951056516, T6g, T6f); |
457 | | T6i = FMA(KP470564281, T6h, T6e); |
458 | | T7J = FNMS(KP126329378, T7v, T7w); |
459 | | T6v = FNMS(KP470564281, T6e, T6h); |
460 | | T7x = FMA(KP126329378, T7w, T7v); |
461 | | } |
462 | | { |
463 | | E T31, T4Y, T38, T4Z, T2X, T34; |
464 | | T2X = FNMS(KP559016994, T2W, T2V); |
465 | | T31 = FMA(KP951056516, T30, T2X); |
466 | | T4Y = FNMS(KP951056516, T30, T2X); |
467 | | T34 = FNMS(KP559016994, T33, T32); |
468 | | T38 = FMA(KP951056516, T37, T34); |
469 | | T4Z = FNMS(KP951056516, T37, T34); |
470 | | T39 = FNMS(KP634619297, T38, T31); |
471 | | T57 = FMA(KP827271945, T4Y, T4Z); |
472 | | T3L = FMA(KP634619297, T31, T38); |
473 | | T50 = FNMS(KP827271945, T4Z, T4Y); |
474 | | } |
475 | | { |
476 | | E T3w, T4O, T3D, T4P, T3s, T3z; |
477 | | T3s = FNMS(KP559016994, T3r, T3q); |
478 | | T3w = FNMS(KP951056516, T3v, T3s); |
479 | | T4O = FMA(KP951056516, T3v, T3s); |
480 | | T3z = FNMS(KP559016994, T3y, T3x); |
481 | | T3D = FNMS(KP951056516, T3C, T3z); |
482 | | T4P = FMA(KP951056516, T3C, T3z); |
483 | | T3E = FMA(KP939062505, T3D, T3w); |
484 | | T59 = FMA(KP126329378, T4O, T4P); |
485 | | T3O = FNMS(KP939062505, T3w, T3D); |
486 | | T4Q = FNMS(KP126329378, T4P, T4O); |
487 | | } |
488 | | { |
489 | | E T5Z, T7r, T62, T7s, T5X, T60; |
490 | | T5X = FMA(KP559016994, T3y, T3x); |
491 | | T5Z = FMA(KP951056516, T5Y, T5X); |
492 | | T7r = FNMS(KP951056516, T5Y, T5X); |
493 | | T60 = FMA(KP559016994, T3r, T3q); |
494 | | T62 = FMA(KP951056516, T61, T60); |
495 | | T7s = FNMS(KP951056516, T61, T60); |
496 | | T63 = FMA(KP062914667, T62, T5Z); |
497 | | T7H = FMA(KP827271945, T7r, T7s); |
498 | | T6y = FNMS(KP062914667, T5Z, T62); |
499 | | T7t = FNMS(KP827271945, T7s, T7r); |
500 | | } |
501 | | { |
502 | | E T3h, T4S, T3o, T4R, T3d, T3k; |
503 | | T3d = FNMS(KP559016994, T3c, T3b); |
504 | | T3h = FNMS(KP951056516, T3g, T3d); |
505 | | T4S = FMA(KP951056516, T3g, T3d); |
506 | | T3k = FNMS(KP559016994, T3j, T3i); |
507 | | T3o = FNMS(KP951056516, T3n, T3k); |
508 | | T4R = FMA(KP951056516, T3n, T3k); |
509 | | T3p = FMA(KP549754652, T3o, T3h); |
510 | | T5a = FMA(KP470564281, T4R, T4S); |
511 | | T3N = FNMS(KP549754652, T3h, T3o); |
512 | | T4T = FNMS(KP470564281, T4S, T4R); |
513 | | } |
514 | | { |
515 | | E T66, T7o, T69, T7p, T64, T67; |
516 | | T64 = FMA(KP559016994, T3j, T3i); |
517 | | T66 = FNMS(KP951056516, T65, T64); |
518 | | T7o = FMA(KP951056516, T65, T64); |
519 | | T67 = FMA(KP559016994, T3c, T3b); |
520 | | T69 = FMA(KP951056516, T68, T67); |
521 | | T7p = FNMS(KP951056516, T68, T67); |
522 | | T6a = FMA(KP634619297, T69, T66); |
523 | | T7G = FNMS(KP062914667, T7o, T7p); |
524 | | T6z = FNMS(KP634619297, T66, T69); |
525 | | T7q = FMA(KP062914667, T7p, T7o); |
526 | | } |
527 | | } |
528 | | cr[0] = T9 + TM; |
529 | | { |
530 | | E T1U, T1X, T2l, T20, T2m, T1F, T2r, T27, T2h; |
531 | | { |
532 | | E T1W, T1V, T1Y, T1Z; |
533 | | T1W = T1S - T1T; |
534 | | T1U = T1S + T1T; |
535 | | T1V = FNMS(KP250000000, T1U, T1R); |
536 | | T1X = FMA(KP559016994, T1W, T1V); |
537 | | T2l = FNMS(KP559016994, T1W, T1V); |
538 | | T1Y = Ti - Tr; |
539 | | T1Z = TB - TK; |
540 | | T20 = FMA(KP618033988, T1Z, T1Y); |
541 | | T2m = FNMS(KP618033988, T1Y, T1Z); |
542 | | { |
543 | | E T1E, T2g, T11, T2f, TZ, T10; |
544 | | T1E = FMA(KP618033988, T1D, T1k); |
545 | | T2g = FNMS(KP618033988, T1k, T1D); |
546 | | TZ = FNMS(KP250000000, TM, T9); |
547 | | T10 = Ts - TL; |
548 | | T11 = FMA(KP559016994, T10, TZ); |
549 | | T2f = FNMS(KP559016994, T10, TZ); |
550 | | T1F = FNMS(KP951056516, T1E, T11); |
551 | | T2r = FNMS(KP951056516, T2g, T2f); |
552 | | T27 = FMA(KP951056516, T1E, T11); |
553 | | T2h = FMA(KP951056516, T2g, T2f); |
554 | | } |
555 | | } |
556 | | { |
557 | | E T2s, T2w, T2v, T1G, T22, T21; |
558 | | ci[0] = T1R + T1U; |
559 | | T2s = T2q * T2r; |
560 | | T2w = T2u * T2r; |
561 | | T2v = FMA(KP951056516, T2m, T2l); |
562 | | cr[WS(rs, 15)] = FNMS(T2u, T2v, T2s); |
563 | | ci[WS(rs, 15)] = FMA(T2q, T2v, T2w); |
564 | | T1G = TY * T1F; |
565 | | T22 = T1I * T1F; |
566 | | T21 = FMA(KP951056516, T20, T1X); |
567 | | cr[WS(rs, 5)] = FNMS(T1I, T21, T1G); |
568 | | ci[WS(rs, 5)] = FMA(TY, T21, T22); |
569 | | { |
570 | | E T28, T2c, T2b, T2i, T2o, T2n; |
571 | | T28 = T26 * T27; |
572 | | T2c = T2a * T27; |
573 | | T2b = FNMS(KP951056516, T20, T1X); |
574 | | cr[WS(rs, 20)] = FNMS(T2a, T2b, T28); |
575 | | ci[WS(rs, 20)] = FMA(T26, T2b, T2c); |
576 | | T2i = T2e * T2h; |
577 | | T2o = T2k * T2h; |
578 | | T2n = FNMS(KP951056516, T2m, T2l); |
579 | | cr[WS(rs, 10)] = FNMS(T2k, T2n, T2i); |
580 | | ci[WS(rs, 10)] = FMA(T2e, T2n, T2o); |
581 | | } |
582 | | } |
583 | | } |
584 | | { |
585 | | E T6B, T73, T6Q, T78, T7j, T6u, T72, T7l, T6N, T77, T7k, T7m; |
586 | | { |
587 | | E T6x, T6A, T6O, T6P; |
588 | | T6x = FMA(KP968479752, T6w, T6v); |
589 | | T6A = FNMS(KP845997307, T6z, T6y); |
590 | | T6B = FNMS(KP681693190, T6A, T6x); |
591 | | T73 = FMA(KP560319534, T6x, T6A); |
592 | | T6O = FNMS(KP968479752, T6p, T6i); |
593 | | T6P = FNMS(KP845997307, T6a, T63); |
594 | | T6Q = FMA(KP681693190, T6P, T6O); |
595 | | T78 = FNMS(KP560319534, T6O, T6P); |
596 | | } |
597 | | { |
598 | | E T6r, T6t, T6b, T6q, T6s; |
599 | | T6b = FMA(KP845997307, T6a, T63); |
600 | | T6q = FMA(KP968479752, T6p, T6i); |
601 | | T6r = FMA(KP906616052, T6q, T6b); |
602 | | T6t = FNMS(KP906616052, T6q, T6b); |
603 | | T7j = FMA(KP998026728, T6r, T5W); |
604 | | T6s = FNMS(KP249506682, T6r, T5W); |
605 | | T6u = FNMS(KP557913902, T6t, T6s); |
606 | | T72 = FMA(KP557913902, T6t, T6s); |
607 | | } |
608 | | { |
609 | | E T6K, T6M, T6I, T6J, T6L; |
610 | | T6I = FMA(KP845997307, T6z, T6y); |
611 | | T6J = FNMS(KP968479752, T6w, T6v); |
612 | | T6K = FNMS(KP906616052, T6J, T6I); |
613 | | T6M = FMA(KP906616052, T6J, T6I); |
614 | | T7l = FMA(KP998026728, T6K, T6H); |
615 | | T6L = FNMS(KP249506682, T6K, T6H); |
616 | | T6N = FNMS(KP557913902, T6M, T6L); |
617 | | T77 = FMA(KP557913902, T6M, T6L); |
618 | | } |
619 | | T7k = T4l * T7j; |
620 | | cr[WS(rs, 2)] = FNMS(T4n, T7l, T7k); |
621 | | T7m = T4l * T7l; |
622 | | ci[WS(rs, 2)] = FMA(T4n, T7j, T7m); |
623 | | { |
624 | | E T6C, T6D, T6R, T6S; |
625 | | T6C = FNMS(KP860541664, T6B, T6u); |
626 | | T6D = T5T * T6C; |
627 | | T6R = FNMS(KP860541664, T6Q, T6N); |
628 | | T6S = T5T * T6R; |
629 | | cr[WS(rs, 7)] = FNMS(T6E, T6R, T6D); |
630 | | ci[WS(rs, 7)] = FMA(T6E, T6C, T6S); |
631 | | } |
632 | | { |
633 | | E T7d, T7e, T7h, T7i; |
634 | | T7d = FMA(KP949179823, T73, T72); |
635 | | T7e = T7c * T7d; |
636 | | T7h = FNMS(KP949179823, T78, T77); |
637 | | T7i = T7c * T7h; |
638 | | cr[WS(rs, 17)] = FNMS(T7g, T7h, T7e); |
639 | | ci[WS(rs, 17)] = FMA(T7g, T7d, T7i); |
640 | | } |
641 | | { |
642 | | E T74, T75, T79, T7a; |
643 | | T74 = FNMS(KP949179823, T73, T72); |
644 | | T75 = T71 * T74; |
645 | | T79 = FMA(KP949179823, T78, T77); |
646 | | T7a = T71 * T79; |
647 | | cr[WS(rs, 12)] = FNMS(T76, T79, T75); |
648 | | ci[WS(rs, 12)] = FMA(T76, T74, T7a); |
649 | | } |
650 | | { |
651 | | E T6V, T6W, T6Z, T70; |
652 | | T6V = FMA(KP860541664, T6B, T6u); |
653 | | T6W = T6U * T6V; |
654 | | T6Z = FMA(KP860541664, T6Q, T6N); |
655 | | T70 = T6U * T6Z; |
656 | | cr[WS(rs, 22)] = FNMS(T6Y, T6Z, T6W); |
657 | | ci[WS(rs, 22)] = FMA(T6Y, T6V, T70); |
658 | | } |
659 | | } |
660 | | { |
661 | | E T7U, T8f, T82, T8k, T7F, T7M, T7X, T7Y, T7D, T7R, T8e, T7I, T7L, T7E, T7O; |
662 | | E T7N; |
663 | | { |
664 | | E T7S, T7T, T80, T81; |
665 | | T7S = FNMS(KP734762448, T7K, T7J); |
666 | | T7T = FNMS(KP772036680, T7H, T7G); |
667 | | T7U = FNMS(KP621716863, T7T, T7S); |
668 | | T8f = FMA(KP614372930, T7S, T7T); |
669 | | T80 = FNMS(KP734762448, T7A, T7x); |
670 | | T81 = FNMS(KP772036680, T7t, T7q); |
671 | | T82 = FNMS(KP621716863, T81, T80); |
672 | | T8k = FMA(KP614372930, T80, T81); |
673 | | } |
674 | | T7F = FMA(KP951056516, T6G, T6F); |
675 | | T7I = FMA(KP772036680, T7H, T7G); |
676 | | T7L = FMA(KP734762448, T7K, T7J); |
677 | | T7M = FMA(KP994076283, T7L, T7I); |
678 | | T7X = FNMS(KP249506682, T7M, T7F); |
679 | | T7Y = FNMS(KP994076283, T7L, T7I); |
680 | | { |
681 | | E T7C, T7Q, T7u, T7B, T7P; |
682 | | T7u = FMA(KP772036680, T7t, T7q); |
683 | | T7B = FMA(KP734762448, T7A, T7x); |
684 | | T7C = FMA(KP994076283, T7B, T7u); |
685 | | T7Q = FNMS(KP994076283, T7B, T7u); |
686 | | T7D = FMA(KP998026728, T7C, T7n); |
687 | | T7P = FNMS(KP249506682, T7C, T7n); |
688 | | T7R = FNMS(KP557913902, T7Q, T7P); |
689 | | T8e = FMA(KP557913902, T7Q, T7P); |
690 | | } |
691 | | T7E = TO * T7D; |
692 | | T7O = TR * T7D; |
693 | | T7N = FMA(KP998026728, T7M, T7F); |
694 | | cr[WS(rs, 3)] = FNMS(TR, T7N, T7E); |
695 | | ci[WS(rs, 3)] = FMA(TO, T7N, T7O); |
696 | | { |
697 | | E T8l, T8t, T8q, T8u, T8h, T8m, T8j, T8p, T8g; |
698 | | T8j = FMA(KP557913902, T7Y, T7X); |
699 | | T8l = FNMS(KP949179823, T8k, T8j); |
700 | | T8t = FMA(KP949179823, T8k, T8j); |
701 | | T8p = FNMS(KP949179823, T8f, T8e); |
702 | | T8q = T8o * T8p; |
703 | | T8u = T8s * T8p; |
704 | | T8g = FMA(KP949179823, T8f, T8e); |
705 | | T8h = T8d * T8g; |
706 | | T8m = T8i * T8g; |
707 | | cr[WS(rs, 13)] = FNMS(T8i, T8l, T8h); |
708 | | ci[WS(rs, 13)] = FMA(T8d, T8l, T8m); |
709 | | cr[WS(rs, 18)] = FNMS(T8s, T8t, T8q); |
710 | | ci[WS(rs, 18)] = FMA(T8o, T8t, T8u); |
711 | | } |
712 | | { |
713 | | E T83, T8b, T88, T8c, T7W, T84, T7Z, T87, T7V; |
714 | | T7Z = FNMS(KP557913902, T7Y, T7X); |
715 | | T83 = FNMS(KP943557151, T82, T7Z); |
716 | | T8b = FMA(KP943557151, T82, T7Z); |
717 | | T87 = FNMS(KP943557151, T7U, T7R); |
718 | | T88 = T86 * T87; |
719 | | T8c = T8a * T87; |
720 | | T7V = FMA(KP943557151, T7U, T7R); |
721 | | T7W = T4z * T7V; |
722 | | T84 = T4B * T7V; |
723 | | cr[WS(rs, 8)] = FNMS(T4B, T83, T7W); |
724 | | ci[WS(rs, 8)] = FMA(T4z, T83, T84); |
725 | | cr[WS(rs, 23)] = FNMS(T8a, T8b, T88); |
726 | | ci[WS(rs, 23)] = FMA(T86, T8b, T8c); |
727 | | } |
728 | | } |
729 | | { |
730 | | E T5c, T5y, T5o, T5E, T5f, T5i, T5j, T5k, T5P, T55, T5x, T5g, T5h, T5Q, T5S; |
731 | | E T5R; |
732 | | { |
733 | | E T58, T5b, T5m, T5n; |
734 | | T58 = FNMS(KP912575812, T57, T56); |
735 | | T5b = FNMS(KP912018591, T5a, T59); |
736 | | T5c = FNMS(KP726211448, T5b, T58); |
737 | | T5y = FMA(KP525970792, T58, T5b); |
738 | | T5m = FNMS(KP912575812, T50, T4X); |
739 | | T5n = FMA(KP912018591, T4T, T4Q); |
740 | | T5o = FNMS(KP726211448, T5n, T5m); |
741 | | T5E = FMA(KP525970792, T5m, T5n); |
742 | | } |
743 | | T5f = FNMS(KP951056516, T40, T3X); |
744 | | T5g = FMA(KP912018591, T5a, T59); |
745 | | T5h = FMA(KP912575812, T57, T56); |
746 | | T5i = FMA(KP851038619, T5h, T5g); |
747 | | T5j = FNMS(KP248028675, T5i, T5f); |
748 | | T5k = FNMS(KP851038619, T5h, T5g); |
749 | | { |
750 | | E T52, T54, T4U, T51, T53; |
751 | | T4U = FNMS(KP912018591, T4T, T4Q); |
752 | | T51 = FMA(KP912575812, T50, T4X); |
753 | | T52 = FMA(KP851038619, T51, T4U); |
754 | | T54 = FNMS(KP851038619, T51, T4U); |
755 | | T5P = FNMS(KP992114701, T52, T4N); |
756 | | T53 = FMA(KP248028675, T52, T4N); |
757 | | T55 = FMA(KP554608978, T54, T53); |
758 | | T5x = FNMS(KP554608978, T54, T53); |
759 | | } |
760 | | T5Q = TS * T5P; |
761 | | T5S = TW * T5P; |
762 | | T5R = FMA(KP992114701, T5i, T5f); |
763 | | cr[WS(rs, 4)] = FNMS(TW, T5R, T5Q); |
764 | | ci[WS(rs, 4)] = FMA(TS, T5R, T5S); |
765 | | { |
766 | | E T5F, T5N, T5K, T5O, T5A, T5G, T5D, T5J, T5z; |
767 | | T5D = FMA(KP554608978, T5k, T5j); |
768 | | T5F = FNMS(KP943557151, T5E, T5D); |
769 | | T5N = FMA(KP943557151, T5E, T5D); |
770 | | T5J = FMA(KP943557151, T5y, T5x); |
771 | | T5K = T5I * T5J; |
772 | | T5O = T5M * T5J; |
773 | | T5z = FNMS(KP943557151, T5y, T5x); |
774 | | T5A = T5w * T5z; |
775 | | T5G = T5C * T5z; |
776 | | cr[WS(rs, 14)] = FNMS(T5C, T5F, T5A); |
777 | | ci[WS(rs, 14)] = FMA(T5w, T5F, T5G); |
778 | | cr[WS(rs, 19)] = FNMS(T5M, T5N, T5K); |
779 | | ci[WS(rs, 19)] = FMA(T5I, T5N, T5O); |
780 | | } |
781 | | { |
782 | | E T5p, T5t, T5s, T5u, T5e, T5q, T5l, T5r, T5d; |
783 | | T5l = FNMS(KP554608978, T5k, T5j); |
784 | | T5p = FNMS(KP803003575, T5o, T5l); |
785 | | T5t = FMA(KP803003575, T5o, T5l); |
786 | | T5r = FMA(KP803003575, T5c, T55); |
787 | | T5s = T23 * T5r; |
788 | | T5u = T25 * T5r; |
789 | | T5d = FNMS(KP803003575, T5c, T55); |
790 | | T5e = TT * T5d; |
791 | | T5q = TX * T5d; |
792 | | cr[WS(rs, 9)] = FNMS(TX, T5p, T5e); |
793 | | ci[WS(rs, 9)] = FMA(TT, T5p, T5q); |
794 | | cr[WS(rs, 24)] = FNMS(T25, T5t, T5s); |
795 | | ci[WS(rs, 24)] = FMA(T23, T5t, T5u); |
796 | | } |
797 | | } |
798 | | { |
799 | | E T3Q, T4q, T4a, T4w, T41, T44, T45, T46, T4J, T3J, T4p, T42, T43, T4K, T4M; |
800 | | E T4L; |
801 | | { |
802 | | E T3M, T3P, T48, T49; |
803 | | T3M = FMA(KP871714437, T3L, T3K); |
804 | | T3P = FNMS(KP831864738, T3O, T3N); |
805 | | T3Q = FNMS(KP559154169, T3P, T3M); |
806 | | T4q = FMA(KP683113946, T3M, T3P); |
807 | | T48 = FNMS(KP871714437, T39, T2U); |
808 | | T49 = FNMS(KP831864738, T3E, T3p); |
809 | | T4a = FMA(KP559154169, T49, T48); |
810 | | T4w = FNMS(KP683113946, T48, T49); |
811 | | } |
812 | | T41 = FMA(KP951056516, T40, T3X); |
813 | | T42 = FNMS(KP871714437, T3L, T3K); |
814 | | T43 = FMA(KP831864738, T3O, T3N); |
815 | | T44 = FNMS(KP904730450, T43, T42); |
816 | | T45 = FNMS(KP242145790, T44, T41); |
817 | | T46 = FMA(KP904730450, T43, T42); |
818 | | { |
819 | | E T3G, T3I, T3a, T3F, T3H; |
820 | | T3a = FMA(KP871714437, T39, T2U); |
821 | | T3F = FMA(KP831864738, T3E, T3p); |
822 | | T3G = FMA(KP904730450, T3F, T3a); |
823 | | T3I = FNMS(KP904730450, T3F, T3a); |
824 | | T4J = FMA(KP968583161, T3G, T2F); |
825 | | T3H = FNMS(KP242145790, T3G, T2F); |
826 | | T3J = FMA(KP541454447, T3I, T3H); |
827 | | T4p = FNMS(KP541454447, T3I, T3H); |
828 | | } |
829 | | T4K = TN * T4J; |
830 | | T4M = TQ * T4J; |
831 | | T4L = FMA(KP968583161, T44, T41); |
832 | | cr[WS(rs, 1)] = FNMS(TQ, T4L, T4K); |
833 | | ci[WS(rs, 1)] = FMA(TN, T4L, T4M); |
834 | | { |
835 | | E T4x, T4H, T4E, T4I, T4s, T4y, T4v, T4D, T4r; |
836 | | T4v = FNMS(KP541454447, T46, T45); |
837 | | T4x = FNMS(KP833417178, T4w, T4v); |
838 | | T4H = FMA(KP833417178, T4w, T4v); |
839 | | T4D = FMA(KP833417178, T4q, T4p); |
840 | | T4E = T4C * T4D; |
841 | | T4I = T4G * T4D; |
842 | | T4r = FNMS(KP833417178, T4q, T4p); |
843 | | T4s = T4o * T4r; |
844 | | T4y = T4u * T4r; |
845 | | cr[WS(rs, 11)] = FNMS(T4u, T4x, T4s); |
846 | | ci[WS(rs, 11)] = FMA(T4o, T4x, T4y); |
847 | | cr[WS(rs, 16)] = FNMS(T4G, T4H, T4E); |
848 | | ci[WS(rs, 16)] = FMA(T4C, T4H, T4I); |
849 | | } |
850 | | { |
851 | | E T4b, T4j, T4g, T4k, T3S, T4c, T47, T4f, T3R; |
852 | | T47 = FMA(KP541454447, T46, T45); |
853 | | T4b = FMA(KP921177326, T4a, T47); |
854 | | T4j = FNMS(KP921177326, T4a, T47); |
855 | | T4f = FMA(KP921177326, T3Q, T3J); |
856 | | T4g = T4e * T4f; |
857 | | T4k = T4i * T4f; |
858 | | T3R = FNMS(KP921177326, T3Q, T3J); |
859 | | T3S = T2y * T3R; |
860 | | T4c = T3U * T3R; |
861 | | cr[WS(rs, 6)] = FNMS(T3U, T4b, T3S); |
862 | | ci[WS(rs, 6)] = FMA(T2y, T4b, T4c); |
863 | | cr[WS(rs, 21)] = FNMS(T4i, T4j, T4g); |
864 | | ci[WS(rs, 21)] = FMA(T4e, T4j, T4k); |
865 | | } |
866 | | } |
867 | | } |
868 | | } |
869 | | } |
870 | | } |
871 | | |
872 | | static const tw_instr twinstr[] = { |
873 | | { TW_CEXP, 1, 1 }, |
874 | | { TW_CEXP, 1, 3 }, |
875 | | { TW_CEXP, 1, 9 }, |
876 | | { TW_CEXP, 1, 24 }, |
877 | | { TW_NEXT, 1, 0 } |
878 | | }; |
879 | | |
880 | | static const hc2hc_desc desc = { 25, "hb2_25", twinstr, &GENUS, { 84, 78, 356, 0 } }; |
881 | | |
882 | | void X(codelet_hb2_25) (planner *p) { |
883 | | X(khc2hc_register) (p, hb2_25, &desc); |
884 | | } |
885 | | #else |
886 | | |
887 | | /* Generated by: ../../../genfft/gen_hc2hc.native -compact -variables 4 -pipeline-latency 4 -sign 1 -twiddle-log3 -precompute-twiddles -n 25 -dif -name hb2_25 -include rdft/scalar/hb.h */ |
888 | | |
889 | | /* |
890 | | * This function contains 440 FP additions, 340 FP multiplications, |
891 | | * (or, 280 additions, 180 multiplications, 160 fused multiply/add), |
892 | | * 155 stack variables, 20 constants, and 100 memory accesses |
893 | | */ |
894 | | #include "rdft/scalar/hb.h" |
895 | | |
896 | | static void hb2_25(R *cr, R *ci, const R *W, stride rs, INT mb, INT me, INT ms) |
897 | 0 | { |
898 | 0 | DK(KP998026728, +0.998026728428271561952336806863450553336905220); |
899 | 0 | DK(KP062790519, +0.062790519529313376076178224565631133122484832); |
900 | 0 | DK(KP992114701, +0.992114701314477831049793042785778521453036709); |
901 | 0 | DK(KP125333233, +0.125333233564304245373118759816508793942918247); |
902 | 0 | DK(KP425779291, +0.425779291565072648862502445744251703979973042); |
903 | 0 | DK(KP904827052, +0.904827052466019527713668647932697593970413911); |
904 | 0 | DK(KP248689887, +0.248689887164854788242283746006447968417567406); |
905 | 0 | DK(KP968583161, +0.968583161128631119490168375464735813836012403); |
906 | 0 | DK(KP770513242, +0.770513242775789230803009636396177847271667672); |
907 | 0 | DK(KP637423989, +0.637423989748689710176712811676016195434917298); |
908 | 0 | DK(KP844327925, +0.844327925502015078548558063966681505381659241); |
909 | 0 | DK(KP535826794, +0.535826794978996618271308767867639978063575346); |
910 | 0 | DK(KP684547105, +0.684547105928688673732283357621209269889519233); |
911 | 0 | DK(KP728968627, +0.728968627421411523146730319055259111372571664); |
912 | 0 | DK(KP481753674, +0.481753674101715274987191502872129653528542010); |
913 | 0 | DK(KP876306680, +0.876306680043863587308115903922062583399064238); |
914 | 0 | DK(KP559016994, +0.559016994374947424102293417182819058860154590); |
915 | 0 | DK(KP250000000, +0.250000000000000000000000000000000000000000000); |
916 | 0 | DK(KP587785252, +0.587785252292473129168705954639072768597652438); |
917 | 0 | DK(KP951056516, +0.951056516295153572116439333379382143405698634); |
918 | 0 | { |
919 | 0 | INT m; |
920 | 0 | for (m = mb, W = W + ((mb - 1) * 8); m < me; m = m + 1, cr = cr + ms, ci = ci - ms, W = W + 8, MAKE_VOLATILE_STRIDE(50, rs)) { |
921 | 0 | E TN, TQ, TO, TR, TT, TY, T2t, T2r, TZ, TU, T4f, T4l, T2d, T4v, T5m; |
922 | 0 | E T2j, T5l, T4X, T2v, T11, T3R, T1L, T5d, T6x, T5h, T6t, T25, T26, T27, T29; |
923 | 0 | E T6D, T7v, T49, T7l, T7p, T7t, T2p, T2n, T4b, T4p, T5n, T6B, T5b, T5p, T6p; |
924 | 0 | E T6r, T59, T4r; |
925 | 0 | { |
926 | 0 | E T2c, T4j, T2h, T4e, T2b, T4k, T2i, T4d; |
927 | 0 | { |
928 | 0 | E TP, TX, TS, TW; |
929 | 0 | TN = W[0]; |
930 | 0 | TQ = W[1]; |
931 | 0 | TO = W[2]; |
932 | 0 | TR = W[3]; |
933 | 0 | TP = TN * TO; |
934 | 0 | TX = TQ * TO; |
935 | 0 | TS = TQ * TR; |
936 | 0 | TW = TN * TR; |
937 | 0 | TT = TP - TS; |
938 | 0 | TY = TW + TX; |
939 | 0 | T2t = TW - TX; |
940 | 0 | T2r = TP + TS; |
941 | 0 | TZ = W[5]; |
942 | 0 | T2c = TQ * TZ; |
943 | 0 | T4j = TO * TZ; |
944 | 0 | T2h = TN * TZ; |
945 | 0 | T4e = TR * TZ; |
946 | 0 | TU = W[4]; |
947 | 0 | T2b = TN * TU; |
948 | 0 | T4k = TR * TU; |
949 | 0 | T2i = TQ * TU; |
950 | 0 | T4d = TO * TU; |
951 | 0 | } |
952 | 0 | T4f = T4d - T4e; |
953 | 0 | T4l = T4j + T4k; |
954 | 0 | { |
955 | 0 | E T2s, T2u, TV, T10, T3P, T3Q, T1J, T1K; |
956 | 0 | T2d = T2b - T2c; |
957 | 0 | T4v = T2b + T2c; |
958 | 0 | T5m = T4j - T4k; |
959 | 0 | T2j = T2h + T2i; |
960 | 0 | T5l = T4d + T4e; |
961 | 0 | T4X = T2h - T2i; |
962 | 0 | T2s = T2r * TU; |
963 | 0 | T2u = T2t * TZ; |
964 | 0 | T2v = T2s + T2u; |
965 | 0 | TV = TT * TU; |
966 | 0 | T10 = TY * TZ; |
967 | 0 | T11 = TV + T10; |
968 | 0 | T3P = T2r * TZ; |
969 | 0 | T3Q = T2t * TU; |
970 | 0 | T3R = T3P - T3Q; |
971 | 0 | T1J = TT * TZ; |
972 | 0 | T1K = TY * TU; |
973 | 0 | T1L = T1J - T1K; |
974 | 0 | T5d = TV - T10; |
975 | 0 | T6x = T3P + T3Q; |
976 | 0 | T5h = T1J + T1K; |
977 | 0 | T6t = T2s - T2u; |
978 | 0 | T25 = W[6]; |
979 | 0 | T26 = W[7]; |
980 | 0 | T27 = FMA(TT, T25, TY * T26); |
981 | 0 | T29 = FNMS(TY, T25, TT * T26); |
982 | 0 | T6D = FNMS(T4X, T25, T4v * T26); |
983 | 0 | T7v = FNMS(T1L, T25, T11 * T26); |
984 | 0 | T49 = FMA(T2r, T25, T2t * T26); |
985 | 0 | T7l = FMA(T2d, T25, T2j * T26); |
986 | 0 | T7p = FNMS(T2j, T25, T2d * T26); |
987 | 0 | T7t = FMA(T11, T25, T1L * T26); |
988 | 0 | T2p = FNMS(TZ, T25, TU * T26); |
989 | 0 | T2n = FMA(TU, T25, TZ * T26); |
990 | 0 | T4b = FNMS(T2t, T25, T2r * T26); |
991 | 0 | T4p = FMA(T2v, T25, T3R * T26); |
992 | 0 | T5n = FMA(T5l, T25, T5m * T26); |
993 | 0 | T6B = FMA(T4v, T25, T4X * T26); |
994 | 0 | T5b = FNMS(TQ, T25, TN * T26); |
995 | 0 | T5p = FNMS(T5m, T25, T5l * T26); |
996 | 0 | T6p = FMA(TO, T25, TR * T26); |
997 | 0 | T6r = FNMS(TR, T25, TO * T26); |
998 | 0 | T59 = FMA(TN, T25, TQ * T26); |
999 | 0 | T4r = FNMS(T3R, T25, T2v * T26); |
1000 | 0 | } |
1001 | 0 | } |
1002 | 0 | { |
1003 | 0 | E T9, T6i, T40, T3z, T5Y, Ti, Tr, Ts, T1d, T1m, T1P, T2K, T4P, T3H, T4y; |
1004 | 0 | E T5G, T71, T65, T6N, T5z, T70, T64, T6K, T2Z, T4Q, T3I, T4B, T20, T5Z, T3C; |
1005 | 0 | E T43, T6j, TB, TK, TL, T1w, T1F, T1Q, T3f, T4S, T3K, T4F, T5V, T74, T68; |
1006 | 0 | E T6U, T5O, T73, T67, T6R, T3u, T4T, T3L, T4I; |
1007 | 0 | { |
1008 | 0 | E T1, T4, T7, T8, T3Z, T3Y, T3x, T3y; |
1009 | 0 | T1 = cr[0]; |
1010 | 0 | { |
1011 | 0 | E T2, T3, T5, T6; |
1012 | 0 | T2 = cr[WS(rs, 5)]; |
1013 | 0 | T3 = ci[WS(rs, 4)]; |
1014 | 0 | T4 = T2 + T3; |
1015 | 0 | T5 = cr[WS(rs, 10)]; |
1016 | 0 | T6 = ci[WS(rs, 9)]; |
1017 | 0 | T7 = T5 + T6; |
1018 | 0 | T8 = T4 + T7; |
1019 | 0 | T3Z = T5 - T6; |
1020 | 0 | T3Y = T2 - T3; |
1021 | 0 | } |
1022 | 0 | T9 = T1 + T8; |
1023 | 0 | T6i = FMA(KP951056516, T3Y, KP587785252 * T3Z); |
1024 | 0 | T40 = FNMS(KP951056516, T3Z, KP587785252 * T3Y); |
1025 | 0 | T3x = FNMS(KP250000000, T8, T1); |
1026 | 0 | T3y = KP559016994 * (T4 - T7); |
1027 | 0 | T3z = T3x - T3y; |
1028 | 0 | T5Y = T3y + T3x; |
1029 | 0 | } |
1030 | 0 | { |
1031 | 0 | E Ta, T2x, T5w, T2F, Th, T2w, T1e, T2P, T5B, T2X, T1l, T2O, Tj, T2N, T5D; |
1032 | 0 | E T2T, Tq, T2S, T15, T2B, T5u, T2H, T1c, T2G; |
1033 | 0 | { |
1034 | 0 | E Tg, T2E, Td, T2D; |
1035 | 0 | Ta = cr[WS(rs, 1)]; |
1036 | 0 | { |
1037 | 0 | E Te, Tf, Tb, Tc; |
1038 | 0 | Te = cr[WS(rs, 11)]; |
1039 | 0 | Tf = ci[WS(rs, 8)]; |
1040 | 0 | Tg = Te + Tf; |
1041 | 0 | T2E = Te - Tf; |
1042 | 0 | Tb = cr[WS(rs, 6)]; |
1043 | 0 | Tc = ci[WS(rs, 3)]; |
1044 | 0 | Td = Tb + Tc; |
1045 | 0 | T2D = Tb - Tc; |
1046 | 0 | } |
1047 | 0 | T2x = KP559016994 * (Td - Tg); |
1048 | 0 | T5w = FMA(KP951056516, T2D, KP587785252 * T2E); |
1049 | 0 | T2F = FNMS(KP951056516, T2E, KP587785252 * T2D); |
1050 | 0 | Th = Td + Tg; |
1051 | 0 | T2w = FNMS(KP250000000, Th, Ta); |
1052 | 0 | } |
1053 | 0 | { |
1054 | 0 | E T1k, T2W, T1h, T2V; |
1055 | 0 | T1e = ci[WS(rs, 20)]; |
1056 | 0 | { |
1057 | 0 | E T1i, T1j, T1f, T1g; |
1058 | 0 | T1i = cr[WS(rs, 14)]; |
1059 | 0 | T1j = cr[WS(rs, 19)]; |
1060 | 0 | T1k = T1i + T1j; |
1061 | 0 | T2W = T1j - T1i; |
1062 | 0 | T1f = ci[WS(rs, 15)]; |
1063 | 0 | T1g = cr[WS(rs, 24)]; |
1064 | 0 | T1h = T1f - T1g; |
1065 | 0 | T2V = T1f + T1g; |
1066 | 0 | } |
1067 | 0 | T2P = KP559016994 * (T1h + T1k); |
1068 | 0 | T5B = FMA(KP951056516, T2V, KP587785252 * T2W); |
1069 | 0 | T2X = FNMS(KP951056516, T2W, KP587785252 * T2V); |
1070 | 0 | T1l = T1h - T1k; |
1071 | 0 | T2O = FNMS(KP250000000, T1l, T1e); |
1072 | 0 | } |
1073 | 0 | { |
1074 | 0 | E Tp, T2M, Tm, T2L; |
1075 | 0 | Tj = cr[WS(rs, 4)]; |
1076 | 0 | { |
1077 | 0 | E Tn, To, Tk, Tl; |
1078 | 0 | Tn = ci[WS(rs, 10)]; |
1079 | 0 | To = ci[WS(rs, 5)]; |
1080 | 0 | Tp = Tn + To; |
1081 | 0 | T2M = Tn - To; |
1082 | 0 | Tk = cr[WS(rs, 9)]; |
1083 | 0 | Tl = ci[0]; |
1084 | 0 | Tm = Tk + Tl; |
1085 | 0 | T2L = Tk - Tl; |
1086 | 0 | } |
1087 | 0 | T2N = FNMS(KP951056516, T2M, KP587785252 * T2L); |
1088 | 0 | T5D = FMA(KP951056516, T2L, KP587785252 * T2M); |
1089 | 0 | T2T = KP559016994 * (Tm - Tp); |
1090 | 0 | Tq = Tm + Tp; |
1091 | 0 | T2S = FNMS(KP250000000, Tq, Tj); |
1092 | 0 | } |
1093 | 0 | { |
1094 | 0 | E T1b, T2A, T18, T2z; |
1095 | 0 | T15 = ci[WS(rs, 23)]; |
1096 | 0 | { |
1097 | 0 | E T19, T1a, T16, T17; |
1098 | 0 | T19 = ci[WS(rs, 13)]; |
1099 | 0 | T1a = cr[WS(rs, 16)]; |
1100 | 0 | T1b = T19 - T1a; |
1101 | 0 | T2A = T19 + T1a; |
1102 | 0 | T16 = ci[WS(rs, 18)]; |
1103 | 0 | T17 = cr[WS(rs, 21)]; |
1104 | 0 | T18 = T16 - T17; |
1105 | 0 | T2z = T16 + T17; |
1106 | 0 | } |
1107 | 0 | T2B = FNMS(KP951056516, T2A, KP587785252 * T2z); |
1108 | 0 | T5u = FMA(KP951056516, T2z, KP587785252 * T2A); |
1109 | 0 | T2H = KP559016994 * (T18 - T1b); |
1110 | 0 | T1c = T18 + T1b; |
1111 | 0 | T2G = FNMS(KP250000000, T1c, T15); |
1112 | 0 | } |
1113 | 0 | Ti = Ta + Th; |
1114 | 0 | Tr = Tj + Tq; |
1115 | 0 | Ts = Ti + Tr; |
1116 | 0 | T1d = T15 + T1c; |
1117 | 0 | T1m = T1e + T1l; |
1118 | 0 | T1P = T1d + T1m; |
1119 | 0 | { |
1120 | 0 | E T2C, T4w, T2J, T4x, T2y, T2I; |
1121 | 0 | T2y = T2w - T2x; |
1122 | 0 | T2C = T2y - T2B; |
1123 | 0 | T4w = T2y + T2B; |
1124 | 0 | T2I = T2G - T2H; |
1125 | 0 | T2J = T2F + T2I; |
1126 | 0 | T4x = T2I - T2F; |
1127 | 0 | T2K = FNMS(KP481753674, T2J, KP876306680 * T2C); |
1128 | 0 | T4P = FMA(KP728968627, T4x, KP684547105 * T4w); |
1129 | 0 | T3H = FMA(KP876306680, T2J, KP481753674 * T2C); |
1130 | 0 | T4y = FNMS(KP684547105, T4x, KP728968627 * T4w); |
1131 | 0 | } |
1132 | 0 | { |
1133 | 0 | E T5C, T6M, T5F, T6L, T5A, T5E; |
1134 | 0 | T5A = T2T + T2S; |
1135 | 0 | T5C = T5A - T5B; |
1136 | 0 | T6M = T5A + T5B; |
1137 | 0 | T5E = T2O + T2P; |
1138 | 0 | T5F = T5D + T5E; |
1139 | 0 | T6L = T5E - T5D; |
1140 | 0 | T5G = FNMS(KP844327925, T5F, KP535826794 * T5C); |
1141 | 0 | T71 = FMA(KP637423989, T6L, KP770513242 * T6M); |
1142 | 0 | T65 = FMA(KP535826794, T5F, KP844327925 * T5C); |
1143 | 0 | T6N = FNMS(KP637423989, T6M, KP770513242 * T6L); |
1144 | 0 | } |
1145 | 0 | { |
1146 | 0 | E T5v, T6I, T5y, T6J, T5t, T5x; |
1147 | 0 | T5t = T2x + T2w; |
1148 | 0 | T5v = T5t - T5u; |
1149 | 0 | T6I = T5t + T5u; |
1150 | 0 | T5x = T2H + T2G; |
1151 | 0 | T5y = T5w + T5x; |
1152 | 0 | T6J = T5x - T5w; |
1153 | 0 | T5z = FNMS(KP248689887, T5y, KP968583161 * T5v); |
1154 | 0 | T70 = FMA(KP535826794, T6J, KP844327925 * T6I); |
1155 | 0 | T64 = FMA(KP968583161, T5y, KP248689887 * T5v); |
1156 | 0 | T6K = FNMS(KP844327925, T6J, KP535826794 * T6I); |
1157 | 0 | } |
1158 | 0 | { |
1159 | 0 | E T2R, T4z, T2Y, T4A, T2Q, T2U; |
1160 | 0 | T2Q = T2O - T2P; |
1161 | 0 | T2R = T2N + T2Q; |
1162 | 0 | T4z = T2Q - T2N; |
1163 | 0 | T2U = T2S - T2T; |
1164 | 0 | T2Y = T2U - T2X; |
1165 | 0 | T4A = T2U + T2X; |
1166 | 0 | T2Z = FMA(KP904827052, T2R, KP425779291 * T2Y); |
1167 | 0 | T4Q = FNMS(KP992114701, T4z, KP125333233 * T4A); |
1168 | 0 | T3I = FNMS(KP425779291, T2R, KP904827052 * T2Y); |
1169 | 0 | T4B = FMA(KP125333233, T4z, KP992114701 * T4A); |
1170 | 0 | } |
1171 | 0 | } |
1172 | 0 | { |
1173 | 0 | E T1S, T1V, T1Y, T1Z, T3B, T3A, T41, T42; |
1174 | 0 | T1S = ci[WS(rs, 24)]; |
1175 | 0 | { |
1176 | 0 | E T1T, T1U, T1W, T1X; |
1177 | 0 | T1T = ci[WS(rs, 19)]; |
1178 | 0 | T1U = cr[WS(rs, 20)]; |
1179 | 0 | T1V = T1T - T1U; |
1180 | 0 | T1W = ci[WS(rs, 14)]; |
1181 | 0 | T1X = cr[WS(rs, 15)]; |
1182 | 0 | T1Y = T1W - T1X; |
1183 | 0 | T1Z = T1V + T1Y; |
1184 | 0 | T3B = T1W + T1X; |
1185 | 0 | T3A = T1T + T1U; |
1186 | 0 | } |
1187 | 0 | T20 = T1S + T1Z; |
1188 | 0 | T5Z = FMA(KP951056516, T3A, KP587785252 * T3B); |
1189 | 0 | T3C = FNMS(KP951056516, T3B, KP587785252 * T3A); |
1190 | 0 | T41 = FNMS(KP250000000, T1Z, T1S); |
1191 | 0 | T42 = KP559016994 * (T1V - T1Y); |
1192 | 0 | T43 = T41 - T42; |
1193 | 0 | T6j = T42 + T41; |
1194 | 0 | } |
1195 | 0 | { |
1196 | 0 | E Tt, T32, T5L, T3a, TA, T31, T1o, T36, T5J, T3c, T1v, T3b, TC, T3h, T5S; |
1197 | 0 | E T3p, TJ, T3g, T1x, T3l, T5Q, T3r, T1E, T3q; |
1198 | 0 | { |
1199 | 0 | E Tw, T38, Tz, T39; |
1200 | 0 | Tt = cr[WS(rs, 2)]; |
1201 | 0 | { |
1202 | 0 | E Tu, Tv, Tx, Ty; |
1203 | 0 | Tu = cr[WS(rs, 7)]; |
1204 | 0 | Tv = ci[WS(rs, 2)]; |
1205 | 0 | Tw = Tu + Tv; |
1206 | 0 | T38 = Tu - Tv; |
1207 | 0 | Tx = cr[WS(rs, 12)]; |
1208 | 0 | Ty = ci[WS(rs, 7)]; |
1209 | 0 | Tz = Tx + Ty; |
1210 | 0 | T39 = Tx - Ty; |
1211 | 0 | } |
1212 | 0 | T32 = KP559016994 * (Tw - Tz); |
1213 | 0 | T5L = FMA(KP951056516, T38, KP587785252 * T39); |
1214 | 0 | T3a = FNMS(KP951056516, T39, KP587785252 * T38); |
1215 | 0 | TA = Tw + Tz; |
1216 | 0 | T31 = FNMS(KP250000000, TA, Tt); |
1217 | 0 | } |
1218 | 0 | { |
1219 | 0 | E T1r, T34, T1u, T35; |
1220 | 0 | T1o = ci[WS(rs, 22)]; |
1221 | 0 | { |
1222 | 0 | E T1p, T1q, T1s, T1t; |
1223 | 0 | T1p = ci[WS(rs, 17)]; |
1224 | 0 | T1q = cr[WS(rs, 22)]; |
1225 | 0 | T1r = T1p - T1q; |
1226 | 0 | T34 = T1p + T1q; |
1227 | 0 | T1s = ci[WS(rs, 12)]; |
1228 | 0 | T1t = cr[WS(rs, 17)]; |
1229 | 0 | T1u = T1s - T1t; |
1230 | 0 | T35 = T1s + T1t; |
1231 | 0 | } |
1232 | 0 | T36 = FNMS(KP951056516, T35, KP587785252 * T34); |
1233 | 0 | T5J = FMA(KP951056516, T34, KP587785252 * T35); |
1234 | 0 | T3c = KP559016994 * (T1r - T1u); |
1235 | 0 | T1v = T1r + T1u; |
1236 | 0 | T3b = FNMS(KP250000000, T1v, T1o); |
1237 | 0 | } |
1238 | 0 | { |
1239 | 0 | E TI, T3o, TF, T3n; |
1240 | 0 | TC = cr[WS(rs, 3)]; |
1241 | 0 | { |
1242 | 0 | E TG, TH, TD, TE; |
1243 | 0 | TG = ci[WS(rs, 11)]; |
1244 | 0 | TH = ci[WS(rs, 6)]; |
1245 | 0 | TI = TG + TH; |
1246 | 0 | T3o = TG - TH; |
1247 | 0 | TD = cr[WS(rs, 8)]; |
1248 | 0 | TE = ci[WS(rs, 1)]; |
1249 | 0 | TF = TD + TE; |
1250 | 0 | T3n = TD - TE; |
1251 | 0 | } |
1252 | 0 | T3h = KP559016994 * (TF - TI); |
1253 | 0 | T5S = FMA(KP951056516, T3n, KP587785252 * T3o); |
1254 | 0 | T3p = FNMS(KP951056516, T3o, KP587785252 * T3n); |
1255 | 0 | TJ = TF + TI; |
1256 | 0 | T3g = FNMS(KP250000000, TJ, TC); |
1257 | 0 | } |
1258 | 0 | { |
1259 | 0 | E T1D, T3k, T1A, T3j; |
1260 | 0 | T1x = ci[WS(rs, 21)]; |
1261 | 0 | { |
1262 | 0 | E T1B, T1C, T1y, T1z; |
1263 | 0 | T1B = cr[WS(rs, 13)]; |
1264 | 0 | T1C = cr[WS(rs, 18)]; |
1265 | 0 | T1D = T1B + T1C; |
1266 | 0 | T3k = T1C - T1B; |
1267 | 0 | T1y = ci[WS(rs, 16)]; |
1268 | 0 | T1z = cr[WS(rs, 23)]; |
1269 | 0 | T1A = T1y - T1z; |
1270 | 0 | T3j = T1y + T1z; |
1271 | 0 | } |
1272 | 0 | T3l = FNMS(KP951056516, T3k, KP587785252 * T3j); |
1273 | 0 | T5Q = FMA(KP951056516, T3j, KP587785252 * T3k); |
1274 | 0 | T3r = KP559016994 * (T1A + T1D); |
1275 | 0 | T1E = T1A - T1D; |
1276 | 0 | T3q = FNMS(KP250000000, T1E, T1x); |
1277 | 0 | } |
1278 | 0 | TB = Tt + TA; |
1279 | 0 | TK = TC + TJ; |
1280 | 0 | TL = TB + TK; |
1281 | 0 | T1w = T1o + T1v; |
1282 | 0 | T1F = T1x + T1E; |
1283 | 0 | T1Q = T1w + T1F; |
1284 | 0 | { |
1285 | 0 | E T37, T4D, T3e, T4E, T33, T3d; |
1286 | 0 | T33 = T31 - T32; |
1287 | 0 | T37 = T33 - T36; |
1288 | 0 | T4D = T33 + T36; |
1289 | 0 | T3d = T3b - T3c; |
1290 | 0 | T3e = T3a + T3d; |
1291 | 0 | T4E = T3d - T3a; |
1292 | 0 | T3f = FNMS(KP844327925, T3e, KP535826794 * T37); |
1293 | 0 | T4S = FMA(KP062790519, T4E, KP998026728 * T4D); |
1294 | 0 | T3K = FMA(KP535826794, T3e, KP844327925 * T37); |
1295 | 0 | T4F = FNMS(KP998026728, T4E, KP062790519 * T4D); |
1296 | 0 | } |
1297 | 0 | { |
1298 | 0 | E T5R, T6T, T5U, T6S, T5P, T5T; |
1299 | 0 | T5P = T3h + T3g; |
1300 | 0 | T5R = T5P - T5Q; |
1301 | 0 | T6T = T5P + T5Q; |
1302 | 0 | T5T = T3q + T3r; |
1303 | 0 | T5U = T5S + T5T; |
1304 | 0 | T6S = T5T - T5S; |
1305 | 0 | T5V = FNMS(KP684547105, T5U, KP728968627 * T5R); |
1306 | 0 | T74 = FNMS(KP992114701, T6S, KP125333233 * T6T); |
1307 | 0 | T68 = FMA(KP728968627, T5U, KP684547105 * T5R); |
1308 | 0 | T6U = FMA(KP125333233, T6S, KP992114701 * T6T); |
1309 | 0 | } |
1310 | 0 | { |
1311 | 0 | E T5K, T6Q, T5N, T6P, T5I, T5M; |
1312 | 0 | T5I = T32 + T31; |
1313 | 0 | T5K = T5I - T5J; |
1314 | 0 | T6Q = T5I + T5J; |
1315 | 0 | T5M = T3c + T3b; |
1316 | 0 | T5N = T5L + T5M; |
1317 | 0 | T6P = T5M - T5L; |
1318 | 0 | T5O = FNMS(KP481753674, T5N, KP876306680 * T5K); |
1319 | 0 | T73 = FNMS(KP425779291, T6P, KP904827052 * T6Q); |
1320 | 0 | T67 = FMA(KP876306680, T5N, KP481753674 * T5K); |
1321 | 0 | T6R = FMA(KP904827052, T6P, KP425779291 * T6Q); |
1322 | 0 | } |
1323 | 0 | { |
1324 | 0 | E T3m, T4H, T3t, T4G, T3i, T3s; |
1325 | 0 | T3i = T3g - T3h; |
1326 | 0 | T3m = T3i - T3l; |
1327 | 0 | T4H = T3i + T3l; |
1328 | 0 | T3s = T3q - T3r; |
1329 | 0 | T3t = T3p + T3s; |
1330 | 0 | T4G = T3s - T3p; |
1331 | 0 | T3u = FNMS(KP998026728, T3t, KP062790519 * T3m); |
1332 | 0 | T4T = FNMS(KP637423989, T4G, KP770513242 * T4H); |
1333 | 0 | T3L = FMA(KP062790519, T3t, KP998026728 * T3m); |
1334 | 0 | T4I = FMA(KP770513242, T4G, KP637423989 * T4H); |
1335 | 0 | } |
1336 | 0 | } |
1337 | 0 | { |
1338 | 0 | E TM, T14, T2e, T21, T23, T2l, T1H, T2f, T1O, T2k; |
1339 | 0 | { |
1340 | 0 | E T12, T13, T1R, T22; |
1341 | 0 | T12 = KP559016994 * (Ts - TL); |
1342 | 0 | TM = Ts + TL; |
1343 | 0 | T13 = FNMS(KP250000000, TM, T9); |
1344 | 0 | T14 = T12 + T13; |
1345 | 0 | T2e = T13 - T12; |
1346 | 0 | T1R = KP559016994 * (T1P - T1Q); |
1347 | 0 | T21 = T1P + T1Q; |
1348 | 0 | T22 = FNMS(KP250000000, T21, T20); |
1349 | 0 | T23 = T1R + T22; |
1350 | 0 | T2l = T22 - T1R; |
1351 | 0 | } |
1352 | 0 | { |
1353 | 0 | E T1n, T1G, T1M, T1N; |
1354 | 0 | T1n = T1d - T1m; |
1355 | 0 | T1G = T1w - T1F; |
1356 | 0 | T1H = FMA(KP951056516, T1n, KP587785252 * T1G); |
1357 | 0 | T2f = FNMS(KP951056516, T1G, KP587785252 * T1n); |
1358 | 0 | T1M = Ti - Tr; |
1359 | 0 | T1N = TB - TK; |
1360 | 0 | T1O = FMA(KP951056516, T1M, KP587785252 * T1N); |
1361 | 0 | T2k = FNMS(KP951056516, T1N, KP587785252 * T1M); |
1362 | 0 | } |
1363 | 0 | { |
1364 | 0 | E T1I, T24, T2o, T2q; |
1365 | 0 | cr[0] = T9 + TM; |
1366 | 0 | ci[0] = T20 + T21; |
1367 | 0 | T1I = T14 - T1H; |
1368 | 0 | T24 = T1O + T23; |
1369 | 0 | cr[WS(rs, 5)] = FNMS(T1L, T24, T11 * T1I); |
1370 | 0 | ci[WS(rs, 5)] = FMA(T1L, T1I, T11 * T24); |
1371 | 0 | T2o = T2e + T2f; |
1372 | 0 | T2q = T2l - T2k; |
1373 | 0 | cr[WS(rs, 15)] = FNMS(T2p, T2q, T2n * T2o); |
1374 | 0 | ci[WS(rs, 15)] = FMA(T2p, T2o, T2n * T2q); |
1375 | 0 | { |
1376 | 0 | E T2g, T2m, T28, T2a; |
1377 | 0 | T2g = T2e - T2f; |
1378 | 0 | T2m = T2k + T2l; |
1379 | 0 | cr[WS(rs, 10)] = FNMS(T2j, T2m, T2d * T2g); |
1380 | 0 | ci[WS(rs, 10)] = FMA(T2j, T2g, T2d * T2m); |
1381 | 0 | T28 = T14 + T1H; |
1382 | 0 | T2a = T23 - T1O; |
1383 | 0 | cr[WS(rs, 20)] = FNMS(T29, T2a, T27 * T28); |
1384 | 0 | ci[WS(rs, 20)] = FMA(T29, T28, T27 * T2a); |
1385 | 0 | } |
1386 | 0 | } |
1387 | 0 | } |
1388 | 0 | { |
1389 | 0 | E T76, T7n, T7a, T7q, T6H, T6W, T6X, T6Y, T7e, T7f, T7d, T7g, T7x, T7y; |
1390 | 0 | { |
1391 | 0 | E T72, T75, T78, T79; |
1392 | 0 | T72 = T70 + T71; |
1393 | 0 | T75 = T73 - T74; |
1394 | 0 | T76 = FMA(KP951056516, T72, KP587785252 * T75); |
1395 | 0 | T7n = FNMS(KP951056516, T75, KP587785252 * T72); |
1396 | 0 | T78 = T6K - T6N; |
1397 | 0 | T79 = T6U - T6R; |
1398 | 0 | T7a = FMA(KP951056516, T78, KP587785252 * T79); |
1399 | 0 | T7q = FNMS(KP951056516, T79, KP587785252 * T78); |
1400 | 0 | } |
1401 | 0 | { |
1402 | 0 | E T6O, T6V, T7b, T7c; |
1403 | 0 | T6H = T5Y + T5Z; |
1404 | 0 | T6O = T6K + T6N; |
1405 | 0 | T6V = T6R + T6U; |
1406 | 0 | T6W = T6O - T6V; |
1407 | 0 | T6X = FNMS(KP250000000, T6W, T6H); |
1408 | 0 | T6Y = KP559016994 * (T6O + T6V); |
1409 | 0 | T7e = T6j - T6i; |
1410 | 0 | T7b = T70 - T71; |
1411 | 0 | T7c = T73 + T74; |
1412 | 0 | T7f = T7b + T7c; |
1413 | 0 | T7d = KP559016994 * (T7b - T7c); |
1414 | 0 | T7g = FNMS(KP250000000, T7f, T7e); |
1415 | 0 | } |
1416 | 0 | T7x = T6H + T6W; |
1417 | 0 | T7y = T7e + T7f; |
1418 | 0 | cr[WS(rs, 4)] = FNMS(TY, T7y, TT * T7x); |
1419 | 0 | ci[WS(rs, 4)] = FMA(TY, T7x, TT * T7y); |
1420 | 0 | { |
1421 | 0 | E T7o, T7u, T7s, T7w, T7m, T7r; |
1422 | 0 | T7m = T6X - T6Y; |
1423 | 0 | T7o = T7m - T7n; |
1424 | 0 | T7u = T7m + T7n; |
1425 | 0 | T7r = T7g - T7d; |
1426 | 0 | T7s = T7q + T7r; |
1427 | 0 | T7w = T7r - T7q; |
1428 | 0 | cr[WS(rs, 14)] = FNMS(T7p, T7s, T7l * T7o); |
1429 | 0 | ci[WS(rs, 14)] = FMA(T7p, T7o, T7l * T7s); |
1430 | 0 | cr[WS(rs, 19)] = FNMS(T7v, T7w, T7t * T7u); |
1431 | 0 | ci[WS(rs, 19)] = FMA(T7v, T7u, T7t * T7w); |
1432 | 0 | } |
1433 | 0 | { |
1434 | 0 | E T77, T7j, T7i, T7k, T6Z, T7h; |
1435 | 0 | T6Z = T6X + T6Y; |
1436 | 0 | T77 = T6Z - T76; |
1437 | 0 | T7j = T6Z + T76; |
1438 | 0 | T7h = T7d + T7g; |
1439 | 0 | T7i = T7a + T7h; |
1440 | 0 | T7k = T7h - T7a; |
1441 | 0 | cr[WS(rs, 9)] = FNMS(TZ, T7i, TU * T77); |
1442 | 0 | ci[WS(rs, 9)] = FMA(TZ, T77, TU * T7i); |
1443 | 0 | cr[WS(rs, 24)] = FNMS(T26, T7k, T25 * T7j); |
1444 | 0 | ci[WS(rs, 24)] = FMA(T26, T7j, T25 * T7k); |
1445 | 0 | } |
1446 | 0 | } |
1447 | 0 | { |
1448 | 0 | E T3N, T4h, T3U, T4m, T3D, T3E, T3w, T3F, T44, T45, T3X, T46, T4t, T4u; |
1449 | 0 | { |
1450 | 0 | E T3J, T3M, T3S, T3T; |
1451 | 0 | T3J = T3H - T3I; |
1452 | 0 | T3M = T3K - T3L; |
1453 | 0 | T3N = FMA(KP951056516, T3J, KP587785252 * T3M); |
1454 | 0 | T4h = FNMS(KP951056516, T3M, KP587785252 * T3J); |
1455 | 0 | T3S = T2K + T2Z; |
1456 | 0 | T3T = T3f - T3u; |
1457 | 0 | T3U = FMA(KP951056516, T3S, KP587785252 * T3T); |
1458 | 0 | T4m = FNMS(KP951056516, T3T, KP587785252 * T3S); |
1459 | 0 | } |
1460 | 0 | { |
1461 | 0 | E T30, T3v, T3V, T3W; |
1462 | 0 | T3D = T3z - T3C; |
1463 | 0 | T30 = T2K - T2Z; |
1464 | 0 | T3v = T3f + T3u; |
1465 | 0 | T3E = T30 + T3v; |
1466 | 0 | T3w = KP559016994 * (T30 - T3v); |
1467 | 0 | T3F = FNMS(KP250000000, T3E, T3D); |
1468 | 0 | T44 = T40 + T43; |
1469 | 0 | T3V = T3H + T3I; |
1470 | 0 | T3W = T3K + T3L; |
1471 | 0 | T45 = T3V + T3W; |
1472 | 0 | T3X = KP559016994 * (T3V - T3W); |
1473 | 0 | T46 = FNMS(KP250000000, T45, T44); |
1474 | 0 | } |
1475 | 0 | T4t = T3D + T3E; |
1476 | 0 | T4u = T44 + T45; |
1477 | 0 | cr[WS(rs, 2)] = FNMS(T2t, T4u, T2r * T4t); |
1478 | 0 | ci[WS(rs, 2)] = FMA(T2t, T4t, T2r * T4u); |
1479 | 0 | { |
1480 | 0 | E T4i, T4q, T4o, T4s, T4g, T4n; |
1481 | 0 | T4g = T3F - T3w; |
1482 | 0 | T4i = T4g - T4h; |
1483 | 0 | T4q = T4g + T4h; |
1484 | 0 | T4n = T46 - T3X; |
1485 | 0 | T4o = T4m + T4n; |
1486 | 0 | T4s = T4n - T4m; |
1487 | 0 | cr[WS(rs, 12)] = FNMS(T4l, T4o, T4f * T4i); |
1488 | 0 | ci[WS(rs, 12)] = FMA(T4l, T4i, T4f * T4o); |
1489 | 0 | cr[WS(rs, 17)] = FNMS(T4r, T4s, T4p * T4q); |
1490 | 0 | ci[WS(rs, 17)] = FMA(T4r, T4q, T4p * T4s); |
1491 | 0 | } |
1492 | 0 | { |
1493 | 0 | E T3O, T4a, T48, T4c, T3G, T47; |
1494 | 0 | T3G = T3w + T3F; |
1495 | 0 | T3O = T3G - T3N; |
1496 | 0 | T4a = T3G + T3N; |
1497 | 0 | T47 = T3X + T46; |
1498 | 0 | T48 = T3U + T47; |
1499 | 0 | T4c = T47 - T3U; |
1500 | 0 | cr[WS(rs, 7)] = FNMS(T3R, T48, T2v * T3O); |
1501 | 0 | ci[WS(rs, 7)] = FMA(T3R, T3O, T2v * T48); |
1502 | 0 | cr[WS(rs, 22)] = FNMS(T4b, T4c, T49 * T4a); |
1503 | 0 | ci[WS(rs, 22)] = FMA(T4b, T4a, T49 * T4c); |
1504 | 0 | } |
1505 | 0 | } |
1506 | 0 | { |
1507 | 0 | E T4V, T5f, T50, T5i, T4L, T4M, T4K, T4N, T54, T55, T53, T56, T5r, T5s; |
1508 | 0 | { |
1509 | 0 | E T4R, T4U, T4Y, T4Z; |
1510 | 0 | T4R = T4P - T4Q; |
1511 | 0 | T4U = T4S - T4T; |
1512 | 0 | T4V = FMA(KP951056516, T4R, KP587785252 * T4U); |
1513 | 0 | T5f = FNMS(KP951056516, T4U, KP587785252 * T4R); |
1514 | 0 | T4Y = T4y + T4B; |
1515 | 0 | T4Z = T4F + T4I; |
1516 | 0 | T50 = FMA(KP951056516, T4Y, KP587785252 * T4Z); |
1517 | 0 | T5i = FNMS(KP951056516, T4Z, KP587785252 * T4Y); |
1518 | 0 | } |
1519 | 0 | { |
1520 | 0 | E T4C, T4J, T51, T52; |
1521 | 0 | T4L = T3z + T3C; |
1522 | 0 | T4C = T4y - T4B; |
1523 | 0 | T4J = T4F - T4I; |
1524 | 0 | T4M = T4C + T4J; |
1525 | 0 | T4K = KP559016994 * (T4C - T4J); |
1526 | 0 | T4N = FNMS(KP250000000, T4M, T4L); |
1527 | 0 | T54 = T43 - T40; |
1528 | 0 | T51 = T4P + T4Q; |
1529 | 0 | T52 = T4S + T4T; |
1530 | 0 | T55 = T51 + T52; |
1531 | 0 | T53 = KP559016994 * (T51 - T52); |
1532 | 0 | T56 = FNMS(KP250000000, T55, T54); |
1533 | 0 | } |
1534 | 0 | T5r = T4L + T4M; |
1535 | 0 | T5s = T54 + T55; |
1536 | 0 | cr[WS(rs, 3)] = FNMS(TR, T5s, TO * T5r); |
1537 | 0 | ci[WS(rs, 3)] = FMA(TR, T5r, TO * T5s); |
1538 | 0 | { |
1539 | 0 | E T5g, T5o, T5k, T5q, T5e, T5j; |
1540 | 0 | T5e = T4N - T4K; |
1541 | 0 | T5g = T5e - T5f; |
1542 | 0 | T5o = T5e + T5f; |
1543 | 0 | T5j = T56 - T53; |
1544 | 0 | T5k = T5i + T5j; |
1545 | 0 | T5q = T5j - T5i; |
1546 | 0 | cr[WS(rs, 13)] = FNMS(T5h, T5k, T5d * T5g); |
1547 | 0 | ci[WS(rs, 13)] = FMA(T5h, T5g, T5d * T5k); |
1548 | 0 | cr[WS(rs, 18)] = FNMS(T5p, T5q, T5n * T5o); |
1549 | 0 | ci[WS(rs, 18)] = FMA(T5p, T5o, T5n * T5q); |
1550 | 0 | } |
1551 | 0 | { |
1552 | 0 | E T4W, T5a, T58, T5c, T4O, T57; |
1553 | 0 | T4O = T4K + T4N; |
1554 | 0 | T4W = T4O - T4V; |
1555 | 0 | T5a = T4O + T4V; |
1556 | 0 | T57 = T53 + T56; |
1557 | 0 | T58 = T50 + T57; |
1558 | 0 | T5c = T57 - T50; |
1559 | 0 | cr[WS(rs, 8)] = FNMS(T4X, T58, T4v * T4W); |
1560 | 0 | ci[WS(rs, 8)] = FMA(T4X, T4W, T4v * T58); |
1561 | 0 | cr[WS(rs, 23)] = FNMS(T5b, T5c, T59 * T5a); |
1562 | 0 | ci[WS(rs, 23)] = FMA(T5b, T5a, T59 * T5c); |
1563 | 0 | } |
1564 | 0 | } |
1565 | 0 | { |
1566 | 0 | E T6a, T6v, T6e, T6y, T60, T61, T5X, T62, T6k, T6l, T6h, T6m, T6F, T6G; |
1567 | 0 | { |
1568 | 0 | E T66, T69, T6c, T6d; |
1569 | 0 | T66 = T64 - T65; |
1570 | 0 | T69 = T67 - T68; |
1571 | 0 | T6a = FMA(KP951056516, T66, KP587785252 * T69); |
1572 | 0 | T6v = FNMS(KP951056516, T69, KP587785252 * T66); |
1573 | 0 | T6c = T5z - T5G; |
1574 | 0 | T6d = T5O - T5V; |
1575 | 0 | T6e = FMA(KP951056516, T6c, KP587785252 * T6d); |
1576 | 0 | T6y = FNMS(KP951056516, T6d, KP587785252 * T6c); |
1577 | 0 | } |
1578 | 0 | { |
1579 | 0 | E T5H, T5W, T6f, T6g; |
1580 | 0 | T60 = T5Y - T5Z; |
1581 | 0 | T5H = T5z + T5G; |
1582 | 0 | T5W = T5O + T5V; |
1583 | 0 | T61 = T5H + T5W; |
1584 | 0 | T5X = KP559016994 * (T5H - T5W); |
1585 | 0 | T62 = FNMS(KP250000000, T61, T60); |
1586 | 0 | T6k = T6i + T6j; |
1587 | 0 | T6f = T64 + T65; |
1588 | 0 | T6g = T67 + T68; |
1589 | 0 | T6l = T6f + T6g; |
1590 | 0 | T6h = KP559016994 * (T6f - T6g); |
1591 | 0 | T6m = FNMS(KP250000000, T6l, T6k); |
1592 | 0 | } |
1593 | 0 | T6F = T60 + T61; |
1594 | 0 | T6G = T6k + T6l; |
1595 | 0 | cr[WS(rs, 1)] = FNMS(TQ, T6G, TN * T6F); |
1596 | 0 | ci[WS(rs, 1)] = FMA(TQ, T6F, TN * T6G); |
1597 | 0 | { |
1598 | 0 | E T6w, T6C, T6A, T6E, T6u, T6z; |
1599 | 0 | T6u = T62 - T5X; |
1600 | 0 | T6w = T6u - T6v; |
1601 | 0 | T6C = T6u + T6v; |
1602 | 0 | T6z = T6m - T6h; |
1603 | 0 | T6A = T6y + T6z; |
1604 | 0 | T6E = T6z - T6y; |
1605 | 0 | cr[WS(rs, 11)] = FNMS(T6x, T6A, T6t * T6w); |
1606 | 0 | ci[WS(rs, 11)] = FMA(T6x, T6w, T6t * T6A); |
1607 | 0 | cr[WS(rs, 16)] = FNMS(T6D, T6E, T6B * T6C); |
1608 | 0 | ci[WS(rs, 16)] = FMA(T6D, T6C, T6B * T6E); |
1609 | 0 | } |
1610 | 0 | { |
1611 | 0 | E T6b, T6q, T6o, T6s, T63, T6n; |
1612 | 0 | T63 = T5X + T62; |
1613 | 0 | T6b = T63 - T6a; |
1614 | 0 | T6q = T63 + T6a; |
1615 | 0 | T6n = T6h + T6m; |
1616 | 0 | T6o = T6e + T6n; |
1617 | 0 | T6s = T6n - T6e; |
1618 | 0 | cr[WS(rs, 6)] = FNMS(T5m, T6o, T5l * T6b); |
1619 | 0 | ci[WS(rs, 6)] = FMA(T5m, T6b, T5l * T6o); |
1620 | 0 | cr[WS(rs, 21)] = FNMS(T6r, T6s, T6p * T6q); |
1621 | 0 | ci[WS(rs, 21)] = FMA(T6r, T6q, T6p * T6s); |
1622 | 0 | } |
1623 | 0 | } |
1624 | 0 | } |
1625 | 0 | } |
1626 | 0 | } |
1627 | 0 | } |
1628 | | |
1629 | | static const tw_instr twinstr[] = { |
1630 | | { TW_CEXP, 1, 1 }, |
1631 | | { TW_CEXP, 1, 3 }, |
1632 | | { TW_CEXP, 1, 9 }, |
1633 | | { TW_CEXP, 1, 24 }, |
1634 | | { TW_NEXT, 1, 0 } |
1635 | | }; |
1636 | | |
1637 | | static const hc2hc_desc desc = { 25, "hb2_25", twinstr, &GENUS, { 280, 180, 160, 0 } }; |
1638 | | |
1639 | 1 | void X(codelet_hb2_25) (planner *p) { |
1640 | 1 | X(khc2hc_register) (p, hb2_25, &desc); |
1641 | 1 | } |
1642 | | #endif |