/src/fftw3/rdft/scalar/r2cf/r2cf_64.c
Line | Count | Source |
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 Sun Dec 14 06:36:16 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_r2cf.native -fma -compact -variables 4 -pipeline-latency 4 -n 64 -name r2cf_64 -include rdft/scalar/r2cf.h */ |
29 | | |
30 | | /* |
31 | | * This function contains 394 FP additions, 196 FP multiplications, |
32 | | * (or, 198 additions, 0 multiplications, 196 fused multiply/add), |
33 | | * 106 stack variables, 15 constants, and 128 memory accesses |
34 | | */ |
35 | | #include "rdft/scalar/r2cf.h" |
36 | | |
37 | | static void r2cf_64(R *R0, R *R1, R *Cr, R *Ci, stride rs, stride csr, stride csi, INT v, INT ivs, INT ovs) |
38 | | { |
39 | | DK(KP773010453, +0.773010453362736960810906609758469800971041293); |
40 | | DK(KP995184726, +0.995184726672196886244836953109479921575474869); |
41 | | DK(KP098491403, +0.098491403357164253077197521291327432293052451); |
42 | | DK(KP820678790, +0.820678790828660330972281985331011598767386482); |
43 | | DK(KP956940335, +0.956940335732208864935797886980269969482849206); |
44 | | DK(KP881921264, +0.881921264348355029712756863660388349508442621); |
45 | | DK(KP534511135, +0.534511135950791641089685961295362908582039528); |
46 | | DK(KP303346683, +0.303346683607342391675883946941299872384187453); |
47 | | DK(KP980785280, +0.980785280403230449126182236134239036973933731); |
48 | | DK(KP198912367, +0.198912367379658006911597622644676228597850501); |
49 | | DK(KP831469612, +0.831469612302545237078788377617905756738560812); |
50 | | DK(KP668178637, +0.668178637919298919997757686523080761552472251); |
51 | | DK(KP923879532, +0.923879532511286756128183189396788286822416626); |
52 | | DK(KP414213562, +0.414213562373095048801688724209698078569671875); |
53 | | DK(KP707106781, +0.707106781186547524400844362104849039284835938); |
54 | | { |
55 | | INT i; |
56 | | for (i = v; i > 0; i = i - 1, R0 = R0 + ivs, R1 = R1 + ivs, Cr = Cr + ovs, Ci = Ci + ovs, MAKE_VOLATILE_STRIDE(256, rs), MAKE_VOLATILE_STRIDE(256, csr), MAKE_VOLATILE_STRIDE(256, csi)) { |
57 | | E T11, T2j, T3D, T5p, T4P, T5P, T7, Te, Tf, T1k, T1H, T4a, T5A, T4l, T5D; |
58 | | E T2U, T3i, T1R, T2e, T4v, T5H, T4G, T5K, T31, T3l, T3Z, T5t, T42, T5s, TZ; |
59 | | E T3f, T1b, T2n, T3Q, T5w, T3T, T5v, TK, T3e, T18, T2m, Tm, Tt, Tu, T4S; |
60 | | E T5q, T14, T2k, T3K, T5Q, T1z, T1I, T4o, T5B, T2X, T3j, T4h, T5E, T26, T2f; |
61 | | E T4J, T5I, T34, T3m, T4C, T5L; |
62 | | { |
63 | | E T3, T3z, Td, T3B, T6, T4N, Ta, T3A, T3C, T4O; |
64 | | { |
65 | | E T1, T2, Tb, Tc; |
66 | | T1 = R0[0]; |
67 | | T2 = R0[WS(rs, 16)]; |
68 | | T3 = T1 + T2; |
69 | | T3z = T1 - T2; |
70 | | Tb = R0[WS(rs, 28)]; |
71 | | Tc = R0[WS(rs, 12)]; |
72 | | Td = Tb + Tc; |
73 | | T3B = Tb - Tc; |
74 | | } |
75 | | { |
76 | | E T4, T5, T8, T9; |
77 | | T4 = R0[WS(rs, 8)]; |
78 | | T5 = R0[WS(rs, 24)]; |
79 | | T6 = T4 + T5; |
80 | | T4N = T4 - T5; |
81 | | T8 = R0[WS(rs, 4)]; |
82 | | T9 = R0[WS(rs, 20)]; |
83 | | Ta = T8 + T9; |
84 | | T3A = T8 - T9; |
85 | | } |
86 | | T11 = T3 - T6; |
87 | | T2j = Td - Ta; |
88 | | T3C = T3A + T3B; |
89 | | T3D = FMA(KP707106781, T3C, T3z); |
90 | | T5p = FNMS(KP707106781, T3C, T3z); |
91 | | T4O = T3B - T3A; |
92 | | T4P = FNMS(KP707106781, T4O, T4N); |
93 | | T5P = FMA(KP707106781, T4O, T4N); |
94 | | T7 = T3 + T6; |
95 | | Te = Ta + Td; |
96 | | Tf = T7 + Te; |
97 | | } |
98 | | { |
99 | | E T1g, T46, T1G, T47, T1j, T4j, T1D, T48; |
100 | | { |
101 | | E T1e, T1f, T1E, T1F; |
102 | | T1e = R1[0]; |
103 | | T1f = R1[WS(rs, 16)]; |
104 | | T1g = T1e + T1f; |
105 | | T46 = T1e - T1f; |
106 | | T1E = R1[WS(rs, 4)]; |
107 | | T1F = R1[WS(rs, 20)]; |
108 | | T1G = T1E + T1F; |
109 | | T47 = T1E - T1F; |
110 | | } |
111 | | { |
112 | | E T1h, T1i, T1B, T1C; |
113 | | T1h = R1[WS(rs, 8)]; |
114 | | T1i = R1[WS(rs, 24)]; |
115 | | T1j = T1h + T1i; |
116 | | T4j = T1h - T1i; |
117 | | T1B = R1[WS(rs, 28)]; |
118 | | T1C = R1[WS(rs, 12)]; |
119 | | T1D = T1B + T1C; |
120 | | T48 = T1B - T1C; |
121 | | } |
122 | | T1k = T1g - T1j; |
123 | | T1H = T1D - T1G; |
124 | | { |
125 | | E T49, T4k, T2S, T2T; |
126 | | T49 = T47 + T48; |
127 | | T4a = FMA(KP707106781, T49, T46); |
128 | | T5A = FNMS(KP707106781, T49, T46); |
129 | | T4k = T47 - T48; |
130 | | T4l = FMA(KP707106781, T4k, T4j); |
131 | | T5D = FNMS(KP707106781, T4k, T4j); |
132 | | T2S = T1g + T1j; |
133 | | T2T = T1G + T1D; |
134 | | T2U = T2S + T2T; |
135 | | T3i = T2S - T2T; |
136 | | } |
137 | | } |
138 | | { |
139 | | E T1N, T4r, T2d, T4s, T1Q, T4E, T2a, T4t; |
140 | | { |
141 | | E T1L, T1M, T2b, T2c; |
142 | | T1L = R1[WS(rs, 31)]; |
143 | | T1M = R1[WS(rs, 15)]; |
144 | | T1N = T1L + T1M; |
145 | | T4r = T1L - T1M; |
146 | | T2b = R1[WS(rs, 3)]; |
147 | | T2c = R1[WS(rs, 19)]; |
148 | | T2d = T2b + T2c; |
149 | | T4s = T2b - T2c; |
150 | | } |
151 | | { |
152 | | E T1O, T1P, T28, T29; |
153 | | T1O = R1[WS(rs, 7)]; |
154 | | T1P = R1[WS(rs, 23)]; |
155 | | T1Q = T1O + T1P; |
156 | | T4E = T1P - T1O; |
157 | | T28 = R1[WS(rs, 27)]; |
158 | | T29 = R1[WS(rs, 11)]; |
159 | | T2a = T28 + T29; |
160 | | T4t = T28 - T29; |
161 | | } |
162 | | T1R = T1N - T1Q; |
163 | | T2e = T2a - T2d; |
164 | | { |
165 | | E T4u, T4F, T2Z, T30; |
166 | | T4u = T4s + T4t; |
167 | | T4v = FMA(KP707106781, T4u, T4r); |
168 | | T5H = FNMS(KP707106781, T4u, T4r); |
169 | | T4F = T4t - T4s; |
170 | | T4G = FMA(KP707106781, T4F, T4E); |
171 | | T5K = FNMS(KP707106781, T4F, T4E); |
172 | | T2Z = T1N + T1Q; |
173 | | T30 = T2d + T2a; |
174 | | T31 = T2Z + T30; |
175 | | T3l = T2Z - T30; |
176 | | } |
177 | | } |
178 | | { |
179 | | E TN, T3V, TX, T3X, TQ, T40, TU, T3W, T3Y, T41; |
180 | | { |
181 | | E TL, TM, TV, TW; |
182 | | TL = R0[WS(rs, 31)]; |
183 | | TM = R0[WS(rs, 15)]; |
184 | | TN = TL + TM; |
185 | | T3V = TL - TM; |
186 | | TV = R0[WS(rs, 27)]; |
187 | | TW = R0[WS(rs, 11)]; |
188 | | TX = TV + TW; |
189 | | T3X = TV - TW; |
190 | | } |
191 | | { |
192 | | E TO, TP, TS, TT; |
193 | | TO = R0[WS(rs, 7)]; |
194 | | TP = R0[WS(rs, 23)]; |
195 | | TQ = TO + TP; |
196 | | T40 = TO - TP; |
197 | | TS = R0[WS(rs, 3)]; |
198 | | TT = R0[WS(rs, 19)]; |
199 | | TU = TS + TT; |
200 | | T3W = TS - TT; |
201 | | } |
202 | | T3Y = T3W + T3X; |
203 | | T3Z = FMA(KP707106781, T3Y, T3V); |
204 | | T5t = FNMS(KP707106781, T3Y, T3V); |
205 | | T41 = T3W - T3X; |
206 | | T42 = FMA(KP707106781, T41, T40); |
207 | | T5s = FNMS(KP707106781, T41, T40); |
208 | | { |
209 | | E TR, TY, T19, T1a; |
210 | | TR = TN + TQ; |
211 | | TY = TU + TX; |
212 | | TZ = TR + TY; |
213 | | T3f = TR - TY; |
214 | | T19 = TN - TQ; |
215 | | T1a = TX - TU; |
216 | | T1b = FNMS(KP414213562, T1a, T19); |
217 | | T2n = FMA(KP414213562, T19, T1a); |
218 | | } |
219 | | } |
220 | | { |
221 | | E Ty, T3M, TI, T3O, TB, T3R, TF, T3N, T3P, T3S; |
222 | | { |
223 | | E Tw, Tx, TG, TH; |
224 | | Tw = R0[WS(rs, 1)]; |
225 | | Tx = R0[WS(rs, 17)]; |
226 | | Ty = Tw + Tx; |
227 | | T3M = Tw - Tx; |
228 | | TG = R0[WS(rs, 29)]; |
229 | | TH = R0[WS(rs, 13)]; |
230 | | TI = TG + TH; |
231 | | T3O = TG - TH; |
232 | | } |
233 | | { |
234 | | E Tz, TA, TD, TE; |
235 | | Tz = R0[WS(rs, 9)]; |
236 | | TA = R0[WS(rs, 25)]; |
237 | | TB = Tz + TA; |
238 | | T3R = Tz - TA; |
239 | | TD = R0[WS(rs, 5)]; |
240 | | TE = R0[WS(rs, 21)]; |
241 | | TF = TD + TE; |
242 | | T3N = TD - TE; |
243 | | } |
244 | | T3P = T3N + T3O; |
245 | | T3Q = FMA(KP707106781, T3P, T3M); |
246 | | T5w = FNMS(KP707106781, T3P, T3M); |
247 | | T3S = T3N - T3O; |
248 | | T3T = FMA(KP707106781, T3S, T3R); |
249 | | T5v = FNMS(KP707106781, T3S, T3R); |
250 | | { |
251 | | E TC, TJ, T16, T17; |
252 | | TC = Ty + TB; |
253 | | TJ = TF + TI; |
254 | | TK = TC + TJ; |
255 | | T3e = TC - TJ; |
256 | | T16 = Ty - TB; |
257 | | T17 = TI - TF; |
258 | | T18 = FMA(KP414213562, T17, T16); |
259 | | T2m = FNMS(KP414213562, T16, T17); |
260 | | } |
261 | | } |
262 | | { |
263 | | E Ti, T3E, Ts, T3I, Tl, T3F, Tp, T3H, T4Q, T4R; |
264 | | { |
265 | | E Tg, Th, Tq, Tr; |
266 | | Tg = R0[WS(rs, 2)]; |
267 | | Th = R0[WS(rs, 18)]; |
268 | | Ti = Tg + Th; |
269 | | T3E = Tg - Th; |
270 | | Tq = R0[WS(rs, 6)]; |
271 | | Tr = R0[WS(rs, 22)]; |
272 | | Ts = Tq + Tr; |
273 | | T3I = Tq - Tr; |
274 | | } |
275 | | { |
276 | | E Tj, Tk, Tn, To; |
277 | | Tj = R0[WS(rs, 10)]; |
278 | | Tk = R0[WS(rs, 26)]; |
279 | | Tl = Tj + Tk; |
280 | | T3F = Tj - Tk; |
281 | | Tn = R0[WS(rs, 30)]; |
282 | | To = R0[WS(rs, 14)]; |
283 | | Tp = Tn + To; |
284 | | T3H = Tn - To; |
285 | | } |
286 | | Tm = Ti + Tl; |
287 | | Tt = Tp + Ts; |
288 | | Tu = Tm + Tt; |
289 | | T4Q = FMA(KP414213562, T3E, T3F); |
290 | | T4R = FNMS(KP414213562, T3H, T3I); |
291 | | T4S = T4Q + T4R; |
292 | | T5q = T4Q - T4R; |
293 | | { |
294 | | E T12, T13, T3G, T3J; |
295 | | T12 = Ti - Tl; |
296 | | T13 = Tp - Ts; |
297 | | T14 = T12 + T13; |
298 | | T2k = T13 - T12; |
299 | | T3G = FNMS(KP414213562, T3F, T3E); |
300 | | T3J = FMA(KP414213562, T3I, T3H); |
301 | | T3K = T3G + T3J; |
302 | | T5Q = T3J - T3G; |
303 | | } |
304 | | } |
305 | | { |
306 | | E T1n, T4b, T1x, T4f, T1q, T4c, T1u, T4e; |
307 | | { |
308 | | E T1l, T1m, T1v, T1w; |
309 | | T1l = R1[WS(rs, 2)]; |
310 | | T1m = R1[WS(rs, 18)]; |
311 | | T1n = T1l + T1m; |
312 | | T4b = T1l - T1m; |
313 | | T1v = R1[WS(rs, 6)]; |
314 | | T1w = R1[WS(rs, 22)]; |
315 | | T1x = T1v + T1w; |
316 | | T4f = T1v - T1w; |
317 | | } |
318 | | { |
319 | | E T1o, T1p, T1s, T1t; |
320 | | T1o = R1[WS(rs, 10)]; |
321 | | T1p = R1[WS(rs, 26)]; |
322 | | T1q = T1o + T1p; |
323 | | T4c = T1o - T1p; |
324 | | T1s = R1[WS(rs, 30)]; |
325 | | T1t = R1[WS(rs, 14)]; |
326 | | T1u = T1s + T1t; |
327 | | T4e = T1s - T1t; |
328 | | } |
329 | | { |
330 | | E T1r, T1y, T4m, T4n; |
331 | | T1r = T1n - T1q; |
332 | | T1y = T1u - T1x; |
333 | | T1z = T1r + T1y; |
334 | | T1I = T1y - T1r; |
335 | | T4m = FMA(KP414213562, T4b, T4c); |
336 | | T4n = FNMS(KP414213562, T4e, T4f); |
337 | | T4o = T4m + T4n; |
338 | | T5B = T4m - T4n; |
339 | | } |
340 | | { |
341 | | E T2V, T2W, T4d, T4g; |
342 | | T2V = T1n + T1q; |
343 | | T2W = T1u + T1x; |
344 | | T2X = T2V + T2W; |
345 | | T3j = T2W - T2V; |
346 | | T4d = FNMS(KP414213562, T4c, T4b); |
347 | | T4g = FMA(KP414213562, T4f, T4e); |
348 | | T4h = T4d + T4g; |
349 | | T5E = T4g - T4d; |
350 | | } |
351 | | } |
352 | | { |
353 | | E T1U, T4w, T24, T4A, T1X, T4x, T21, T4z; |
354 | | { |
355 | | E T1S, T1T, T22, T23; |
356 | | T1S = R1[WS(rs, 1)]; |
357 | | T1T = R1[WS(rs, 17)]; |
358 | | T1U = T1S + T1T; |
359 | | T4w = T1S - T1T; |
360 | | T22 = R1[WS(rs, 5)]; |
361 | | T23 = R1[WS(rs, 21)]; |
362 | | T24 = T22 + T23; |
363 | | T4A = T23 - T22; |
364 | | } |
365 | | { |
366 | | E T1V, T1W, T1Z, T20; |
367 | | T1V = R1[WS(rs, 9)]; |
368 | | T1W = R1[WS(rs, 25)]; |
369 | | T1X = T1V + T1W; |
370 | | T4x = T1W - T1V; |
371 | | T1Z = R1[WS(rs, 29)]; |
372 | | T20 = R1[WS(rs, 13)]; |
373 | | T21 = T1Z + T20; |
374 | | T4z = T1Z - T20; |
375 | | } |
376 | | { |
377 | | E T1Y, T25, T4H, T4I; |
378 | | T1Y = T1U - T1X; |
379 | | T25 = T21 - T24; |
380 | | T26 = T1Y + T25; |
381 | | T2f = T25 - T1Y; |
382 | | T4H = FNMS(KP414213562, T4w, T4x); |
383 | | T4I = FMA(KP414213562, T4z, T4A); |
384 | | T4J = T4H + T4I; |
385 | | T5I = T4I - T4H; |
386 | | } |
387 | | { |
388 | | E T32, T33, T4y, T4B; |
389 | | T32 = T1U + T1X; |
390 | | T33 = T21 + T24; |
391 | | T34 = T32 + T33; |
392 | | T3m = T33 - T32; |
393 | | T4y = FMA(KP414213562, T4x, T4w); |
394 | | T4B = FNMS(KP414213562, T4A, T4z); |
395 | | T4C = T4y + T4B; |
396 | | T5L = T4B - T4y; |
397 | | } |
398 | | } |
399 | | { |
400 | | E Tv, T10, T39, T3a, T3b, T3c; |
401 | | Tv = Tf + Tu; |
402 | | T10 = TK + TZ; |
403 | | T39 = Tv + T10; |
404 | | T3a = T2U + T2X; |
405 | | T3b = T31 + T34; |
406 | | T3c = T3a + T3b; |
407 | | Cr[WS(csr, 16)] = Tv - T10; |
408 | | Ci[WS(csi, 16)] = T3b - T3a; |
409 | | Cr[WS(csr, 32)] = T39 - T3c; |
410 | | Cr[0] = T39 + T3c; |
411 | | } |
412 | | { |
413 | | E T2R, T37, T36, T38, T2Y, T35; |
414 | | T2R = Tf - Tu; |
415 | | T37 = TZ - TK; |
416 | | T2Y = T2U - T2X; |
417 | | T35 = T31 - T34; |
418 | | T36 = T2Y + T35; |
419 | | T38 = T35 - T2Y; |
420 | | Cr[WS(csr, 24)] = FNMS(KP707106781, T36, T2R); |
421 | | Ci[WS(csi, 24)] = FMS(KP707106781, T38, T37); |
422 | | Cr[WS(csr, 8)] = FMA(KP707106781, T36, T2R); |
423 | | Ci[WS(csi, 8)] = FMA(KP707106781, T38, T37); |
424 | | } |
425 | | { |
426 | | E T3h, T3x, T3w, T3y, T3o, T3s, T3r, T3t; |
427 | | { |
428 | | E T3d, T3g, T3u, T3v; |
429 | | T3d = T7 - Te; |
430 | | T3g = T3e + T3f; |
431 | | T3h = FMA(KP707106781, T3g, T3d); |
432 | | T3x = FNMS(KP707106781, T3g, T3d); |
433 | | T3u = FNMS(KP414213562, T3i, T3j); |
434 | | T3v = FMA(KP414213562, T3l, T3m); |
435 | | T3w = T3u + T3v; |
436 | | T3y = T3v - T3u; |
437 | | } |
438 | | { |
439 | | E T3k, T3n, T3p, T3q; |
440 | | T3k = FMA(KP414213562, T3j, T3i); |
441 | | T3n = FNMS(KP414213562, T3m, T3l); |
442 | | T3o = T3k + T3n; |
443 | | T3s = T3n - T3k; |
444 | | T3p = Tt - Tm; |
445 | | T3q = T3f - T3e; |
446 | | T3r = FNMS(KP707106781, T3q, T3p); |
447 | | T3t = FMA(KP707106781, T3q, T3p); |
448 | | } |
449 | | Cr[WS(csr, 28)] = FNMS(KP923879532, T3o, T3h); |
450 | | Ci[WS(csi, 28)] = FMS(KP923879532, T3w, T3t); |
451 | | Cr[WS(csr, 4)] = FMA(KP923879532, T3o, T3h); |
452 | | Ci[WS(csi, 4)] = FMA(KP923879532, T3w, T3t); |
453 | | Ci[WS(csi, 12)] = FMS(KP923879532, T3s, T3r); |
454 | | Cr[WS(csr, 12)] = FMA(KP923879532, T3y, T3x); |
455 | | Ci[WS(csi, 20)] = FMA(KP923879532, T3s, T3r); |
456 | | Cr[WS(csr, 20)] = FNMS(KP923879532, T3y, T3x); |
457 | | } |
458 | | { |
459 | | E T2z, T2P, T2J, T2L, T2C, T2M, T2F, T2N; |
460 | | { |
461 | | E T2x, T2y, T2H, T2I; |
462 | | T2x = FNMS(KP707106781, T14, T11); |
463 | | T2y = T2n - T2m; |
464 | | T2z = FMA(KP923879532, T2y, T2x); |
465 | | T2P = FNMS(KP923879532, T2y, T2x); |
466 | | T2H = FNMS(KP707106781, T2k, T2j); |
467 | | T2I = T1b - T18; |
468 | | T2J = FMA(KP923879532, T2I, T2H); |
469 | | T2L = FNMS(KP923879532, T2I, T2H); |
470 | | } |
471 | | { |
472 | | E T2A, T2B, T2D, T2E; |
473 | | T2A = FNMS(KP707106781, T1z, T1k); |
474 | | T2B = FNMS(KP707106781, T1I, T1H); |
475 | | T2C = FNMS(KP668178637, T2B, T2A); |
476 | | T2M = FMA(KP668178637, T2A, T2B); |
477 | | T2D = FNMS(KP707106781, T26, T1R); |
478 | | T2E = FNMS(KP707106781, T2f, T2e); |
479 | | T2F = FMA(KP668178637, T2E, T2D); |
480 | | T2N = FNMS(KP668178637, T2D, T2E); |
481 | | } |
482 | | { |
483 | | E T2G, T2O, T2K, T2Q; |
484 | | T2G = T2C + T2F; |
485 | | Cr[WS(csr, 26)] = FNMS(KP831469612, T2G, T2z); |
486 | | Cr[WS(csr, 6)] = FMA(KP831469612, T2G, T2z); |
487 | | T2O = T2M + T2N; |
488 | | Ci[WS(csi, 6)] = -(FMA(KP831469612, T2O, T2L)); |
489 | | Ci[WS(csi, 26)] = FNMS(KP831469612, T2O, T2L); |
490 | | T2K = T2F - T2C; |
491 | | Ci[WS(csi, 10)] = FMA(KP831469612, T2K, T2J); |
492 | | Ci[WS(csi, 22)] = FMS(KP831469612, T2K, T2J); |
493 | | T2Q = T2M - T2N; |
494 | | Cr[WS(csr, 22)] = FNMS(KP831469612, T2Q, T2P); |
495 | | Cr[WS(csr, 10)] = FMA(KP831469612, T2Q, T2P); |
496 | | } |
497 | | } |
498 | | { |
499 | | E T1d, T2v, T2p, T2r, T1K, T2s, T2h, T2t; |
500 | | { |
501 | | E T15, T1c, T2l, T2o; |
502 | | T15 = FMA(KP707106781, T14, T11); |
503 | | T1c = T18 + T1b; |
504 | | T1d = FMA(KP923879532, T1c, T15); |
505 | | T2v = FNMS(KP923879532, T1c, T15); |
506 | | T2l = FMA(KP707106781, T2k, T2j); |
507 | | T2o = T2m + T2n; |
508 | | T2p = FNMS(KP923879532, T2o, T2l); |
509 | | T2r = FMA(KP923879532, T2o, T2l); |
510 | | } |
511 | | { |
512 | | E T1A, T1J, T27, T2g; |
513 | | T1A = FMA(KP707106781, T1z, T1k); |
514 | | T1J = FMA(KP707106781, T1I, T1H); |
515 | | T1K = FMA(KP198912367, T1J, T1A); |
516 | | T2s = FNMS(KP198912367, T1A, T1J); |
517 | | T27 = FMA(KP707106781, T26, T1R); |
518 | | T2g = FMA(KP707106781, T2f, T2e); |
519 | | T2h = FNMS(KP198912367, T2g, T27); |
520 | | T2t = FMA(KP198912367, T27, T2g); |
521 | | } |
522 | | { |
523 | | E T2i, T2u, T2q, T2w; |
524 | | T2i = T1K + T2h; |
525 | | Cr[WS(csr, 30)] = FNMS(KP980785280, T2i, T1d); |
526 | | Cr[WS(csr, 2)] = FMA(KP980785280, T2i, T1d); |
527 | | T2u = T2s + T2t; |
528 | | Ci[WS(csi, 2)] = FMA(KP980785280, T2u, T2r); |
529 | | Ci[WS(csi, 30)] = FMS(KP980785280, T2u, T2r); |
530 | | T2q = T2h - T1K; |
531 | | Ci[WS(csi, 14)] = FMS(KP980785280, T2q, T2p); |
532 | | Ci[WS(csi, 18)] = FMA(KP980785280, T2q, T2p); |
533 | | T2w = T2t - T2s; |
534 | | Cr[WS(csr, 18)] = FNMS(KP980785280, T2w, T2v); |
535 | | Cr[WS(csr, 14)] = FMA(KP980785280, T2w, T2v); |
536 | | } |
537 | | } |
538 | | { |
539 | | E T5r, T63, T6d, T5R, T5y, T6e, T6b, T6j, T5U, T64, T5G, T5Z, T68, T6i, T5N; |
540 | | E T5Y; |
541 | | { |
542 | | E T5u, T5x, T5C, T5F; |
543 | | T5r = FNMS(KP923879532, T5q, T5p); |
544 | | T63 = FMA(KP923879532, T5q, T5p); |
545 | | T6d = FMA(KP923879532, T5Q, T5P); |
546 | | T5R = FNMS(KP923879532, T5Q, T5P); |
547 | | T5u = FMA(KP668178637, T5t, T5s); |
548 | | T5x = FNMS(KP668178637, T5w, T5v); |
549 | | T5y = T5u - T5x; |
550 | | T6e = T5x + T5u; |
551 | | { |
552 | | E T69, T6a, T5S, T5T; |
553 | | T69 = FMA(KP923879532, T5I, T5H); |
554 | | T6a = FNMS(KP923879532, T5L, T5K); |
555 | | T6b = FMA(KP303346683, T6a, T69); |
556 | | T6j = FNMS(KP303346683, T69, T6a); |
557 | | T5S = FNMS(KP668178637, T5s, T5t); |
558 | | T5T = FMA(KP668178637, T5v, T5w); |
559 | | T5U = T5S - T5T; |
560 | | T64 = T5T + T5S; |
561 | | } |
562 | | T5C = FNMS(KP923879532, T5B, T5A); |
563 | | T5F = FNMS(KP923879532, T5E, T5D); |
564 | | T5G = FNMS(KP534511135, T5F, T5C); |
565 | | T5Z = FMA(KP534511135, T5C, T5F); |
566 | | { |
567 | | E T66, T67, T5J, T5M; |
568 | | T66 = FMA(KP923879532, T5B, T5A); |
569 | | T67 = FMA(KP923879532, T5E, T5D); |
570 | | T68 = FMA(KP303346683, T67, T66); |
571 | | T6i = FNMS(KP303346683, T66, T67); |
572 | | T5J = FNMS(KP923879532, T5I, T5H); |
573 | | T5M = FMA(KP923879532, T5L, T5K); |
574 | | T5N = FNMS(KP534511135, T5M, T5J); |
575 | | T5Y = FMA(KP534511135, T5J, T5M); |
576 | | } |
577 | | } |
578 | | { |
579 | | E T5z, T5O, T5X, T60; |
580 | | T5z = FMA(KP831469612, T5y, T5r); |
581 | | T5O = T5G + T5N; |
582 | | Cr[WS(csr, 27)] = FNMS(KP881921264, T5O, T5z); |
583 | | Cr[WS(csr, 5)] = FMA(KP881921264, T5O, T5z); |
584 | | T5X = FNMS(KP831469612, T5U, T5R); |
585 | | T60 = T5Y - T5Z; |
586 | | Ci[WS(csi, 5)] = FMS(KP881921264, T60, T5X); |
587 | | Ci[WS(csi, 27)] = FMA(KP881921264, T60, T5X); |
588 | | } |
589 | | { |
590 | | E T5V, T5W, T61, T62; |
591 | | T5V = FMA(KP831469612, T5U, T5R); |
592 | | T5W = T5N - T5G; |
593 | | Ci[WS(csi, 11)] = FMA(KP881921264, T5W, T5V); |
594 | | Ci[WS(csi, 21)] = FMS(KP881921264, T5W, T5V); |
595 | | T61 = FNMS(KP831469612, T5y, T5r); |
596 | | T62 = T5Z + T5Y; |
597 | | Cr[WS(csr, 21)] = FNMS(KP881921264, T62, T61); |
598 | | Cr[WS(csr, 11)] = FMA(KP881921264, T62, T61); |
599 | | } |
600 | | { |
601 | | E T65, T6c, T6h, T6k; |
602 | | T65 = FMA(KP831469612, T64, T63); |
603 | | T6c = T68 + T6b; |
604 | | Cr[WS(csr, 29)] = FNMS(KP956940335, T6c, T65); |
605 | | Cr[WS(csr, 3)] = FMA(KP956940335, T6c, T65); |
606 | | T6h = FMA(KP831469612, T6e, T6d); |
607 | | T6k = T6i - T6j; |
608 | | Ci[WS(csi, 3)] = FMA(KP956940335, T6k, T6h); |
609 | | Ci[WS(csi, 29)] = FMS(KP956940335, T6k, T6h); |
610 | | } |
611 | | { |
612 | | E T6f, T6g, T6l, T6m; |
613 | | T6f = FNMS(KP831469612, T6e, T6d); |
614 | | T6g = T6b - T68; |
615 | | Ci[WS(csi, 13)] = FMS(KP956940335, T6g, T6f); |
616 | | Ci[WS(csi, 19)] = FMA(KP956940335, T6g, T6f); |
617 | | T6l = FNMS(KP831469612, T64, T63); |
618 | | T6m = T6i + T6j; |
619 | | Cr[WS(csr, 19)] = FMA(KP956940335, T6m, T6l); |
620 | | Cr[WS(csr, 13)] = FNMS(KP956940335, T6m, T6l); |
621 | | } |
622 | | } |
623 | | { |
624 | | E T3L, T55, T5f, T4T, T44, T5g, T5d, T5l, T4W, T56, T4q, T51, T5a, T5k, T4L; |
625 | | E T50; |
626 | | { |
627 | | E T3U, T43, T4i, T4p; |
628 | | T3L = FMA(KP923879532, T3K, T3D); |
629 | | T55 = FNMS(KP923879532, T3K, T3D); |
630 | | T5f = FNMS(KP923879532, T4S, T4P); |
631 | | T4T = FMA(KP923879532, T4S, T4P); |
632 | | T3U = FNMS(KP198912367, T3T, T3Q); |
633 | | T43 = FMA(KP198912367, T42, T3Z); |
634 | | T44 = T3U + T43; |
635 | | T5g = T43 - T3U; |
636 | | { |
637 | | E T5b, T5c, T4U, T4V; |
638 | | T5b = FNMS(KP923879532, T4C, T4v); |
639 | | T5c = FNMS(KP923879532, T4J, T4G); |
640 | | T5d = FMA(KP820678790, T5c, T5b); |
641 | | T5l = FNMS(KP820678790, T5b, T5c); |
642 | | T4U = FMA(KP198912367, T3Q, T3T); |
643 | | T4V = FNMS(KP198912367, T3Z, T42); |
644 | | T4W = T4U + T4V; |
645 | | T56 = T4U - T4V; |
646 | | } |
647 | | T4i = FMA(KP923879532, T4h, T4a); |
648 | | T4p = FMA(KP923879532, T4o, T4l); |
649 | | T4q = FNMS(KP098491403, T4p, T4i); |
650 | | T51 = FMA(KP098491403, T4i, T4p); |
651 | | { |
652 | | E T58, T59, T4D, T4K; |
653 | | T58 = FNMS(KP923879532, T4h, T4a); |
654 | | T59 = FNMS(KP923879532, T4o, T4l); |
655 | | T5a = FMA(KP820678790, T59, T58); |
656 | | T5k = FNMS(KP820678790, T58, T59); |
657 | | T4D = FMA(KP923879532, T4C, T4v); |
658 | | T4K = FMA(KP923879532, T4J, T4G); |
659 | | T4L = FNMS(KP098491403, T4K, T4D); |
660 | | T50 = FMA(KP098491403, T4D, T4K); |
661 | | } |
662 | | } |
663 | | { |
664 | | E T45, T4M, T4Z, T52; |
665 | | T45 = FMA(KP980785280, T44, T3L); |
666 | | T4M = T4q + T4L; |
667 | | Cr[WS(csr, 31)] = FNMS(KP995184726, T4M, T45); |
668 | | Cr[WS(csr, 1)] = FMA(KP995184726, T4M, T45); |
669 | | T4Z = FMA(KP980785280, T4W, T4T); |
670 | | T52 = T50 - T51; |
671 | | Ci[WS(csi, 1)] = FMS(KP995184726, T52, T4Z); |
672 | | Ci[WS(csi, 31)] = FMA(KP995184726, T52, T4Z); |
673 | | } |
674 | | { |
675 | | E T4X, T4Y, T53, T54; |
676 | | T4X = FNMS(KP980785280, T4W, T4T); |
677 | | T4Y = T4L - T4q; |
678 | | Ci[WS(csi, 15)] = FMA(KP995184726, T4Y, T4X); |
679 | | Ci[WS(csi, 17)] = FMS(KP995184726, T4Y, T4X); |
680 | | T53 = FNMS(KP980785280, T44, T3L); |
681 | | T54 = T51 + T50; |
682 | | Cr[WS(csr, 17)] = FNMS(KP995184726, T54, T53); |
683 | | Cr[WS(csr, 15)] = FMA(KP995184726, T54, T53); |
684 | | } |
685 | | { |
686 | | E T57, T5e, T5j, T5m; |
687 | | T57 = FMA(KP980785280, T56, T55); |
688 | | T5e = T5a + T5d; |
689 | | Cr[WS(csr, 25)] = FNMS(KP773010453, T5e, T57); |
690 | | Cr[WS(csr, 7)] = FMA(KP773010453, T5e, T57); |
691 | | T5j = FMA(KP980785280, T5g, T5f); |
692 | | T5m = T5k - T5l; |
693 | | Ci[WS(csi, 7)] = FMA(KP773010453, T5m, T5j); |
694 | | Ci[WS(csi, 25)] = FMS(KP773010453, T5m, T5j); |
695 | | } |
696 | | { |
697 | | E T5h, T5i, T5n, T5o; |
698 | | T5h = FNMS(KP980785280, T5g, T5f); |
699 | | T5i = T5d - T5a; |
700 | | Ci[WS(csi, 9)] = FMS(KP773010453, T5i, T5h); |
701 | | Ci[WS(csi, 23)] = FMA(KP773010453, T5i, T5h); |
702 | | T5n = FNMS(KP980785280, T56, T55); |
703 | | T5o = T5k + T5l; |
704 | | Cr[WS(csr, 23)] = FMA(KP773010453, T5o, T5n); |
705 | | Cr[WS(csr, 9)] = FNMS(KP773010453, T5o, T5n); |
706 | | } |
707 | | } |
708 | | } |
709 | | } |
710 | | } |
711 | | |
712 | | static const kr2c_desc desc = { 64, "r2cf_64", { 198, 0, 196, 0 }, &GENUS }; |
713 | | |
714 | | void X(codelet_r2cf_64) (planner *p) { X(kr2c_register) (p, r2cf_64, &desc); |
715 | | } |
716 | | |
717 | | #else |
718 | | |
719 | | /* Generated by: ../../../genfft/gen_r2cf.native -compact -variables 4 -pipeline-latency 4 -n 64 -name r2cf_64 -include rdft/scalar/r2cf.h */ |
720 | | |
721 | | /* |
722 | | * This function contains 394 FP additions, 124 FP multiplications, |
723 | | * (or, 342 additions, 72 multiplications, 52 fused multiply/add), |
724 | | * 106 stack variables, 15 constants, and 128 memory accesses |
725 | | */ |
726 | | #include "rdft/scalar/r2cf.h" |
727 | | |
728 | | static void r2cf_64(R *R0, R *R1, R *Cr, R *Ci, stride rs, stride csr, stride csi, INT v, INT ivs, INT ovs) |
729 | 0 | { |
730 | 0 | DK(KP773010453, +0.773010453362736960810906609758469800971041293); |
731 | 0 | DK(KP634393284, +0.634393284163645498215171613225493370675687095); |
732 | 0 | DK(KP098017140, +0.098017140329560601994195563888641845861136673); |
733 | 0 | DK(KP995184726, +0.995184726672196886244836953109479921575474869); |
734 | 0 | DK(KP290284677, +0.290284677254462367636192375817395274691476278); |
735 | 0 | DK(KP956940335, +0.956940335732208864935797886980269969482849206); |
736 | 0 | DK(KP471396736, +0.471396736825997648556387625905254377657460319); |
737 | 0 | DK(KP881921264, +0.881921264348355029712756863660388349508442621); |
738 | 0 | DK(KP195090322, +0.195090322016128267848284868477022240927691618); |
739 | 0 | DK(KP980785280, +0.980785280403230449126182236134239036973933731); |
740 | 0 | DK(KP555570233, +0.555570233019602224742830813948532874374937191); |
741 | 0 | DK(KP831469612, +0.831469612302545237078788377617905756738560812); |
742 | 0 | DK(KP382683432, +0.382683432365089771728459984030398866761344562); |
743 | 0 | DK(KP923879532, +0.923879532511286756128183189396788286822416626); |
744 | 0 | DK(KP707106781, +0.707106781186547524400844362104849039284835938); |
745 | 0 | { |
746 | 0 | INT i; |
747 | 0 | for (i = v; i > 0; i = i - 1, R0 = R0 + ivs, R1 = R1 + ivs, Cr = Cr + ovs, Ci = Ci + ovs, MAKE_VOLATILE_STRIDE(256, rs), MAKE_VOLATILE_STRIDE(256, csr), MAKE_VOLATILE_STRIDE(256, csi)) { |
748 | 0 | E T4l, T5a, T15, T3n, T2T, T3Q, T7, Te, Tf, T4A, T4L, T1X, T3B, T23, T3y; |
749 | 0 | E T5I, T66, T4R, T52, T2j, T3F, T2H, T3I, T5P, T69, T1i, T3t, T1l, T3u, TZ; |
750 | 0 | E T63, T4v, T58, T1r, T3r, T1u, T3q, TK, T62, T4s, T57, Tm, Tt, Tu, T4o; |
751 | 0 | E T5b, T1c, T3R, T2Q, T3o, T1M, T3z, T5L, T67, T26, T3C, T4H, T4M, T2y, T3J; |
752 | 0 | E T5S, T6a, T2C, T3G, T4Y, T53; |
753 | 0 | { |
754 | 0 | E T3, T11, Td, T13, T6, T2S, Ta, T12, T14, T2R; |
755 | 0 | { |
756 | 0 | E T1, T2, Tb, Tc; |
757 | 0 | T1 = R0[0]; |
758 | 0 | T2 = R0[WS(rs, 16)]; |
759 | 0 | T3 = T1 + T2; |
760 | 0 | T11 = T1 - T2; |
761 | 0 | Tb = R0[WS(rs, 28)]; |
762 | 0 | Tc = R0[WS(rs, 12)]; |
763 | 0 | Td = Tb + Tc; |
764 | 0 | T13 = Tb - Tc; |
765 | 0 | } |
766 | 0 | { |
767 | 0 | E T4, T5, T8, T9; |
768 | 0 | T4 = R0[WS(rs, 8)]; |
769 | 0 | T5 = R0[WS(rs, 24)]; |
770 | 0 | T6 = T4 + T5; |
771 | 0 | T2S = T4 - T5; |
772 | 0 | T8 = R0[WS(rs, 4)]; |
773 | 0 | T9 = R0[WS(rs, 20)]; |
774 | 0 | Ta = T8 + T9; |
775 | 0 | T12 = T8 - T9; |
776 | 0 | } |
777 | 0 | T4l = T3 - T6; |
778 | 0 | T5a = Td - Ta; |
779 | 0 | T14 = KP707106781 * (T12 + T13); |
780 | 0 | T15 = T11 + T14; |
781 | 0 | T3n = T11 - T14; |
782 | 0 | T2R = KP707106781 * (T13 - T12); |
783 | 0 | T2T = T2R - T2S; |
784 | 0 | T3Q = T2S + T2R; |
785 | 0 | T7 = T3 + T6; |
786 | 0 | Te = Ta + Td; |
787 | 0 | Tf = T7 + Te; |
788 | 0 | } |
789 | 0 | { |
790 | 0 | E T1P, T4J, T21, T4y, T1S, T4K, T1W, T4z; |
791 | 0 | { |
792 | 0 | E T1N, T1O, T1Z, T20; |
793 | 0 | T1N = R1[WS(rs, 28)]; |
794 | 0 | T1O = R1[WS(rs, 12)]; |
795 | 0 | T1P = T1N - T1O; |
796 | 0 | T4J = T1N + T1O; |
797 | 0 | T1Z = R1[0]; |
798 | 0 | T20 = R1[WS(rs, 16)]; |
799 | 0 | T21 = T1Z - T20; |
800 | 0 | T4y = T1Z + T20; |
801 | 0 | } |
802 | 0 | { |
803 | 0 | E T1Q, T1R, T1U, T1V; |
804 | 0 | T1Q = R1[WS(rs, 4)]; |
805 | 0 | T1R = R1[WS(rs, 20)]; |
806 | 0 | T1S = T1Q - T1R; |
807 | 0 | T4K = T1Q + T1R; |
808 | 0 | T1U = R1[WS(rs, 8)]; |
809 | 0 | T1V = R1[WS(rs, 24)]; |
810 | 0 | T1W = T1U - T1V; |
811 | 0 | T4z = T1U + T1V; |
812 | 0 | } |
813 | 0 | T4A = T4y - T4z; |
814 | 0 | T4L = T4J - T4K; |
815 | 0 | { |
816 | 0 | E T1T, T22, T5G, T5H; |
817 | 0 | T1T = KP707106781 * (T1P - T1S); |
818 | 0 | T1X = T1T - T1W; |
819 | 0 | T3B = T1W + T1T; |
820 | 0 | T22 = KP707106781 * (T1S + T1P); |
821 | 0 | T23 = T21 + T22; |
822 | 0 | T3y = T21 - T22; |
823 | 0 | T5G = T4y + T4z; |
824 | 0 | T5H = T4K + T4J; |
825 | 0 | T5I = T5G + T5H; |
826 | 0 | T66 = T5G - T5H; |
827 | 0 | } |
828 | 0 | } |
829 | 0 | { |
830 | 0 | E T2b, T4P, T2G, T4Q, T2e, T51, T2h, T50; |
831 | 0 | { |
832 | 0 | E T29, T2a, T2E, T2F; |
833 | 0 | T29 = R1[WS(rs, 31)]; |
834 | 0 | T2a = R1[WS(rs, 15)]; |
835 | 0 | T2b = T29 - T2a; |
836 | 0 | T4P = T29 + T2a; |
837 | 0 | T2E = R1[WS(rs, 7)]; |
838 | 0 | T2F = R1[WS(rs, 23)]; |
839 | 0 | T2G = T2E - T2F; |
840 | 0 | T4Q = T2E + T2F; |
841 | 0 | } |
842 | 0 | { |
843 | 0 | E T2c, T2d, T2f, T2g; |
844 | 0 | T2c = R1[WS(rs, 3)]; |
845 | 0 | T2d = R1[WS(rs, 19)]; |
846 | 0 | T2e = T2c - T2d; |
847 | 0 | T51 = T2c + T2d; |
848 | 0 | T2f = R1[WS(rs, 27)]; |
849 | 0 | T2g = R1[WS(rs, 11)]; |
850 | 0 | T2h = T2f - T2g; |
851 | 0 | T50 = T2f + T2g; |
852 | 0 | } |
853 | 0 | T4R = T4P - T4Q; |
854 | 0 | T52 = T50 - T51; |
855 | 0 | { |
856 | 0 | E T2i, T2D, T5N, T5O; |
857 | 0 | T2i = KP707106781 * (T2e + T2h); |
858 | 0 | T2j = T2b + T2i; |
859 | 0 | T3F = T2b - T2i; |
860 | 0 | T2D = KP707106781 * (T2h - T2e); |
861 | 0 | T2H = T2D - T2G; |
862 | 0 | T3I = T2G + T2D; |
863 | 0 | T5N = T4P + T4Q; |
864 | 0 | T5O = T51 + T50; |
865 | 0 | T5P = T5N + T5O; |
866 | 0 | T69 = T5N - T5O; |
867 | 0 | } |
868 | 0 | } |
869 | 0 | { |
870 | 0 | E TN, T1e, TX, T1g, TQ, T1k, TU, T1f, T1h, T1j; |
871 | 0 | { |
872 | 0 | E TL, TM, TV, TW; |
873 | 0 | TL = R0[WS(rs, 31)]; |
874 | 0 | TM = R0[WS(rs, 15)]; |
875 | 0 | TN = TL + TM; |
876 | 0 | T1e = TL - TM; |
877 | 0 | TV = R0[WS(rs, 27)]; |
878 | 0 | TW = R0[WS(rs, 11)]; |
879 | 0 | TX = TV + TW; |
880 | 0 | T1g = TV - TW; |
881 | 0 | } |
882 | 0 | { |
883 | 0 | E TO, TP, TS, TT; |
884 | 0 | TO = R0[WS(rs, 7)]; |
885 | 0 | TP = R0[WS(rs, 23)]; |
886 | 0 | TQ = TO + TP; |
887 | 0 | T1k = TO - TP; |
888 | 0 | TS = R0[WS(rs, 3)]; |
889 | 0 | TT = R0[WS(rs, 19)]; |
890 | 0 | TU = TS + TT; |
891 | 0 | T1f = TS - TT; |
892 | 0 | } |
893 | 0 | T1h = KP707106781 * (T1f + T1g); |
894 | 0 | T1i = T1e + T1h; |
895 | 0 | T3t = T1e - T1h; |
896 | 0 | T1j = KP707106781 * (T1g - T1f); |
897 | 0 | T1l = T1j - T1k; |
898 | 0 | T3u = T1k + T1j; |
899 | 0 | { |
900 | 0 | E TR, TY, T4t, T4u; |
901 | 0 | TR = TN + TQ; |
902 | 0 | TY = TU + TX; |
903 | 0 | TZ = TR + TY; |
904 | 0 | T63 = TR - TY; |
905 | 0 | T4t = TN - TQ; |
906 | 0 | T4u = TX - TU; |
907 | 0 | T4v = FNMS(KP382683432, T4u, KP923879532 * T4t); |
908 | 0 | T58 = FMA(KP382683432, T4t, KP923879532 * T4u); |
909 | 0 | } |
910 | 0 | } |
911 | 0 | { |
912 | 0 | E Ty, T1s, TI, T1n, TB, T1q, TF, T1o, T1p, T1t; |
913 | 0 | { |
914 | 0 | E Tw, Tx, TG, TH; |
915 | 0 | Tw = R0[WS(rs, 1)]; |
916 | 0 | Tx = R0[WS(rs, 17)]; |
917 | 0 | Ty = Tw + Tx; |
918 | 0 | T1s = Tw - Tx; |
919 | 0 | TG = R0[WS(rs, 29)]; |
920 | 0 | TH = R0[WS(rs, 13)]; |
921 | 0 | TI = TG + TH; |
922 | 0 | T1n = TG - TH; |
923 | 0 | } |
924 | 0 | { |
925 | 0 | E Tz, TA, TD, TE; |
926 | 0 | Tz = R0[WS(rs, 9)]; |
927 | 0 | TA = R0[WS(rs, 25)]; |
928 | 0 | TB = Tz + TA; |
929 | 0 | T1q = Tz - TA; |
930 | 0 | TD = R0[WS(rs, 5)]; |
931 | 0 | TE = R0[WS(rs, 21)]; |
932 | 0 | TF = TD + TE; |
933 | 0 | T1o = TD - TE; |
934 | 0 | } |
935 | 0 | T1p = KP707106781 * (T1n - T1o); |
936 | 0 | T1r = T1p - T1q; |
937 | 0 | T3r = T1q + T1p; |
938 | 0 | T1t = KP707106781 * (T1o + T1n); |
939 | 0 | T1u = T1s + T1t; |
940 | 0 | T3q = T1s - T1t; |
941 | 0 | { |
942 | 0 | E TC, TJ, T4q, T4r; |
943 | 0 | TC = Ty + TB; |
944 | 0 | TJ = TF + TI; |
945 | 0 | TK = TC + TJ; |
946 | 0 | T62 = TC - TJ; |
947 | 0 | T4q = Ty - TB; |
948 | 0 | T4r = TI - TF; |
949 | 0 | T4s = FMA(KP923879532, T4q, KP382683432 * T4r); |
950 | 0 | T57 = FNMS(KP382683432, T4q, KP923879532 * T4r); |
951 | 0 | } |
952 | 0 | } |
953 | 0 | { |
954 | 0 | E Ti, T16, Ts, T1a, Tl, T17, Tp, T19, T4m, T4n; |
955 | 0 | { |
956 | 0 | E Tg, Th, Tq, Tr; |
957 | 0 | Tg = R0[WS(rs, 2)]; |
958 | 0 | Th = R0[WS(rs, 18)]; |
959 | 0 | Ti = Tg + Th; |
960 | 0 | T16 = Tg - Th; |
961 | 0 | Tq = R0[WS(rs, 6)]; |
962 | 0 | Tr = R0[WS(rs, 22)]; |
963 | 0 | Ts = Tq + Tr; |
964 | 0 | T1a = Tq - Tr; |
965 | 0 | } |
966 | 0 | { |
967 | 0 | E Tj, Tk, Tn, To; |
968 | 0 | Tj = R0[WS(rs, 10)]; |
969 | 0 | Tk = R0[WS(rs, 26)]; |
970 | 0 | Tl = Tj + Tk; |
971 | 0 | T17 = Tj - Tk; |
972 | 0 | Tn = R0[WS(rs, 30)]; |
973 | 0 | To = R0[WS(rs, 14)]; |
974 | 0 | Tp = Tn + To; |
975 | 0 | T19 = Tn - To; |
976 | 0 | } |
977 | 0 | Tm = Ti + Tl; |
978 | 0 | Tt = Tp + Ts; |
979 | 0 | Tu = Tm + Tt; |
980 | 0 | T4m = Ti - Tl; |
981 | 0 | T4n = Tp - Ts; |
982 | 0 | T4o = KP707106781 * (T4m + T4n); |
983 | 0 | T5b = KP707106781 * (T4n - T4m); |
984 | 0 | { |
985 | 0 | E T18, T1b, T2O, T2P; |
986 | 0 | T18 = FNMS(KP382683432, T17, KP923879532 * T16); |
987 | 0 | T1b = FMA(KP923879532, T19, KP382683432 * T1a); |
988 | 0 | T1c = T18 + T1b; |
989 | 0 | T3R = T1b - T18; |
990 | 0 | T2O = FNMS(KP923879532, T1a, KP382683432 * T19); |
991 | 0 | T2P = FMA(KP382683432, T16, KP923879532 * T17); |
992 | 0 | T2Q = T2O - T2P; |
993 | 0 | T3o = T2P + T2O; |
994 | 0 | } |
995 | 0 | } |
996 | 0 | { |
997 | 0 | E T1A, T4E, T1K, T4C, T1D, T4F, T1H, T4B; |
998 | 0 | { |
999 | 0 | E T1y, T1z, T1I, T1J; |
1000 | 0 | T1y = R1[WS(rs, 30)]; |
1001 | 0 | T1z = R1[WS(rs, 14)]; |
1002 | 0 | T1A = T1y - T1z; |
1003 | 0 | T4E = T1y + T1z; |
1004 | 0 | T1I = R1[WS(rs, 10)]; |
1005 | 0 | T1J = R1[WS(rs, 26)]; |
1006 | 0 | T1K = T1I - T1J; |
1007 | 0 | T4C = T1I + T1J; |
1008 | 0 | } |
1009 | 0 | { |
1010 | 0 | E T1B, T1C, T1F, T1G; |
1011 | 0 | T1B = R1[WS(rs, 6)]; |
1012 | 0 | T1C = R1[WS(rs, 22)]; |
1013 | 0 | T1D = T1B - T1C; |
1014 | 0 | T4F = T1B + T1C; |
1015 | 0 | T1F = R1[WS(rs, 2)]; |
1016 | 0 | T1G = R1[WS(rs, 18)]; |
1017 | 0 | T1H = T1F - T1G; |
1018 | 0 | T4B = T1F + T1G; |
1019 | 0 | } |
1020 | 0 | { |
1021 | 0 | E T1E, T1L, T5J, T5K; |
1022 | 0 | T1E = FNMS(KP923879532, T1D, KP382683432 * T1A); |
1023 | 0 | T1L = FMA(KP382683432, T1H, KP923879532 * T1K); |
1024 | 0 | T1M = T1E - T1L; |
1025 | 0 | T3z = T1L + T1E; |
1026 | 0 | T5J = T4B + T4C; |
1027 | 0 | T5K = T4E + T4F; |
1028 | 0 | T5L = T5J + T5K; |
1029 | 0 | T67 = T5K - T5J; |
1030 | 0 | } |
1031 | 0 | { |
1032 | 0 | E T24, T25, T4D, T4G; |
1033 | 0 | T24 = FNMS(KP382683432, T1K, KP923879532 * T1H); |
1034 | 0 | T25 = FMA(KP923879532, T1A, KP382683432 * T1D); |
1035 | 0 | T26 = T24 + T25; |
1036 | 0 | T3C = T25 - T24; |
1037 | 0 | T4D = T4B - T4C; |
1038 | 0 | T4G = T4E - T4F; |
1039 | 0 | T4H = KP707106781 * (T4D + T4G); |
1040 | 0 | T4M = KP707106781 * (T4G - T4D); |
1041 | 0 | } |
1042 | 0 | } |
1043 | 0 | { |
1044 | 0 | E T2m, T4S, T2w, T4W, T2p, T4T, T2t, T4V; |
1045 | 0 | { |
1046 | 0 | E T2k, T2l, T2u, T2v; |
1047 | 0 | T2k = R1[WS(rs, 1)]; |
1048 | 0 | T2l = R1[WS(rs, 17)]; |
1049 | 0 | T2m = T2k - T2l; |
1050 | 0 | T4S = T2k + T2l; |
1051 | 0 | T2u = R1[WS(rs, 5)]; |
1052 | 0 | T2v = R1[WS(rs, 21)]; |
1053 | 0 | T2w = T2u - T2v; |
1054 | 0 | T4W = T2u + T2v; |
1055 | 0 | } |
1056 | 0 | { |
1057 | 0 | E T2n, T2o, T2r, T2s; |
1058 | 0 | T2n = R1[WS(rs, 9)]; |
1059 | 0 | T2o = R1[WS(rs, 25)]; |
1060 | 0 | T2p = T2n - T2o; |
1061 | 0 | T4T = T2n + T2o; |
1062 | 0 | T2r = R1[WS(rs, 29)]; |
1063 | 0 | T2s = R1[WS(rs, 13)]; |
1064 | 0 | T2t = T2r - T2s; |
1065 | 0 | T4V = T2r + T2s; |
1066 | 0 | } |
1067 | 0 | { |
1068 | 0 | E T2q, T2x, T5Q, T5R; |
1069 | 0 | T2q = FNMS(KP382683432, T2p, KP923879532 * T2m); |
1070 | 0 | T2x = FMA(KP923879532, T2t, KP382683432 * T2w); |
1071 | 0 | T2y = T2q + T2x; |
1072 | 0 | T3J = T2x - T2q; |
1073 | 0 | T5Q = T4S + T4T; |
1074 | 0 | T5R = T4V + T4W; |
1075 | 0 | T5S = T5Q + T5R; |
1076 | 0 | T6a = T5R - T5Q; |
1077 | 0 | } |
1078 | 0 | { |
1079 | 0 | E T2A, T2B, T4U, T4X; |
1080 | 0 | T2A = FNMS(KP923879532, T2w, KP382683432 * T2t); |
1081 | 0 | T2B = FMA(KP382683432, T2m, KP923879532 * T2p); |
1082 | 0 | T2C = T2A - T2B; |
1083 | 0 | T3G = T2B + T2A; |
1084 | 0 | T4U = T4S - T4T; |
1085 | 0 | T4X = T4V - T4W; |
1086 | 0 | T4Y = KP707106781 * (T4U + T4X); |
1087 | 0 | T53 = KP707106781 * (T4X - T4U); |
1088 | 0 | } |
1089 | 0 | } |
1090 | 0 | { |
1091 | 0 | E Tv, T10, T5X, T5Y, T5Z, T60; |
1092 | 0 | Tv = Tf + Tu; |
1093 | 0 | T10 = TK + TZ; |
1094 | 0 | T5X = Tv + T10; |
1095 | 0 | T5Y = T5I + T5L; |
1096 | 0 | T5Z = T5P + T5S; |
1097 | 0 | T60 = T5Y + T5Z; |
1098 | 0 | Cr[WS(csr, 16)] = Tv - T10; |
1099 | 0 | Ci[WS(csi, 16)] = T5Z - T5Y; |
1100 | 0 | Cr[WS(csr, 32)] = T5X - T60; |
1101 | 0 | Cr[0] = T5X + T60; |
1102 | 0 | } |
1103 | 0 | { |
1104 | 0 | E T5F, T5V, T5U, T5W, T5M, T5T; |
1105 | 0 | T5F = Tf - Tu; |
1106 | 0 | T5V = TZ - TK; |
1107 | 0 | T5M = T5I - T5L; |
1108 | 0 | T5T = T5P - T5S; |
1109 | 0 | T5U = KP707106781 * (T5M + T5T); |
1110 | 0 | T5W = KP707106781 * (T5T - T5M); |
1111 | 0 | Cr[WS(csr, 24)] = T5F - T5U; |
1112 | 0 | Ci[WS(csi, 24)] = T5W - T5V; |
1113 | 0 | Cr[WS(csr, 8)] = T5F + T5U; |
1114 | 0 | Ci[WS(csi, 8)] = T5V + T5W; |
1115 | 0 | } |
1116 | 0 | { |
1117 | 0 | E T65, T6l, T6k, T6m, T6c, T6g, T6f, T6h; |
1118 | 0 | { |
1119 | 0 | E T61, T64, T6i, T6j; |
1120 | 0 | T61 = T7 - Te; |
1121 | 0 | T64 = KP707106781 * (T62 + T63); |
1122 | 0 | T65 = T61 + T64; |
1123 | 0 | T6l = T61 - T64; |
1124 | 0 | T6i = FNMS(KP382683432, T66, KP923879532 * T67); |
1125 | 0 | T6j = FMA(KP382683432, T69, KP923879532 * T6a); |
1126 | 0 | T6k = T6i + T6j; |
1127 | 0 | T6m = T6j - T6i; |
1128 | 0 | } |
1129 | 0 | { |
1130 | 0 | E T68, T6b, T6d, T6e; |
1131 | 0 | T68 = FMA(KP923879532, T66, KP382683432 * T67); |
1132 | 0 | T6b = FNMS(KP382683432, T6a, KP923879532 * T69); |
1133 | 0 | T6c = T68 + T6b; |
1134 | 0 | T6g = T6b - T68; |
1135 | 0 | T6d = KP707106781 * (T63 - T62); |
1136 | 0 | T6e = Tt - Tm; |
1137 | 0 | T6f = T6d - T6e; |
1138 | 0 | T6h = T6e + T6d; |
1139 | 0 | } |
1140 | 0 | Cr[WS(csr, 28)] = T65 - T6c; |
1141 | 0 | Ci[WS(csi, 28)] = T6k - T6h; |
1142 | 0 | Cr[WS(csr, 4)] = T65 + T6c; |
1143 | 0 | Ci[WS(csi, 4)] = T6h + T6k; |
1144 | 0 | Ci[WS(csi, 12)] = T6f + T6g; |
1145 | 0 | Cr[WS(csr, 12)] = T6l + T6m; |
1146 | 0 | Ci[WS(csi, 20)] = T6g - T6f; |
1147 | 0 | Cr[WS(csr, 20)] = T6l - T6m; |
1148 | 0 | } |
1149 | 0 | { |
1150 | 0 | E T5n, T5D, T5x, T5z, T5q, T5A, T5t, T5B; |
1151 | 0 | { |
1152 | 0 | E T5l, T5m, T5v, T5w; |
1153 | 0 | T5l = T4l - T4o; |
1154 | 0 | T5m = T58 - T57; |
1155 | 0 | T5n = T5l + T5m; |
1156 | 0 | T5D = T5l - T5m; |
1157 | 0 | T5v = T4v - T4s; |
1158 | 0 | T5w = T5b - T5a; |
1159 | 0 | T5x = T5v - T5w; |
1160 | 0 | T5z = T5w + T5v; |
1161 | 0 | } |
1162 | 0 | { |
1163 | 0 | E T5o, T5p, T5r, T5s; |
1164 | 0 | T5o = T4A - T4H; |
1165 | 0 | T5p = T4M - T4L; |
1166 | 0 | T5q = FMA(KP831469612, T5o, KP555570233 * T5p); |
1167 | 0 | T5A = FNMS(KP555570233, T5o, KP831469612 * T5p); |
1168 | 0 | T5r = T4R - T4Y; |
1169 | 0 | T5s = T53 - T52; |
1170 | 0 | T5t = FNMS(KP555570233, T5s, KP831469612 * T5r); |
1171 | 0 | T5B = FMA(KP555570233, T5r, KP831469612 * T5s); |
1172 | 0 | } |
1173 | 0 | { |
1174 | 0 | E T5u, T5C, T5y, T5E; |
1175 | 0 | T5u = T5q + T5t; |
1176 | 0 | Cr[WS(csr, 26)] = T5n - T5u; |
1177 | 0 | Cr[WS(csr, 6)] = T5n + T5u; |
1178 | 0 | T5C = T5A + T5B; |
1179 | 0 | Ci[WS(csi, 6)] = T5z + T5C; |
1180 | 0 | Ci[WS(csi, 26)] = T5C - T5z; |
1181 | 0 | T5y = T5t - T5q; |
1182 | 0 | Ci[WS(csi, 10)] = T5x + T5y; |
1183 | 0 | Ci[WS(csi, 22)] = T5y - T5x; |
1184 | 0 | T5E = T5B - T5A; |
1185 | 0 | Cr[WS(csr, 22)] = T5D - T5E; |
1186 | 0 | Cr[WS(csr, 10)] = T5D + T5E; |
1187 | 0 | } |
1188 | 0 | } |
1189 | 0 | { |
1190 | 0 | E T4x, T5j, T5d, T5f, T4O, T5g, T55, T5h; |
1191 | 0 | { |
1192 | 0 | E T4p, T4w, T59, T5c; |
1193 | 0 | T4p = T4l + T4o; |
1194 | 0 | T4w = T4s + T4v; |
1195 | 0 | T4x = T4p + T4w; |
1196 | 0 | T5j = T4p - T4w; |
1197 | 0 | T59 = T57 + T58; |
1198 | 0 | T5c = T5a + T5b; |
1199 | 0 | T5d = T59 - T5c; |
1200 | 0 | T5f = T5c + T59; |
1201 | 0 | } |
1202 | 0 | { |
1203 | 0 | E T4I, T4N, T4Z, T54; |
1204 | 0 | T4I = T4A + T4H; |
1205 | 0 | T4N = T4L + T4M; |
1206 | 0 | T4O = FMA(KP980785280, T4I, KP195090322 * T4N); |
1207 | 0 | T5g = FNMS(KP195090322, T4I, KP980785280 * T4N); |
1208 | 0 | T4Z = T4R + T4Y; |
1209 | 0 | T54 = T52 + T53; |
1210 | 0 | T55 = FNMS(KP195090322, T54, KP980785280 * T4Z); |
1211 | 0 | T5h = FMA(KP195090322, T4Z, KP980785280 * T54); |
1212 | 0 | } |
1213 | 0 | { |
1214 | 0 | E T56, T5i, T5e, T5k; |
1215 | 0 | T56 = T4O + T55; |
1216 | 0 | Cr[WS(csr, 30)] = T4x - T56; |
1217 | 0 | Cr[WS(csr, 2)] = T4x + T56; |
1218 | 0 | T5i = T5g + T5h; |
1219 | 0 | Ci[WS(csi, 2)] = T5f + T5i; |
1220 | 0 | Ci[WS(csi, 30)] = T5i - T5f; |
1221 | 0 | T5e = T55 - T4O; |
1222 | 0 | Ci[WS(csi, 14)] = T5d + T5e; |
1223 | 0 | Ci[WS(csi, 18)] = T5e - T5d; |
1224 | 0 | T5k = T5h - T5g; |
1225 | 0 | Cr[WS(csr, 18)] = T5j - T5k; |
1226 | 0 | Cr[WS(csr, 14)] = T5j + T5k; |
1227 | 0 | } |
1228 | 0 | } |
1229 | 0 | { |
1230 | 0 | E T3p, T41, T4c, T3S, T3w, T4b, T49, T4h, T3P, T42, T3E, T3W, T46, T4g, T3L; |
1231 | 0 | E T3X; |
1232 | 0 | { |
1233 | 0 | E T3s, T3v, T3A, T3D; |
1234 | 0 | T3p = T3n + T3o; |
1235 | 0 | T41 = T3n - T3o; |
1236 | 0 | T4c = T3R - T3Q; |
1237 | 0 | T3S = T3Q + T3R; |
1238 | 0 | T3s = FMA(KP831469612, T3q, KP555570233 * T3r); |
1239 | 0 | T3v = FNMS(KP555570233, T3u, KP831469612 * T3t); |
1240 | 0 | T3w = T3s + T3v; |
1241 | 0 | T4b = T3v - T3s; |
1242 | 0 | { |
1243 | 0 | E T47, T48, T3N, T3O; |
1244 | 0 | T47 = T3F - T3G; |
1245 | 0 | T48 = T3J - T3I; |
1246 | 0 | T49 = FNMS(KP471396736, T48, KP881921264 * T47); |
1247 | 0 | T4h = FMA(KP471396736, T47, KP881921264 * T48); |
1248 | 0 | T3N = FNMS(KP555570233, T3q, KP831469612 * T3r); |
1249 | 0 | T3O = FMA(KP555570233, T3t, KP831469612 * T3u); |
1250 | 0 | T3P = T3N + T3O; |
1251 | 0 | T42 = T3O - T3N; |
1252 | 0 | } |
1253 | 0 | T3A = T3y + T3z; |
1254 | 0 | T3D = T3B + T3C; |
1255 | 0 | T3E = FMA(KP956940335, T3A, KP290284677 * T3D); |
1256 | 0 | T3W = FNMS(KP290284677, T3A, KP956940335 * T3D); |
1257 | 0 | { |
1258 | 0 | E T44, T45, T3H, T3K; |
1259 | 0 | T44 = T3y - T3z; |
1260 | 0 | T45 = T3C - T3B; |
1261 | 0 | T46 = FMA(KP881921264, T44, KP471396736 * T45); |
1262 | 0 | T4g = FNMS(KP471396736, T44, KP881921264 * T45); |
1263 | 0 | T3H = T3F + T3G; |
1264 | 0 | T3K = T3I + T3J; |
1265 | 0 | T3L = FNMS(KP290284677, T3K, KP956940335 * T3H); |
1266 | 0 | T3X = FMA(KP290284677, T3H, KP956940335 * T3K); |
1267 | 0 | } |
1268 | 0 | } |
1269 | 0 | { |
1270 | 0 | E T3x, T3M, T3V, T3Y; |
1271 | 0 | T3x = T3p + T3w; |
1272 | 0 | T3M = T3E + T3L; |
1273 | 0 | Cr[WS(csr, 29)] = T3x - T3M; |
1274 | 0 | Cr[WS(csr, 3)] = T3x + T3M; |
1275 | 0 | T3V = T3S + T3P; |
1276 | 0 | T3Y = T3W + T3X; |
1277 | 0 | Ci[WS(csi, 3)] = T3V + T3Y; |
1278 | 0 | Ci[WS(csi, 29)] = T3Y - T3V; |
1279 | 0 | } |
1280 | 0 | { |
1281 | 0 | E T3T, T3U, T3Z, T40; |
1282 | 0 | T3T = T3P - T3S; |
1283 | 0 | T3U = T3L - T3E; |
1284 | 0 | Ci[WS(csi, 13)] = T3T + T3U; |
1285 | 0 | Ci[WS(csi, 19)] = T3U - T3T; |
1286 | 0 | T3Z = T3p - T3w; |
1287 | 0 | T40 = T3X - T3W; |
1288 | 0 | Cr[WS(csr, 19)] = T3Z - T40; |
1289 | 0 | Cr[WS(csr, 13)] = T3Z + T40; |
1290 | 0 | } |
1291 | 0 | { |
1292 | 0 | E T43, T4a, T4f, T4i; |
1293 | 0 | T43 = T41 + T42; |
1294 | 0 | T4a = T46 + T49; |
1295 | 0 | Cr[WS(csr, 27)] = T43 - T4a; |
1296 | 0 | Cr[WS(csr, 5)] = T43 + T4a; |
1297 | 0 | T4f = T4c + T4b; |
1298 | 0 | T4i = T4g + T4h; |
1299 | 0 | Ci[WS(csi, 5)] = T4f + T4i; |
1300 | 0 | Ci[WS(csi, 27)] = T4i - T4f; |
1301 | 0 | } |
1302 | 0 | { |
1303 | 0 | E T4d, T4e, T4j, T4k; |
1304 | 0 | T4d = T4b - T4c; |
1305 | 0 | T4e = T49 - T46; |
1306 | 0 | Ci[WS(csi, 11)] = T4d + T4e; |
1307 | 0 | Ci[WS(csi, 21)] = T4e - T4d; |
1308 | 0 | T4j = T41 - T42; |
1309 | 0 | T4k = T4h - T4g; |
1310 | 0 | Cr[WS(csr, 21)] = T4j - T4k; |
1311 | 0 | Cr[WS(csr, 11)] = T4j + T4k; |
1312 | 0 | } |
1313 | 0 | } |
1314 | 0 | { |
1315 | 0 | E T1d, T33, T3e, T2U, T1w, T3d, T3b, T3j, T2N, T34, T28, T2Y, T38, T3i, T2J; |
1316 | 0 | E T2Z; |
1317 | 0 | { |
1318 | 0 | E T1m, T1v, T1Y, T27; |
1319 | 0 | T1d = T15 - T1c; |
1320 | 0 | T33 = T15 + T1c; |
1321 | 0 | T3e = T2T + T2Q; |
1322 | 0 | T2U = T2Q - T2T; |
1323 | 0 | T1m = FMA(KP195090322, T1i, KP980785280 * T1l); |
1324 | 0 | T1v = FNMS(KP195090322, T1u, KP980785280 * T1r); |
1325 | 0 | T1w = T1m - T1v; |
1326 | 0 | T3d = T1v + T1m; |
1327 | 0 | { |
1328 | 0 | E T39, T3a, T2L, T2M; |
1329 | 0 | T39 = T2j + T2y; |
1330 | 0 | T3a = T2H + T2C; |
1331 | 0 | T3b = FNMS(KP098017140, T3a, KP995184726 * T39); |
1332 | 0 | T3j = FMA(KP995184726, T3a, KP098017140 * T39); |
1333 | 0 | T2L = FNMS(KP195090322, T1l, KP980785280 * T1i); |
1334 | 0 | T2M = FMA(KP980785280, T1u, KP195090322 * T1r); |
1335 | 0 | T2N = T2L - T2M; |
1336 | 0 | T34 = T2M + T2L; |
1337 | 0 | } |
1338 | 0 | T1Y = T1M - T1X; |
1339 | 0 | T27 = T23 - T26; |
1340 | 0 | T28 = FMA(KP634393284, T1Y, KP773010453 * T27); |
1341 | 0 | T2Y = FNMS(KP634393284, T27, KP773010453 * T1Y); |
1342 | 0 | { |
1343 | 0 | E T36, T37, T2z, T2I; |
1344 | 0 | T36 = T1X + T1M; |
1345 | 0 | T37 = T23 + T26; |
1346 | 0 | T38 = FMA(KP098017140, T36, KP995184726 * T37); |
1347 | 0 | T3i = FNMS(KP098017140, T37, KP995184726 * T36); |
1348 | 0 | T2z = T2j - T2y; |
1349 | 0 | T2I = T2C - T2H; |
1350 | 0 | T2J = FNMS(KP634393284, T2I, KP773010453 * T2z); |
1351 | 0 | T2Z = FMA(KP773010453, T2I, KP634393284 * T2z); |
1352 | 0 | } |
1353 | 0 | } |
1354 | 0 | { |
1355 | 0 | E T1x, T2K, T2X, T30; |
1356 | 0 | T1x = T1d + T1w; |
1357 | 0 | T2K = T28 + T2J; |
1358 | 0 | Cr[WS(csr, 25)] = T1x - T2K; |
1359 | 0 | Cr[WS(csr, 7)] = T1x + T2K; |
1360 | 0 | T2X = T2U + T2N; |
1361 | 0 | T30 = T2Y + T2Z; |
1362 | 0 | Ci[WS(csi, 7)] = T2X + T30; |
1363 | 0 | Ci[WS(csi, 25)] = T30 - T2X; |
1364 | 0 | } |
1365 | 0 | { |
1366 | 0 | E T2V, T2W, T31, T32; |
1367 | 0 | T2V = T2N - T2U; |
1368 | 0 | T2W = T2J - T28; |
1369 | 0 | Ci[WS(csi, 9)] = T2V + T2W; |
1370 | 0 | Ci[WS(csi, 23)] = T2W - T2V; |
1371 | 0 | T31 = T1d - T1w; |
1372 | 0 | T32 = T2Z - T2Y; |
1373 | 0 | Cr[WS(csr, 23)] = T31 - T32; |
1374 | 0 | Cr[WS(csr, 9)] = T31 + T32; |
1375 | 0 | } |
1376 | 0 | { |
1377 | 0 | E T35, T3c, T3h, T3k; |
1378 | 0 | T35 = T33 + T34; |
1379 | 0 | T3c = T38 + T3b; |
1380 | 0 | Cr[WS(csr, 31)] = T35 - T3c; |
1381 | 0 | Cr[WS(csr, 1)] = T35 + T3c; |
1382 | 0 | T3h = T3e + T3d; |
1383 | 0 | T3k = T3i + T3j; |
1384 | 0 | Ci[WS(csi, 1)] = T3h + T3k; |
1385 | 0 | Ci[WS(csi, 31)] = T3k - T3h; |
1386 | 0 | } |
1387 | 0 | { |
1388 | 0 | E T3f, T3g, T3l, T3m; |
1389 | 0 | T3f = T3d - T3e; |
1390 | 0 | T3g = T3b - T38; |
1391 | 0 | Ci[WS(csi, 15)] = T3f + T3g; |
1392 | 0 | Ci[WS(csi, 17)] = T3g - T3f; |
1393 | 0 | T3l = T33 - T34; |
1394 | 0 | T3m = T3j - T3i; |
1395 | 0 | Cr[WS(csr, 17)] = T3l - T3m; |
1396 | 0 | Cr[WS(csr, 15)] = T3l + T3m; |
1397 | 0 | } |
1398 | 0 | } |
1399 | 0 | } |
1400 | 0 | } |
1401 | 0 | } |
1402 | | |
1403 | | static const kr2c_desc desc = { 64, "r2cf_64", { 342, 72, 52, 0 }, &GENUS }; |
1404 | | |
1405 | 1 | void X(codelet_r2cf_64) (planner *p) { X(kr2c_register) (p, r2cf_64, &desc); |
1406 | 1 | } |
1407 | | |
1408 | | #endif |