/src/fftw3/rdft/scalar/r2cf/hf_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 Sat Jul 25 06:29:19 UTC 2026 */ |
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 -n 64 -dit -name hf_64 -include rdft/scalar/hf.h */ |
29 | | |
30 | | /* |
31 | | * This function contains 1038 FP additions, 644 FP multiplications, |
32 | | * (or, 520 additions, 126 multiplications, 518 fused multiply/add), |
33 | | * 190 stack variables, 15 constants, and 256 memory accesses |
34 | | */ |
35 | | #include "rdft/scalar/hf.h" |
36 | | |
37 | | static void hf_64(R *cr, R *ci, const R *W, stride rs, INT mb, INT me, INT ms) |
38 | | { |
39 | | DK(KP881921264, +0.881921264348355029712756863660388349508442621); |
40 | | DK(KP956940335, +0.956940335732208864935797886980269969482849206); |
41 | | DK(KP831469612, +0.831469612302545237078788377617905756738560812); |
42 | | DK(KP773010453, +0.773010453362736960810906609758469800971041293); |
43 | | DK(KP995184726, +0.995184726672196886244836953109479921575474869); |
44 | | DK(KP980785280, +0.980785280403230449126182236134239036973933731); |
45 | | DK(KP668178637, +0.668178637919298919997757686523080761552472251); |
46 | | DK(KP303346683, +0.303346683607342391675883946941299872384187453); |
47 | | DK(KP534511135, +0.534511135950791641089685961295362908582039528); |
48 | | DK(KP820678790, +0.820678790828660330972281985331011598767386482); |
49 | | DK(KP098491403, +0.098491403357164253077197521291327432293052451); |
50 | | DK(KP923879532, +0.923879532511286756128183189396788286822416626); |
51 | | DK(KP198912367, +0.198912367379658006911597622644676228597850501); |
52 | | DK(KP707106781, +0.707106781186547524400844362104849039284835938); |
53 | | DK(KP414213562, +0.414213562373095048801688724209698078569671875); |
54 | | { |
55 | | INT m; |
56 | | for (m = mb, W = W + ((mb - 1) * 126); m < me; m = m + 1, cr = cr + ms, ci = ci - ms, W = W + 126, MAKE_VOLATILE_STRIDE(128, rs)) { |
57 | | E Tm, TeM, TjR, Tkl, T7e, TcA, TiV, Tjm, T1G, TeW, TeZ, Thr, T7Q, TcI, T7X; |
58 | | E TcJ, T29, Tf8, Tf5, Thw, T87, TcN, T8u, TcQ, T5K, TfS, Tgc, ThX, Taq, Tdm; |
59 | | E Tbj, Tdx, TN, Tjl, TeP, TiP, T7l, TcB, T7s, TcC, T1f, TeR, TeU, Ths, T7B; |
60 | | E TcF, T7I, TcG, T32, Tfj, Tfg, ThB, T8G, TcU, T93, TcX, T3X, Tfr, TfK, ThM; |
61 | | E T9h, Td3, Taa, Tde, T2A, Tf6, Tfb, Thx, T8m, TcR, T8x, TcO, T3t, Tfh, Tfm; |
62 | | E ThC, T8V, TcY, T96, TcV, T4o, TfL, Tfu, ThN, T9w, Tdf, Tad, Td4, T6b, Tg9; |
63 | | E TfV, ThY, TaF, Tdy, Tbm, Tdn, T4Q, ThJ, TfA, TfN, T9M, Tdh, Taf, Td8, T5h; |
64 | | E ThI, TfF, TfO, Ta1, Tdi, Tag, Tdb, T6D, ThU, Tg1, Tgf, TaV, TdA, Tbo, Tdr; |
65 | | E T74, ThT, Tg6, Tge, Tba, TdB, Tbp, Tdu; |
66 | | { |
67 | | E T1, TiT, T7, TiS, Te, T7a, Tk, T7c; |
68 | | T1 = cr[0]; |
69 | | TiT = ci[0]; |
70 | | { |
71 | | E T3, T6, T4, TiR, T2, T5; |
72 | | T3 = cr[WS(rs, 32)]; |
73 | | T6 = ci[WS(rs, 32)]; |
74 | | T2 = W[62]; |
75 | | T4 = T2 * T3; |
76 | | TiR = T2 * T6; |
77 | | T5 = W[63]; |
78 | | T7 = FMA(T5, T6, T4); |
79 | | TiS = FNMS(T5, T3, TiR); |
80 | | } |
81 | | { |
82 | | E Ta, Td, Tb, T79, T9, Tc; |
83 | | Ta = cr[WS(rs, 16)]; |
84 | | Td = ci[WS(rs, 16)]; |
85 | | T9 = W[30]; |
86 | | Tb = T9 * Ta; |
87 | | T79 = T9 * Td; |
88 | | Tc = W[31]; |
89 | | Te = FMA(Tc, Td, Tb); |
90 | | T7a = FNMS(Tc, Ta, T79); |
91 | | } |
92 | | { |
93 | | E Tg, Tj, Th, T7b, Tf, Ti; |
94 | | Tg = cr[WS(rs, 48)]; |
95 | | Tj = ci[WS(rs, 48)]; |
96 | | Tf = W[94]; |
97 | | Th = Tf * Tg; |
98 | | T7b = Tf * Tj; |
99 | | Ti = W[95]; |
100 | | Tk = FMA(Ti, Tj, Th); |
101 | | T7c = FNMS(Ti, Tg, T7b); |
102 | | } |
103 | | { |
104 | | E T8, Tl, TjP, TjQ; |
105 | | T8 = T1 + T7; |
106 | | Tl = Te + Tk; |
107 | | Tm = T8 + Tl; |
108 | | TeM = T8 - Tl; |
109 | | TjP = Te - Tk; |
110 | | TjQ = TiT - TiS; |
111 | | TjR = TjP + TjQ; |
112 | | Tkl = TjQ - TjP; |
113 | | } |
114 | | { |
115 | | E T78, T7d, TiQ, TiU; |
116 | | T78 = T1 - T7; |
117 | | T7d = T7a - T7c; |
118 | | T7e = T78 - T7d; |
119 | | TcA = T78 + T7d; |
120 | | TiQ = T7a + T7c; |
121 | | TiU = TiS + TiT; |
122 | | TiV = TiQ + TiU; |
123 | | Tjm = TiU - TiQ; |
124 | | } |
125 | | } |
126 | | { |
127 | | E T1l, T7S, T1E, T7O, T1r, T7U, T1y, T7M; |
128 | | { |
129 | | E T1h, T1k, T1i, T7R, T1g, T1j; |
130 | | T1h = cr[WS(rs, 60)]; |
131 | | T1k = ci[WS(rs, 60)]; |
132 | | T1g = W[118]; |
133 | | T1i = T1g * T1h; |
134 | | T7R = T1g * T1k; |
135 | | T1j = W[119]; |
136 | | T1l = FMA(T1j, T1k, T1i); |
137 | | T7S = FNMS(T1j, T1h, T7R); |
138 | | } |
139 | | { |
140 | | E T1A, T1D, T1B, T7N, T1z, T1C; |
141 | | T1A = cr[WS(rs, 44)]; |
142 | | T1D = ci[WS(rs, 44)]; |
143 | | T1z = W[86]; |
144 | | T1B = T1z * T1A; |
145 | | T7N = T1z * T1D; |
146 | | T1C = W[87]; |
147 | | T1E = FMA(T1C, T1D, T1B); |
148 | | T7O = FNMS(T1C, T1A, T7N); |
149 | | } |
150 | | { |
151 | | E T1n, T1q, T1o, T7T, T1m, T1p; |
152 | | T1n = cr[WS(rs, 28)]; |
153 | | T1q = ci[WS(rs, 28)]; |
154 | | T1m = W[54]; |
155 | | T1o = T1m * T1n; |
156 | | T7T = T1m * T1q; |
157 | | T1p = W[55]; |
158 | | T1r = FMA(T1p, T1q, T1o); |
159 | | T7U = FNMS(T1p, T1n, T7T); |
160 | | } |
161 | | { |
162 | | E T1u, T1x, T1v, T7L, T1t, T1w; |
163 | | T1u = cr[WS(rs, 12)]; |
164 | | T1x = ci[WS(rs, 12)]; |
165 | | T1t = W[22]; |
166 | | T1v = T1t * T1u; |
167 | | T7L = T1t * T1x; |
168 | | T1w = W[23]; |
169 | | T1y = FMA(T1w, T1x, T1v); |
170 | | T7M = FNMS(T1w, T1u, T7L); |
171 | | } |
172 | | { |
173 | | E T1s, T1F, TeX, TeY; |
174 | | T1s = T1l + T1r; |
175 | | T1F = T1y + T1E; |
176 | | T1G = T1s + T1F; |
177 | | TeW = T1s - T1F; |
178 | | TeX = T7S + T7U; |
179 | | TeY = T7M + T7O; |
180 | | TeZ = TeX - TeY; |
181 | | Thr = TeX + TeY; |
182 | | } |
183 | | { |
184 | | E T7K, T7P, T7V, T7W; |
185 | | T7K = T1l - T1r; |
186 | | T7P = T7M - T7O; |
187 | | T7Q = T7K - T7P; |
188 | | TcI = T7K + T7P; |
189 | | T7V = T7S - T7U; |
190 | | T7W = T1y - T1E; |
191 | | T7X = T7V + T7W; |
192 | | TcJ = T7V - T7W; |
193 | | } |
194 | | } |
195 | | { |
196 | | E T1O, T8p, T27, T85, T1U, T8r, T21, T83; |
197 | | { |
198 | | E T1K, T1N, T1L, T8o, T1J, T1M; |
199 | | T1K = cr[WS(rs, 2)]; |
200 | | T1N = ci[WS(rs, 2)]; |
201 | | T1J = W[2]; |
202 | | T1L = T1J * T1K; |
203 | | T8o = T1J * T1N; |
204 | | T1M = W[3]; |
205 | | T1O = FMA(T1M, T1N, T1L); |
206 | | T8p = FNMS(T1M, T1K, T8o); |
207 | | } |
208 | | { |
209 | | E T23, T26, T24, T84, T22, T25; |
210 | | T23 = cr[WS(rs, 50)]; |
211 | | T26 = ci[WS(rs, 50)]; |
212 | | T22 = W[98]; |
213 | | T24 = T22 * T23; |
214 | | T84 = T22 * T26; |
215 | | T25 = W[99]; |
216 | | T27 = FMA(T25, T26, T24); |
217 | | T85 = FNMS(T25, T23, T84); |
218 | | } |
219 | | { |
220 | | E T1Q, T1T, T1R, T8q, T1P, T1S; |
221 | | T1Q = cr[WS(rs, 34)]; |
222 | | T1T = ci[WS(rs, 34)]; |
223 | | T1P = W[66]; |
224 | | T1R = T1P * T1Q; |
225 | | T8q = T1P * T1T; |
226 | | T1S = W[67]; |
227 | | T1U = FMA(T1S, T1T, T1R); |
228 | | T8r = FNMS(T1S, T1Q, T8q); |
229 | | } |
230 | | { |
231 | | E T1X, T20, T1Y, T82, T1W, T1Z; |
232 | | T1X = cr[WS(rs, 18)]; |
233 | | T20 = ci[WS(rs, 18)]; |
234 | | T1W = W[34]; |
235 | | T1Y = T1W * T1X; |
236 | | T82 = T1W * T20; |
237 | | T1Z = W[35]; |
238 | | T21 = FMA(T1Z, T20, T1Y); |
239 | | T83 = FNMS(T1Z, T1X, T82); |
240 | | } |
241 | | { |
242 | | E T1V, T28, Tf3, Tf4; |
243 | | T1V = T1O + T1U; |
244 | | T28 = T21 + T27; |
245 | | T29 = T1V + T28; |
246 | | Tf8 = T1V - T28; |
247 | | Tf3 = T8p + T8r; |
248 | | Tf4 = T83 + T85; |
249 | | Tf5 = Tf3 - Tf4; |
250 | | Thw = Tf3 + Tf4; |
251 | | } |
252 | | { |
253 | | E T81, T86, T8s, T8t; |
254 | | T81 = T1O - T1U; |
255 | | T86 = T83 - T85; |
256 | | T87 = T81 - T86; |
257 | | TcN = T81 + T86; |
258 | | T8s = T8p - T8r; |
259 | | T8t = T21 - T27; |
260 | | T8u = T8s + T8t; |
261 | | TcQ = T8s - T8t; |
262 | | } |
263 | | } |
264 | | { |
265 | | E T5p, Tbf, T5I, Tao, T5v, Tbh, T5C, Tam; |
266 | | { |
267 | | E T5l, T5o, T5m, Tbe, T5k, T5n; |
268 | | T5l = cr[WS(rs, 63)]; |
269 | | T5o = ci[WS(rs, 63)]; |
270 | | T5k = W[124]; |
271 | | T5m = T5k * T5l; |
272 | | Tbe = T5k * T5o; |
273 | | T5n = W[125]; |
274 | | T5p = FMA(T5n, T5o, T5m); |
275 | | Tbf = FNMS(T5n, T5l, Tbe); |
276 | | } |
277 | | { |
278 | | E T5E, T5H, T5F, Tan, T5D, T5G; |
279 | | T5E = cr[WS(rs, 47)]; |
280 | | T5H = ci[WS(rs, 47)]; |
281 | | T5D = W[92]; |
282 | | T5F = T5D * T5E; |
283 | | Tan = T5D * T5H; |
284 | | T5G = W[93]; |
285 | | T5I = FMA(T5G, T5H, T5F); |
286 | | Tao = FNMS(T5G, T5E, Tan); |
287 | | } |
288 | | { |
289 | | E T5r, T5u, T5s, Tbg, T5q, T5t; |
290 | | T5r = cr[WS(rs, 31)]; |
291 | | T5u = ci[WS(rs, 31)]; |
292 | | T5q = W[60]; |
293 | | T5s = T5q * T5r; |
294 | | Tbg = T5q * T5u; |
295 | | T5t = W[61]; |
296 | | T5v = FMA(T5t, T5u, T5s); |
297 | | Tbh = FNMS(T5t, T5r, Tbg); |
298 | | } |
299 | | { |
300 | | E T5y, T5B, T5z, Tal, T5x, T5A; |
301 | | T5y = cr[WS(rs, 15)]; |
302 | | T5B = ci[WS(rs, 15)]; |
303 | | T5x = W[28]; |
304 | | T5z = T5x * T5y; |
305 | | Tal = T5x * T5B; |
306 | | T5A = W[29]; |
307 | | T5C = FMA(T5A, T5B, T5z); |
308 | | Tam = FNMS(T5A, T5y, Tal); |
309 | | } |
310 | | { |
311 | | E T5w, T5J, Tga, Tgb; |
312 | | T5w = T5p + T5v; |
313 | | T5J = T5C + T5I; |
314 | | T5K = T5w + T5J; |
315 | | TfS = T5w - T5J; |
316 | | Tga = Tbf + Tbh; |
317 | | Tgb = Tam + Tao; |
318 | | Tgc = Tga - Tgb; |
319 | | ThX = Tga + Tgb; |
320 | | } |
321 | | { |
322 | | E Tak, Tap, Tbd, Tbi; |
323 | | Tak = T5p - T5v; |
324 | | Tap = Tam - Tao; |
325 | | Taq = Tak - Tap; |
326 | | Tdm = Tak + Tap; |
327 | | Tbd = T5I - T5C; |
328 | | Tbi = Tbf - Tbh; |
329 | | Tbj = Tbd - Tbi; |
330 | | Tdx = Tbi + Tbd; |
331 | | } |
332 | | } |
333 | | { |
334 | | E Ts, T7h, TL, T7q, Ty, T7j, TF, T7o; |
335 | | { |
336 | | E To, Tr, Tp, T7g, Tn, Tq; |
337 | | To = cr[WS(rs, 8)]; |
338 | | Tr = ci[WS(rs, 8)]; |
339 | | Tn = W[14]; |
340 | | Tp = Tn * To; |
341 | | T7g = Tn * Tr; |
342 | | Tq = W[15]; |
343 | | Ts = FMA(Tq, Tr, Tp); |
344 | | T7h = FNMS(Tq, To, T7g); |
345 | | } |
346 | | { |
347 | | E TH, TK, TI, T7p, TG, TJ; |
348 | | TH = cr[WS(rs, 24)]; |
349 | | TK = ci[WS(rs, 24)]; |
350 | | TG = W[46]; |
351 | | TI = TG * TH; |
352 | | T7p = TG * TK; |
353 | | TJ = W[47]; |
354 | | TL = FMA(TJ, TK, TI); |
355 | | T7q = FNMS(TJ, TH, T7p); |
356 | | } |
357 | | { |
358 | | E Tu, Tx, Tv, T7i, Tt, Tw; |
359 | | Tu = cr[WS(rs, 40)]; |
360 | | Tx = ci[WS(rs, 40)]; |
361 | | Tt = W[78]; |
362 | | Tv = Tt * Tu; |
363 | | T7i = Tt * Tx; |
364 | | Tw = W[79]; |
365 | | Ty = FMA(Tw, Tx, Tv); |
366 | | T7j = FNMS(Tw, Tu, T7i); |
367 | | } |
368 | | { |
369 | | E TB, TE, TC, T7n, TA, TD; |
370 | | TB = cr[WS(rs, 56)]; |
371 | | TE = ci[WS(rs, 56)]; |
372 | | TA = W[110]; |
373 | | TC = TA * TB; |
374 | | T7n = TA * TE; |
375 | | TD = W[111]; |
376 | | TF = FMA(TD, TE, TC); |
377 | | T7o = FNMS(TD, TB, T7n); |
378 | | } |
379 | | { |
380 | | E Tz, TM, TeN, TeO; |
381 | | Tz = Ts + Ty; |
382 | | TM = TF + TL; |
383 | | TN = Tz + TM; |
384 | | Tjl = Tz - TM; |
385 | | TeN = T7o + T7q; |
386 | | TeO = T7h + T7j; |
387 | | TeP = TeN - TeO; |
388 | | TiP = TeO + TeN; |
389 | | } |
390 | | { |
391 | | E T7f, T7k, T7m, T7r; |
392 | | T7f = Ts - Ty; |
393 | | T7k = T7h - T7j; |
394 | | T7l = T7f - T7k; |
395 | | TcB = T7f + T7k; |
396 | | T7m = TF - TL; |
397 | | T7r = T7o - T7q; |
398 | | T7s = T7m + T7r; |
399 | | TcC = T7m - T7r; |
400 | | } |
401 | | } |
402 | | { |
403 | | E TU, T7D, T1d, T7z, T10, T7F, T17, T7x; |
404 | | { |
405 | | E TQ, TT, TR, T7C, TP, TS; |
406 | | TQ = cr[WS(rs, 4)]; |
407 | | TT = ci[WS(rs, 4)]; |
408 | | TP = W[6]; |
409 | | TR = TP * TQ; |
410 | | T7C = TP * TT; |
411 | | TS = W[7]; |
412 | | TU = FMA(TS, TT, TR); |
413 | | T7D = FNMS(TS, TQ, T7C); |
414 | | } |
415 | | { |
416 | | E T19, T1c, T1a, T7y, T18, T1b; |
417 | | T19 = cr[WS(rs, 52)]; |
418 | | T1c = ci[WS(rs, 52)]; |
419 | | T18 = W[102]; |
420 | | T1a = T18 * T19; |
421 | | T7y = T18 * T1c; |
422 | | T1b = W[103]; |
423 | | T1d = FMA(T1b, T1c, T1a); |
424 | | T7z = FNMS(T1b, T19, T7y); |
425 | | } |
426 | | { |
427 | | E TW, TZ, TX, T7E, TV, TY; |
428 | | TW = cr[WS(rs, 36)]; |
429 | | TZ = ci[WS(rs, 36)]; |
430 | | TV = W[70]; |
431 | | TX = TV * TW; |
432 | | T7E = TV * TZ; |
433 | | TY = W[71]; |
434 | | T10 = FMA(TY, TZ, TX); |
435 | | T7F = FNMS(TY, TW, T7E); |
436 | | } |
437 | | { |
438 | | E T13, T16, T14, T7w, T12, T15; |
439 | | T13 = cr[WS(rs, 20)]; |
440 | | T16 = ci[WS(rs, 20)]; |
441 | | T12 = W[38]; |
442 | | T14 = T12 * T13; |
443 | | T7w = T12 * T16; |
444 | | T15 = W[39]; |
445 | | T17 = FMA(T15, T16, T14); |
446 | | T7x = FNMS(T15, T13, T7w); |
447 | | } |
448 | | { |
449 | | E T11, T1e, TeS, TeT; |
450 | | T11 = TU + T10; |
451 | | T1e = T17 + T1d; |
452 | | T1f = T11 + T1e; |
453 | | TeR = T11 - T1e; |
454 | | TeS = T7D + T7F; |
455 | | TeT = T7x + T7z; |
456 | | TeU = TeS - TeT; |
457 | | Ths = TeS + TeT; |
458 | | } |
459 | | { |
460 | | E T7v, T7A, T7G, T7H; |
461 | | T7v = TU - T10; |
462 | | T7A = T7x - T7z; |
463 | | T7B = T7v - T7A; |
464 | | TcF = T7v + T7A; |
465 | | T7G = T7D - T7F; |
466 | | T7H = T17 - T1d; |
467 | | T7I = T7G + T7H; |
468 | | TcG = T7G - T7H; |
469 | | } |
470 | | } |
471 | | { |
472 | | E T2H, T8Y, T30, T8E, T2N, T90, T2U, T8C; |
473 | | { |
474 | | E T2D, T2G, T2E, T8X, T2C, T2F; |
475 | | T2D = cr[WS(rs, 62)]; |
476 | | T2G = ci[WS(rs, 62)]; |
477 | | T2C = W[122]; |
478 | | T2E = T2C * T2D; |
479 | | T8X = T2C * T2G; |
480 | | T2F = W[123]; |
481 | | T2H = FMA(T2F, T2G, T2E); |
482 | | T8Y = FNMS(T2F, T2D, T8X); |
483 | | } |
484 | | { |
485 | | E T2W, T2Z, T2X, T8D, T2V, T2Y; |
486 | | T2W = cr[WS(rs, 46)]; |
487 | | T2Z = ci[WS(rs, 46)]; |
488 | | T2V = W[90]; |
489 | | T2X = T2V * T2W; |
490 | | T8D = T2V * T2Z; |
491 | | T2Y = W[91]; |
492 | | T30 = FMA(T2Y, T2Z, T2X); |
493 | | T8E = FNMS(T2Y, T2W, T8D); |
494 | | } |
495 | | { |
496 | | E T2J, T2M, T2K, T8Z, T2I, T2L; |
497 | | T2J = cr[WS(rs, 30)]; |
498 | | T2M = ci[WS(rs, 30)]; |
499 | | T2I = W[58]; |
500 | | T2K = T2I * T2J; |
501 | | T8Z = T2I * T2M; |
502 | | T2L = W[59]; |
503 | | T2N = FMA(T2L, T2M, T2K); |
504 | | T90 = FNMS(T2L, T2J, T8Z); |
505 | | } |
506 | | { |
507 | | E T2Q, T2T, T2R, T8B, T2P, T2S; |
508 | | T2Q = cr[WS(rs, 14)]; |
509 | | T2T = ci[WS(rs, 14)]; |
510 | | T2P = W[26]; |
511 | | T2R = T2P * T2Q; |
512 | | T8B = T2P * T2T; |
513 | | T2S = W[27]; |
514 | | T2U = FMA(T2S, T2T, T2R); |
515 | | T8C = FNMS(T2S, T2Q, T8B); |
516 | | } |
517 | | { |
518 | | E T2O, T31, Tfe, Tff; |
519 | | T2O = T2H + T2N; |
520 | | T31 = T2U + T30; |
521 | | T32 = T2O + T31; |
522 | | Tfj = T2O - T31; |
523 | | Tfe = T8Y + T90; |
524 | | Tff = T8C + T8E; |
525 | | Tfg = Tfe - Tff; |
526 | | ThB = Tfe + Tff; |
527 | | } |
528 | | { |
529 | | E T8A, T8F, T91, T92; |
530 | | T8A = T2H - T2N; |
531 | | T8F = T8C - T8E; |
532 | | T8G = T8A - T8F; |
533 | | TcU = T8A + T8F; |
534 | | T91 = T8Y - T90; |
535 | | T92 = T2U - T30; |
536 | | T93 = T91 + T92; |
537 | | TcX = T91 - T92; |
538 | | } |
539 | | } |
540 | | { |
541 | | E T3C, Ta5, T3V, T9f, T3I, Ta7, T3P, T9d; |
542 | | { |
543 | | E T3y, T3B, T3z, Ta4, T3x, T3A; |
544 | | T3y = cr[WS(rs, 1)]; |
545 | | T3B = ci[WS(rs, 1)]; |
546 | | T3x = W[0]; |
547 | | T3z = T3x * T3y; |
548 | | Ta4 = T3x * T3B; |
549 | | T3A = W[1]; |
550 | | T3C = FMA(T3A, T3B, T3z); |
551 | | Ta5 = FNMS(T3A, T3y, Ta4); |
552 | | } |
553 | | { |
554 | | E T3R, T3U, T3S, T9e, T3Q, T3T; |
555 | | T3R = cr[WS(rs, 49)]; |
556 | | T3U = ci[WS(rs, 49)]; |
557 | | T3Q = W[96]; |
558 | | T3S = T3Q * T3R; |
559 | | T9e = T3Q * T3U; |
560 | | T3T = W[97]; |
561 | | T3V = FMA(T3T, T3U, T3S); |
562 | | T9f = FNMS(T3T, T3R, T9e); |
563 | | } |
564 | | { |
565 | | E T3E, T3H, T3F, Ta6, T3D, T3G; |
566 | | T3E = cr[WS(rs, 33)]; |
567 | | T3H = ci[WS(rs, 33)]; |
568 | | T3D = W[64]; |
569 | | T3F = T3D * T3E; |
570 | | Ta6 = T3D * T3H; |
571 | | T3G = W[65]; |
572 | | T3I = FMA(T3G, T3H, T3F); |
573 | | Ta7 = FNMS(T3G, T3E, Ta6); |
574 | | } |
575 | | { |
576 | | E T3L, T3O, T3M, T9c, T3K, T3N; |
577 | | T3L = cr[WS(rs, 17)]; |
578 | | T3O = ci[WS(rs, 17)]; |
579 | | T3K = W[32]; |
580 | | T3M = T3K * T3L; |
581 | | T9c = T3K * T3O; |
582 | | T3N = W[33]; |
583 | | T3P = FMA(T3N, T3O, T3M); |
584 | | T9d = FNMS(T3N, T3L, T9c); |
585 | | } |
586 | | { |
587 | | E T3J, T3W, TfI, TfJ; |
588 | | T3J = T3C + T3I; |
589 | | T3W = T3P + T3V; |
590 | | T3X = T3J + T3W; |
591 | | Tfr = T3J - T3W; |
592 | | TfI = Ta5 + Ta7; |
593 | | TfJ = T9d + T9f; |
594 | | TfK = TfI - TfJ; |
595 | | ThM = TfI + TfJ; |
596 | | } |
597 | | { |
598 | | E T9b, T9g, Ta8, Ta9; |
599 | | T9b = T3C - T3I; |
600 | | T9g = T9d - T9f; |
601 | | T9h = T9b - T9g; |
602 | | Td3 = T9b + T9g; |
603 | | Ta8 = Ta5 - Ta7; |
604 | | Ta9 = T3P - T3V; |
605 | | Taa = Ta8 + Ta9; |
606 | | Tde = Ta8 - Ta9; |
607 | | } |
608 | | } |
609 | | { |
610 | | E T2f, T8a, T2y, T8j, T2l, T8c, T2s, T8h; |
611 | | { |
612 | | E T2b, T2e, T2c, T89, T2a, T2d; |
613 | | T2b = cr[WS(rs, 10)]; |
614 | | T2e = ci[WS(rs, 10)]; |
615 | | T2a = W[18]; |
616 | | T2c = T2a * T2b; |
617 | | T89 = T2a * T2e; |
618 | | T2d = W[19]; |
619 | | T2f = FMA(T2d, T2e, T2c); |
620 | | T8a = FNMS(T2d, T2b, T89); |
621 | | } |
622 | | { |
623 | | E T2u, T2x, T2v, T8i, T2t, T2w; |
624 | | T2u = cr[WS(rs, 26)]; |
625 | | T2x = ci[WS(rs, 26)]; |
626 | | T2t = W[50]; |
627 | | T2v = T2t * T2u; |
628 | | T8i = T2t * T2x; |
629 | | T2w = W[51]; |
630 | | T2y = FMA(T2w, T2x, T2v); |
631 | | T8j = FNMS(T2w, T2u, T8i); |
632 | | } |
633 | | { |
634 | | E T2h, T2k, T2i, T8b, T2g, T2j; |
635 | | T2h = cr[WS(rs, 42)]; |
636 | | T2k = ci[WS(rs, 42)]; |
637 | | T2g = W[82]; |
638 | | T2i = T2g * T2h; |
639 | | T8b = T2g * T2k; |
640 | | T2j = W[83]; |
641 | | T2l = FMA(T2j, T2k, T2i); |
642 | | T8c = FNMS(T2j, T2h, T8b); |
643 | | } |
644 | | { |
645 | | E T2o, T2r, T2p, T8g, T2n, T2q; |
646 | | T2o = cr[WS(rs, 58)]; |
647 | | T2r = ci[WS(rs, 58)]; |
648 | | T2n = W[114]; |
649 | | T2p = T2n * T2o; |
650 | | T8g = T2n * T2r; |
651 | | T2q = W[115]; |
652 | | T2s = FMA(T2q, T2r, T2p); |
653 | | T8h = FNMS(T2q, T2o, T8g); |
654 | | } |
655 | | { |
656 | | E T2m, T2z, Tf9, Tfa; |
657 | | T2m = T2f + T2l; |
658 | | T2z = T2s + T2y; |
659 | | T2A = T2m + T2z; |
660 | | Tf6 = T2m - T2z; |
661 | | Tf9 = T8h + T8j; |
662 | | Tfa = T8a + T8c; |
663 | | Tfb = Tf9 - Tfa; |
664 | | Thx = Tfa + Tf9; |
665 | | { |
666 | | E T8e, T8v, T8l, T8w; |
667 | | { |
668 | | E T88, T8d, T8f, T8k; |
669 | | T88 = T2f - T2l; |
670 | | T8d = T8a - T8c; |
671 | | T8e = T88 - T8d; |
672 | | T8v = T88 + T8d; |
673 | | T8f = T2s - T2y; |
674 | | T8k = T8h - T8j; |
675 | | T8l = T8f + T8k; |
676 | | T8w = T8k - T8f; |
677 | | } |
678 | | T8m = T8e + T8l; |
679 | | TcR = T8l - T8e; |
680 | | T8x = T8v + T8w; |
681 | | TcO = T8v - T8w; |
682 | | } |
683 | | } |
684 | | } |
685 | | { |
686 | | E T38, T8J, T3r, T8S, T3e, T8L, T3l, T8Q; |
687 | | { |
688 | | E T34, T37, T35, T8I, T33, T36; |
689 | | T34 = cr[WS(rs, 6)]; |
690 | | T37 = ci[WS(rs, 6)]; |
691 | | T33 = W[10]; |
692 | | T35 = T33 * T34; |
693 | | T8I = T33 * T37; |
694 | | T36 = W[11]; |
695 | | T38 = FMA(T36, T37, T35); |
696 | | T8J = FNMS(T36, T34, T8I); |
697 | | } |
698 | | { |
699 | | E T3n, T3q, T3o, T8R, T3m, T3p; |
700 | | T3n = cr[WS(rs, 22)]; |
701 | | T3q = ci[WS(rs, 22)]; |
702 | | T3m = W[42]; |
703 | | T3o = T3m * T3n; |
704 | | T8R = T3m * T3q; |
705 | | T3p = W[43]; |
706 | | T3r = FMA(T3p, T3q, T3o); |
707 | | T8S = FNMS(T3p, T3n, T8R); |
708 | | } |
709 | | { |
710 | | E T3a, T3d, T3b, T8K, T39, T3c; |
711 | | T3a = cr[WS(rs, 38)]; |
712 | | T3d = ci[WS(rs, 38)]; |
713 | | T39 = W[74]; |
714 | | T3b = T39 * T3a; |
715 | | T8K = T39 * T3d; |
716 | | T3c = W[75]; |
717 | | T3e = FMA(T3c, T3d, T3b); |
718 | | T8L = FNMS(T3c, T3a, T8K); |
719 | | } |
720 | | { |
721 | | E T3h, T3k, T3i, T8P, T3g, T3j; |
722 | | T3h = cr[WS(rs, 54)]; |
723 | | T3k = ci[WS(rs, 54)]; |
724 | | T3g = W[106]; |
725 | | T3i = T3g * T3h; |
726 | | T8P = T3g * T3k; |
727 | | T3j = W[107]; |
728 | | T3l = FMA(T3j, T3k, T3i); |
729 | | T8Q = FNMS(T3j, T3h, T8P); |
730 | | } |
731 | | { |
732 | | E T3f, T3s, Tfk, Tfl; |
733 | | T3f = T38 + T3e; |
734 | | T3s = T3l + T3r; |
735 | | T3t = T3f + T3s; |
736 | | Tfh = T3f - T3s; |
737 | | Tfk = T8Q + T8S; |
738 | | Tfl = T8J + T8L; |
739 | | Tfm = Tfk - Tfl; |
740 | | ThC = Tfl + Tfk; |
741 | | { |
742 | | E T8N, T94, T8U, T95; |
743 | | { |
744 | | E T8H, T8M, T8O, T8T; |
745 | | T8H = T38 - T3e; |
746 | | T8M = T8J - T8L; |
747 | | T8N = T8H - T8M; |
748 | | T94 = T8H + T8M; |
749 | | T8O = T3l - T3r; |
750 | | T8T = T8Q - T8S; |
751 | | T8U = T8O + T8T; |
752 | | T95 = T8T - T8O; |
753 | | } |
754 | | T8V = T8N + T8U; |
755 | | TcY = T8U - T8N; |
756 | | T96 = T94 + T95; |
757 | | TcV = T94 - T95; |
758 | | } |
759 | | } |
760 | | } |
761 | | { |
762 | | E T43, T9k, T4m, T9t, T49, T9m, T4g, T9r; |
763 | | { |
764 | | E T3Z, T42, T40, T9j, T3Y, T41; |
765 | | T3Z = cr[WS(rs, 9)]; |
766 | | T42 = ci[WS(rs, 9)]; |
767 | | T3Y = W[16]; |
768 | | T40 = T3Y * T3Z; |
769 | | T9j = T3Y * T42; |
770 | | T41 = W[17]; |
771 | | T43 = FMA(T41, T42, T40); |
772 | | T9k = FNMS(T41, T3Z, T9j); |
773 | | } |
774 | | { |
775 | | E T4i, T4l, T4j, T9s, T4h, T4k; |
776 | | T4i = cr[WS(rs, 25)]; |
777 | | T4l = ci[WS(rs, 25)]; |
778 | | T4h = W[48]; |
779 | | T4j = T4h * T4i; |
780 | | T9s = T4h * T4l; |
781 | | T4k = W[49]; |
782 | | T4m = FMA(T4k, T4l, T4j); |
783 | | T9t = FNMS(T4k, T4i, T9s); |
784 | | } |
785 | | { |
786 | | E T45, T48, T46, T9l, T44, T47; |
787 | | T45 = cr[WS(rs, 41)]; |
788 | | T48 = ci[WS(rs, 41)]; |
789 | | T44 = W[80]; |
790 | | T46 = T44 * T45; |
791 | | T9l = T44 * T48; |
792 | | T47 = W[81]; |
793 | | T49 = FMA(T47, T48, T46); |
794 | | T9m = FNMS(T47, T45, T9l); |
795 | | } |
796 | | { |
797 | | E T4c, T4f, T4d, T9q, T4b, T4e; |
798 | | T4c = cr[WS(rs, 57)]; |
799 | | T4f = ci[WS(rs, 57)]; |
800 | | T4b = W[112]; |
801 | | T4d = T4b * T4c; |
802 | | T9q = T4b * T4f; |
803 | | T4e = W[113]; |
804 | | T4g = FMA(T4e, T4f, T4d); |
805 | | T9r = FNMS(T4e, T4c, T9q); |
806 | | } |
807 | | { |
808 | | E T4a, T4n, Tfs, Tft; |
809 | | T4a = T43 + T49; |
810 | | T4n = T4g + T4m; |
811 | | T4o = T4a + T4n; |
812 | | TfL = T4a - T4n; |
813 | | Tfs = T9r + T9t; |
814 | | Tft = T9k + T9m; |
815 | | Tfu = Tfs - Tft; |
816 | | ThN = Tft + Tfs; |
817 | | { |
818 | | E T9o, Tab, T9v, Tac; |
819 | | { |
820 | | E T9i, T9n, T9p, T9u; |
821 | | T9i = T43 - T49; |
822 | | T9n = T9k - T9m; |
823 | | T9o = T9i - T9n; |
824 | | Tab = T9i + T9n; |
825 | | T9p = T4g - T4m; |
826 | | T9u = T9r - T9t; |
827 | | T9v = T9p + T9u; |
828 | | Tac = T9u - T9p; |
829 | | } |
830 | | T9w = T9o + T9v; |
831 | | Tdf = T9v - T9o; |
832 | | Tad = Tab + Tac; |
833 | | Td4 = Tab - Tac; |
834 | | } |
835 | | } |
836 | | } |
837 | | { |
838 | | E T5Q, Tat, T69, TaC, T5W, Tav, T63, TaA; |
839 | | { |
840 | | E T5M, T5P, T5N, Tas, T5L, T5O; |
841 | | T5M = cr[WS(rs, 7)]; |
842 | | T5P = ci[WS(rs, 7)]; |
843 | | T5L = W[12]; |
844 | | T5N = T5L * T5M; |
845 | | Tas = T5L * T5P; |
846 | | T5O = W[13]; |
847 | | T5Q = FMA(T5O, T5P, T5N); |
848 | | Tat = FNMS(T5O, T5M, Tas); |
849 | | } |
850 | | { |
851 | | E T65, T68, T66, TaB, T64, T67; |
852 | | T65 = cr[WS(rs, 23)]; |
853 | | T68 = ci[WS(rs, 23)]; |
854 | | T64 = W[44]; |
855 | | T66 = T64 * T65; |
856 | | TaB = T64 * T68; |
857 | | T67 = W[45]; |
858 | | T69 = FMA(T67, T68, T66); |
859 | | TaC = FNMS(T67, T65, TaB); |
860 | | } |
861 | | { |
862 | | E T5S, T5V, T5T, Tau, T5R, T5U; |
863 | | T5S = cr[WS(rs, 39)]; |
864 | | T5V = ci[WS(rs, 39)]; |
865 | | T5R = W[76]; |
866 | | T5T = T5R * T5S; |
867 | | Tau = T5R * T5V; |
868 | | T5U = W[77]; |
869 | | T5W = FMA(T5U, T5V, T5T); |
870 | | Tav = FNMS(T5U, T5S, Tau); |
871 | | } |
872 | | { |
873 | | E T5Z, T62, T60, Taz, T5Y, T61; |
874 | | T5Z = cr[WS(rs, 55)]; |
875 | | T62 = ci[WS(rs, 55)]; |
876 | | T5Y = W[108]; |
877 | | T60 = T5Y * T5Z; |
878 | | Taz = T5Y * T62; |
879 | | T61 = W[109]; |
880 | | T63 = FMA(T61, T62, T60); |
881 | | TaA = FNMS(T61, T5Z, Taz); |
882 | | } |
883 | | { |
884 | | E T5X, T6a, TfT, TfU; |
885 | | T5X = T5Q + T5W; |
886 | | T6a = T63 + T69; |
887 | | T6b = T5X + T6a; |
888 | | Tg9 = T6a - T5X; |
889 | | TfT = TaA + TaC; |
890 | | TfU = Tat + Tav; |
891 | | TfV = TfT - TfU; |
892 | | ThY = TfU + TfT; |
893 | | { |
894 | | E Tax, Tbl, TaE, Tbk; |
895 | | { |
896 | | E Tar, Taw, Tay, TaD; |
897 | | Tar = T5Q - T5W; |
898 | | Taw = Tat - Tav; |
899 | | Tax = Tar - Taw; |
900 | | Tbl = Tar + Taw; |
901 | | Tay = T63 - T69; |
902 | | TaD = TaA - TaC; |
903 | | TaE = Tay + TaD; |
904 | | Tbk = Tay - TaD; |
905 | | } |
906 | | TaF = Tax + TaE; |
907 | | Tdy = TaE - Tax; |
908 | | Tbm = Tbk - Tbl; |
909 | | Tdn = Tbl + Tbk; |
910 | | } |
911 | | } |
912 | | } |
913 | | { |
914 | | E T4v, T9G, T4O, T9C, T4B, T9I, T4I, T9A; |
915 | | { |
916 | | E T4r, T4u, T4s, T9F, T4q, T4t; |
917 | | T4r = cr[WS(rs, 5)]; |
918 | | T4u = ci[WS(rs, 5)]; |
919 | | T4q = W[8]; |
920 | | T4s = T4q * T4r; |
921 | | T9F = T4q * T4u; |
922 | | T4t = W[9]; |
923 | | T4v = FMA(T4t, T4u, T4s); |
924 | | T9G = FNMS(T4t, T4r, T9F); |
925 | | } |
926 | | { |
927 | | E T4K, T4N, T4L, T9B, T4J, T4M; |
928 | | T4K = cr[WS(rs, 53)]; |
929 | | T4N = ci[WS(rs, 53)]; |
930 | | T4J = W[104]; |
931 | | T4L = T4J * T4K; |
932 | | T9B = T4J * T4N; |
933 | | T4M = W[105]; |
934 | | T4O = FMA(T4M, T4N, T4L); |
935 | | T9C = FNMS(T4M, T4K, T9B); |
936 | | } |
937 | | { |
938 | | E T4x, T4A, T4y, T9H, T4w, T4z; |
939 | | T4x = cr[WS(rs, 37)]; |
940 | | T4A = ci[WS(rs, 37)]; |
941 | | T4w = W[72]; |
942 | | T4y = T4w * T4x; |
943 | | T9H = T4w * T4A; |
944 | | T4z = W[73]; |
945 | | T4B = FMA(T4z, T4A, T4y); |
946 | | T9I = FNMS(T4z, T4x, T9H); |
947 | | } |
948 | | { |
949 | | E T4E, T4H, T4F, T9z, T4D, T4G; |
950 | | T4E = cr[WS(rs, 21)]; |
951 | | T4H = ci[WS(rs, 21)]; |
952 | | T4D = W[40]; |
953 | | T4F = T4D * T4E; |
954 | | T9z = T4D * T4H; |
955 | | T4G = W[41]; |
956 | | T4I = FMA(T4G, T4H, T4F); |
957 | | T9A = FNMS(T4G, T4E, T9z); |
958 | | } |
959 | | { |
960 | | E T4C, T4P, Tfw, Tfx, Tfy, Tfz; |
961 | | T4C = T4v + T4B; |
962 | | T4P = T4I + T4O; |
963 | | Tfw = T4C - T4P; |
964 | | Tfx = T9G + T9I; |
965 | | Tfy = T9A + T9C; |
966 | | Tfz = Tfx - Tfy; |
967 | | T4Q = T4C + T4P; |
968 | | ThJ = Tfx + Tfy; |
969 | | TfA = Tfw - Tfz; |
970 | | TfN = Tfw + Tfz; |
971 | | } |
972 | | { |
973 | | E T9E, Td6, T9L, Td7; |
974 | | { |
975 | | E T9y, T9D, T9J, T9K; |
976 | | T9y = T4v - T4B; |
977 | | T9D = T9A - T9C; |
978 | | T9E = T9y - T9D; |
979 | | Td6 = T9y + T9D; |
980 | | T9J = T9G - T9I; |
981 | | T9K = T4I - T4O; |
982 | | T9L = T9J + T9K; |
983 | | Td7 = T9J - T9K; |
984 | | } |
985 | | T9M = FNMS(KP414213562, T9L, T9E); |
986 | | Tdh = FNMS(KP414213562, Td6, Td7); |
987 | | Taf = FMA(KP414213562, T9E, T9L); |
988 | | Td8 = FMA(KP414213562, Td7, Td6); |
989 | | } |
990 | | } |
991 | | { |
992 | | E T4W, T9V, T5f, T9R, T52, T9X, T59, T9P; |
993 | | { |
994 | | E T4S, T4V, T4T, T9U, T4R, T4U; |
995 | | T4S = cr[WS(rs, 61)]; |
996 | | T4V = ci[WS(rs, 61)]; |
997 | | T4R = W[120]; |
998 | | T4T = T4R * T4S; |
999 | | T9U = T4R * T4V; |
1000 | | T4U = W[121]; |
1001 | | T4W = FMA(T4U, T4V, T4T); |
1002 | | T9V = FNMS(T4U, T4S, T9U); |
1003 | | } |
1004 | | { |
1005 | | E T5b, T5e, T5c, T9Q, T5a, T5d; |
1006 | | T5b = cr[WS(rs, 45)]; |
1007 | | T5e = ci[WS(rs, 45)]; |
1008 | | T5a = W[88]; |
1009 | | T5c = T5a * T5b; |
1010 | | T9Q = T5a * T5e; |
1011 | | T5d = W[89]; |
1012 | | T5f = FMA(T5d, T5e, T5c); |
1013 | | T9R = FNMS(T5d, T5b, T9Q); |
1014 | | } |
1015 | | { |
1016 | | E T4Y, T51, T4Z, T9W, T4X, T50; |
1017 | | T4Y = cr[WS(rs, 29)]; |
1018 | | T51 = ci[WS(rs, 29)]; |
1019 | | T4X = W[56]; |
1020 | | T4Z = T4X * T4Y; |
1021 | | T9W = T4X * T51; |
1022 | | T50 = W[57]; |
1023 | | T52 = FMA(T50, T51, T4Z); |
1024 | | T9X = FNMS(T50, T4Y, T9W); |
1025 | | } |
1026 | | { |
1027 | | E T55, T58, T56, T9O, T54, T57; |
1028 | | T55 = cr[WS(rs, 13)]; |
1029 | | T58 = ci[WS(rs, 13)]; |
1030 | | T54 = W[24]; |
1031 | | T56 = T54 * T55; |
1032 | | T9O = T54 * T58; |
1033 | | T57 = W[25]; |
1034 | | T59 = FMA(T57, T58, T56); |
1035 | | T9P = FNMS(T57, T55, T9O); |
1036 | | } |
1037 | | { |
1038 | | E T53, T5g, TfB, TfC, TfD, TfE; |
1039 | | T53 = T4W + T52; |
1040 | | T5g = T59 + T5f; |
1041 | | TfB = T53 - T5g; |
1042 | | TfC = T9V + T9X; |
1043 | | TfD = T9P + T9R; |
1044 | | TfE = TfC - TfD; |
1045 | | T5h = T53 + T5g; |
1046 | | ThI = TfC + TfD; |
1047 | | TfF = TfB + TfE; |
1048 | | TfO = TfE - TfB; |
1049 | | } |
1050 | | { |
1051 | | E T9T, Td9, Ta0, Tda; |
1052 | | { |
1053 | | E T9N, T9S, T9Y, T9Z; |
1054 | | T9N = T4W - T52; |
1055 | | T9S = T9P - T9R; |
1056 | | T9T = T9N - T9S; |
1057 | | Td9 = T9N + T9S; |
1058 | | T9Y = T9V - T9X; |
1059 | | T9Z = T59 - T5f; |
1060 | | Ta0 = T9Y + T9Z; |
1061 | | Tda = T9Y - T9Z; |
1062 | | } |
1063 | | Ta1 = FMA(KP414213562, Ta0, T9T); |
1064 | | Tdi = FMA(KP414213562, Td9, Tda); |
1065 | | Tag = FNMS(KP414213562, T9T, Ta0); |
1066 | | Tdb = FNMS(KP414213562, Tda, Td9); |
1067 | | } |
1068 | | } |
1069 | | { |
1070 | | E T6i, TaQ, T6B, TaL, T6o, TaS, T6v, TaJ; |
1071 | | { |
1072 | | E T6e, T6h, T6f, TaP, T6d, T6g; |
1073 | | T6e = cr[WS(rs, 3)]; |
1074 | | T6h = ci[WS(rs, 3)]; |
1075 | | T6d = W[4]; |
1076 | | T6f = T6d * T6e; |
1077 | | TaP = T6d * T6h; |
1078 | | T6g = W[5]; |
1079 | | T6i = FMA(T6g, T6h, T6f); |
1080 | | TaQ = FNMS(T6g, T6e, TaP); |
1081 | | } |
1082 | | { |
1083 | | E T6x, T6A, T6y, TaK, T6w, T6z; |
1084 | | T6x = cr[WS(rs, 51)]; |
1085 | | T6A = ci[WS(rs, 51)]; |
1086 | | T6w = W[100]; |
1087 | | T6y = T6w * T6x; |
1088 | | TaK = T6w * T6A; |
1089 | | T6z = W[101]; |
1090 | | T6B = FMA(T6z, T6A, T6y); |
1091 | | TaL = FNMS(T6z, T6x, TaK); |
1092 | | } |
1093 | | { |
1094 | | E T6k, T6n, T6l, TaR, T6j, T6m; |
1095 | | T6k = cr[WS(rs, 35)]; |
1096 | | T6n = ci[WS(rs, 35)]; |
1097 | | T6j = W[68]; |
1098 | | T6l = T6j * T6k; |
1099 | | TaR = T6j * T6n; |
1100 | | T6m = W[69]; |
1101 | | T6o = FMA(T6m, T6n, T6l); |
1102 | | TaS = FNMS(T6m, T6k, TaR); |
1103 | | } |
1104 | | { |
1105 | | E T6r, T6u, T6s, TaI, T6q, T6t; |
1106 | | T6r = cr[WS(rs, 19)]; |
1107 | | T6u = ci[WS(rs, 19)]; |
1108 | | T6q = W[36]; |
1109 | | T6s = T6q * T6r; |
1110 | | TaI = T6q * T6u; |
1111 | | T6t = W[37]; |
1112 | | T6v = FMA(T6t, T6u, T6s); |
1113 | | TaJ = FNMS(T6t, T6r, TaI); |
1114 | | } |
1115 | | { |
1116 | | E T6p, T6C, TfX, TfY, TfZ, Tg0; |
1117 | | T6p = T6i + T6o; |
1118 | | T6C = T6v + T6B; |
1119 | | TfX = T6p - T6C; |
1120 | | TfY = TaQ + TaS; |
1121 | | TfZ = TaJ + TaL; |
1122 | | Tg0 = TfY - TfZ; |
1123 | | T6D = T6p + T6C; |
1124 | | ThU = TfY + TfZ; |
1125 | | Tg1 = TfX - Tg0; |
1126 | | Tgf = TfX + Tg0; |
1127 | | } |
1128 | | { |
1129 | | E TaN, Tdp, TaU, Tdq; |
1130 | | { |
1131 | | E TaH, TaM, TaO, TaT; |
1132 | | TaH = T6i - T6o; |
1133 | | TaM = TaJ - TaL; |
1134 | | TaN = TaH - TaM; |
1135 | | Tdp = TaH + TaM; |
1136 | | TaO = T6B - T6v; |
1137 | | TaT = TaQ - TaS; |
1138 | | TaU = TaO - TaT; |
1139 | | Tdq = TaT + TaO; |
1140 | | } |
1141 | | TaV = FMA(KP414213562, TaU, TaN); |
1142 | | TdA = FNMS(KP414213562, Tdp, Tdq); |
1143 | | Tbo = FNMS(KP414213562, TaN, TaU); |
1144 | | Tdr = FMA(KP414213562, Tdq, Tdp); |
1145 | | } |
1146 | | } |
1147 | | { |
1148 | | E T6J, Tb5, T72, Tb0, T6P, Tb7, T6W, TaY; |
1149 | | { |
1150 | | E T6F, T6I, T6G, Tb4, T6E, T6H; |
1151 | | T6F = cr[WS(rs, 59)]; |
1152 | | T6I = ci[WS(rs, 59)]; |
1153 | | T6E = W[116]; |
1154 | | T6G = T6E * T6F; |
1155 | | Tb4 = T6E * T6I; |
1156 | | T6H = W[117]; |
1157 | | T6J = FMA(T6H, T6I, T6G); |
1158 | | Tb5 = FNMS(T6H, T6F, Tb4); |
1159 | | } |
1160 | | { |
1161 | | E T6Y, T71, T6Z, TaZ, T6X, T70; |
1162 | | T6Y = cr[WS(rs, 43)]; |
1163 | | T71 = ci[WS(rs, 43)]; |
1164 | | T6X = W[84]; |
1165 | | T6Z = T6X * T6Y; |
1166 | | TaZ = T6X * T71; |
1167 | | T70 = W[85]; |
1168 | | T72 = FMA(T70, T71, T6Z); |
1169 | | Tb0 = FNMS(T70, T6Y, TaZ); |
1170 | | } |
1171 | | { |
1172 | | E T6L, T6O, T6M, Tb6, T6K, T6N; |
1173 | | T6L = cr[WS(rs, 27)]; |
1174 | | T6O = ci[WS(rs, 27)]; |
1175 | | T6K = W[52]; |
1176 | | T6M = T6K * T6L; |
1177 | | Tb6 = T6K * T6O; |
1178 | | T6N = W[53]; |
1179 | | T6P = FMA(T6N, T6O, T6M); |
1180 | | Tb7 = FNMS(T6N, T6L, Tb6); |
1181 | | } |
1182 | | { |
1183 | | E T6S, T6V, T6T, TaX, T6R, T6U; |
1184 | | T6S = cr[WS(rs, 11)]; |
1185 | | T6V = ci[WS(rs, 11)]; |
1186 | | T6R = W[20]; |
1187 | | T6T = T6R * T6S; |
1188 | | TaX = T6R * T6V; |
1189 | | T6U = W[21]; |
1190 | | T6W = FMA(T6U, T6V, T6T); |
1191 | | TaY = FNMS(T6U, T6S, TaX); |
1192 | | } |
1193 | | { |
1194 | | E T6Q, T73, Tg2, Tg3, Tg4, Tg5; |
1195 | | T6Q = T6J + T6P; |
1196 | | T73 = T6W + T72; |
1197 | | Tg2 = T6Q - T73; |
1198 | | Tg3 = Tb5 + Tb7; |
1199 | | Tg4 = TaY + Tb0; |
1200 | | Tg5 = Tg3 - Tg4; |
1201 | | T74 = T6Q + T73; |
1202 | | ThT = Tg3 + Tg4; |
1203 | | Tg6 = Tg2 + Tg5; |
1204 | | Tge = Tg2 - Tg5; |
1205 | | } |
1206 | | { |
1207 | | E Tb2, Tds, Tb9, Tdt; |
1208 | | { |
1209 | | E TaW, Tb1, Tb3, Tb8; |
1210 | | TaW = T6J - T6P; |
1211 | | Tb1 = TaY - Tb0; |
1212 | | Tb2 = TaW - Tb1; |
1213 | | Tds = TaW + Tb1; |
1214 | | Tb3 = T72 - T6W; |
1215 | | Tb8 = Tb5 - Tb7; |
1216 | | Tb9 = Tb3 - Tb8; |
1217 | | Tdt = Tb8 + Tb3; |
1218 | | } |
1219 | | Tba = FNMS(KP414213562, Tb9, Tb2); |
1220 | | TdB = FMA(KP414213562, Tds, Tdt); |
1221 | | Tbp = FMA(KP414213562, Tb2, Tb9); |
1222 | | Tdu = FNMS(KP414213562, Tdt, Tds); |
1223 | | } |
1224 | | } |
1225 | | { |
1226 | | E T1I, Tio, T3v, Tj1, TiX, Tj2, Tir, TiN, T76, TiJ, TiC, TiG, T5j, TiK, Tix; |
1227 | | E TiF; |
1228 | | { |
1229 | | E TO, T1H, Tip, Tiq; |
1230 | | TO = Tm + TN; |
1231 | | T1H = T1f + T1G; |
1232 | | T1I = TO + T1H; |
1233 | | Tio = TO - T1H; |
1234 | | { |
1235 | | E T2B, T3u, TiO, TiW; |
1236 | | T2B = T29 + T2A; |
1237 | | T3u = T32 + T3t; |
1238 | | T3v = T2B + T3u; |
1239 | | Tj1 = T2B - T3u; |
1240 | | TiO = Ths + Thr; |
1241 | | TiW = TiP + TiV; |
1242 | | TiX = TiO + TiW; |
1243 | | Tj2 = TiW - TiO; |
1244 | | } |
1245 | | Tip = ThB + ThC; |
1246 | | Tiq = Thw + Thx; |
1247 | | Tir = Tip - Tiq; |
1248 | | TiN = Tiq + Tip; |
1249 | | { |
1250 | | E T6c, T75, Tiy, Tiz, TiA, TiB; |
1251 | | T6c = T5K + T6b; |
1252 | | T75 = T6D + T74; |
1253 | | Tiy = T6c - T75; |
1254 | | Tiz = ThX + ThY; |
1255 | | TiA = ThU + ThT; |
1256 | | TiB = Tiz - TiA; |
1257 | | T76 = T6c + T75; |
1258 | | TiJ = Tiz + TiA; |
1259 | | TiC = Tiy - TiB; |
1260 | | TiG = Tiy + TiB; |
1261 | | } |
1262 | | { |
1263 | | E T4p, T5i, Tit, Tiu, Tiv, Tiw; |
1264 | | T4p = T3X + T4o; |
1265 | | T5i = T4Q + T5h; |
1266 | | Tit = T4p - T5i; |
1267 | | Tiu = ThM + ThN; |
1268 | | Tiv = ThJ + ThI; |
1269 | | Tiw = Tiu - Tiv; |
1270 | | T5j = T4p + T5i; |
1271 | | TiK = Tiu + Tiv; |
1272 | | Tix = Tit + Tiw; |
1273 | | TiF = Tit - Tiw; |
1274 | | } |
1275 | | } |
1276 | | { |
1277 | | E T3w, T77, Tj3, Tj4; |
1278 | | T3w = T1I + T3v; |
1279 | | T77 = T5j + T76; |
1280 | | ci[WS(rs, 31)] = T3w - T77; |
1281 | | cr[0] = T3w + T77; |
1282 | | Tj3 = Tj1 + Tj2; |
1283 | | Tj4 = TiC - Tix; |
1284 | | cr[WS(rs, 56)] = FMS(KP707106781, Tj4, Tj3); |
1285 | | ci[WS(rs, 39)] = FMA(KP707106781, Tj4, Tj3); |
1286 | | } |
1287 | | { |
1288 | | E Tj5, Tj6, Tis, TiD; |
1289 | | Tj5 = Tj2 - Tj1; |
1290 | | Tj6 = TiG - TiF; |
1291 | | cr[WS(rs, 40)] = FMS(KP707106781, Tj6, Tj5); |
1292 | | ci[WS(rs, 55)] = FMA(KP707106781, Tj6, Tj5); |
1293 | | Tis = Tio - Tir; |
1294 | | TiD = Tix + TiC; |
1295 | | ci[WS(rs, 23)] = FNMS(KP707106781, TiD, Tis); |
1296 | | cr[WS(rs, 8)] = FMA(KP707106781, TiD, Tis); |
1297 | | } |
1298 | | { |
1299 | | E TiE, TiH, TiM, TiY; |
1300 | | TiE = Tio + Tir; |
1301 | | TiH = TiF + TiG; |
1302 | | cr[WS(rs, 24)] = FNMS(KP707106781, TiH, TiE); |
1303 | | ci[WS(rs, 7)] = FMA(KP707106781, TiH, TiE); |
1304 | | TiM = TiK + TiJ; |
1305 | | TiY = TiN + TiX; |
1306 | | cr[WS(rs, 32)] = TiM - TiY; |
1307 | | ci[WS(rs, 63)] = TiM + TiY; |
1308 | | } |
1309 | | { |
1310 | | E TiZ, Tj0, TiI, TiL; |
1311 | | TiZ = T76 - T5j; |
1312 | | Tj0 = TiX - TiN; |
1313 | | cr[WS(rs, 48)] = TiZ - Tj0; |
1314 | | ci[WS(rs, 47)] = TiZ + Tj0; |
1315 | | TiI = T1I - T3v; |
1316 | | TiL = TiJ - TiK; |
1317 | | cr[WS(rs, 16)] = TiI - TiL; |
1318 | | ci[WS(rs, 15)] = TiI + TiL; |
1319 | | } |
1320 | | } |
1321 | | { |
1322 | | E T99, Tk2, TbB, TjW, Taj, TbL, Tbv, TbF, Tce, Tcy, Tci, Tcu, Tc7, Tcx, Tch; |
1323 | | E Tcr, TbZ, Tkg, Tcn, Tka, Tbs, TbM, Tbw, TbI, T80, Tk9, Tkf, Tby, TbS, TjV; |
1324 | | E Tk1, Tck; |
1325 | | { |
1326 | | E T8z, Tbz, T98, TbA; |
1327 | | { |
1328 | | E T8n, T8y, T8W, T97; |
1329 | | T8n = FMA(KP707106781, T8m, T87); |
1330 | | T8y = FMA(KP707106781, T8x, T8u); |
1331 | | T8z = FNMS(KP198912367, T8y, T8n); |
1332 | | Tbz = FMA(KP198912367, T8n, T8y); |
1333 | | T8W = FMA(KP707106781, T8V, T8G); |
1334 | | T97 = FMA(KP707106781, T96, T93); |
1335 | | T98 = FMA(KP198912367, T97, T8W); |
1336 | | TbA = FNMS(KP198912367, T8W, T97); |
1337 | | } |
1338 | | T99 = T8z + T98; |
1339 | | Tk2 = T98 - T8z; |
1340 | | TbB = Tbz - TbA; |
1341 | | TjW = Tbz + TbA; |
1342 | | } |
1343 | | { |
1344 | | E Ta3, TbD, Tai, TbE; |
1345 | | { |
1346 | | E T9x, Ta2, Tae, Tah; |
1347 | | T9x = FMA(KP707106781, T9w, T9h); |
1348 | | Ta2 = T9M + Ta1; |
1349 | | Ta3 = FMA(KP923879532, Ta2, T9x); |
1350 | | TbD = FNMS(KP923879532, Ta2, T9x); |
1351 | | Tae = FMA(KP707106781, Tad, Taa); |
1352 | | Tah = Taf + Tag; |
1353 | | Tai = FMA(KP923879532, Tah, Tae); |
1354 | | TbE = FNMS(KP923879532, Tah, Tae); |
1355 | | } |
1356 | | Taj = FNMS(KP098491403, Tai, Ta3); |
1357 | | TbL = FNMS(KP820678790, TbD, TbE); |
1358 | | Tbv = FMA(KP098491403, Ta3, Tai); |
1359 | | TbF = FMA(KP820678790, TbE, TbD); |
1360 | | } |
1361 | | { |
1362 | | E Tca, Tcs, Tcd, Tct; |
1363 | | { |
1364 | | E Tc8, Tc9, Tcb, Tcc; |
1365 | | Tc8 = FNMS(KP707106781, TaF, Taq); |
1366 | | Tc9 = Tbp - Tbo; |
1367 | | Tca = FNMS(KP923879532, Tc9, Tc8); |
1368 | | Tcs = FMA(KP923879532, Tc9, Tc8); |
1369 | | Tcb = FNMS(KP707106781, Tbm, Tbj); |
1370 | | Tcc = Tba - TaV; |
1371 | | Tcd = FMA(KP923879532, Tcc, Tcb); |
1372 | | Tct = FNMS(KP923879532, Tcc, Tcb); |
1373 | | } |
1374 | | Tce = FNMS(KP534511135, Tcd, Tca); |
1375 | | Tcy = FNMS(KP303346683, Tcs, Tct); |
1376 | | Tci = FMA(KP534511135, Tca, Tcd); |
1377 | | Tcu = FMA(KP303346683, Tct, Tcs); |
1378 | | } |
1379 | | { |
1380 | | E Tc3, Tcp, Tc6, Tcq; |
1381 | | { |
1382 | | E Tc1, Tc2, Tc4, Tc5; |
1383 | | Tc1 = FNMS(KP707106781, T9w, T9h); |
1384 | | Tc2 = Taf - Tag; |
1385 | | Tc3 = FNMS(KP923879532, Tc2, Tc1); |
1386 | | Tcp = FMA(KP923879532, Tc2, Tc1); |
1387 | | Tc4 = FNMS(KP707106781, Tad, Taa); |
1388 | | Tc5 = Ta1 - T9M; |
1389 | | Tc6 = FNMS(KP923879532, Tc5, Tc4); |
1390 | | Tcq = FMA(KP923879532, Tc5, Tc4); |
1391 | | } |
1392 | | Tc7 = FNMS(KP534511135, Tc6, Tc3); |
1393 | | Tcx = FNMS(KP303346683, Tcp, Tcq); |
1394 | | Tch = FMA(KP534511135, Tc3, Tc6); |
1395 | | Tcr = FMA(KP303346683, Tcq, Tcp); |
1396 | | } |
1397 | | { |
1398 | | E TbV, Tcm, TbY, Tcl; |
1399 | | { |
1400 | | E TbT, TbU, TbW, TbX; |
1401 | | TbT = FNMS(KP707106781, T96, T93); |
1402 | | TbU = FNMS(KP707106781, T8V, T8G); |
1403 | | TbV = FMA(KP668178637, TbU, TbT); |
1404 | | Tcm = FNMS(KP668178637, TbT, TbU); |
1405 | | TbW = FNMS(KP707106781, T8x, T8u); |
1406 | | TbX = FNMS(KP707106781, T8m, T87); |
1407 | | TbY = FNMS(KP668178637, TbX, TbW); |
1408 | | Tcl = FMA(KP668178637, TbW, TbX); |
1409 | | } |
1410 | | TbZ = TbV - TbY; |
1411 | | Tkg = Tcl - Tcm; |
1412 | | Tcn = Tcl + Tcm; |
1413 | | Tka = TbY + TbV; |
1414 | | } |
1415 | | { |
1416 | | E Tbc, TbG, Tbr, TbH; |
1417 | | { |
1418 | | E TaG, Tbb, Tbn, Tbq; |
1419 | | TaG = FMA(KP707106781, TaF, Taq); |
1420 | | Tbb = TaV + Tba; |
1421 | | Tbc = FMA(KP923879532, Tbb, TaG); |
1422 | | TbG = FNMS(KP923879532, Tbb, TaG); |
1423 | | Tbn = FMA(KP707106781, Tbm, Tbj); |
1424 | | Tbq = Tbo + Tbp; |
1425 | | Tbr = FMA(KP923879532, Tbq, Tbn); |
1426 | | TbH = FNMS(KP923879532, Tbq, Tbn); |
1427 | | } |
1428 | | Tbs = FNMS(KP098491403, Tbr, Tbc); |
1429 | | TbM = FNMS(KP820678790, TbG, TbH); |
1430 | | Tbw = FMA(KP098491403, Tbc, Tbr); |
1431 | | TbI = FMA(KP820678790, TbH, TbG); |
1432 | | } |
1433 | | { |
1434 | | E T7u, TbO, TjT, Tk7, T7Z, Tk8, TbR, TjU, T7t, TjS; |
1435 | | T7t = T7l + T7s; |
1436 | | T7u = FMA(KP707106781, T7t, T7e); |
1437 | | TbO = FNMS(KP707106781, T7t, T7e); |
1438 | | TjS = TcB - TcC; |
1439 | | TjT = FMA(KP707106781, TjS, TjR); |
1440 | | Tk7 = FNMS(KP707106781, TjS, TjR); |
1441 | | { |
1442 | | E T7J, T7Y, TbP, TbQ; |
1443 | | T7J = FNMS(KP414213562, T7I, T7B); |
1444 | | T7Y = FMA(KP414213562, T7X, T7Q); |
1445 | | T7Z = T7J + T7Y; |
1446 | | Tk8 = T7Y - T7J; |
1447 | | TbP = FMA(KP414213562, T7B, T7I); |
1448 | | TbQ = FNMS(KP414213562, T7Q, T7X); |
1449 | | TbR = TbP - TbQ; |
1450 | | TjU = TbP + TbQ; |
1451 | | } |
1452 | | T80 = FMA(KP923879532, T7Z, T7u); |
1453 | | Tk9 = FMA(KP923879532, Tk8, Tk7); |
1454 | | Tkf = FNMS(KP923879532, Tk8, Tk7); |
1455 | | Tby = FNMS(KP923879532, T7Z, T7u); |
1456 | | TbS = FNMS(KP923879532, TbR, TbO); |
1457 | | TjV = FMA(KP923879532, TjU, TjT); |
1458 | | Tk1 = FNMS(KP923879532, TjU, TjT); |
1459 | | Tck = FMA(KP923879532, TbR, TbO); |
1460 | | } |
1461 | | { |
1462 | | E T9a, Tbt, TbK, TbN; |
1463 | | T9a = FMA(KP980785280, T99, T80); |
1464 | | Tbt = Taj + Tbs; |
1465 | | cr[WS(rs, 31)] = FNMS(KP995184726, Tbt, T9a); |
1466 | | ci[0] = FMA(KP995184726, Tbt, T9a); |
1467 | | TbK = FNMS(KP980785280, TbB, Tby); |
1468 | | TbN = TbL + TbM; |
1469 | | cr[WS(rs, 23)] = FMA(KP773010453, TbN, TbK); |
1470 | | ci[WS(rs, 8)] = FNMS(KP773010453, TbN, TbK); |
1471 | | } |
1472 | | { |
1473 | | E Tkb, Tkc, Tkj, Tkk; |
1474 | | Tkb = FMA(KP831469612, Tka, Tk9); |
1475 | | Tkc = Tcx - Tcy; |
1476 | | cr[WS(rs, 35)] = FMS(KP956940335, Tkc, Tkb); |
1477 | | ci[WS(rs, 60)] = FMA(KP956940335, Tkc, Tkb); |
1478 | | Tkj = FNMS(KP831469612, Tkg, Tkf); |
1479 | | Tkk = Tce - Tc7; |
1480 | | cr[WS(rs, 43)] = FMS(KP881921264, Tkk, Tkj); |
1481 | | ci[WS(rs, 52)] = FMA(KP881921264, Tkk, Tkj); |
1482 | | } |
1483 | | { |
1484 | | E Tbu, Tbx, TbC, TbJ; |
1485 | | Tbu = FNMS(KP980785280, T99, T80); |
1486 | | Tbx = Tbv + Tbw; |
1487 | | ci[WS(rs, 16)] = FNMS(KP995184726, Tbx, Tbu); |
1488 | | cr[WS(rs, 15)] = FMA(KP995184726, Tbx, Tbu); |
1489 | | TbC = FMA(KP980785280, TbB, Tby); |
1490 | | TbJ = TbF + TbI; |
1491 | | ci[WS(rs, 24)] = FNMS(KP773010453, TbJ, TbC); |
1492 | | cr[WS(rs, 7)] = FMA(KP773010453, TbJ, TbC); |
1493 | | } |
1494 | | { |
1495 | | E Tkd, Tke, Tkh, Tki; |
1496 | | Tkd = FNMS(KP831469612, Tka, Tk9); |
1497 | | Tke = Tcu - Tcr; |
1498 | | cr[WS(rs, 51)] = FMS(KP956940335, Tke, Tkd); |
1499 | | ci[WS(rs, 44)] = FMA(KP956940335, Tke, Tkd); |
1500 | | Tkh = FMA(KP831469612, Tkg, Tkf); |
1501 | | Tki = Tci - Tch; |
1502 | | cr[WS(rs, 59)] = FMS(KP881921264, Tki, Tkh); |
1503 | | ci[WS(rs, 36)] = FMA(KP881921264, Tki, Tkh); |
1504 | | } |
1505 | | { |
1506 | | E Tc0, Tcf, Tcw, Tcz; |
1507 | | Tc0 = FMA(KP831469612, TbZ, TbS); |
1508 | | Tcf = Tc7 + Tce; |
1509 | | cr[WS(rs, 27)] = FNMS(KP881921264, Tcf, Tc0); |
1510 | | ci[WS(rs, 4)] = FMA(KP881921264, Tcf, Tc0); |
1511 | | Tcw = FNMS(KP831469612, Tcn, Tck); |
1512 | | Tcz = Tcx + Tcy; |
1513 | | cr[WS(rs, 19)] = FMA(KP956940335, Tcz, Tcw); |
1514 | | ci[WS(rs, 12)] = FNMS(KP956940335, Tcz, Tcw); |
1515 | | } |
1516 | | { |
1517 | | E TjX, TjY, Tk5, Tk6; |
1518 | | TjX = FMA(KP980785280, TjW, TjV); |
1519 | | TjY = Tbw - Tbv; |
1520 | | cr[WS(rs, 63)] = FMS(KP995184726, TjY, TjX); |
1521 | | ci[WS(rs, 32)] = FMA(KP995184726, TjY, TjX); |
1522 | | Tk5 = FNMS(KP980785280, Tk2, Tk1); |
1523 | | Tk6 = TbI - TbF; |
1524 | | cr[WS(rs, 55)] = FMS(KP773010453, Tk6, Tk5); |
1525 | | ci[WS(rs, 40)] = FMA(KP773010453, Tk6, Tk5); |
1526 | | } |
1527 | | { |
1528 | | E Tcg, Tcj, Tco, Tcv; |
1529 | | Tcg = FNMS(KP831469612, TbZ, TbS); |
1530 | | Tcj = Tch + Tci; |
1531 | | ci[WS(rs, 20)] = FNMS(KP881921264, Tcj, Tcg); |
1532 | | cr[WS(rs, 11)] = FMA(KP881921264, Tcj, Tcg); |
1533 | | Tco = FMA(KP831469612, Tcn, Tck); |
1534 | | Tcv = Tcr + Tcu; |
1535 | | ci[WS(rs, 28)] = FNMS(KP956940335, Tcv, Tco); |
1536 | | cr[WS(rs, 3)] = FMA(KP956940335, Tcv, Tco); |
1537 | | } |
1538 | | { |
1539 | | E TjZ, Tk0, Tk3, Tk4; |
1540 | | TjZ = FNMS(KP980785280, TjW, TjV); |
1541 | | Tk0 = Tbs - Taj; |
1542 | | cr[WS(rs, 47)] = FMS(KP995184726, Tk0, TjZ); |
1543 | | ci[WS(rs, 48)] = FMA(KP995184726, Tk0, TjZ); |
1544 | | Tk3 = FMA(KP980785280, Tk2, Tk1); |
1545 | | Tk4 = TbL - TbM; |
1546 | | cr[WS(rs, 39)] = FMS(KP773010453, Tk4, Tk3); |
1547 | | ci[WS(rs, 56)] = FMA(KP773010453, Tk4, Tk3); |
1548 | | } |
1549 | | } |
1550 | | { |
1551 | | E Thu, Ti8, Tj9, Tjf, ThF, Tjg, Tib, Tja, ThR, Til, Ti6, Tif, Ti2, Tim, Ti5; |
1552 | | E Tii; |
1553 | | { |
1554 | | E Thq, Tht, Tj7, Tj8; |
1555 | | Thq = Tm - TN; |
1556 | | Tht = Thr - Ths; |
1557 | | Thu = Thq - Tht; |
1558 | | Ti8 = Thq + Tht; |
1559 | | Tj7 = T1f - T1G; |
1560 | | Tj8 = TiV - TiP; |
1561 | | Tj9 = Tj7 + Tj8; |
1562 | | Tjf = Tj8 - Tj7; |
1563 | | } |
1564 | | { |
1565 | | E Thz, Ti9, ThE, Tia; |
1566 | | { |
1567 | | E Thv, Thy, ThA, ThD; |
1568 | | Thv = T29 - T2A; |
1569 | | Thy = Thw - Thx; |
1570 | | Thz = Thv + Thy; |
1571 | | Ti9 = Thv - Thy; |
1572 | | ThA = T32 - T3t; |
1573 | | ThD = ThB - ThC; |
1574 | | ThE = ThA - ThD; |
1575 | | Tia = ThA + ThD; |
1576 | | } |
1577 | | ThF = Thz + ThE; |
1578 | | Tjg = Tia - Ti9; |
1579 | | Tib = Ti9 + Tia; |
1580 | | Tja = Thz - ThE; |
1581 | | } |
1582 | | { |
1583 | | E ThL, Tid, ThQ, Tie; |
1584 | | { |
1585 | | E ThH, ThK, ThO, ThP; |
1586 | | ThH = T3X - T4o; |
1587 | | ThK = ThI - ThJ; |
1588 | | ThL = ThH - ThK; |
1589 | | Tid = ThH + ThK; |
1590 | | ThO = ThM - ThN; |
1591 | | ThP = T4Q - T5h; |
1592 | | ThQ = ThO - ThP; |
1593 | | Tie = ThO + ThP; |
1594 | | } |
1595 | | ThR = FMA(KP414213562, ThQ, ThL); |
1596 | | Til = FMA(KP414213562, Tid, Tie); |
1597 | | Ti6 = FNMS(KP414213562, ThL, ThQ); |
1598 | | Tif = FNMS(KP414213562, Tie, Tid); |
1599 | | } |
1600 | | { |
1601 | | E ThW, Tig, Ti1, Tih; |
1602 | | { |
1603 | | E ThS, ThV, ThZ, Ti0; |
1604 | | ThS = T5K - T6b; |
1605 | | ThV = ThT - ThU; |
1606 | | ThW = ThS - ThV; |
1607 | | Tig = ThS + ThV; |
1608 | | ThZ = ThX - ThY; |
1609 | | Ti0 = T74 - T6D; |
1610 | | Ti1 = ThZ + Ti0; |
1611 | | Tih = Ti0 - ThZ; |
1612 | | } |
1613 | | Ti2 = FNMS(KP414213562, Ti1, ThW); |
1614 | | Tim = FMA(KP414213562, Tig, Tih); |
1615 | | Ti5 = FMA(KP414213562, ThW, Ti1); |
1616 | | Tii = FNMS(KP414213562, Tih, Tig); |
1617 | | } |
1618 | | { |
1619 | | E ThG, Ti3, Tjh, Tji; |
1620 | | ThG = FMA(KP707106781, ThF, Thu); |
1621 | | Ti3 = ThR + Ti2; |
1622 | | ci[WS(rs, 27)] = FNMS(KP923879532, Ti3, ThG); |
1623 | | cr[WS(rs, 4)] = FMA(KP923879532, Ti3, ThG); |
1624 | | Tjh = FMA(KP707106781, Tjg, Tjf); |
1625 | | Tji = Ti6 + Ti5; |
1626 | | cr[WS(rs, 36)] = FMS(KP923879532, Tji, Tjh); |
1627 | | ci[WS(rs, 59)] = FMA(KP923879532, Tji, Tjh); |
1628 | | } |
1629 | | { |
1630 | | E Tjj, Tjk, Ti4, Ti7; |
1631 | | Tjj = FNMS(KP707106781, Tjg, Tjf); |
1632 | | Tjk = Ti2 - ThR; |
1633 | | cr[WS(rs, 52)] = FMS(KP923879532, Tjk, Tjj); |
1634 | | ci[WS(rs, 43)] = FMA(KP923879532, Tjk, Tjj); |
1635 | | Ti4 = FNMS(KP707106781, ThF, Thu); |
1636 | | Ti7 = Ti5 - Ti6; |
1637 | | cr[WS(rs, 20)] = FNMS(KP923879532, Ti7, Ti4); |
1638 | | ci[WS(rs, 11)] = FMA(KP923879532, Ti7, Ti4); |
1639 | | } |
1640 | | { |
1641 | | E Tic, Tij, Tjb, Tjc; |
1642 | | Tic = FMA(KP707106781, Tib, Ti8); |
1643 | | Tij = Tif + Tii; |
1644 | | cr[WS(rs, 28)] = FNMS(KP923879532, Tij, Tic); |
1645 | | ci[WS(rs, 3)] = FMA(KP923879532, Tij, Tic); |
1646 | | Tjb = FMA(KP707106781, Tja, Tj9); |
1647 | | Tjc = Tim - Til; |
1648 | | cr[WS(rs, 60)] = FMS(KP923879532, Tjc, Tjb); |
1649 | | ci[WS(rs, 35)] = FMA(KP923879532, Tjc, Tjb); |
1650 | | } |
1651 | | { |
1652 | | E Tjd, Tje, Tik, Tin; |
1653 | | Tjd = FNMS(KP707106781, Tja, Tj9); |
1654 | | Tje = Tii - Tif; |
1655 | | cr[WS(rs, 44)] = FMS(KP923879532, Tje, Tjd); |
1656 | | ci[WS(rs, 51)] = FMA(KP923879532, Tje, Tjd); |
1657 | | Tik = FNMS(KP707106781, Tib, Ti8); |
1658 | | Tin = Til + Tim; |
1659 | | ci[WS(rs, 19)] = FNMS(KP923879532, Tin, Tik); |
1660 | | cr[WS(rs, 12)] = FMA(KP923879532, Tin, Tik); |
1661 | | } |
1662 | | } |
1663 | | { |
1664 | | E Tf2, TjJ, Tgo, TjD, TgI, Tjv, Tha, Tjp, Tfp, Tjw, Tgr, Tjq, Th4, Tho, Th7; |
1665 | | E Thk, TfR, TgB, Tgl, Tgv, TgP, TjK, Thd, TjE, TgX, Thn, Th8, Thh, Tgi, TgC; |
1666 | | E Tgm, Tgy; |
1667 | | { |
1668 | | E TeQ, TjB, Tf1, TjC, TeV, Tf0; |
1669 | | TeQ = TeM + TeP; |
1670 | | TjB = Tjm - Tjl; |
1671 | | TeV = TeR - TeU; |
1672 | | Tf0 = TeW + TeZ; |
1673 | | Tf1 = TeV + Tf0; |
1674 | | TjC = Tf0 - TeV; |
1675 | | Tf2 = FNMS(KP707106781, Tf1, TeQ); |
1676 | | TjJ = FNMS(KP707106781, TjC, TjB); |
1677 | | Tgo = FMA(KP707106781, Tf1, TeQ); |
1678 | | TjD = FMA(KP707106781, TjC, TjB); |
1679 | | } |
1680 | | { |
1681 | | E TgE, Tjn, TgH, Tjo, TgF, TgG; |
1682 | | TgE = TeM - TeP; |
1683 | | Tjn = Tjl + Tjm; |
1684 | | TgF = TeR + TeU; |
1685 | | TgG = TeW - TeZ; |
1686 | | TgH = TgF + TgG; |
1687 | | Tjo = TgF - TgG; |
1688 | | TgI = FMA(KP707106781, TgH, TgE); |
1689 | | Tjv = FNMS(KP707106781, Tjo, Tjn); |
1690 | | Tha = FNMS(KP707106781, TgH, TgE); |
1691 | | Tjp = FMA(KP707106781, Tjo, Tjn); |
1692 | | } |
1693 | | { |
1694 | | E Tfd, Tgp, Tfo, Tgq; |
1695 | | { |
1696 | | E Tf7, Tfc, Tfi, Tfn; |
1697 | | Tf7 = Tf5 + Tf6; |
1698 | | Tfc = Tf8 + Tfb; |
1699 | | Tfd = FMA(KP414213562, Tfc, Tf7); |
1700 | | Tgp = FNMS(KP414213562, Tf7, Tfc); |
1701 | | Tfi = Tfg + Tfh; |
1702 | | Tfn = Tfj + Tfm; |
1703 | | Tfo = FNMS(KP414213562, Tfn, Tfi); |
1704 | | Tgq = FMA(KP414213562, Tfi, Tfn); |
1705 | | } |
1706 | | Tfp = Tfd - Tfo; |
1707 | | Tjw = Tgq - Tgp; |
1708 | | Tgr = Tgp + Tgq; |
1709 | | Tjq = Tfd + Tfo; |
1710 | | } |
1711 | | { |
1712 | | E Th0, Thi, Th3, Thj; |
1713 | | { |
1714 | | E TgY, TgZ, Th1, Th2; |
1715 | | TgY = TfS - TfV; |
1716 | | TgZ = Tgf + Tge; |
1717 | | Th0 = FMA(KP707106781, TgZ, TgY); |
1718 | | Thi = FNMS(KP707106781, TgZ, TgY); |
1719 | | Th1 = Tgc + Tg9; |
1720 | | Th2 = Tg6 - Tg1; |
1721 | | Th3 = FMA(KP707106781, Th2, Th1); |
1722 | | Thj = FNMS(KP707106781, Th2, Th1); |
1723 | | } |
1724 | | Th4 = FNMS(KP198912367, Th3, Th0); |
1725 | | Tho = FNMS(KP668178637, Thi, Thj); |
1726 | | Th7 = FMA(KP198912367, Th0, Th3); |
1727 | | Thk = FMA(KP668178637, Thj, Thi); |
1728 | | } |
1729 | | { |
1730 | | E TfH, Tgt, TfQ, Tgu; |
1731 | | { |
1732 | | E Tfv, TfG, TfM, TfP; |
1733 | | Tfv = Tfr + Tfu; |
1734 | | TfG = TfA + TfF; |
1735 | | TfH = FNMS(KP707106781, TfG, Tfv); |
1736 | | Tgt = FMA(KP707106781, TfG, Tfv); |
1737 | | TfM = TfK + TfL; |
1738 | | TfP = TfN + TfO; |
1739 | | TfQ = FNMS(KP707106781, TfP, TfM); |
1740 | | Tgu = FMA(KP707106781, TfP, TfM); |
1741 | | } |
1742 | | TfR = FMA(KP668178637, TfQ, TfH); |
1743 | | TgB = FMA(KP198912367, Tgt, Tgu); |
1744 | | Tgl = FNMS(KP668178637, TfH, TfQ); |
1745 | | Tgv = FNMS(KP198912367, Tgu, Tgt); |
1746 | | } |
1747 | | { |
1748 | | E TgL, Thc, TgO, Thb; |
1749 | | { |
1750 | | E TgJ, TgK, TgM, TgN; |
1751 | | TgJ = Tf8 - Tfb; |
1752 | | TgK = Tf5 - Tf6; |
1753 | | TgL = FMA(KP414213562, TgK, TgJ); |
1754 | | Thc = FNMS(KP414213562, TgJ, TgK); |
1755 | | TgM = Tfj - Tfm; |
1756 | | TgN = Tfg - Tfh; |
1757 | | TgO = FNMS(KP414213562, TgN, TgM); |
1758 | | Thb = FMA(KP414213562, TgM, TgN); |
1759 | | } |
1760 | | TgP = TgL + TgO; |
1761 | | TjK = TgL - TgO; |
1762 | | Thd = Thb - Thc; |
1763 | | TjE = Thc + Thb; |
1764 | | } |
1765 | | { |
1766 | | E TgT, Thf, TgW, Thg; |
1767 | | { |
1768 | | E TgR, TgS, TgU, TgV; |
1769 | | TgR = Tfr - Tfu; |
1770 | | TgS = TfN - TfO; |
1771 | | TgT = FMA(KP707106781, TgS, TgR); |
1772 | | Thf = FNMS(KP707106781, TgS, TgR); |
1773 | | TgU = TfK - TfL; |
1774 | | TgV = TfF - TfA; |
1775 | | TgW = FMA(KP707106781, TgV, TgU); |
1776 | | Thg = FNMS(KP707106781, TgV, TgU); |
1777 | | } |
1778 | | TgX = FMA(KP198912367, TgW, TgT); |
1779 | | Thn = FMA(KP668178637, Thf, Thg); |
1780 | | Th8 = FNMS(KP198912367, TgT, TgW); |
1781 | | Thh = FNMS(KP668178637, Thg, Thf); |
1782 | | } |
1783 | | { |
1784 | | E Tg8, Tgw, Tgh, Tgx; |
1785 | | { |
1786 | | E TfW, Tg7, Tgd, Tgg; |
1787 | | TfW = TfS + TfV; |
1788 | | Tg7 = Tg1 + Tg6; |
1789 | | Tg8 = FNMS(KP707106781, Tg7, TfW); |
1790 | | Tgw = FMA(KP707106781, Tg7, TfW); |
1791 | | Tgd = Tg9 - Tgc; |
1792 | | Tgg = Tge - Tgf; |
1793 | | Tgh = FNMS(KP707106781, Tgg, Tgd); |
1794 | | Tgx = FMA(KP707106781, Tgg, Tgd); |
1795 | | } |
1796 | | Tgi = FMA(KP668178637, Tgh, Tg8); |
1797 | | TgC = FMA(KP198912367, Tgw, Tgx); |
1798 | | Tgm = FNMS(KP668178637, Tg8, Tgh); |
1799 | | Tgy = FNMS(KP198912367, Tgx, Tgw); |
1800 | | } |
1801 | | { |
1802 | | E Tfq, Tgj, TgA, TgD; |
1803 | | Tfq = FMA(KP923879532, Tfp, Tf2); |
1804 | | Tgj = TfR + Tgi; |
1805 | | ci[WS(rs, 25)] = FNMS(KP831469612, Tgj, Tfq); |
1806 | | cr[WS(rs, 6)] = FMA(KP831469612, Tgj, Tfq); |
1807 | | TgA = FNMS(KP923879532, Tgr, Tgo); |
1808 | | TgD = TgB + TgC; |
1809 | | ci[WS(rs, 17)] = FNMS(KP980785280, TgD, TgA); |
1810 | | cr[WS(rs, 14)] = FMA(KP980785280, TgD, TgA); |
1811 | | } |
1812 | | { |
1813 | | E TjF, TjG, TjN, TjO; |
1814 | | TjF = FMA(KP923879532, TjE, TjD); |
1815 | | TjG = Th8 + Th7; |
1816 | | cr[WS(rs, 34)] = FMS(KP980785280, TjG, TjF); |
1817 | | ci[WS(rs, 61)] = FMA(KP980785280, TjG, TjF); |
1818 | | TjN = FNMS(KP923879532, TjK, TjJ); |
1819 | | TjO = Thk - Thh; |
1820 | | cr[WS(rs, 42)] = FMS(KP831469612, TjO, TjN); |
1821 | | ci[WS(rs, 53)] = FMA(KP831469612, TjO, TjN); |
1822 | | } |
1823 | | { |
1824 | | E Tgk, Tgn, Tgs, Tgz; |
1825 | | Tgk = FNMS(KP923879532, Tfp, Tf2); |
1826 | | Tgn = Tgl + Tgm; |
1827 | | cr[WS(rs, 22)] = FMA(KP831469612, Tgn, Tgk); |
1828 | | ci[WS(rs, 9)] = FNMS(KP831469612, Tgn, Tgk); |
1829 | | Tgs = FMA(KP923879532, Tgr, Tgo); |
1830 | | Tgz = Tgv + Tgy; |
1831 | | cr[WS(rs, 30)] = FNMS(KP980785280, Tgz, Tgs); |
1832 | | ci[WS(rs, 1)] = FMA(KP980785280, Tgz, Tgs); |
1833 | | } |
1834 | | { |
1835 | | E TjH, TjI, TjL, TjM; |
1836 | | TjH = FNMS(KP923879532, TjE, TjD); |
1837 | | TjI = Th4 - TgX; |
1838 | | cr[WS(rs, 50)] = FMS(KP980785280, TjI, TjH); |
1839 | | ci[WS(rs, 45)] = FMA(KP980785280, TjI, TjH); |
1840 | | TjL = FMA(KP923879532, TjK, TjJ); |
1841 | | TjM = Thn + Tho; |
1842 | | cr[WS(rs, 58)] = -(FMA(KP831469612, TjM, TjL)); |
1843 | | ci[WS(rs, 37)] = FNMS(KP831469612, TjM, TjL); |
1844 | | } |
1845 | | { |
1846 | | E TgQ, Th5, Thm, Thp; |
1847 | | TgQ = FMA(KP923879532, TgP, TgI); |
1848 | | Th5 = TgX + Th4; |
1849 | | ci[WS(rs, 29)] = FNMS(KP980785280, Th5, TgQ); |
1850 | | cr[WS(rs, 2)] = FMA(KP980785280, Th5, TgQ); |
1851 | | Thm = FNMS(KP923879532, Thd, Tha); |
1852 | | Thp = Thn - Tho; |
1853 | | ci[WS(rs, 21)] = FNMS(KP831469612, Thp, Thm); |
1854 | | cr[WS(rs, 10)] = FMA(KP831469612, Thp, Thm); |
1855 | | } |
1856 | | { |
1857 | | E Tjr, Tjs, Tjz, TjA; |
1858 | | Tjr = FMA(KP923879532, Tjq, Tjp); |
1859 | | Tjs = TgC - TgB; |
1860 | | cr[WS(rs, 62)] = FMS(KP980785280, Tjs, Tjr); |
1861 | | ci[WS(rs, 33)] = FMA(KP980785280, Tjs, Tjr); |
1862 | | Tjz = FNMS(KP923879532, Tjw, Tjv); |
1863 | | TjA = Tgi - TfR; |
1864 | | cr[WS(rs, 54)] = FMS(KP831469612, TjA, Tjz); |
1865 | | ci[WS(rs, 41)] = FMA(KP831469612, TjA, Tjz); |
1866 | | } |
1867 | | { |
1868 | | E Th6, Th9, The, Thl; |
1869 | | Th6 = FNMS(KP923879532, TgP, TgI); |
1870 | | Th9 = Th7 - Th8; |
1871 | | cr[WS(rs, 18)] = FNMS(KP980785280, Th9, Th6); |
1872 | | ci[WS(rs, 13)] = FMA(KP980785280, Th9, Th6); |
1873 | | The = FMA(KP923879532, Thd, Tha); |
1874 | | Thl = Thh + Thk; |
1875 | | cr[WS(rs, 26)] = FNMS(KP831469612, Thl, The); |
1876 | | ci[WS(rs, 5)] = FMA(KP831469612, Thl, The); |
1877 | | } |
1878 | | { |
1879 | | E Tjt, Tju, Tjx, Tjy; |
1880 | | Tjt = FNMS(KP923879532, Tjq, Tjp); |
1881 | | Tju = Tgy - Tgv; |
1882 | | cr[WS(rs, 46)] = FMS(KP980785280, Tju, Tjt); |
1883 | | ci[WS(rs, 49)] = FMA(KP980785280, Tju, Tjt); |
1884 | | Tjx = FMA(KP923879532, Tjw, Tjv); |
1885 | | Tjy = Tgl - Tgm; |
1886 | | cr[WS(rs, 38)] = FMS(KP831469612, Tjy, Tjx); |
1887 | | ci[WS(rs, 57)] = FMA(KP831469612, Tjy, Tjx); |
1888 | | } |
1889 | | } |
1890 | | { |
1891 | | E Td1, Tkw, TdN, Tkq, Tdl, TdX, TdI, TdR, Teq, TeK, Tet, TeG, Tej, TeJ, Teu; |
1892 | | E TeD, Teb, TkK, Tez, TkE, TdE, TdY, TdH, TdU, TcM, TkD, TkJ, TdK, Te4, Tkp; |
1893 | | E Tkv, Tew; |
1894 | | { |
1895 | | E TcT, TdM, Td0, TdL; |
1896 | | { |
1897 | | E TcP, TcS, TcW, TcZ; |
1898 | | TcP = FMA(KP707106781, TcO, TcN); |
1899 | | TcS = FMA(KP707106781, TcR, TcQ); |
1900 | | TcT = FMA(KP198912367, TcS, TcP); |
1901 | | TdM = FNMS(KP198912367, TcP, TcS); |
1902 | | TcW = FMA(KP707106781, TcV, TcU); |
1903 | | TcZ = FMA(KP707106781, TcY, TcX); |
1904 | | Td0 = FNMS(KP198912367, TcZ, TcW); |
1905 | | TdL = FMA(KP198912367, TcW, TcZ); |
1906 | | } |
1907 | | Td1 = TcT + Td0; |
1908 | | Tkw = TcT - Td0; |
1909 | | TdN = TdL - TdM; |
1910 | | Tkq = TdM + TdL; |
1911 | | } |
1912 | | { |
1913 | | E Tdd, TdP, Tdk, TdQ; |
1914 | | { |
1915 | | E Td5, Tdc, Tdg, Tdj; |
1916 | | Td5 = FMA(KP707106781, Td4, Td3); |
1917 | | Tdc = Td8 + Tdb; |
1918 | | Tdd = FMA(KP923879532, Tdc, Td5); |
1919 | | TdP = FNMS(KP923879532, Tdc, Td5); |
1920 | | Tdg = FMA(KP707106781, Tdf, Tde); |
1921 | | Tdj = Tdh + Tdi; |
1922 | | Tdk = FMA(KP923879532, Tdj, Tdg); |
1923 | | TdQ = FNMS(KP923879532, Tdj, Tdg); |
1924 | | } |
1925 | | Tdl = FMA(KP098491403, Tdk, Tdd); |
1926 | | TdX = FMA(KP820678790, TdP, TdQ); |
1927 | | TdI = FNMS(KP098491403, Tdd, Tdk); |
1928 | | TdR = FNMS(KP820678790, TdQ, TdP); |
1929 | | } |
1930 | | { |
1931 | | E Tem, TeE, Tep, TeF; |
1932 | | { |
1933 | | E Tek, Tel, Ten, Teo; |
1934 | | Tek = FNMS(KP707106781, Tdn, Tdm); |
1935 | | Tel = TdB - TdA; |
1936 | | Tem = FNMS(KP923879532, Tel, Tek); |
1937 | | TeE = FMA(KP923879532, Tel, Tek); |
1938 | | Ten = FNMS(KP707106781, Tdy, Tdx); |
1939 | | Teo = Tdu - Tdr; |
1940 | | Tep = FMA(KP923879532, Teo, Ten); |
1941 | | TeF = FNMS(KP923879532, Teo, Ten); |
1942 | | } |
1943 | | Teq = FNMS(KP534511135, Tep, Tem); |
1944 | | TeK = FNMS(KP303346683, TeE, TeF); |
1945 | | Tet = FMA(KP534511135, Tem, Tep); |
1946 | | TeG = FMA(KP303346683, TeF, TeE); |
1947 | | } |
1948 | | { |
1949 | | E Tef, TeB, Tei, TeC; |
1950 | | { |
1951 | | E Ted, Tee, Teg, Teh; |
1952 | | Ted = FNMS(KP707106781, Td4, Td3); |
1953 | | Tee = Tdi - Tdh; |
1954 | | Tef = FNMS(KP923879532, Tee, Ted); |
1955 | | TeB = FMA(KP923879532, Tee, Ted); |
1956 | | Teg = FNMS(KP707106781, Tdf, Tde); |
1957 | | Teh = Td8 - Tdb; |
1958 | | Tei = FNMS(KP923879532, Teh, Teg); |
1959 | | TeC = FMA(KP923879532, Teh, Teg); |
1960 | | } |
1961 | | Tej = FMA(KP534511135, Tei, Tef); |
1962 | | TeJ = FMA(KP303346683, TeB, TeC); |
1963 | | Teu = FNMS(KP534511135, Tef, Tei); |
1964 | | TeD = FNMS(KP303346683, TeC, TeB); |
1965 | | } |
1966 | | { |
1967 | | E Te7, Tex, Tea, Tey; |
1968 | | { |
1969 | | E Te5, Te6, Te8, Te9; |
1970 | | Te5 = FNMS(KP707106781, TcR, TcQ); |
1971 | | Te6 = FNMS(KP707106781, TcO, TcN); |
1972 | | Te7 = FMA(KP668178637, Te6, Te5); |
1973 | | Tex = FNMS(KP668178637, Te5, Te6); |
1974 | | Te8 = FNMS(KP707106781, TcY, TcX); |
1975 | | Te9 = FNMS(KP707106781, TcV, TcU); |
1976 | | Tea = FNMS(KP668178637, Te9, Te8); |
1977 | | Tey = FMA(KP668178637, Te8, Te9); |
1978 | | } |
1979 | | Teb = Te7 - Tea; |
1980 | | TkK = Tey - Tex; |
1981 | | Tez = Tex + Tey; |
1982 | | TkE = Te7 + Tea; |
1983 | | } |
1984 | | { |
1985 | | E Tdw, TdS, TdD, TdT; |
1986 | | { |
1987 | | E Tdo, Tdv, Tdz, TdC; |
1988 | | Tdo = FMA(KP707106781, Tdn, Tdm); |
1989 | | Tdv = Tdr + Tdu; |
1990 | | Tdw = FMA(KP923879532, Tdv, Tdo); |
1991 | | TdS = FNMS(KP923879532, Tdv, Tdo); |
1992 | | Tdz = FMA(KP707106781, Tdy, Tdx); |
1993 | | TdC = TdA + TdB; |
1994 | | TdD = FMA(KP923879532, TdC, Tdz); |
1995 | | TdT = FNMS(KP923879532, TdC, Tdz); |
1996 | | } |
1997 | | TdE = FNMS(KP098491403, TdD, Tdw); |
1998 | | TdY = FNMS(KP820678790, TdS, TdT); |
1999 | | TdH = FMA(KP098491403, Tdw, TdD); |
2000 | | TdU = FMA(KP820678790, TdT, TdS); |
2001 | | } |
2002 | | { |
2003 | | E TcE, Te0, Tkn, TkB, TcL, TkC, Te3, Tko, TcD, Tkm; |
2004 | | TcD = TcB + TcC; |
2005 | | TcE = FMA(KP707106781, TcD, TcA); |
2006 | | Te0 = FNMS(KP707106781, TcD, TcA); |
2007 | | Tkm = T7s - T7l; |
2008 | | Tkn = FMA(KP707106781, Tkm, Tkl); |
2009 | | TkB = FNMS(KP707106781, Tkm, Tkl); |
2010 | | { |
2011 | | E TcH, TcK, Te1, Te2; |
2012 | | TcH = FMA(KP414213562, TcG, TcF); |
2013 | | TcK = FNMS(KP414213562, TcJ, TcI); |
2014 | | TcL = TcH + TcK; |
2015 | | TkC = TcH - TcK; |
2016 | | Te1 = FMA(KP414213562, TcI, TcJ); |
2017 | | Te2 = FNMS(KP414213562, TcF, TcG); |
2018 | | Te3 = Te1 - Te2; |
2019 | | Tko = Te2 + Te1; |
2020 | | } |
2021 | | TcM = FMA(KP923879532, TcL, TcE); |
2022 | | TkD = FMA(KP923879532, TkC, TkB); |
2023 | | TkJ = FNMS(KP923879532, TkC, TkB); |
2024 | | TdK = FNMS(KP923879532, TcL, TcE); |
2025 | | Te4 = FNMS(KP923879532, Te3, Te0); |
2026 | | Tkp = FMA(KP923879532, Tko, Tkn); |
2027 | | Tkv = FNMS(KP923879532, Tko, Tkn); |
2028 | | Tew = FMA(KP923879532, Te3, Te0); |
2029 | | } |
2030 | | { |
2031 | | E Td2, TdF, TdW, TdZ; |
2032 | | Td2 = FMA(KP980785280, Td1, TcM); |
2033 | | TdF = Tdl + TdE; |
2034 | | ci[WS(rs, 30)] = FNMS(KP995184726, TdF, Td2); |
2035 | | cr[WS(rs, 1)] = FMA(KP995184726, TdF, Td2); |
2036 | | TdW = FNMS(KP980785280, TdN, TdK); |
2037 | | TdZ = TdX - TdY; |
2038 | | ci[WS(rs, 22)] = FNMS(KP773010453, TdZ, TdW); |
2039 | | cr[WS(rs, 9)] = FMA(KP773010453, TdZ, TdW); |
2040 | | } |
2041 | | { |
2042 | | E TkF, TkG, TkN, TkO; |
2043 | | TkF = FMA(KP831469612, TkE, TkD); |
2044 | | TkG = TeJ + TeK; |
2045 | | cr[WS(rs, 61)] = -(FMA(KP956940335, TkG, TkF)); |
2046 | | ci[WS(rs, 34)] = FNMS(KP956940335, TkG, TkF); |
2047 | | TkN = FNMS(KP831469612, TkK, TkJ); |
2048 | | TkO = Teq - Tej; |
2049 | | cr[WS(rs, 53)] = FMS(KP881921264, TkO, TkN); |
2050 | | ci[WS(rs, 42)] = FMA(KP881921264, TkO, TkN); |
2051 | | } |
2052 | | { |
2053 | | E TdG, TdJ, TdO, TdV; |
2054 | | TdG = FNMS(KP980785280, Td1, TcM); |
2055 | | TdJ = TdH - TdI; |
2056 | | cr[WS(rs, 17)] = FNMS(KP995184726, TdJ, TdG); |
2057 | | ci[WS(rs, 14)] = FMA(KP995184726, TdJ, TdG); |
2058 | | TdO = FMA(KP980785280, TdN, TdK); |
2059 | | TdV = TdR + TdU; |
2060 | | cr[WS(rs, 25)] = FNMS(KP773010453, TdV, TdO); |
2061 | | ci[WS(rs, 6)] = FMA(KP773010453, TdV, TdO); |
2062 | | } |
2063 | | { |
2064 | | E TkH, TkI, TkL, TkM; |
2065 | | TkH = FNMS(KP831469612, TkE, TkD); |
2066 | | TkI = TeG - TeD; |
2067 | | cr[WS(rs, 45)] = FMS(KP956940335, TkI, TkH); |
2068 | | ci[WS(rs, 50)] = FMA(KP956940335, TkI, TkH); |
2069 | | TkL = FMA(KP831469612, TkK, TkJ); |
2070 | | TkM = Teu + Tet; |
2071 | | cr[WS(rs, 37)] = FMS(KP881921264, TkM, TkL); |
2072 | | ci[WS(rs, 58)] = FMA(KP881921264, TkM, TkL); |
2073 | | } |
2074 | | { |
2075 | | E Tec, Ter, TeI, TeL; |
2076 | | Tec = FMA(KP831469612, Teb, Te4); |
2077 | | Ter = Tej + Teq; |
2078 | | ci[WS(rs, 26)] = FNMS(KP881921264, Ter, Tec); |
2079 | | cr[WS(rs, 5)] = FMA(KP881921264, Ter, Tec); |
2080 | | TeI = FNMS(KP831469612, Tez, Tew); |
2081 | | TeL = TeJ - TeK; |
2082 | | ci[WS(rs, 18)] = FNMS(KP956940335, TeL, TeI); |
2083 | | cr[WS(rs, 13)] = FMA(KP956940335, TeL, TeI); |
2084 | | } |
2085 | | { |
2086 | | E Tkr, Tks, Tkz, TkA; |
2087 | | Tkr = FMA(KP980785280, Tkq, Tkp); |
2088 | | Tks = TdI + TdH; |
2089 | | cr[WS(rs, 33)] = FMS(KP995184726, Tks, Tkr); |
2090 | | ci[WS(rs, 62)] = FMA(KP995184726, Tks, Tkr); |
2091 | | Tkz = FNMS(KP980785280, Tkw, Tkv); |
2092 | | TkA = TdU - TdR; |
2093 | | cr[WS(rs, 41)] = FMS(KP773010453, TkA, Tkz); |
2094 | | ci[WS(rs, 54)] = FMA(KP773010453, TkA, Tkz); |
2095 | | } |
2096 | | { |
2097 | | E Tes, Tev, TeA, TeH; |
2098 | | Tes = FNMS(KP831469612, Teb, Te4); |
2099 | | Tev = Tet - Teu; |
2100 | | cr[WS(rs, 21)] = FNMS(KP881921264, Tev, Tes); |
2101 | | ci[WS(rs, 10)] = FMA(KP881921264, Tev, Tes); |
2102 | | TeA = FMA(KP831469612, Tez, Tew); |
2103 | | TeH = TeD + TeG; |
2104 | | cr[WS(rs, 29)] = FNMS(KP956940335, TeH, TeA); |
2105 | | ci[WS(rs, 2)] = FMA(KP956940335, TeH, TeA); |
2106 | | } |
2107 | | { |
2108 | | E Tkt, Tku, Tkx, Tky; |
2109 | | Tkt = FNMS(KP980785280, Tkq, Tkp); |
2110 | | Tku = TdE - Tdl; |
2111 | | cr[WS(rs, 49)] = FMS(KP995184726, Tku, Tkt); |
2112 | | ci[WS(rs, 46)] = FMA(KP995184726, Tku, Tkt); |
2113 | | Tkx = FMA(KP980785280, Tkw, Tkv); |
2114 | | Tky = TdX + TdY; |
2115 | | cr[WS(rs, 57)] = -(FMA(KP773010453, Tky, Tkx)); |
2116 | | ci[WS(rs, 38)] = FNMS(KP773010453, Tky, Tkx); |
2117 | | } |
2118 | | } |
2119 | | } |
2120 | | } |
2121 | | } |
2122 | | |
2123 | | static const tw_instr twinstr[] = { |
2124 | | { TW_FULL, 1, 64 }, |
2125 | | { TW_NEXT, 1, 0 } |
2126 | | }; |
2127 | | |
2128 | | static const hc2hc_desc desc = { 64, "hf_64", twinstr, &GENUS, { 520, 126, 518, 0 } }; |
2129 | | |
2130 | | void X(codelet_hf_64) (planner *p) { |
2131 | | X(khc2hc_register) (p, hf_64, &desc); |
2132 | | } |
2133 | | #else |
2134 | | |
2135 | | /* Generated by: ../../../genfft/gen_hc2hc.native -compact -variables 4 -pipeline-latency 4 -n 64 -dit -name hf_64 -include rdft/scalar/hf.h */ |
2136 | | |
2137 | | /* |
2138 | | * This function contains 1038 FP additions, 500 FP multiplications, |
2139 | | * (or, 808 additions, 270 multiplications, 230 fused multiply/add), |
2140 | | * 176 stack variables, 15 constants, and 256 memory accesses |
2141 | | */ |
2142 | | #include "rdft/scalar/hf.h" |
2143 | | |
2144 | | static void hf_64(R *cr, R *ci, const R *W, stride rs, INT mb, INT me, INT ms) |
2145 | 0 | { |
2146 | 0 | DK(KP290284677, +0.290284677254462367636192375817395274691476278); |
2147 | 0 | DK(KP956940335, +0.956940335732208864935797886980269969482849206); |
2148 | 0 | DK(KP881921264, +0.881921264348355029712756863660388349508442621); |
2149 | 0 | DK(KP471396736, +0.471396736825997648556387625905254377657460319); |
2150 | 0 | DK(KP555570233, +0.555570233019602224742830813948532874374937191); |
2151 | 0 | DK(KP831469612, +0.831469612302545237078788377617905756738560812); |
2152 | 0 | DK(KP098017140, +0.098017140329560601994195563888641845861136673); |
2153 | 0 | DK(KP995184726, +0.995184726672196886244836953109479921575474869); |
2154 | 0 | DK(KP773010453, +0.773010453362736960810906609758469800971041293); |
2155 | 0 | DK(KP634393284, +0.634393284163645498215171613225493370675687095); |
2156 | 0 | DK(KP980785280, +0.980785280403230449126182236134239036973933731); |
2157 | 0 | DK(KP195090322, +0.195090322016128267848284868477022240927691618); |
2158 | 0 | DK(KP382683432, +0.382683432365089771728459984030398866761344562); |
2159 | 0 | DK(KP923879532, +0.923879532511286756128183189396788286822416626); |
2160 | 0 | DK(KP707106781, +0.707106781186547524400844362104849039284835938); |
2161 | 0 | { |
2162 | 0 | INT m; |
2163 | 0 | for (m = mb, W = W + ((mb - 1) * 126); m < me; m = m + 1, cr = cr + ms, ci = ci - ms, W = W + 126, MAKE_VOLATILE_STRIDE(128, rs)) { |
2164 | 0 | E Tj, TcL, ThT, Tin, T6b, Taz, TgT, Thn, TG, Thm, TcO, TgO, T6m, Tim, TaC; |
2165 | 0 | E ThQ, T14, Tfr, T6y, T9O, TaG, Tc0, TcU, TeE, T1r, Tfq, T6J, T9P, TaJ, Tc1; |
2166 | 0 | E TcZ, TeF, T1Q, T2d, Tfu, Tfv, Tfw, Tfx, T6Q, TaM, Tdb, TeI, T71, TaQ, T7a; |
2167 | 0 | E TaN, Td6, TeJ, T77, TaP, T2B, T2Y, Tfz, TfA, TfB, TfC, T7h, TaW, Tdm, TeL; |
2168 | 0 | E T7s, TaU, T7B, TaX, Tdh, TeM, T7y, TaT, T5j, TfR, Tec, TeX, TfY, Tgy, T8D; |
2169 | 0 | E Tbl, T8O, Tbx, T9l, Tbm, TdV, Tf0, T9i, Tbw, T3M, TfL, TdL, TeT, TfI, Tgt; |
2170 | 0 | E T7K, Tbd, T7V, Tb3, T8s, Tbe, Tdu, TeQ, T8p, Tb2, T4x, TfJ, TdE, TdM, TfO; |
2171 | 0 | E Tgu, T87, T8u, T8i, T8v, Tba, Tbh, Tdz, TdN, Tb7, Tbg, T64, TfZ, Te5, Ted; |
2172 | 0 | E TfU, Tgz, T90, T9n, T9b, T9o, Tbt, TbA, Te0, Tee, Tbq, Tbz; |
2173 | 0 | { |
2174 | 0 | E T1, TgR, T6, TgQ, Tc, T68, Th, T69; |
2175 | 0 | T1 = cr[0]; |
2176 | 0 | TgR = ci[0]; |
2177 | 0 | { |
2178 | 0 | E T3, T5, T2, T4; |
2179 | 0 | T3 = cr[WS(rs, 32)]; |
2180 | 0 | T5 = ci[WS(rs, 32)]; |
2181 | 0 | T2 = W[62]; |
2182 | 0 | T4 = W[63]; |
2183 | 0 | T6 = FMA(T2, T3, T4 * T5); |
2184 | 0 | TgQ = FNMS(T4, T3, T2 * T5); |
2185 | 0 | } |
2186 | 0 | { |
2187 | 0 | E T9, Tb, T8, Ta; |
2188 | 0 | T9 = cr[WS(rs, 16)]; |
2189 | 0 | Tb = ci[WS(rs, 16)]; |
2190 | 0 | T8 = W[30]; |
2191 | 0 | Ta = W[31]; |
2192 | 0 | Tc = FMA(T8, T9, Ta * Tb); |
2193 | 0 | T68 = FNMS(Ta, T9, T8 * Tb); |
2194 | 0 | } |
2195 | 0 | { |
2196 | 0 | E Te, Tg, Td, Tf; |
2197 | 0 | Te = cr[WS(rs, 48)]; |
2198 | 0 | Tg = ci[WS(rs, 48)]; |
2199 | 0 | Td = W[94]; |
2200 | 0 | Tf = W[95]; |
2201 | 0 | Th = FMA(Td, Te, Tf * Tg); |
2202 | 0 | T69 = FNMS(Tf, Te, Td * Tg); |
2203 | 0 | } |
2204 | 0 | { |
2205 | 0 | E T7, Ti, ThR, ThS; |
2206 | 0 | T7 = T1 + T6; |
2207 | 0 | Ti = Tc + Th; |
2208 | 0 | Tj = T7 + Ti; |
2209 | 0 | TcL = T7 - Ti; |
2210 | 0 | ThR = Tc - Th; |
2211 | 0 | ThS = TgR - TgQ; |
2212 | 0 | ThT = ThR + ThS; |
2213 | 0 | Tin = ThS - ThR; |
2214 | 0 | } |
2215 | 0 | { |
2216 | 0 | E T67, T6a, TgP, TgS; |
2217 | 0 | T67 = T1 - T6; |
2218 | 0 | T6a = T68 - T69; |
2219 | 0 | T6b = T67 - T6a; |
2220 | 0 | Taz = T67 + T6a; |
2221 | 0 | TgP = T68 + T69; |
2222 | 0 | TgS = TgQ + TgR; |
2223 | 0 | TgT = TgP + TgS; |
2224 | 0 | Thn = TgS - TgP; |
2225 | 0 | } |
2226 | 0 | } |
2227 | 0 | { |
2228 | 0 | E To, T6d, Tt, T6e, T6c, T6f, Tz, T6i, TE, T6j, T6h, T6k; |
2229 | 0 | { |
2230 | 0 | E Tl, Tn, Tk, Tm; |
2231 | 0 | Tl = cr[WS(rs, 8)]; |
2232 | 0 | Tn = ci[WS(rs, 8)]; |
2233 | 0 | Tk = W[14]; |
2234 | 0 | Tm = W[15]; |
2235 | 0 | To = FMA(Tk, Tl, Tm * Tn); |
2236 | 0 | T6d = FNMS(Tm, Tl, Tk * Tn); |
2237 | 0 | } |
2238 | 0 | { |
2239 | 0 | E Tq, Ts, Tp, Tr; |
2240 | 0 | Tq = cr[WS(rs, 40)]; |
2241 | 0 | Ts = ci[WS(rs, 40)]; |
2242 | 0 | Tp = W[78]; |
2243 | 0 | Tr = W[79]; |
2244 | 0 | Tt = FMA(Tp, Tq, Tr * Ts); |
2245 | 0 | T6e = FNMS(Tr, Tq, Tp * Ts); |
2246 | 0 | } |
2247 | 0 | T6c = To - Tt; |
2248 | 0 | T6f = T6d - T6e; |
2249 | 0 | { |
2250 | 0 | E Tw, Ty, Tv, Tx; |
2251 | 0 | Tw = cr[WS(rs, 56)]; |
2252 | 0 | Ty = ci[WS(rs, 56)]; |
2253 | 0 | Tv = W[110]; |
2254 | 0 | Tx = W[111]; |
2255 | 0 | Tz = FMA(Tv, Tw, Tx * Ty); |
2256 | 0 | T6i = FNMS(Tx, Tw, Tv * Ty); |
2257 | 0 | } |
2258 | 0 | { |
2259 | 0 | E TB, TD, TA, TC; |
2260 | 0 | TB = cr[WS(rs, 24)]; |
2261 | 0 | TD = ci[WS(rs, 24)]; |
2262 | 0 | TA = W[46]; |
2263 | 0 | TC = W[47]; |
2264 | 0 | TE = FMA(TA, TB, TC * TD); |
2265 | 0 | T6j = FNMS(TC, TB, TA * TD); |
2266 | 0 | } |
2267 | 0 | T6h = Tz - TE; |
2268 | 0 | T6k = T6i - T6j; |
2269 | 0 | { |
2270 | 0 | E Tu, TF, TcM, TcN; |
2271 | 0 | Tu = To + Tt; |
2272 | 0 | TF = Tz + TE; |
2273 | 0 | TG = Tu + TF; |
2274 | 0 | Thm = Tu - TF; |
2275 | 0 | TcM = T6i + T6j; |
2276 | 0 | TcN = T6d + T6e; |
2277 | 0 | TcO = TcM - TcN; |
2278 | 0 | TgO = TcN + TcM; |
2279 | 0 | } |
2280 | 0 | { |
2281 | 0 | E T6g, T6l, TaA, TaB; |
2282 | 0 | T6g = T6c - T6f; |
2283 | 0 | T6l = T6h + T6k; |
2284 | 0 | T6m = KP707106781 * (T6g + T6l); |
2285 | 0 | Tim = KP707106781 * (T6l - T6g); |
2286 | 0 | TaA = T6c + T6f; |
2287 | 0 | TaB = T6h - T6k; |
2288 | 0 | TaC = KP707106781 * (TaA + TaB); |
2289 | 0 | ThQ = KP707106781 * (TaA - TaB); |
2290 | 0 | } |
2291 | 0 | } |
2292 | 0 | { |
2293 | 0 | E TS, TcR, T6o, T6v, T13, TcS, T6r, T6w, T6s, T6x; |
2294 | 0 | { |
2295 | 0 | E TM, T6t, TR, T6u; |
2296 | 0 | { |
2297 | 0 | E TJ, TL, TI, TK; |
2298 | 0 | TJ = cr[WS(rs, 4)]; |
2299 | 0 | TL = ci[WS(rs, 4)]; |
2300 | 0 | TI = W[6]; |
2301 | 0 | TK = W[7]; |
2302 | 0 | TM = FMA(TI, TJ, TK * TL); |
2303 | 0 | T6t = FNMS(TK, TJ, TI * TL); |
2304 | 0 | } |
2305 | 0 | { |
2306 | 0 | E TO, TQ, TN, TP; |
2307 | 0 | TO = cr[WS(rs, 36)]; |
2308 | 0 | TQ = ci[WS(rs, 36)]; |
2309 | 0 | TN = W[70]; |
2310 | 0 | TP = W[71]; |
2311 | 0 | TR = FMA(TN, TO, TP * TQ); |
2312 | 0 | T6u = FNMS(TP, TO, TN * TQ); |
2313 | 0 | } |
2314 | 0 | TS = TM + TR; |
2315 | 0 | TcR = T6t + T6u; |
2316 | 0 | T6o = TM - TR; |
2317 | 0 | T6v = T6t - T6u; |
2318 | 0 | } |
2319 | 0 | { |
2320 | 0 | E TX, T6p, T12, T6q; |
2321 | 0 | { |
2322 | 0 | E TU, TW, TT, TV; |
2323 | 0 | TU = cr[WS(rs, 20)]; |
2324 | 0 | TW = ci[WS(rs, 20)]; |
2325 | 0 | TT = W[38]; |
2326 | 0 | TV = W[39]; |
2327 | 0 | TX = FMA(TT, TU, TV * TW); |
2328 | 0 | T6p = FNMS(TV, TU, TT * TW); |
2329 | 0 | } |
2330 | 0 | { |
2331 | 0 | E TZ, T11, TY, T10; |
2332 | 0 | TZ = cr[WS(rs, 52)]; |
2333 | 0 | T11 = ci[WS(rs, 52)]; |
2334 | 0 | TY = W[102]; |
2335 | 0 | T10 = W[103]; |
2336 | 0 | T12 = FMA(TY, TZ, T10 * T11); |
2337 | 0 | T6q = FNMS(T10, TZ, TY * T11); |
2338 | 0 | } |
2339 | 0 | T13 = TX + T12; |
2340 | 0 | TcS = T6p + T6q; |
2341 | 0 | T6r = T6p - T6q; |
2342 | 0 | T6w = TX - T12; |
2343 | 0 | } |
2344 | 0 | T14 = TS + T13; |
2345 | 0 | Tfr = TcR + TcS; |
2346 | 0 | T6s = T6o - T6r; |
2347 | 0 | T6x = T6v + T6w; |
2348 | 0 | T6y = FNMS(KP382683432, T6x, KP923879532 * T6s); |
2349 | 0 | T9O = FMA(KP923879532, T6x, KP382683432 * T6s); |
2350 | 0 | { |
2351 | 0 | E TaE, TaF, TcQ, TcT; |
2352 | 0 | TaE = T6v - T6w; |
2353 | 0 | TaF = T6o + T6r; |
2354 | 0 | TaG = FMA(KP382683432, TaE, KP923879532 * TaF); |
2355 | 0 | Tc0 = FNMS(KP923879532, TaE, KP382683432 * TaF); |
2356 | 0 | TcQ = TS - T13; |
2357 | 0 | TcT = TcR - TcS; |
2358 | 0 | TcU = TcQ + TcT; |
2359 | 0 | TeE = TcQ - TcT; |
2360 | 0 | } |
2361 | 0 | } |
2362 | 0 | { |
2363 | 0 | E T1f, TcW, T6B, T6E, T1q, TcX, T6C, T6H, T6D, T6I; |
2364 | 0 | { |
2365 | 0 | E T19, T6z, T1e, T6A; |
2366 | 0 | { |
2367 | 0 | E T16, T18, T15, T17; |
2368 | 0 | T16 = cr[WS(rs, 60)]; |
2369 | 0 | T18 = ci[WS(rs, 60)]; |
2370 | 0 | T15 = W[118]; |
2371 | 0 | T17 = W[119]; |
2372 | 0 | T19 = FMA(T15, T16, T17 * T18); |
2373 | 0 | T6z = FNMS(T17, T16, T15 * T18); |
2374 | 0 | } |
2375 | 0 | { |
2376 | 0 | E T1b, T1d, T1a, T1c; |
2377 | 0 | T1b = cr[WS(rs, 28)]; |
2378 | 0 | T1d = ci[WS(rs, 28)]; |
2379 | 0 | T1a = W[54]; |
2380 | 0 | T1c = W[55]; |
2381 | 0 | T1e = FMA(T1a, T1b, T1c * T1d); |
2382 | 0 | T6A = FNMS(T1c, T1b, T1a * T1d); |
2383 | 0 | } |
2384 | 0 | T1f = T19 + T1e; |
2385 | 0 | TcW = T6z + T6A; |
2386 | 0 | T6B = T6z - T6A; |
2387 | 0 | T6E = T19 - T1e; |
2388 | 0 | } |
2389 | 0 | { |
2390 | 0 | E T1k, T6F, T1p, T6G; |
2391 | 0 | { |
2392 | 0 | E T1h, T1j, T1g, T1i; |
2393 | 0 | T1h = cr[WS(rs, 12)]; |
2394 | 0 | T1j = ci[WS(rs, 12)]; |
2395 | 0 | T1g = W[22]; |
2396 | 0 | T1i = W[23]; |
2397 | 0 | T1k = FMA(T1g, T1h, T1i * T1j); |
2398 | 0 | T6F = FNMS(T1i, T1h, T1g * T1j); |
2399 | 0 | } |
2400 | 0 | { |
2401 | 0 | E T1m, T1o, T1l, T1n; |
2402 | 0 | T1m = cr[WS(rs, 44)]; |
2403 | 0 | T1o = ci[WS(rs, 44)]; |
2404 | 0 | T1l = W[86]; |
2405 | 0 | T1n = W[87]; |
2406 | 0 | T1p = FMA(T1l, T1m, T1n * T1o); |
2407 | 0 | T6G = FNMS(T1n, T1m, T1l * T1o); |
2408 | 0 | } |
2409 | 0 | T1q = T1k + T1p; |
2410 | 0 | TcX = T6F + T6G; |
2411 | 0 | T6C = T1k - T1p; |
2412 | 0 | T6H = T6F - T6G; |
2413 | 0 | } |
2414 | 0 | T1r = T1f + T1q; |
2415 | 0 | Tfq = TcW + TcX; |
2416 | 0 | T6D = T6B + T6C; |
2417 | 0 | T6I = T6E - T6H; |
2418 | 0 | T6J = FMA(KP382683432, T6D, KP923879532 * T6I); |
2419 | 0 | T9P = FNMS(KP923879532, T6D, KP382683432 * T6I); |
2420 | 0 | { |
2421 | 0 | E TaH, TaI, TcV, TcY; |
2422 | 0 | TaH = T6E + T6H; |
2423 | 0 | TaI = T6B - T6C; |
2424 | 0 | TaJ = FNMS(KP382683432, TaI, KP923879532 * TaH); |
2425 | 0 | Tc1 = FMA(KP923879532, TaI, KP382683432 * TaH); |
2426 | 0 | TcV = T1f - T1q; |
2427 | 0 | TcY = TcW - TcX; |
2428 | 0 | TcZ = TcV - TcY; |
2429 | 0 | TeF = TcV + TcY; |
2430 | 0 | } |
2431 | 0 | } |
2432 | 0 | { |
2433 | 0 | E T1y, T73, T1D, T74, T1E, Td7, T1J, T6N, T1O, T6O, T1P, Td8, T21, Td4, T6R; |
2434 | 0 | E T6U, T2c, Td3, T6W, T6Z; |
2435 | 0 | { |
2436 | 0 | E T1v, T1x, T1u, T1w; |
2437 | 0 | T1v = cr[WS(rs, 2)]; |
2438 | 0 | T1x = ci[WS(rs, 2)]; |
2439 | 0 | T1u = W[2]; |
2440 | 0 | T1w = W[3]; |
2441 | 0 | T1y = FMA(T1u, T1v, T1w * T1x); |
2442 | 0 | T73 = FNMS(T1w, T1v, T1u * T1x); |
2443 | 0 | } |
2444 | 0 | { |
2445 | 0 | E T1A, T1C, T1z, T1B; |
2446 | 0 | T1A = cr[WS(rs, 34)]; |
2447 | 0 | T1C = ci[WS(rs, 34)]; |
2448 | 0 | T1z = W[66]; |
2449 | 0 | T1B = W[67]; |
2450 | 0 | T1D = FMA(T1z, T1A, T1B * T1C); |
2451 | 0 | T74 = FNMS(T1B, T1A, T1z * T1C); |
2452 | 0 | } |
2453 | 0 | T1E = T1y + T1D; |
2454 | 0 | Td7 = T73 + T74; |
2455 | 0 | { |
2456 | 0 | E T1G, T1I, T1F, T1H; |
2457 | 0 | T1G = cr[WS(rs, 18)]; |
2458 | 0 | T1I = ci[WS(rs, 18)]; |
2459 | 0 | T1F = W[34]; |
2460 | 0 | T1H = W[35]; |
2461 | 0 | T1J = FMA(T1F, T1G, T1H * T1I); |
2462 | 0 | T6N = FNMS(T1H, T1G, T1F * T1I); |
2463 | 0 | } |
2464 | 0 | { |
2465 | 0 | E T1L, T1N, T1K, T1M; |
2466 | 0 | T1L = cr[WS(rs, 50)]; |
2467 | 0 | T1N = ci[WS(rs, 50)]; |
2468 | 0 | T1K = W[98]; |
2469 | 0 | T1M = W[99]; |
2470 | 0 | T1O = FMA(T1K, T1L, T1M * T1N); |
2471 | 0 | T6O = FNMS(T1M, T1L, T1K * T1N); |
2472 | 0 | } |
2473 | 0 | T1P = T1J + T1O; |
2474 | 0 | Td8 = T6N + T6O; |
2475 | 0 | { |
2476 | 0 | E T1V, T6S, T20, T6T; |
2477 | 0 | { |
2478 | 0 | E T1S, T1U, T1R, T1T; |
2479 | 0 | T1S = cr[WS(rs, 10)]; |
2480 | 0 | T1U = ci[WS(rs, 10)]; |
2481 | 0 | T1R = W[18]; |
2482 | 0 | T1T = W[19]; |
2483 | 0 | T1V = FMA(T1R, T1S, T1T * T1U); |
2484 | 0 | T6S = FNMS(T1T, T1S, T1R * T1U); |
2485 | 0 | } |
2486 | 0 | { |
2487 | 0 | E T1X, T1Z, T1W, T1Y; |
2488 | 0 | T1X = cr[WS(rs, 42)]; |
2489 | 0 | T1Z = ci[WS(rs, 42)]; |
2490 | 0 | T1W = W[82]; |
2491 | 0 | T1Y = W[83]; |
2492 | 0 | T20 = FMA(T1W, T1X, T1Y * T1Z); |
2493 | 0 | T6T = FNMS(T1Y, T1X, T1W * T1Z); |
2494 | 0 | } |
2495 | 0 | T21 = T1V + T20; |
2496 | 0 | Td4 = T6S + T6T; |
2497 | 0 | T6R = T1V - T20; |
2498 | 0 | T6U = T6S - T6T; |
2499 | 0 | } |
2500 | 0 | { |
2501 | 0 | E T26, T6X, T2b, T6Y; |
2502 | 0 | { |
2503 | 0 | E T23, T25, T22, T24; |
2504 | 0 | T23 = cr[WS(rs, 58)]; |
2505 | 0 | T25 = ci[WS(rs, 58)]; |
2506 | 0 | T22 = W[114]; |
2507 | 0 | T24 = W[115]; |
2508 | 0 | T26 = FMA(T22, T23, T24 * T25); |
2509 | 0 | T6X = FNMS(T24, T23, T22 * T25); |
2510 | 0 | } |
2511 | 0 | { |
2512 | 0 | E T28, T2a, T27, T29; |
2513 | 0 | T28 = cr[WS(rs, 26)]; |
2514 | 0 | T2a = ci[WS(rs, 26)]; |
2515 | 0 | T27 = W[50]; |
2516 | 0 | T29 = W[51]; |
2517 | 0 | T2b = FMA(T27, T28, T29 * T2a); |
2518 | 0 | T6Y = FNMS(T29, T28, T27 * T2a); |
2519 | 0 | } |
2520 | 0 | T2c = T26 + T2b; |
2521 | 0 | Td3 = T6X + T6Y; |
2522 | 0 | T6W = T26 - T2b; |
2523 | 0 | T6Z = T6X - T6Y; |
2524 | 0 | } |
2525 | 0 | T1Q = T1E + T1P; |
2526 | 0 | T2d = T21 + T2c; |
2527 | 0 | Tfu = T1Q - T2d; |
2528 | 0 | Tfv = Td7 + Td8; |
2529 | 0 | Tfw = Td4 + Td3; |
2530 | 0 | Tfx = Tfv - Tfw; |
2531 | 0 | { |
2532 | 0 | E T6M, T6P, Td9, Tda; |
2533 | 0 | T6M = T1y - T1D; |
2534 | 0 | T6P = T6N - T6O; |
2535 | 0 | T6Q = T6M - T6P; |
2536 | 0 | TaM = T6M + T6P; |
2537 | 0 | Td9 = Td7 - Td8; |
2538 | 0 | Tda = T21 - T2c; |
2539 | 0 | Tdb = Td9 - Tda; |
2540 | 0 | TeI = Td9 + Tda; |
2541 | 0 | } |
2542 | 0 | { |
2543 | 0 | E T6V, T70, T78, T79; |
2544 | 0 | T6V = T6R - T6U; |
2545 | 0 | T70 = T6W + T6Z; |
2546 | 0 | T71 = KP707106781 * (T6V + T70); |
2547 | 0 | TaQ = KP707106781 * (T70 - T6V); |
2548 | 0 | T78 = T6R + T6U; |
2549 | 0 | T79 = T6Z - T6W; |
2550 | 0 | T7a = KP707106781 * (T78 + T79); |
2551 | 0 | TaN = KP707106781 * (T78 - T79); |
2552 | 0 | } |
2553 | 0 | { |
2554 | 0 | E Td2, Td5, T75, T76; |
2555 | 0 | Td2 = T1E - T1P; |
2556 | 0 | Td5 = Td3 - Td4; |
2557 | 0 | Td6 = Td2 - Td5; |
2558 | 0 | TeJ = Td2 + Td5; |
2559 | 0 | T75 = T73 - T74; |
2560 | 0 | T76 = T1J - T1O; |
2561 | 0 | T77 = T75 + T76; |
2562 | 0 | TaP = T75 - T76; |
2563 | 0 | } |
2564 | 0 | } |
2565 | 0 | { |
2566 | 0 | E T2j, T7u, T2o, T7v, T2p, Tdd, T2u, T7e, T2z, T7f, T2A, Tde, T2M, Tdk, T7i; |
2567 | 0 | E T7l, T2X, Tdj, T7n, T7q; |
2568 | 0 | { |
2569 | 0 | E T2g, T2i, T2f, T2h; |
2570 | 0 | T2g = cr[WS(rs, 62)]; |
2571 | 0 | T2i = ci[WS(rs, 62)]; |
2572 | 0 | T2f = W[122]; |
2573 | 0 | T2h = W[123]; |
2574 | 0 | T2j = FMA(T2f, T2g, T2h * T2i); |
2575 | 0 | T7u = FNMS(T2h, T2g, T2f * T2i); |
2576 | 0 | } |
2577 | 0 | { |
2578 | 0 | E T2l, T2n, T2k, T2m; |
2579 | 0 | T2l = cr[WS(rs, 30)]; |
2580 | 0 | T2n = ci[WS(rs, 30)]; |
2581 | 0 | T2k = W[58]; |
2582 | 0 | T2m = W[59]; |
2583 | 0 | T2o = FMA(T2k, T2l, T2m * T2n); |
2584 | 0 | T7v = FNMS(T2m, T2l, T2k * T2n); |
2585 | 0 | } |
2586 | 0 | T2p = T2j + T2o; |
2587 | 0 | Tdd = T7u + T7v; |
2588 | 0 | { |
2589 | 0 | E T2r, T2t, T2q, T2s; |
2590 | 0 | T2r = cr[WS(rs, 14)]; |
2591 | 0 | T2t = ci[WS(rs, 14)]; |
2592 | 0 | T2q = W[26]; |
2593 | 0 | T2s = W[27]; |
2594 | 0 | T2u = FMA(T2q, T2r, T2s * T2t); |
2595 | 0 | T7e = FNMS(T2s, T2r, T2q * T2t); |
2596 | 0 | } |
2597 | 0 | { |
2598 | 0 | E T2w, T2y, T2v, T2x; |
2599 | 0 | T2w = cr[WS(rs, 46)]; |
2600 | 0 | T2y = ci[WS(rs, 46)]; |
2601 | 0 | T2v = W[90]; |
2602 | 0 | T2x = W[91]; |
2603 | 0 | T2z = FMA(T2v, T2w, T2x * T2y); |
2604 | 0 | T7f = FNMS(T2x, T2w, T2v * T2y); |
2605 | 0 | } |
2606 | 0 | T2A = T2u + T2z; |
2607 | 0 | Tde = T7e + T7f; |
2608 | 0 | { |
2609 | 0 | E T2G, T7j, T2L, T7k; |
2610 | 0 | { |
2611 | 0 | E T2D, T2F, T2C, T2E; |
2612 | 0 | T2D = cr[WS(rs, 6)]; |
2613 | 0 | T2F = ci[WS(rs, 6)]; |
2614 | 0 | T2C = W[10]; |
2615 | 0 | T2E = W[11]; |
2616 | 0 | T2G = FMA(T2C, T2D, T2E * T2F); |
2617 | 0 | T7j = FNMS(T2E, T2D, T2C * T2F); |
2618 | 0 | } |
2619 | 0 | { |
2620 | 0 | E T2I, T2K, T2H, T2J; |
2621 | 0 | T2I = cr[WS(rs, 38)]; |
2622 | 0 | T2K = ci[WS(rs, 38)]; |
2623 | 0 | T2H = W[74]; |
2624 | 0 | T2J = W[75]; |
2625 | 0 | T2L = FMA(T2H, T2I, T2J * T2K); |
2626 | 0 | T7k = FNMS(T2J, T2I, T2H * T2K); |
2627 | 0 | } |
2628 | 0 | T2M = T2G + T2L; |
2629 | 0 | Tdk = T7j + T7k; |
2630 | 0 | T7i = T2G - T2L; |
2631 | 0 | T7l = T7j - T7k; |
2632 | 0 | } |
2633 | 0 | { |
2634 | 0 | E T2R, T7o, T2W, T7p; |
2635 | 0 | { |
2636 | 0 | E T2O, T2Q, T2N, T2P; |
2637 | 0 | T2O = cr[WS(rs, 54)]; |
2638 | 0 | T2Q = ci[WS(rs, 54)]; |
2639 | 0 | T2N = W[106]; |
2640 | 0 | T2P = W[107]; |
2641 | 0 | T2R = FMA(T2N, T2O, T2P * T2Q); |
2642 | 0 | T7o = FNMS(T2P, T2O, T2N * T2Q); |
2643 | 0 | } |
2644 | 0 | { |
2645 | 0 | E T2T, T2V, T2S, T2U; |
2646 | 0 | T2T = cr[WS(rs, 22)]; |
2647 | 0 | T2V = ci[WS(rs, 22)]; |
2648 | 0 | T2S = W[42]; |
2649 | 0 | T2U = W[43]; |
2650 | 0 | T2W = FMA(T2S, T2T, T2U * T2V); |
2651 | 0 | T7p = FNMS(T2U, T2T, T2S * T2V); |
2652 | 0 | } |
2653 | 0 | T2X = T2R + T2W; |
2654 | 0 | Tdj = T7o + T7p; |
2655 | 0 | T7n = T2R - T2W; |
2656 | 0 | T7q = T7o - T7p; |
2657 | 0 | } |
2658 | 0 | T2B = T2p + T2A; |
2659 | 0 | T2Y = T2M + T2X; |
2660 | 0 | Tfz = T2B - T2Y; |
2661 | 0 | TfA = Tdd + Tde; |
2662 | 0 | TfB = Tdk + Tdj; |
2663 | 0 | TfC = TfA - TfB; |
2664 | 0 | { |
2665 | 0 | E T7d, T7g, Tdi, Tdl; |
2666 | 0 | T7d = T2j - T2o; |
2667 | 0 | T7g = T7e - T7f; |
2668 | 0 | T7h = T7d - T7g; |
2669 | 0 | TaW = T7d + T7g; |
2670 | 0 | Tdi = T2p - T2A; |
2671 | 0 | Tdl = Tdj - Tdk; |
2672 | 0 | Tdm = Tdi - Tdl; |
2673 | 0 | TeL = Tdi + Tdl; |
2674 | 0 | } |
2675 | 0 | { |
2676 | 0 | E T7m, T7r, T7z, T7A; |
2677 | 0 | T7m = T7i - T7l; |
2678 | 0 | T7r = T7n + T7q; |
2679 | 0 | T7s = KP707106781 * (T7m + T7r); |
2680 | 0 | TaU = KP707106781 * (T7r - T7m); |
2681 | 0 | T7z = T7i + T7l; |
2682 | 0 | T7A = T7q - T7n; |
2683 | 0 | T7B = KP707106781 * (T7z + T7A); |
2684 | 0 | TaX = KP707106781 * (T7z - T7A); |
2685 | 0 | } |
2686 | 0 | { |
2687 | 0 | E Tdf, Tdg, T7w, T7x; |
2688 | 0 | Tdf = Tdd - Tde; |
2689 | 0 | Tdg = T2M - T2X; |
2690 | 0 | Tdh = Tdf - Tdg; |
2691 | 0 | TeM = Tdf + Tdg; |
2692 | 0 | T7w = T7u - T7v; |
2693 | 0 | T7x = T2u - T2z; |
2694 | 0 | T7y = T7w + T7x; |
2695 | 0 | TaT = T7w - T7x; |
2696 | 0 | } |
2697 | 0 | } |
2698 | 0 | { |
2699 | 0 | E T4D, T9e, T4I, T9f, T4J, TdR, T4O, T8A, T4T, T8B, T4U, TdS, T56, Tea, T8E; |
2700 | 0 | E T8H, T5h, Te9, T8J, T8M; |
2701 | 0 | { |
2702 | 0 | E T4A, T4C, T4z, T4B; |
2703 | 0 | T4A = cr[WS(rs, 63)]; |
2704 | 0 | T4C = ci[WS(rs, 63)]; |
2705 | 0 | T4z = W[124]; |
2706 | 0 | T4B = W[125]; |
2707 | 0 | T4D = FMA(T4z, T4A, T4B * T4C); |
2708 | 0 | T9e = FNMS(T4B, T4A, T4z * T4C); |
2709 | 0 | } |
2710 | 0 | { |
2711 | 0 | E T4F, T4H, T4E, T4G; |
2712 | 0 | T4F = cr[WS(rs, 31)]; |
2713 | 0 | T4H = ci[WS(rs, 31)]; |
2714 | 0 | T4E = W[60]; |
2715 | 0 | T4G = W[61]; |
2716 | 0 | T4I = FMA(T4E, T4F, T4G * T4H); |
2717 | 0 | T9f = FNMS(T4G, T4F, T4E * T4H); |
2718 | 0 | } |
2719 | 0 | T4J = T4D + T4I; |
2720 | 0 | TdR = T9e + T9f; |
2721 | 0 | { |
2722 | 0 | E T4L, T4N, T4K, T4M; |
2723 | 0 | T4L = cr[WS(rs, 15)]; |
2724 | 0 | T4N = ci[WS(rs, 15)]; |
2725 | 0 | T4K = W[28]; |
2726 | 0 | T4M = W[29]; |
2727 | 0 | T4O = FMA(T4K, T4L, T4M * T4N); |
2728 | 0 | T8A = FNMS(T4M, T4L, T4K * T4N); |
2729 | 0 | } |
2730 | 0 | { |
2731 | 0 | E T4Q, T4S, T4P, T4R; |
2732 | 0 | T4Q = cr[WS(rs, 47)]; |
2733 | 0 | T4S = ci[WS(rs, 47)]; |
2734 | 0 | T4P = W[92]; |
2735 | 0 | T4R = W[93]; |
2736 | 0 | T4T = FMA(T4P, T4Q, T4R * T4S); |
2737 | 0 | T8B = FNMS(T4R, T4Q, T4P * T4S); |
2738 | 0 | } |
2739 | 0 | T4U = T4O + T4T; |
2740 | 0 | TdS = T8A + T8B; |
2741 | 0 | { |
2742 | 0 | E T50, T8F, T55, T8G; |
2743 | 0 | { |
2744 | 0 | E T4X, T4Z, T4W, T4Y; |
2745 | 0 | T4X = cr[WS(rs, 7)]; |
2746 | 0 | T4Z = ci[WS(rs, 7)]; |
2747 | 0 | T4W = W[12]; |
2748 | 0 | T4Y = W[13]; |
2749 | 0 | T50 = FMA(T4W, T4X, T4Y * T4Z); |
2750 | 0 | T8F = FNMS(T4Y, T4X, T4W * T4Z); |
2751 | 0 | } |
2752 | 0 | { |
2753 | 0 | E T52, T54, T51, T53; |
2754 | 0 | T52 = cr[WS(rs, 39)]; |
2755 | 0 | T54 = ci[WS(rs, 39)]; |
2756 | 0 | T51 = W[76]; |
2757 | 0 | T53 = W[77]; |
2758 | 0 | T55 = FMA(T51, T52, T53 * T54); |
2759 | 0 | T8G = FNMS(T53, T52, T51 * T54); |
2760 | 0 | } |
2761 | 0 | T56 = T50 + T55; |
2762 | 0 | Tea = T8F + T8G; |
2763 | 0 | T8E = T50 - T55; |
2764 | 0 | T8H = T8F - T8G; |
2765 | 0 | } |
2766 | 0 | { |
2767 | 0 | E T5b, T8K, T5g, T8L; |
2768 | 0 | { |
2769 | 0 | E T58, T5a, T57, T59; |
2770 | 0 | T58 = cr[WS(rs, 55)]; |
2771 | 0 | T5a = ci[WS(rs, 55)]; |
2772 | 0 | T57 = W[108]; |
2773 | 0 | T59 = W[109]; |
2774 | 0 | T5b = FMA(T57, T58, T59 * T5a); |
2775 | 0 | T8K = FNMS(T59, T58, T57 * T5a); |
2776 | 0 | } |
2777 | 0 | { |
2778 | 0 | E T5d, T5f, T5c, T5e; |
2779 | 0 | T5d = cr[WS(rs, 23)]; |
2780 | 0 | T5f = ci[WS(rs, 23)]; |
2781 | 0 | T5c = W[44]; |
2782 | 0 | T5e = W[45]; |
2783 | 0 | T5g = FMA(T5c, T5d, T5e * T5f); |
2784 | 0 | T8L = FNMS(T5e, T5d, T5c * T5f); |
2785 | 0 | } |
2786 | 0 | T5h = T5b + T5g; |
2787 | 0 | Te9 = T8K + T8L; |
2788 | 0 | T8J = T5b - T5g; |
2789 | 0 | T8M = T8K - T8L; |
2790 | 0 | } |
2791 | 0 | { |
2792 | 0 | E T4V, T5i, Te8, Teb; |
2793 | 0 | T4V = T4J + T4U; |
2794 | 0 | T5i = T56 + T5h; |
2795 | 0 | T5j = T4V + T5i; |
2796 | 0 | TfR = T4V - T5i; |
2797 | 0 | Te8 = T4J - T4U; |
2798 | 0 | Teb = Te9 - Tea; |
2799 | 0 | Tec = Te8 - Teb; |
2800 | 0 | TeX = Te8 + Teb; |
2801 | 0 | } |
2802 | 0 | { |
2803 | 0 | E TfW, TfX, T8z, T8C; |
2804 | 0 | TfW = TdR + TdS; |
2805 | 0 | TfX = Tea + Te9; |
2806 | 0 | TfY = TfW - TfX; |
2807 | 0 | Tgy = TfW + TfX; |
2808 | 0 | T8z = T4D - T4I; |
2809 | 0 | T8C = T8A - T8B; |
2810 | 0 | T8D = T8z - T8C; |
2811 | 0 | Tbl = T8z + T8C; |
2812 | 0 | } |
2813 | 0 | { |
2814 | 0 | E T8I, T8N, T9j, T9k; |
2815 | 0 | T8I = T8E - T8H; |
2816 | 0 | T8N = T8J + T8M; |
2817 | 0 | T8O = KP707106781 * (T8I + T8N); |
2818 | 0 | Tbx = KP707106781 * (T8N - T8I); |
2819 | 0 | T9j = T8E + T8H; |
2820 | 0 | T9k = T8M - T8J; |
2821 | 0 | T9l = KP707106781 * (T9j + T9k); |
2822 | 0 | Tbm = KP707106781 * (T9j - T9k); |
2823 | 0 | } |
2824 | 0 | { |
2825 | 0 | E TdT, TdU, T9g, T9h; |
2826 | 0 | TdT = TdR - TdS; |
2827 | 0 | TdU = T56 - T5h; |
2828 | 0 | TdV = TdT - TdU; |
2829 | 0 | Tf0 = TdT + TdU; |
2830 | 0 | T9g = T9e - T9f; |
2831 | 0 | T9h = T4O - T4T; |
2832 | 0 | T9i = T9g + T9h; |
2833 | 0 | Tbw = T9g - T9h; |
2834 | 0 | } |
2835 | 0 | } |
2836 | 0 | { |
2837 | 0 | E T36, T7G, T3b, T7H, T3c, TdH, T3h, T8m, T3m, T8n, T3n, TdI, T3z, Tds, T7L; |
2838 | 0 | E T7O, T3K, Tdr, T7S, T7T; |
2839 | 0 | { |
2840 | 0 | E T33, T35, T32, T34; |
2841 | 0 | T33 = cr[WS(rs, 1)]; |
2842 | 0 | T35 = ci[WS(rs, 1)]; |
2843 | 0 | T32 = W[0]; |
2844 | 0 | T34 = W[1]; |
2845 | 0 | T36 = FMA(T32, T33, T34 * T35); |
2846 | 0 | T7G = FNMS(T34, T33, T32 * T35); |
2847 | 0 | } |
2848 | 0 | { |
2849 | 0 | E T38, T3a, T37, T39; |
2850 | 0 | T38 = cr[WS(rs, 33)]; |
2851 | 0 | T3a = ci[WS(rs, 33)]; |
2852 | 0 | T37 = W[64]; |
2853 | 0 | T39 = W[65]; |
2854 | 0 | T3b = FMA(T37, T38, T39 * T3a); |
2855 | 0 | T7H = FNMS(T39, T38, T37 * T3a); |
2856 | 0 | } |
2857 | 0 | T3c = T36 + T3b; |
2858 | 0 | TdH = T7G + T7H; |
2859 | 0 | { |
2860 | 0 | E T3e, T3g, T3d, T3f; |
2861 | 0 | T3e = cr[WS(rs, 17)]; |
2862 | 0 | T3g = ci[WS(rs, 17)]; |
2863 | 0 | T3d = W[32]; |
2864 | 0 | T3f = W[33]; |
2865 | 0 | T3h = FMA(T3d, T3e, T3f * T3g); |
2866 | 0 | T8m = FNMS(T3f, T3e, T3d * T3g); |
2867 | 0 | } |
2868 | 0 | { |
2869 | 0 | E T3j, T3l, T3i, T3k; |
2870 | 0 | T3j = cr[WS(rs, 49)]; |
2871 | 0 | T3l = ci[WS(rs, 49)]; |
2872 | 0 | T3i = W[96]; |
2873 | 0 | T3k = W[97]; |
2874 | 0 | T3m = FMA(T3i, T3j, T3k * T3l); |
2875 | 0 | T8n = FNMS(T3k, T3j, T3i * T3l); |
2876 | 0 | } |
2877 | 0 | T3n = T3h + T3m; |
2878 | 0 | TdI = T8m + T8n; |
2879 | 0 | { |
2880 | 0 | E T3t, T7M, T3y, T7N; |
2881 | 0 | { |
2882 | 0 | E T3q, T3s, T3p, T3r; |
2883 | 0 | T3q = cr[WS(rs, 9)]; |
2884 | 0 | T3s = ci[WS(rs, 9)]; |
2885 | 0 | T3p = W[16]; |
2886 | 0 | T3r = W[17]; |
2887 | 0 | T3t = FMA(T3p, T3q, T3r * T3s); |
2888 | 0 | T7M = FNMS(T3r, T3q, T3p * T3s); |
2889 | 0 | } |
2890 | 0 | { |
2891 | 0 | E T3v, T3x, T3u, T3w; |
2892 | 0 | T3v = cr[WS(rs, 41)]; |
2893 | 0 | T3x = ci[WS(rs, 41)]; |
2894 | 0 | T3u = W[80]; |
2895 | 0 | T3w = W[81]; |
2896 | 0 | T3y = FMA(T3u, T3v, T3w * T3x); |
2897 | 0 | T7N = FNMS(T3w, T3v, T3u * T3x); |
2898 | 0 | } |
2899 | 0 | T3z = T3t + T3y; |
2900 | 0 | Tds = T7M + T7N; |
2901 | 0 | T7L = T3t - T3y; |
2902 | 0 | T7O = T7M - T7N; |
2903 | 0 | } |
2904 | 0 | { |
2905 | 0 | E T3E, T7Q, T3J, T7R; |
2906 | 0 | { |
2907 | 0 | E T3B, T3D, T3A, T3C; |
2908 | 0 | T3B = cr[WS(rs, 57)]; |
2909 | 0 | T3D = ci[WS(rs, 57)]; |
2910 | 0 | T3A = W[112]; |
2911 | 0 | T3C = W[113]; |
2912 | 0 | T3E = FMA(T3A, T3B, T3C * T3D); |
2913 | 0 | T7Q = FNMS(T3C, T3B, T3A * T3D); |
2914 | 0 | } |
2915 | 0 | { |
2916 | 0 | E T3G, T3I, T3F, T3H; |
2917 | 0 | T3G = cr[WS(rs, 25)]; |
2918 | 0 | T3I = ci[WS(rs, 25)]; |
2919 | 0 | T3F = W[48]; |
2920 | 0 | T3H = W[49]; |
2921 | 0 | T3J = FMA(T3F, T3G, T3H * T3I); |
2922 | 0 | T7R = FNMS(T3H, T3G, T3F * T3I); |
2923 | 0 | } |
2924 | 0 | T3K = T3E + T3J; |
2925 | 0 | Tdr = T7Q + T7R; |
2926 | 0 | T7S = T7Q - T7R; |
2927 | 0 | T7T = T3E - T3J; |
2928 | 0 | } |
2929 | 0 | { |
2930 | 0 | E T3o, T3L, TdJ, TdK; |
2931 | 0 | T3o = T3c + T3n; |
2932 | 0 | T3L = T3z + T3K; |
2933 | 0 | T3M = T3o + T3L; |
2934 | 0 | TfL = T3o - T3L; |
2935 | 0 | TdJ = TdH - TdI; |
2936 | 0 | TdK = T3z - T3K; |
2937 | 0 | TdL = TdJ - TdK; |
2938 | 0 | TeT = TdJ + TdK; |
2939 | 0 | } |
2940 | 0 | { |
2941 | 0 | E TfG, TfH, T7I, T7J; |
2942 | 0 | TfG = TdH + TdI; |
2943 | 0 | TfH = Tds + Tdr; |
2944 | 0 | TfI = TfG - TfH; |
2945 | 0 | Tgt = TfG + TfH; |
2946 | 0 | T7I = T7G - T7H; |
2947 | 0 | T7J = T3h - T3m; |
2948 | 0 | T7K = T7I + T7J; |
2949 | 0 | Tbd = T7I - T7J; |
2950 | 0 | } |
2951 | 0 | { |
2952 | 0 | E T7P, T7U, T8q, T8r; |
2953 | 0 | T7P = T7L + T7O; |
2954 | 0 | T7U = T7S - T7T; |
2955 | 0 | T7V = KP707106781 * (T7P + T7U); |
2956 | 0 | Tb3 = KP707106781 * (T7P - T7U); |
2957 | 0 | T8q = T7L - T7O; |
2958 | 0 | T8r = T7T + T7S; |
2959 | 0 | T8s = KP707106781 * (T8q + T8r); |
2960 | 0 | Tbe = KP707106781 * (T8r - T8q); |
2961 | 0 | } |
2962 | 0 | { |
2963 | 0 | E Tdq, Tdt, T8l, T8o; |
2964 | 0 | Tdq = T3c - T3n; |
2965 | 0 | Tdt = Tdr - Tds; |
2966 | 0 | Tdu = Tdq - Tdt; |
2967 | 0 | TeQ = Tdq + Tdt; |
2968 | 0 | T8l = T36 - T3b; |
2969 | 0 | T8o = T8m - T8n; |
2970 | 0 | T8p = T8l - T8o; |
2971 | 0 | Tb2 = T8l + T8o; |
2972 | 0 | } |
2973 | 0 | } |
2974 | 0 | { |
2975 | 0 | E T3X, Tdw, T7Z, T82, T4v, TdB, T8b, T8g, T48, Tdx, T80, T85, T4k, TdA, T8a; |
2976 | 0 | E T8d; |
2977 | 0 | { |
2978 | 0 | E T3R, T7X, T3W, T7Y; |
2979 | 0 | { |
2980 | 0 | E T3O, T3Q, T3N, T3P; |
2981 | 0 | T3O = cr[WS(rs, 5)]; |
2982 | 0 | T3Q = ci[WS(rs, 5)]; |
2983 | 0 | T3N = W[8]; |
2984 | 0 | T3P = W[9]; |
2985 | 0 | T3R = FMA(T3N, T3O, T3P * T3Q); |
2986 | 0 | T7X = FNMS(T3P, T3O, T3N * T3Q); |
2987 | 0 | } |
2988 | 0 | { |
2989 | 0 | E T3T, T3V, T3S, T3U; |
2990 | 0 | T3T = cr[WS(rs, 37)]; |
2991 | 0 | T3V = ci[WS(rs, 37)]; |
2992 | 0 | T3S = W[72]; |
2993 | 0 | T3U = W[73]; |
2994 | 0 | T3W = FMA(T3S, T3T, T3U * T3V); |
2995 | 0 | T7Y = FNMS(T3U, T3T, T3S * T3V); |
2996 | 0 | } |
2997 | 0 | T3X = T3R + T3W; |
2998 | 0 | Tdw = T7X + T7Y; |
2999 | 0 | T7Z = T7X - T7Y; |
3000 | 0 | T82 = T3R - T3W; |
3001 | 0 | } |
3002 | 0 | { |
3003 | 0 | E T4p, T8e, T4u, T8f; |
3004 | 0 | { |
3005 | 0 | E T4m, T4o, T4l, T4n; |
3006 | 0 | T4m = cr[WS(rs, 13)]; |
3007 | 0 | T4o = ci[WS(rs, 13)]; |
3008 | 0 | T4l = W[24]; |
3009 | 0 | T4n = W[25]; |
3010 | 0 | T4p = FMA(T4l, T4m, T4n * T4o); |
3011 | 0 | T8e = FNMS(T4n, T4m, T4l * T4o); |
3012 | 0 | } |
3013 | 0 | { |
3014 | 0 | E T4r, T4t, T4q, T4s; |
3015 | 0 | T4r = cr[WS(rs, 45)]; |
3016 | 0 | T4t = ci[WS(rs, 45)]; |
3017 | 0 | T4q = W[88]; |
3018 | 0 | T4s = W[89]; |
3019 | 0 | T4u = FMA(T4q, T4r, T4s * T4t); |
3020 | 0 | T8f = FNMS(T4s, T4r, T4q * T4t); |
3021 | 0 | } |
3022 | 0 | T4v = T4p + T4u; |
3023 | 0 | TdB = T8e + T8f; |
3024 | 0 | T8b = T4p - T4u; |
3025 | 0 | T8g = T8e - T8f; |
3026 | 0 | } |
3027 | 0 | { |
3028 | 0 | E T42, T83, T47, T84; |
3029 | 0 | { |
3030 | 0 | E T3Z, T41, T3Y, T40; |
3031 | 0 | T3Z = cr[WS(rs, 21)]; |
3032 | 0 | T41 = ci[WS(rs, 21)]; |
3033 | 0 | T3Y = W[40]; |
3034 | 0 | T40 = W[41]; |
3035 | 0 | T42 = FMA(T3Y, T3Z, T40 * T41); |
3036 | 0 | T83 = FNMS(T40, T3Z, T3Y * T41); |
3037 | 0 | } |
3038 | 0 | { |
3039 | 0 | E T44, T46, T43, T45; |
3040 | 0 | T44 = cr[WS(rs, 53)]; |
3041 | 0 | T46 = ci[WS(rs, 53)]; |
3042 | 0 | T43 = W[104]; |
3043 | 0 | T45 = W[105]; |
3044 | 0 | T47 = FMA(T43, T44, T45 * T46); |
3045 | 0 | T84 = FNMS(T45, T44, T43 * T46); |
3046 | 0 | } |
3047 | 0 | T48 = T42 + T47; |
3048 | 0 | Tdx = T83 + T84; |
3049 | 0 | T80 = T42 - T47; |
3050 | 0 | T85 = T83 - T84; |
3051 | 0 | } |
3052 | 0 | { |
3053 | 0 | E T4e, T88, T4j, T89; |
3054 | 0 | { |
3055 | 0 | E T4b, T4d, T4a, T4c; |
3056 | 0 | T4b = cr[WS(rs, 61)]; |
3057 | 0 | T4d = ci[WS(rs, 61)]; |
3058 | 0 | T4a = W[120]; |
3059 | 0 | T4c = W[121]; |
3060 | 0 | T4e = FMA(T4a, T4b, T4c * T4d); |
3061 | 0 | T88 = FNMS(T4c, T4b, T4a * T4d); |
3062 | 0 | } |
3063 | 0 | { |
3064 | 0 | E T4g, T4i, T4f, T4h; |
3065 | 0 | T4g = cr[WS(rs, 29)]; |
3066 | 0 | T4i = ci[WS(rs, 29)]; |
3067 | 0 | T4f = W[56]; |
3068 | 0 | T4h = W[57]; |
3069 | 0 | T4j = FMA(T4f, T4g, T4h * T4i); |
3070 | 0 | T89 = FNMS(T4h, T4g, T4f * T4i); |
3071 | 0 | } |
3072 | 0 | T4k = T4e + T4j; |
3073 | 0 | TdA = T88 + T89; |
3074 | 0 | T8a = T88 - T89; |
3075 | 0 | T8d = T4e - T4j; |
3076 | 0 | } |
3077 | 0 | { |
3078 | 0 | E T49, T4w, TdC, TdD; |
3079 | 0 | T49 = T3X + T48; |
3080 | 0 | T4w = T4k + T4v; |
3081 | 0 | T4x = T49 + T4w; |
3082 | 0 | TfJ = T49 - T4w; |
3083 | 0 | TdC = TdA - TdB; |
3084 | 0 | TdD = T4k - T4v; |
3085 | 0 | TdE = TdC - TdD; |
3086 | 0 | TdM = TdD + TdC; |
3087 | 0 | } |
3088 | 0 | { |
3089 | 0 | E TfM, TfN, T81, T86; |
3090 | 0 | TfM = TdA + TdB; |
3091 | 0 | TfN = Tdw + Tdx; |
3092 | 0 | TfO = TfM - TfN; |
3093 | 0 | Tgu = TfN + TfM; |
3094 | 0 | T81 = T7Z + T80; |
3095 | 0 | T86 = T82 - T85; |
3096 | 0 | T87 = FMA(KP923879532, T81, KP382683432 * T86); |
3097 | 0 | T8u = FNMS(KP382683432, T81, KP923879532 * T86); |
3098 | 0 | } |
3099 | 0 | { |
3100 | 0 | E T8c, T8h, Tb8, Tb9; |
3101 | 0 | T8c = T8a + T8b; |
3102 | 0 | T8h = T8d - T8g; |
3103 | 0 | T8i = FNMS(KP382683432, T8h, KP923879532 * T8c); |
3104 | 0 | T8v = FMA(KP382683432, T8c, KP923879532 * T8h); |
3105 | 0 | Tb8 = T8d + T8g; |
3106 | 0 | Tb9 = T8a - T8b; |
3107 | 0 | Tba = FNMS(KP382683432, Tb9, KP923879532 * Tb8); |
3108 | 0 | Tbh = FMA(KP923879532, Tb9, KP382683432 * Tb8); |
3109 | 0 | } |
3110 | 0 | { |
3111 | 0 | E Tdv, Tdy, Tb5, Tb6; |
3112 | 0 | Tdv = T3X - T48; |
3113 | 0 | Tdy = Tdw - Tdx; |
3114 | 0 | Tdz = Tdv + Tdy; |
3115 | 0 | TdN = Tdv - Tdy; |
3116 | 0 | Tb5 = T7Z - T80; |
3117 | 0 | Tb6 = T82 + T85; |
3118 | 0 | Tb7 = FMA(KP382683432, Tb5, KP923879532 * Tb6); |
3119 | 0 | Tbg = FNMS(KP382683432, Tb6, KP923879532 * Tb5); |
3120 | 0 | } |
3121 | 0 | } |
3122 | 0 | { |
3123 | 0 | E T5u, Te2, T8Q, T8X, T62, TdY, T94, T99, T5F, Te3, T8T, T8Y, T5R, TdX, T93; |
3124 | 0 | E T96; |
3125 | 0 | { |
3126 | 0 | E T5o, T8V, T5t, T8W; |
3127 | 0 | { |
3128 | 0 | E T5l, T5n, T5k, T5m; |
3129 | 0 | T5l = cr[WS(rs, 3)]; |
3130 | 0 | T5n = ci[WS(rs, 3)]; |
3131 | 0 | T5k = W[4]; |
3132 | 0 | T5m = W[5]; |
3133 | 0 | T5o = FMA(T5k, T5l, T5m * T5n); |
3134 | 0 | T8V = FNMS(T5m, T5l, T5k * T5n); |
3135 | 0 | } |
3136 | 0 | { |
3137 | 0 | E T5q, T5s, T5p, T5r; |
3138 | 0 | T5q = cr[WS(rs, 35)]; |
3139 | 0 | T5s = ci[WS(rs, 35)]; |
3140 | 0 | T5p = W[68]; |
3141 | 0 | T5r = W[69]; |
3142 | 0 | T5t = FMA(T5p, T5q, T5r * T5s); |
3143 | 0 | T8W = FNMS(T5r, T5q, T5p * T5s); |
3144 | 0 | } |
3145 | 0 | T5u = T5o + T5t; |
3146 | 0 | Te2 = T8V + T8W; |
3147 | 0 | T8Q = T5o - T5t; |
3148 | 0 | T8X = T8V - T8W; |
3149 | 0 | } |
3150 | 0 | { |
3151 | 0 | E T5W, T97, T61, T98; |
3152 | 0 | { |
3153 | 0 | E T5T, T5V, T5S, T5U; |
3154 | 0 | T5T = cr[WS(rs, 11)]; |
3155 | 0 | T5V = ci[WS(rs, 11)]; |
3156 | 0 | T5S = W[20]; |
3157 | 0 | T5U = W[21]; |
3158 | 0 | T5W = FMA(T5S, T5T, T5U * T5V); |
3159 | 0 | T97 = FNMS(T5U, T5T, T5S * T5V); |
3160 | 0 | } |
3161 | 0 | { |
3162 | 0 | E T5Y, T60, T5X, T5Z; |
3163 | 0 | T5Y = cr[WS(rs, 43)]; |
3164 | 0 | T60 = ci[WS(rs, 43)]; |
3165 | 0 | T5X = W[84]; |
3166 | 0 | T5Z = W[85]; |
3167 | 0 | T61 = FMA(T5X, T5Y, T5Z * T60); |
3168 | 0 | T98 = FNMS(T5Z, T5Y, T5X * T60); |
3169 | 0 | } |
3170 | 0 | T62 = T5W + T61; |
3171 | 0 | TdY = T97 + T98; |
3172 | 0 | T94 = T5W - T61; |
3173 | 0 | T99 = T97 - T98; |
3174 | 0 | } |
3175 | 0 | { |
3176 | 0 | E T5z, T8R, T5E, T8S; |
3177 | 0 | { |
3178 | 0 | E T5w, T5y, T5v, T5x; |
3179 | 0 | T5w = cr[WS(rs, 19)]; |
3180 | 0 | T5y = ci[WS(rs, 19)]; |
3181 | 0 | T5v = W[36]; |
3182 | 0 | T5x = W[37]; |
3183 | 0 | T5z = FMA(T5v, T5w, T5x * T5y); |
3184 | 0 | T8R = FNMS(T5x, T5w, T5v * T5y); |
3185 | 0 | } |
3186 | 0 | { |
3187 | 0 | E T5B, T5D, T5A, T5C; |
3188 | 0 | T5B = cr[WS(rs, 51)]; |
3189 | 0 | T5D = ci[WS(rs, 51)]; |
3190 | 0 | T5A = W[100]; |
3191 | 0 | T5C = W[101]; |
3192 | 0 | T5E = FMA(T5A, T5B, T5C * T5D); |
3193 | 0 | T8S = FNMS(T5C, T5B, T5A * T5D); |
3194 | 0 | } |
3195 | 0 | T5F = T5z + T5E; |
3196 | 0 | Te3 = T8R + T8S; |
3197 | 0 | T8T = T8R - T8S; |
3198 | 0 | T8Y = T5z - T5E; |
3199 | 0 | } |
3200 | 0 | { |
3201 | 0 | E T5L, T91, T5Q, T92; |
3202 | 0 | { |
3203 | 0 | E T5I, T5K, T5H, T5J; |
3204 | 0 | T5I = cr[WS(rs, 59)]; |
3205 | 0 | T5K = ci[WS(rs, 59)]; |
3206 | 0 | T5H = W[116]; |
3207 | 0 | T5J = W[117]; |
3208 | 0 | T5L = FMA(T5H, T5I, T5J * T5K); |
3209 | 0 | T91 = FNMS(T5J, T5I, T5H * T5K); |
3210 | 0 | } |
3211 | 0 | { |
3212 | 0 | E T5N, T5P, T5M, T5O; |
3213 | 0 | T5N = cr[WS(rs, 27)]; |
3214 | 0 | T5P = ci[WS(rs, 27)]; |
3215 | 0 | T5M = W[52]; |
3216 | 0 | T5O = W[53]; |
3217 | 0 | T5Q = FMA(T5M, T5N, T5O * T5P); |
3218 | 0 | T92 = FNMS(T5O, T5N, T5M * T5P); |
3219 | 0 | } |
3220 | 0 | T5R = T5L + T5Q; |
3221 | 0 | TdX = T91 + T92; |
3222 | 0 | T93 = T91 - T92; |
3223 | 0 | T96 = T5L - T5Q; |
3224 | 0 | } |
3225 | 0 | { |
3226 | 0 | E T5G, T63, Te1, Te4; |
3227 | 0 | T5G = T5u + T5F; |
3228 | 0 | T63 = T5R + T62; |
3229 | 0 | T64 = T5G + T63; |
3230 | 0 | TfZ = T5G - T63; |
3231 | 0 | Te1 = T5u - T5F; |
3232 | 0 | Te4 = Te2 - Te3; |
3233 | 0 | Te5 = Te1 - Te4; |
3234 | 0 | Ted = Te1 + Te4; |
3235 | 0 | } |
3236 | 0 | { |
3237 | 0 | E TfS, TfT, T8U, T8Z; |
3238 | 0 | TfS = TdX + TdY; |
3239 | 0 | TfT = Te2 + Te3; |
3240 | 0 | TfU = TfS - TfT; |
3241 | 0 | Tgz = TfT + TfS; |
3242 | 0 | T8U = T8Q - T8T; |
3243 | 0 | T8Z = T8X + T8Y; |
3244 | 0 | T90 = FNMS(KP382683432, T8Z, KP923879532 * T8U); |
3245 | 0 | T9n = FMA(KP923879532, T8Z, KP382683432 * T8U); |
3246 | 0 | } |
3247 | 0 | { |
3248 | 0 | E T95, T9a, Tbr, Tbs; |
3249 | 0 | T95 = T93 + T94; |
3250 | 0 | T9a = T96 - T99; |
3251 | 0 | T9b = FMA(KP382683432, T95, KP923879532 * T9a); |
3252 | 0 | T9o = FNMS(KP382683432, T9a, KP923879532 * T95); |
3253 | 0 | Tbr = T96 + T99; |
3254 | 0 | Tbs = T93 - T94; |
3255 | 0 | Tbt = FNMS(KP382683432, Tbs, KP923879532 * Tbr); |
3256 | 0 | TbA = FMA(KP923879532, Tbs, KP382683432 * Tbr); |
3257 | 0 | } |
3258 | 0 | { |
3259 | 0 | E TdW, TdZ, Tbo, Tbp; |
3260 | 0 | TdW = T5R - T62; |
3261 | 0 | TdZ = TdX - TdY; |
3262 | 0 | Te0 = TdW + TdZ; |
3263 | 0 | Tee = TdZ - TdW; |
3264 | 0 | Tbo = T8X - T8Y; |
3265 | 0 | Tbp = T8Q + T8T; |
3266 | 0 | Tbq = FMA(KP382683432, Tbo, KP923879532 * Tbp); |
3267 | 0 | Tbz = FNMS(KP382683432, Tbp, KP923879532 * Tbo); |
3268 | 0 | } |
3269 | 0 | } |
3270 | 0 | { |
3271 | 0 | E T1t, Tgn, TgK, TgL, TgV, Th1, T30, Th0, T66, TgX, Tgw, TgE, TgB, TgF, Tgq; |
3272 | 0 | E TgM; |
3273 | 0 | { |
3274 | 0 | E TH, T1s, TgI, TgJ; |
3275 | 0 | TH = Tj + TG; |
3276 | 0 | T1s = T14 + T1r; |
3277 | 0 | T1t = TH + T1s; |
3278 | 0 | Tgn = TH - T1s; |
3279 | 0 | TgI = Tgy + Tgz; |
3280 | 0 | TgJ = Tgt + Tgu; |
3281 | 0 | TgK = TgI - TgJ; |
3282 | 0 | TgL = TgJ + TgI; |
3283 | 0 | } |
3284 | 0 | { |
3285 | 0 | E TgN, TgU, T2e, T2Z; |
3286 | 0 | TgN = Tfr + Tfq; |
3287 | 0 | TgU = TgO + TgT; |
3288 | 0 | TgV = TgN + TgU; |
3289 | 0 | Th1 = TgU - TgN; |
3290 | 0 | T2e = T1Q + T2d; |
3291 | 0 | T2Z = T2B + T2Y; |
3292 | 0 | T30 = T2e + T2Z; |
3293 | 0 | Th0 = T2e - T2Z; |
3294 | 0 | } |
3295 | 0 | { |
3296 | 0 | E T4y, T65, Tgs, Tgv; |
3297 | 0 | T4y = T3M + T4x; |
3298 | 0 | T65 = T5j + T64; |
3299 | 0 | T66 = T4y + T65; |
3300 | 0 | TgX = T65 - T4y; |
3301 | 0 | Tgs = T3M - T4x; |
3302 | 0 | Tgv = Tgt - Tgu; |
3303 | 0 | Tgw = Tgs + Tgv; |
3304 | 0 | TgE = Tgs - Tgv; |
3305 | 0 | } |
3306 | 0 | { |
3307 | 0 | E Tgx, TgA, Tgo, Tgp; |
3308 | 0 | Tgx = T5j - T64; |
3309 | 0 | TgA = Tgy - Tgz; |
3310 | 0 | TgB = Tgx - TgA; |
3311 | 0 | TgF = Tgx + TgA; |
3312 | 0 | Tgo = TfA + TfB; |
3313 | 0 | Tgp = Tfv + Tfw; |
3314 | 0 | Tgq = Tgo - Tgp; |
3315 | 0 | TgM = Tgp + Tgo; |
3316 | 0 | } |
3317 | 0 | { |
3318 | 0 | E T31, TgW, TgY, TgH; |
3319 | 0 | T31 = T1t + T30; |
3320 | 0 | ci[WS(rs, 31)] = T31 - T66; |
3321 | 0 | cr[0] = T31 + T66; |
3322 | 0 | TgW = TgM + TgV; |
3323 | 0 | cr[WS(rs, 32)] = TgL - TgW; |
3324 | 0 | ci[WS(rs, 63)] = TgL + TgW; |
3325 | 0 | TgY = TgV - TgM; |
3326 | 0 | cr[WS(rs, 48)] = TgX - TgY; |
3327 | 0 | ci[WS(rs, 47)] = TgX + TgY; |
3328 | 0 | TgH = T1t - T30; |
3329 | 0 | cr[WS(rs, 16)] = TgH - TgK; |
3330 | 0 | ci[WS(rs, 15)] = TgH + TgK; |
3331 | 0 | } |
3332 | 0 | { |
3333 | 0 | E Tgr, TgC, TgZ, Th2; |
3334 | 0 | Tgr = Tgn - Tgq; |
3335 | 0 | TgC = KP707106781 * (Tgw + TgB); |
3336 | 0 | ci[WS(rs, 23)] = Tgr - TgC; |
3337 | 0 | cr[WS(rs, 8)] = Tgr + TgC; |
3338 | 0 | TgZ = KP707106781 * (TgB - Tgw); |
3339 | 0 | Th2 = Th0 + Th1; |
3340 | 0 | cr[WS(rs, 56)] = TgZ - Th2; |
3341 | 0 | ci[WS(rs, 39)] = TgZ + Th2; |
3342 | 0 | } |
3343 | 0 | { |
3344 | 0 | E Th3, Th4, TgD, TgG; |
3345 | 0 | Th3 = KP707106781 * (TgF - TgE); |
3346 | 0 | Th4 = Th1 - Th0; |
3347 | 0 | cr[WS(rs, 40)] = Th3 - Th4; |
3348 | 0 | ci[WS(rs, 55)] = Th3 + Th4; |
3349 | 0 | TgD = Tgn + Tgq; |
3350 | 0 | TgG = KP707106781 * (TgE + TgF); |
3351 | 0 | cr[WS(rs, 24)] = TgD - TgG; |
3352 | 0 | ci[WS(rs, 7)] = TgD + TgG; |
3353 | 0 | } |
3354 | 0 | } |
3355 | 0 | { |
3356 | 0 | E T6L, T9x, ThV, Ti1, T7E, Ti0, T9A, ThO, T8y, T9K, T9u, T9E, T9r, T9L, T9v; |
3357 | 0 | E T9H; |
3358 | 0 | { |
3359 | 0 | E T6n, T6K, ThP, ThU; |
3360 | 0 | T6n = T6b + T6m; |
3361 | 0 | T6K = T6y + T6J; |
3362 | 0 | T6L = T6n - T6K; |
3363 | 0 | T9x = T6n + T6K; |
3364 | 0 | ThP = T9O - T9P; |
3365 | 0 | ThU = ThQ + ThT; |
3366 | 0 | ThV = ThP + ThU; |
3367 | 0 | Ti1 = ThU - ThP; |
3368 | 0 | } |
3369 | 0 | { |
3370 | 0 | E T7c, T9y, T7D, T9z; |
3371 | 0 | { |
3372 | 0 | E T72, T7b, T7t, T7C; |
3373 | 0 | T72 = T6Q + T71; |
3374 | 0 | T7b = T77 + T7a; |
3375 | 0 | T7c = FMA(KP195090322, T72, KP980785280 * T7b); |
3376 | 0 | T9y = FNMS(KP195090322, T7b, KP980785280 * T72); |
3377 | 0 | T7t = T7h + T7s; |
3378 | 0 | T7C = T7y + T7B; |
3379 | 0 | T7D = FNMS(KP980785280, T7C, KP195090322 * T7t); |
3380 | 0 | T9z = FMA(KP980785280, T7t, KP195090322 * T7C); |
3381 | 0 | } |
3382 | 0 | T7E = T7c + T7D; |
3383 | 0 | Ti0 = T9z - T9y; |
3384 | 0 | T9A = T9y + T9z; |
3385 | 0 | ThO = T7c - T7D; |
3386 | 0 | } |
3387 | 0 | { |
3388 | 0 | E T8k, T9D, T8x, T9C; |
3389 | 0 | { |
3390 | 0 | E T7W, T8j, T8t, T8w; |
3391 | 0 | T7W = T7K + T7V; |
3392 | 0 | T8j = T87 + T8i; |
3393 | 0 | T8k = T7W - T8j; |
3394 | 0 | T9D = T7W + T8j; |
3395 | 0 | T8t = T8p + T8s; |
3396 | 0 | T8w = T8u + T8v; |
3397 | 0 | T8x = T8t - T8w; |
3398 | 0 | T9C = T8t + T8w; |
3399 | 0 | } |
3400 | 0 | T8y = FMA(KP634393284, T8k, KP773010453 * T8x); |
3401 | 0 | T9K = FMA(KP995184726, T9D, KP098017140 * T9C); |
3402 | 0 | T9u = FNMS(KP773010453, T8k, KP634393284 * T8x); |
3403 | 0 | T9E = FNMS(KP098017140, T9D, KP995184726 * T9C); |
3404 | 0 | } |
3405 | 0 | { |
3406 | 0 | E T9d, T9G, T9q, T9F; |
3407 | 0 | { |
3408 | 0 | E T8P, T9c, T9m, T9p; |
3409 | 0 | T8P = T8D + T8O; |
3410 | 0 | T9c = T90 + T9b; |
3411 | 0 | T9d = T8P - T9c; |
3412 | 0 | T9G = T8P + T9c; |
3413 | 0 | T9m = T9i + T9l; |
3414 | 0 | T9p = T9n + T9o; |
3415 | 0 | T9q = T9m - T9p; |
3416 | 0 | T9F = T9m + T9p; |
3417 | 0 | } |
3418 | 0 | T9r = FNMS(KP634393284, T9q, KP773010453 * T9d); |
3419 | 0 | T9L = FNMS(KP995184726, T9F, KP098017140 * T9G); |
3420 | 0 | T9v = FMA(KP773010453, T9q, KP634393284 * T9d); |
3421 | 0 | T9H = FMA(KP098017140, T9F, KP995184726 * T9G); |
3422 | 0 | } |
3423 | 0 | { |
3424 | 0 | E T7F, T9s, ThZ, Ti2; |
3425 | 0 | T7F = T6L + T7E; |
3426 | 0 | T9s = T8y + T9r; |
3427 | 0 | ci[WS(rs, 24)] = T7F - T9s; |
3428 | 0 | cr[WS(rs, 7)] = T7F + T9s; |
3429 | 0 | ThZ = T9v - T9u; |
3430 | 0 | Ti2 = Ti0 + Ti1; |
3431 | 0 | cr[WS(rs, 39)] = ThZ - Ti2; |
3432 | 0 | ci[WS(rs, 56)] = ThZ + Ti2; |
3433 | 0 | } |
3434 | 0 | { |
3435 | 0 | E Ti3, Ti4, T9t, T9w; |
3436 | 0 | Ti3 = T9r - T8y; |
3437 | 0 | Ti4 = Ti1 - Ti0; |
3438 | 0 | cr[WS(rs, 55)] = Ti3 - Ti4; |
3439 | 0 | ci[WS(rs, 40)] = Ti3 + Ti4; |
3440 | 0 | T9t = T6L - T7E; |
3441 | 0 | T9w = T9u + T9v; |
3442 | 0 | cr[WS(rs, 23)] = T9t - T9w; |
3443 | 0 | ci[WS(rs, 8)] = T9t + T9w; |
3444 | 0 | } |
3445 | 0 | { |
3446 | 0 | E T9B, T9I, ThN, ThW; |
3447 | 0 | T9B = T9x + T9A; |
3448 | 0 | T9I = T9E + T9H; |
3449 | 0 | cr[WS(rs, 31)] = T9B - T9I; |
3450 | 0 | ci[0] = T9B + T9I; |
3451 | 0 | ThN = T9L - T9K; |
3452 | 0 | ThW = ThO + ThV; |
3453 | 0 | cr[WS(rs, 63)] = ThN - ThW; |
3454 | 0 | ci[WS(rs, 32)] = ThN + ThW; |
3455 | 0 | } |
3456 | 0 | { |
3457 | 0 | E ThX, ThY, T9J, T9M; |
3458 | 0 | ThX = T9H - T9E; |
3459 | 0 | ThY = ThV - ThO; |
3460 | 0 | cr[WS(rs, 47)] = ThX - ThY; |
3461 | 0 | ci[WS(rs, 48)] = ThX + ThY; |
3462 | 0 | T9J = T9x - T9A; |
3463 | 0 | T9M = T9K + T9L; |
3464 | 0 | ci[WS(rs, 16)] = T9J - T9M; |
3465 | 0 | cr[WS(rs, 15)] = T9J + T9M; |
3466 | 0 | } |
3467 | 0 | } |
3468 | 0 | { |
3469 | 0 | E Tft, Tg7, Tgh, Tgl, Th9, Thf, TfE, Th6, TfQ, Tg4, Tga, The, Tge, Tgk, Tg1; |
3470 | 0 | E Tg5; |
3471 | 0 | { |
3472 | 0 | E Tfp, Tfs, Tgf, Tgg; |
3473 | 0 | Tfp = Tj - TG; |
3474 | 0 | Tfs = Tfq - Tfr; |
3475 | 0 | Tft = Tfp - Tfs; |
3476 | 0 | Tg7 = Tfp + Tfs; |
3477 | 0 | Tgf = TfY + TfZ; |
3478 | 0 | Tgg = TfR + TfU; |
3479 | 0 | Tgh = FMA(KP382683432, Tgf, KP923879532 * Tgg); |
3480 | 0 | Tgl = FNMS(KP923879532, Tgf, KP382683432 * Tgg); |
3481 | 0 | } |
3482 | 0 | { |
3483 | 0 | E Th7, Th8, Tfy, TfD; |
3484 | 0 | Th7 = T14 - T1r; |
3485 | 0 | Th8 = TgT - TgO; |
3486 | 0 | Th9 = Th7 + Th8; |
3487 | 0 | Thf = Th8 - Th7; |
3488 | 0 | Tfy = Tfu + Tfx; |
3489 | 0 | TfD = Tfz - TfC; |
3490 | 0 | TfE = KP707106781 * (Tfy + TfD); |
3491 | 0 | Th6 = KP707106781 * (Tfy - TfD); |
3492 | 0 | } |
3493 | 0 | { |
3494 | 0 | E TfK, TfP, Tg8, Tg9; |
3495 | 0 | TfK = TfI - TfJ; |
3496 | 0 | TfP = TfL - TfO; |
3497 | 0 | TfQ = FMA(KP382683432, TfK, KP923879532 * TfP); |
3498 | 0 | Tg4 = FNMS(KP923879532, TfK, KP382683432 * TfP); |
3499 | 0 | Tg8 = Tfu - Tfx; |
3500 | 0 | Tg9 = Tfz + TfC; |
3501 | 0 | Tga = KP707106781 * (Tg8 + Tg9); |
3502 | 0 | The = KP707106781 * (Tg9 - Tg8); |
3503 | 0 | } |
3504 | 0 | { |
3505 | 0 | E Tgc, Tgd, TfV, Tg0; |
3506 | 0 | Tgc = TfL + TfO; |
3507 | 0 | Tgd = TfI + TfJ; |
3508 | 0 | Tge = FNMS(KP382683432, Tgd, KP923879532 * Tgc); |
3509 | 0 | Tgk = FMA(KP923879532, Tgd, KP382683432 * Tgc); |
3510 | 0 | TfV = TfR - TfU; |
3511 | 0 | Tg0 = TfY - TfZ; |
3512 | 0 | Tg1 = FNMS(KP382683432, Tg0, KP923879532 * TfV); |
3513 | 0 | Tg5 = FMA(KP923879532, Tg0, KP382683432 * TfV); |
3514 | 0 | } |
3515 | 0 | { |
3516 | 0 | E TfF, Tg2, Thd, Thg; |
3517 | 0 | TfF = Tft + TfE; |
3518 | 0 | Tg2 = TfQ + Tg1; |
3519 | 0 | ci[WS(rs, 27)] = TfF - Tg2; |
3520 | 0 | cr[WS(rs, 4)] = TfF + Tg2; |
3521 | 0 | Thd = Tg5 - Tg4; |
3522 | 0 | Thg = The + Thf; |
3523 | 0 | cr[WS(rs, 36)] = Thd - Thg; |
3524 | 0 | ci[WS(rs, 59)] = Thd + Thg; |
3525 | 0 | } |
3526 | 0 | { |
3527 | 0 | E Thh, Thi, Tg3, Tg6; |
3528 | 0 | Thh = Tg1 - TfQ; |
3529 | 0 | Thi = Thf - The; |
3530 | 0 | cr[WS(rs, 52)] = Thh - Thi; |
3531 | 0 | ci[WS(rs, 43)] = Thh + Thi; |
3532 | 0 | Tg3 = Tft - TfE; |
3533 | 0 | Tg6 = Tg4 + Tg5; |
3534 | 0 | cr[WS(rs, 20)] = Tg3 - Tg6; |
3535 | 0 | ci[WS(rs, 11)] = Tg3 + Tg6; |
3536 | 0 | } |
3537 | 0 | { |
3538 | 0 | E Tgb, Tgi, Th5, Tha; |
3539 | 0 | Tgb = Tg7 + Tga; |
3540 | 0 | Tgi = Tge + Tgh; |
3541 | 0 | cr[WS(rs, 28)] = Tgb - Tgi; |
3542 | 0 | ci[WS(rs, 3)] = Tgb + Tgi; |
3543 | 0 | Th5 = Tgl - Tgk; |
3544 | 0 | Tha = Th6 + Th9; |
3545 | 0 | cr[WS(rs, 60)] = Th5 - Tha; |
3546 | 0 | ci[WS(rs, 35)] = Th5 + Tha; |
3547 | 0 | } |
3548 | 0 | { |
3549 | 0 | E Thb, Thc, Tgj, Tgm; |
3550 | 0 | Thb = Tgh - Tge; |
3551 | 0 | Thc = Th9 - Th6; |
3552 | 0 | cr[WS(rs, 44)] = Thb - Thc; |
3553 | 0 | ci[WS(rs, 51)] = Thb + Thc; |
3554 | 0 | Tgj = Tg7 - Tga; |
3555 | 0 | Tgm = Tgk + Tgl; |
3556 | 0 | ci[WS(rs, 19)] = Tgj - Tgm; |
3557 | 0 | cr[WS(rs, 12)] = Tgj + Tgm; |
3558 | 0 | } |
3559 | 0 | } |
3560 | 0 | { |
3561 | 0 | E TeH, Tf9, TeO, Thk, Thp, Thv, Tfc, Thu, Tf3, Tfn, Tf7, Tfj, TeW, Tfm, Tf6; |
3562 | 0 | E Tfg; |
3563 | 0 | { |
3564 | 0 | E TeD, TeG, Tfa, Tfb; |
3565 | 0 | TeD = TcL + TcO; |
3566 | 0 | TeG = KP707106781 * (TeE + TeF); |
3567 | 0 | TeH = TeD - TeG; |
3568 | 0 | Tf9 = TeD + TeG; |
3569 | 0 | { |
3570 | 0 | E TeK, TeN, Thl, Tho; |
3571 | 0 | TeK = FMA(KP923879532, TeI, KP382683432 * TeJ); |
3572 | 0 | TeN = FNMS(KP923879532, TeM, KP382683432 * TeL); |
3573 | 0 | TeO = TeK + TeN; |
3574 | 0 | Thk = TeK - TeN; |
3575 | 0 | Thl = KP707106781 * (TcU - TcZ); |
3576 | 0 | Tho = Thm + Thn; |
3577 | 0 | Thp = Thl + Tho; |
3578 | 0 | Thv = Tho - Thl; |
3579 | 0 | } |
3580 | 0 | Tfa = FNMS(KP382683432, TeI, KP923879532 * TeJ); |
3581 | 0 | Tfb = FMA(KP382683432, TeM, KP923879532 * TeL); |
3582 | 0 | Tfc = Tfa + Tfb; |
3583 | 0 | Thu = Tfb - Tfa; |
3584 | 0 | { |
3585 | 0 | E TeZ, Tfh, Tf2, Tfi, TeY, Tf1; |
3586 | 0 | TeY = KP707106781 * (Te5 + Te0); |
3587 | 0 | TeZ = TeX - TeY; |
3588 | 0 | Tfh = TeX + TeY; |
3589 | 0 | Tf1 = KP707106781 * (Ted + Tee); |
3590 | 0 | Tf2 = Tf0 - Tf1; |
3591 | 0 | Tfi = Tf0 + Tf1; |
3592 | 0 | Tf3 = FNMS(KP555570233, Tf2, KP831469612 * TeZ); |
3593 | 0 | Tfn = FMA(KP980785280, Tfh, KP195090322 * Tfi); |
3594 | 0 | Tf7 = FMA(KP555570233, TeZ, KP831469612 * Tf2); |
3595 | 0 | Tfj = FNMS(KP980785280, Tfi, KP195090322 * Tfh); |
3596 | 0 | } |
3597 | 0 | { |
3598 | 0 | E TeS, Tfe, TeV, Tff, TeR, TeU; |
3599 | 0 | TeR = KP707106781 * (TdN + TdM); |
3600 | 0 | TeS = TeQ - TeR; |
3601 | 0 | Tfe = TeQ + TeR; |
3602 | 0 | TeU = KP707106781 * (Tdz + TdE); |
3603 | 0 | TeV = TeT - TeU; |
3604 | 0 | Tff = TeT + TeU; |
3605 | 0 | TeW = FMA(KP831469612, TeS, KP555570233 * TeV); |
3606 | 0 | Tfm = FNMS(KP195090322, Tff, KP980785280 * Tfe); |
3607 | 0 | Tf6 = FNMS(KP831469612, TeV, KP555570233 * TeS); |
3608 | 0 | Tfg = FMA(KP195090322, Tfe, KP980785280 * Tff); |
3609 | 0 | } |
3610 | 0 | } |
3611 | 0 | { |
3612 | 0 | E TeP, Tf4, Tht, Thw; |
3613 | 0 | TeP = TeH + TeO; |
3614 | 0 | Tf4 = TeW + Tf3; |
3615 | 0 | ci[WS(rs, 25)] = TeP - Tf4; |
3616 | 0 | cr[WS(rs, 6)] = TeP + Tf4; |
3617 | 0 | Tht = Tf7 - Tf6; |
3618 | 0 | Thw = Thu + Thv; |
3619 | 0 | cr[WS(rs, 38)] = Tht - Thw; |
3620 | 0 | ci[WS(rs, 57)] = Tht + Thw; |
3621 | 0 | } |
3622 | 0 | { |
3623 | 0 | E Thx, Thy, Tf5, Tf8; |
3624 | 0 | Thx = Tf3 - TeW; |
3625 | 0 | Thy = Thv - Thu; |
3626 | 0 | cr[WS(rs, 54)] = Thx - Thy; |
3627 | 0 | ci[WS(rs, 41)] = Thx + Thy; |
3628 | 0 | Tf5 = TeH - TeO; |
3629 | 0 | Tf8 = Tf6 + Tf7; |
3630 | 0 | cr[WS(rs, 22)] = Tf5 - Tf8; |
3631 | 0 | ci[WS(rs, 9)] = Tf5 + Tf8; |
3632 | 0 | } |
3633 | 0 | { |
3634 | 0 | E Tfd, Tfk, Thj, Thq; |
3635 | 0 | Tfd = Tf9 - Tfc; |
3636 | 0 | Tfk = Tfg + Tfj; |
3637 | 0 | ci[WS(rs, 17)] = Tfd - Tfk; |
3638 | 0 | cr[WS(rs, 14)] = Tfd + Tfk; |
3639 | 0 | Thj = Tfj - Tfg; |
3640 | 0 | Thq = Thk + Thp; |
3641 | 0 | cr[WS(rs, 62)] = Thj - Thq; |
3642 | 0 | ci[WS(rs, 33)] = Thj + Thq; |
3643 | 0 | } |
3644 | 0 | { |
3645 | 0 | E Thr, Ths, Tfl, Tfo; |
3646 | 0 | Thr = Tfn - Tfm; |
3647 | 0 | Ths = Thp - Thk; |
3648 | 0 | cr[WS(rs, 46)] = Thr - Ths; |
3649 | 0 | ci[WS(rs, 49)] = Thr + Ths; |
3650 | 0 | Tfl = Tf9 + Tfc; |
3651 | 0 | Tfo = Tfm + Tfn; |
3652 | 0 | cr[WS(rs, 30)] = Tfl - Tfo; |
3653 | 0 | ci[WS(rs, 1)] = Tfl + Tfo; |
3654 | 0 | } |
3655 | 0 | } |
3656 | 0 | { |
3657 | 0 | E Td1, Ten, Tdo, ThA, ThD, ThJ, Teq, ThI, Teh, TeB, Tel, Tex, TdQ, TeA, Tek; |
3658 | 0 | E Teu; |
3659 | 0 | { |
3660 | 0 | E TcP, Td0, Teo, Tep; |
3661 | 0 | TcP = TcL - TcO; |
3662 | 0 | Td0 = KP707106781 * (TcU + TcZ); |
3663 | 0 | Td1 = TcP - Td0; |
3664 | 0 | Ten = TcP + Td0; |
3665 | 0 | { |
3666 | 0 | E Tdc, Tdn, ThB, ThC; |
3667 | 0 | Tdc = FNMS(KP923879532, Tdb, KP382683432 * Td6); |
3668 | 0 | Tdn = FMA(KP923879532, Tdh, KP382683432 * Tdm); |
3669 | 0 | Tdo = Tdc + Tdn; |
3670 | 0 | ThA = Tdn - Tdc; |
3671 | 0 | ThB = KP707106781 * (TeF - TeE); |
3672 | 0 | ThC = Thn - Thm; |
3673 | 0 | ThD = ThB + ThC; |
3674 | 0 | ThJ = ThC - ThB; |
3675 | 0 | } |
3676 | 0 | Teo = FMA(KP382683432, Tdb, KP923879532 * Td6); |
3677 | 0 | Tep = FNMS(KP382683432, Tdh, KP923879532 * Tdm); |
3678 | 0 | Teq = Teo + Tep; |
3679 | 0 | ThI = Teo - Tep; |
3680 | 0 | { |
3681 | 0 | E Te7, Tew, Teg, Tev, Te6, Tef; |
3682 | 0 | Te6 = KP707106781 * (Te0 - Te5); |
3683 | 0 | Te7 = TdV - Te6; |
3684 | 0 | Tew = TdV + Te6; |
3685 | 0 | Tef = KP707106781 * (Ted - Tee); |
3686 | 0 | Teg = Tec - Tef; |
3687 | 0 | Tev = Tec + Tef; |
3688 | 0 | Teh = FMA(KP555570233, Te7, KP831469612 * Teg); |
3689 | 0 | TeB = FMA(KP980785280, Tew, KP195090322 * Tev); |
3690 | 0 | Tel = FNMS(KP831469612, Te7, KP555570233 * Teg); |
3691 | 0 | Tex = FNMS(KP195090322, Tew, KP980785280 * Tev); |
3692 | 0 | } |
3693 | 0 | { |
3694 | 0 | E TdG, Tet, TdP, Tes, TdF, TdO; |
3695 | 0 | TdF = KP707106781 * (Tdz - TdE); |
3696 | 0 | TdG = Tdu - TdF; |
3697 | 0 | Tet = Tdu + TdF; |
3698 | 0 | TdO = KP707106781 * (TdM - TdN); |
3699 | 0 | TdP = TdL - TdO; |
3700 | 0 | Tes = TdL + TdO; |
3701 | 0 | TdQ = FNMS(KP555570233, TdP, KP831469612 * TdG); |
3702 | 0 | TeA = FNMS(KP980785280, Tes, KP195090322 * Tet); |
3703 | 0 | Tek = FMA(KP831469612, TdP, KP555570233 * TdG); |
3704 | 0 | Teu = FMA(KP195090322, Tes, KP980785280 * Tet); |
3705 | 0 | } |
3706 | 0 | } |
3707 | 0 | { |
3708 | 0 | E Tdp, Tei, ThH, ThK; |
3709 | 0 | Tdp = Td1 + Tdo; |
3710 | 0 | Tei = TdQ + Teh; |
3711 | 0 | cr[WS(rs, 26)] = Tdp - Tei; |
3712 | 0 | ci[WS(rs, 5)] = Tdp + Tei; |
3713 | 0 | ThH = Tel - Tek; |
3714 | 0 | ThK = ThI + ThJ; |
3715 | 0 | cr[WS(rs, 58)] = ThH - ThK; |
3716 | 0 | ci[WS(rs, 37)] = ThH + ThK; |
3717 | 0 | } |
3718 | 0 | { |
3719 | 0 | E ThL, ThM, Tej, Tem; |
3720 | 0 | ThL = Teh - TdQ; |
3721 | 0 | ThM = ThJ - ThI; |
3722 | 0 | cr[WS(rs, 42)] = ThL - ThM; |
3723 | 0 | ci[WS(rs, 53)] = ThL + ThM; |
3724 | 0 | Tej = Td1 - Tdo; |
3725 | 0 | Tem = Tek + Tel; |
3726 | 0 | ci[WS(rs, 21)] = Tej - Tem; |
3727 | 0 | cr[WS(rs, 10)] = Tej + Tem; |
3728 | 0 | } |
3729 | 0 | { |
3730 | 0 | E Ter, Tey, Thz, ThE; |
3731 | 0 | Ter = Ten + Teq; |
3732 | 0 | Tey = Teu + Tex; |
3733 | 0 | ci[WS(rs, 29)] = Ter - Tey; |
3734 | 0 | cr[WS(rs, 2)] = Ter + Tey; |
3735 | 0 | Thz = TeB - TeA; |
3736 | 0 | ThE = ThA + ThD; |
3737 | 0 | cr[WS(rs, 34)] = Thz - ThE; |
3738 | 0 | ci[WS(rs, 61)] = Thz + ThE; |
3739 | 0 | } |
3740 | 0 | { |
3741 | 0 | E ThF, ThG, Tez, TeC; |
3742 | 0 | ThF = Tex - Teu; |
3743 | 0 | ThG = ThD - ThA; |
3744 | 0 | cr[WS(rs, 50)] = ThF - ThG; |
3745 | 0 | ci[WS(rs, 45)] = ThF + ThG; |
3746 | 0 | Tez = Ten - Teq; |
3747 | 0 | TeC = TeA + TeB; |
3748 | 0 | cr[WS(rs, 18)] = Tez - TeC; |
3749 | 0 | ci[WS(rs, 13)] = Tez + TeC; |
3750 | 0 | } |
3751 | 0 | } |
3752 | 0 | { |
3753 | 0 | E Tc3, Tcv, TiD, TiJ, Tca, TiI, Tcy, TiA, Tci, TcI, Tcs, TcC, Tcp, TcJ, Tct; |
3754 | 0 | E TcF; |
3755 | 0 | { |
3756 | 0 | E TbZ, Tc2, TiB, TiC; |
3757 | 0 | TbZ = Taz - TaC; |
3758 | 0 | Tc2 = Tc0 + Tc1; |
3759 | 0 | Tc3 = TbZ - Tc2; |
3760 | 0 | Tcv = TbZ + Tc2; |
3761 | 0 | TiB = TaG - TaJ; |
3762 | 0 | TiC = Tin - Tim; |
3763 | 0 | TiD = TiB + TiC; |
3764 | 0 | TiJ = TiC - TiB; |
3765 | 0 | } |
3766 | 0 | { |
3767 | 0 | E Tc6, Tcw, Tc9, Tcx; |
3768 | 0 | { |
3769 | 0 | E Tc4, Tc5, Tc7, Tc8; |
3770 | 0 | Tc4 = TaP - TaQ; |
3771 | 0 | Tc5 = TaM - TaN; |
3772 | 0 | Tc6 = FMA(KP831469612, Tc4, KP555570233 * Tc5); |
3773 | 0 | Tcw = FNMS(KP555570233, Tc4, KP831469612 * Tc5); |
3774 | 0 | Tc7 = TaW - TaX; |
3775 | 0 | Tc8 = TaT - TaU; |
3776 | 0 | Tc9 = FNMS(KP831469612, Tc8, KP555570233 * Tc7); |
3777 | 0 | Tcx = FMA(KP555570233, Tc8, KP831469612 * Tc7); |
3778 | 0 | } |
3779 | 0 | Tca = Tc6 + Tc9; |
3780 | 0 | TiI = Tcx - Tcw; |
3781 | 0 | Tcy = Tcw + Tcx; |
3782 | 0 | TiA = Tc6 - Tc9; |
3783 | 0 | } |
3784 | 0 | { |
3785 | 0 | E Tce, TcB, Tch, TcA; |
3786 | 0 | { |
3787 | 0 | E Tcc, Tcd, Tcf, Tcg; |
3788 | 0 | Tcc = Tbd - Tbe; |
3789 | 0 | Tcd = Tb7 - Tba; |
3790 | 0 | Tce = Tcc - Tcd; |
3791 | 0 | TcB = Tcc + Tcd; |
3792 | 0 | Tcf = Tb2 - Tb3; |
3793 | 0 | Tcg = Tbh - Tbg; |
3794 | 0 | Tch = Tcf - Tcg; |
3795 | 0 | TcA = Tcf + Tcg; |
3796 | 0 | } |
3797 | 0 | Tci = FMA(KP471396736, Tce, KP881921264 * Tch); |
3798 | 0 | TcI = FMA(KP956940335, TcB, KP290284677 * TcA); |
3799 | 0 | Tcs = FNMS(KP881921264, Tce, KP471396736 * Tch); |
3800 | 0 | TcC = FNMS(KP290284677, TcB, KP956940335 * TcA); |
3801 | 0 | } |
3802 | 0 | { |
3803 | 0 | E Tcl, TcE, Tco, TcD; |
3804 | 0 | { |
3805 | 0 | E Tcj, Tck, Tcm, Tcn; |
3806 | 0 | Tcj = Tbl - Tbm; |
3807 | 0 | Tck = TbA - Tbz; |
3808 | 0 | Tcl = Tcj - Tck; |
3809 | 0 | TcE = Tcj + Tck; |
3810 | 0 | Tcm = Tbw - Tbx; |
3811 | 0 | Tcn = Tbq - Tbt; |
3812 | 0 | Tco = Tcm - Tcn; |
3813 | 0 | TcD = Tcm + Tcn; |
3814 | 0 | } |
3815 | 0 | Tcp = FNMS(KP471396736, Tco, KP881921264 * Tcl); |
3816 | 0 | TcJ = FNMS(KP956940335, TcD, KP290284677 * TcE); |
3817 | 0 | Tct = FMA(KP881921264, Tco, KP471396736 * Tcl); |
3818 | 0 | TcF = FMA(KP290284677, TcD, KP956940335 * TcE); |
3819 | 0 | } |
3820 | 0 | { |
3821 | 0 | E Tcb, Tcq, TiH, TiK; |
3822 | 0 | Tcb = Tc3 + Tca; |
3823 | 0 | Tcq = Tci + Tcp; |
3824 | 0 | ci[WS(rs, 26)] = Tcb - Tcq; |
3825 | 0 | cr[WS(rs, 5)] = Tcb + Tcq; |
3826 | 0 | TiH = Tct - Tcs; |
3827 | 0 | TiK = TiI + TiJ; |
3828 | 0 | cr[WS(rs, 37)] = TiH - TiK; |
3829 | 0 | ci[WS(rs, 58)] = TiH + TiK; |
3830 | 0 | } |
3831 | 0 | { |
3832 | 0 | E TiL, TiM, Tcr, Tcu; |
3833 | 0 | TiL = Tcp - Tci; |
3834 | 0 | TiM = TiJ - TiI; |
3835 | 0 | cr[WS(rs, 53)] = TiL - TiM; |
3836 | 0 | ci[WS(rs, 42)] = TiL + TiM; |
3837 | 0 | Tcr = Tc3 - Tca; |
3838 | 0 | Tcu = Tcs + Tct; |
3839 | 0 | cr[WS(rs, 21)] = Tcr - Tcu; |
3840 | 0 | ci[WS(rs, 10)] = Tcr + Tcu; |
3841 | 0 | } |
3842 | 0 | { |
3843 | 0 | E Tcz, TcG, Tiz, TiE; |
3844 | 0 | Tcz = Tcv + Tcy; |
3845 | 0 | TcG = TcC + TcF; |
3846 | 0 | cr[WS(rs, 29)] = Tcz - TcG; |
3847 | 0 | ci[WS(rs, 2)] = Tcz + TcG; |
3848 | 0 | Tiz = TcJ - TcI; |
3849 | 0 | TiE = TiA + TiD; |
3850 | 0 | cr[WS(rs, 61)] = Tiz - TiE; |
3851 | 0 | ci[WS(rs, 34)] = Tiz + TiE; |
3852 | 0 | } |
3853 | 0 | { |
3854 | 0 | E TiF, TiG, TcH, TcK; |
3855 | 0 | TiF = TcF - TcC; |
3856 | 0 | TiG = TiD - TiA; |
3857 | 0 | cr[WS(rs, 45)] = TiF - TiG; |
3858 | 0 | ci[WS(rs, 50)] = TiF + TiG; |
3859 | 0 | TcH = Tcv - Tcy; |
3860 | 0 | TcK = TcI + TcJ; |
3861 | 0 | ci[WS(rs, 18)] = TcH - TcK; |
3862 | 0 | cr[WS(rs, 13)] = TcH + TcK; |
3863 | 0 | } |
3864 | 0 | } |
3865 | 0 | { |
3866 | 0 | E TaL, TbJ, Tip, Tiv, Tb0, Tiu, TbM, Tik, Tbk, TbW, TbG, TbQ, TbD, TbX, TbH; |
3867 | 0 | E TbT; |
3868 | 0 | { |
3869 | 0 | E TaD, TaK, Til, Tio; |
3870 | 0 | TaD = Taz + TaC; |
3871 | 0 | TaK = TaG + TaJ; |
3872 | 0 | TaL = TaD - TaK; |
3873 | 0 | TbJ = TaD + TaK; |
3874 | 0 | Til = Tc1 - Tc0; |
3875 | 0 | Tio = Tim + Tin; |
3876 | 0 | Tip = Til + Tio; |
3877 | 0 | Tiv = Tio - Til; |
3878 | 0 | } |
3879 | 0 | { |
3880 | 0 | E TaS, TbK, TaZ, TbL; |
3881 | 0 | { |
3882 | 0 | E TaO, TaR, TaV, TaY; |
3883 | 0 | TaO = TaM + TaN; |
3884 | 0 | TaR = TaP + TaQ; |
3885 | 0 | TaS = FNMS(KP980785280, TaR, KP195090322 * TaO); |
3886 | 0 | TbK = FMA(KP195090322, TaR, KP980785280 * TaO); |
3887 | 0 | TaV = TaT + TaU; |
3888 | 0 | TaY = TaW + TaX; |
3889 | 0 | TaZ = FMA(KP980785280, TaV, KP195090322 * TaY); |
3890 | 0 | TbL = FNMS(KP195090322, TaV, KP980785280 * TaY); |
3891 | 0 | } |
3892 | 0 | Tb0 = TaS + TaZ; |
3893 | 0 | Tiu = TbK - TbL; |
3894 | 0 | TbM = TbK + TbL; |
3895 | 0 | Tik = TaZ - TaS; |
3896 | 0 | } |
3897 | 0 | { |
3898 | 0 | E Tbc, TbO, Tbj, TbP; |
3899 | 0 | { |
3900 | 0 | E Tb4, Tbb, Tbf, Tbi; |
3901 | 0 | Tb4 = Tb2 + Tb3; |
3902 | 0 | Tbb = Tb7 + Tba; |
3903 | 0 | Tbc = Tb4 - Tbb; |
3904 | 0 | TbO = Tb4 + Tbb; |
3905 | 0 | Tbf = Tbd + Tbe; |
3906 | 0 | Tbi = Tbg + Tbh; |
3907 | 0 | Tbj = Tbf - Tbi; |
3908 | 0 | TbP = Tbf + Tbi; |
3909 | 0 | } |
3910 | 0 | Tbk = FMA(KP634393284, Tbc, KP773010453 * Tbj); |
3911 | 0 | TbW = FNMS(KP995184726, TbP, KP098017140 * TbO); |
3912 | 0 | TbG = FNMS(KP634393284, Tbj, KP773010453 * Tbc); |
3913 | 0 | TbQ = FMA(KP995184726, TbO, KP098017140 * TbP); |
3914 | 0 | } |
3915 | 0 | { |
3916 | 0 | E Tbv, TbR, TbC, TbS; |
3917 | 0 | { |
3918 | 0 | E Tbn, Tbu, Tby, TbB; |
3919 | 0 | Tbn = Tbl + Tbm; |
3920 | 0 | Tbu = Tbq + Tbt; |
3921 | 0 | Tbv = Tbn - Tbu; |
3922 | 0 | TbR = Tbn + Tbu; |
3923 | 0 | Tby = Tbw + Tbx; |
3924 | 0 | TbB = Tbz + TbA; |
3925 | 0 | TbC = Tby - TbB; |
3926 | 0 | TbS = Tby + TbB; |
3927 | 0 | } |
3928 | 0 | TbD = FNMS(KP773010453, TbC, KP634393284 * Tbv); |
3929 | 0 | TbX = FMA(KP098017140, TbR, KP995184726 * TbS); |
3930 | 0 | TbH = FMA(KP773010453, Tbv, KP634393284 * TbC); |
3931 | 0 | TbT = FNMS(KP098017140, TbS, KP995184726 * TbR); |
3932 | 0 | } |
3933 | 0 | { |
3934 | 0 | E Tb1, TbE, Tit, Tiw; |
3935 | 0 | Tb1 = TaL - Tb0; |
3936 | 0 | TbE = Tbk + TbD; |
3937 | 0 | ci[WS(rs, 22)] = Tb1 - TbE; |
3938 | 0 | cr[WS(rs, 9)] = Tb1 + TbE; |
3939 | 0 | Tit = TbD - Tbk; |
3940 | 0 | Tiw = Tiu + Tiv; |
3941 | 0 | cr[WS(rs, 57)] = Tit - Tiw; |
3942 | 0 | ci[WS(rs, 38)] = Tit + Tiw; |
3943 | 0 | } |
3944 | 0 | { |
3945 | 0 | E Tix, Tiy, TbF, TbI; |
3946 | 0 | Tix = TbH - TbG; |
3947 | 0 | Tiy = Tiv - Tiu; |
3948 | 0 | cr[WS(rs, 41)] = Tix - Tiy; |
3949 | 0 | ci[WS(rs, 54)] = Tix + Tiy; |
3950 | 0 | TbF = TaL + Tb0; |
3951 | 0 | TbI = TbG + TbH; |
3952 | 0 | cr[WS(rs, 25)] = TbF - TbI; |
3953 | 0 | ci[WS(rs, 6)] = TbF + TbI; |
3954 | 0 | } |
3955 | 0 | { |
3956 | 0 | E TbN, TbU, Tij, Tiq; |
3957 | 0 | TbN = TbJ + TbM; |
3958 | 0 | TbU = TbQ + TbT; |
3959 | 0 | ci[WS(rs, 30)] = TbN - TbU; |
3960 | 0 | cr[WS(rs, 1)] = TbN + TbU; |
3961 | 0 | Tij = TbX - TbW; |
3962 | 0 | Tiq = Tik + Tip; |
3963 | 0 | cr[WS(rs, 33)] = Tij - Tiq; |
3964 | 0 | ci[WS(rs, 62)] = Tij + Tiq; |
3965 | 0 | } |
3966 | 0 | { |
3967 | 0 | E Tir, Tis, TbV, TbY; |
3968 | 0 | Tir = TbT - TbQ; |
3969 | 0 | Tis = Tip - Tik; |
3970 | 0 | cr[WS(rs, 49)] = Tir - Tis; |
3971 | 0 | ci[WS(rs, 46)] = Tir + Tis; |
3972 | 0 | TbV = TbJ - TbM; |
3973 | 0 | TbY = TbW + TbX; |
3974 | 0 | cr[WS(rs, 17)] = TbV - TbY; |
3975 | 0 | ci[WS(rs, 14)] = TbV + TbY; |
3976 | 0 | } |
3977 | 0 | } |
3978 | 0 | { |
3979 | 0 | E T9R, Taj, Ti9, Tif, T9Y, Tie, Tam, Ti6, Ta6, Taw, Tag, Taq, Tad, Tax, Tah; |
3980 | 0 | E Tat; |
3981 | 0 | { |
3982 | 0 | E T9N, T9Q, Ti7, Ti8; |
3983 | 0 | T9N = T6b - T6m; |
3984 | 0 | T9Q = T9O + T9P; |
3985 | 0 | T9R = T9N - T9Q; |
3986 | 0 | Taj = T9N + T9Q; |
3987 | 0 | Ti7 = T6J - T6y; |
3988 | 0 | Ti8 = ThT - ThQ; |
3989 | 0 | Ti9 = Ti7 + Ti8; |
3990 | 0 | Tif = Ti8 - Ti7; |
3991 | 0 | } |
3992 | 0 | { |
3993 | 0 | E T9U, Tak, T9X, Tal; |
3994 | 0 | { |
3995 | 0 | E T9S, T9T, T9V, T9W; |
3996 | 0 | T9S = T6Q - T71; |
3997 | 0 | T9T = T77 - T7a; |
3998 | 0 | T9U = FNMS(KP831469612, T9T, KP555570233 * T9S); |
3999 | 0 | Tak = FMA(KP831469612, T9S, KP555570233 * T9T); |
4000 | 0 | T9V = T7h - T7s; |
4001 | 0 | T9W = T7y - T7B; |
4002 | 0 | T9X = FMA(KP555570233, T9V, KP831469612 * T9W); |
4003 | 0 | Tal = FNMS(KP555570233, T9W, KP831469612 * T9V); |
4004 | 0 | } |
4005 | 0 | T9Y = T9U + T9X; |
4006 | 0 | Tie = Tak - Tal; |
4007 | 0 | Tam = Tak + Tal; |
4008 | 0 | Ti6 = T9X - T9U; |
4009 | 0 | } |
4010 | 0 | { |
4011 | 0 | E Ta2, Tao, Ta5, Tap; |
4012 | 0 | { |
4013 | 0 | E Ta0, Ta1, Ta3, Ta4; |
4014 | 0 | Ta0 = T8p - T8s; |
4015 | 0 | Ta1 = T87 - T8i; |
4016 | 0 | Ta2 = Ta0 - Ta1; |
4017 | 0 | Tao = Ta0 + Ta1; |
4018 | 0 | Ta3 = T7K - T7V; |
4019 | 0 | Ta4 = T8v - T8u; |
4020 | 0 | Ta5 = Ta3 - Ta4; |
4021 | 0 | Tap = Ta3 + Ta4; |
4022 | 0 | } |
4023 | 0 | Ta6 = FMA(KP471396736, Ta2, KP881921264 * Ta5); |
4024 | 0 | Taw = FNMS(KP956940335, Tap, KP290284677 * Tao); |
4025 | 0 | Tag = FNMS(KP471396736, Ta5, KP881921264 * Ta2); |
4026 | 0 | Taq = FMA(KP956940335, Tao, KP290284677 * Tap); |
4027 | 0 | } |
4028 | 0 | { |
4029 | 0 | E Ta9, Tar, Tac, Tas; |
4030 | 0 | { |
4031 | 0 | E Ta7, Ta8, Taa, Tab; |
4032 | 0 | Ta7 = T8D - T8O; |
4033 | 0 | Ta8 = T9n - T9o; |
4034 | 0 | Ta9 = Ta7 - Ta8; |
4035 | 0 | Tar = Ta7 + Ta8; |
4036 | 0 | Taa = T9i - T9l; |
4037 | 0 | Tab = T9b - T90; |
4038 | 0 | Tac = Taa - Tab; |
4039 | 0 | Tas = Taa + Tab; |
4040 | 0 | } |
4041 | 0 | Tad = FNMS(KP881921264, Tac, KP471396736 * Ta9); |
4042 | 0 | Tax = FMA(KP290284677, Tar, KP956940335 * Tas); |
4043 | 0 | Tah = FMA(KP881921264, Ta9, KP471396736 * Tac); |
4044 | 0 | Tat = FNMS(KP290284677, Tas, KP956940335 * Tar); |
4045 | 0 | } |
4046 | 0 | { |
4047 | 0 | E T9Z, Tae, Tid, Tig; |
4048 | 0 | T9Z = T9R - T9Y; |
4049 | 0 | Tae = Ta6 + Tad; |
4050 | 0 | ci[WS(rs, 20)] = T9Z - Tae; |
4051 | 0 | cr[WS(rs, 11)] = T9Z + Tae; |
4052 | 0 | Tid = Tad - Ta6; |
4053 | 0 | Tig = Tie + Tif; |
4054 | 0 | cr[WS(rs, 59)] = Tid - Tig; |
4055 | 0 | ci[WS(rs, 36)] = Tid + Tig; |
4056 | 0 | } |
4057 | 0 | { |
4058 | 0 | E Tih, Tii, Taf, Tai; |
4059 | 0 | Tih = Tah - Tag; |
4060 | 0 | Tii = Tif - Tie; |
4061 | 0 | cr[WS(rs, 43)] = Tih - Tii; |
4062 | 0 | ci[WS(rs, 52)] = Tih + Tii; |
4063 | 0 | Taf = T9R + T9Y; |
4064 | 0 | Tai = Tag + Tah; |
4065 | 0 | cr[WS(rs, 27)] = Taf - Tai; |
4066 | 0 | ci[WS(rs, 4)] = Taf + Tai; |
4067 | 0 | } |
4068 | 0 | { |
4069 | 0 | E Tan, Tau, Ti5, Tia; |
4070 | 0 | Tan = Taj + Tam; |
4071 | 0 | Tau = Taq + Tat; |
4072 | 0 | ci[WS(rs, 28)] = Tan - Tau; |
4073 | 0 | cr[WS(rs, 3)] = Tan + Tau; |
4074 | 0 | Ti5 = Tax - Taw; |
4075 | 0 | Tia = Ti6 + Ti9; |
4076 | 0 | cr[WS(rs, 35)] = Ti5 - Tia; |
4077 | 0 | ci[WS(rs, 60)] = Ti5 + Tia; |
4078 | 0 | } |
4079 | 0 | { |
4080 | 0 | E Tib, Tic, Tav, Tay; |
4081 | 0 | Tib = Tat - Taq; |
4082 | 0 | Tic = Ti9 - Ti6; |
4083 | 0 | cr[WS(rs, 51)] = Tib - Tic; |
4084 | 0 | ci[WS(rs, 44)] = Tib + Tic; |
4085 | 0 | Tav = Taj - Tam; |
4086 | 0 | Tay = Taw + Tax; |
4087 | 0 | cr[WS(rs, 19)] = Tav - Tay; |
4088 | 0 | ci[WS(rs, 12)] = Tav + Tay; |
4089 | 0 | } |
4090 | 0 | } |
4091 | 0 | } |
4092 | 0 | } |
4093 | 0 | } |
4094 | | |
4095 | | static const tw_instr twinstr[] = { |
4096 | | { TW_FULL, 1, 64 }, |
4097 | | { TW_NEXT, 1, 0 } |
4098 | | }; |
4099 | | |
4100 | | static const hc2hc_desc desc = { 64, "hf_64", twinstr, &GENUS, { 808, 270, 230, 0 } }; |
4101 | | |
4102 | 1 | void X(codelet_hf_64) (planner *p) { |
4103 | 1 | X(khc2hc_register) (p, hf_64, &desc); |
4104 | 1 | } |
4105 | | #endif |