diff options
Diffstat (limited to 'src/fftw3/rdft/codelets/hc2r/hb_64.c')
| -rw-r--r-- | src/fftw3/rdft/codelets/hc2r/hb_64.c | 1972 |
1 files changed, 1972 insertions, 0 deletions
diff --git a/src/fftw3/rdft/codelets/hc2r/hb_64.c b/src/fftw3/rdft/codelets/hc2r/hb_64.c new file mode 100644 index 0000000..d3f9afc --- /dev/null +++ b/src/fftw3/rdft/codelets/hc2r/hb_64.c @@ -0,0 +1,1972 @@ +/* + * Copyright (c) 2003 Matteo Frigo + * Copyright (c) 2003 Massachusetts Institute of Technology + * + * This program is free software; you can redistribute it and/or modify + * it under the terms of the GNU General Public License as published by + * the Free Software Foundation; either version 2 of the License, or + * (at your option) any later version. + * + * This program is distributed in the hope that it will be useful, + * but WITHOUT ANY WARRANTY; without even the implied warranty of + * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the + * GNU General Public License for more details. + * + * You should have received a copy of the GNU General Public License + * along with this program; if not, write to the Free Software + * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA + * + */ + +/* This file was automatically generated --- DO NOT EDIT */ +/* Generated on Sat Jul 5 22:11:42 EDT 2003 */ + +#include "codelet-rdft.h" + +/* Generated by: /homee/stevenj/cvs/fftw3.0.1/genfft/gen_hc2hc -compact -variables 4 -sign 1 -n 64 -dif -name hb_64 -include hb.h */ + +/* + * This function contains 1038 FP additions, 500 FP multiplications, + * (or, 808 additions, 270 multiplications, 230 fused multiply/add), + * 196 stack variables, and 256 memory accesses + */ +/* + * Generator Id's : + * $Id: hb_64.c,v 1.1 2008/10/17 06:12:08 scuri Exp $ + * $Id: hb_64.c,v 1.1 2008/10/17 06:12:08 scuri Exp $ + * $Id: hb_64.c,v 1.1 2008/10/17 06:12:08 scuri Exp $ + */ + +#include "hb.h" + +static const R *hb_64(R *rio, R *iio, const R *W, stride ios, int m, int dist) +{ + DK(KP634393284, +0.634393284163645498215171613225493370675687095); + DK(KP773010453, +0.773010453362736960810906609758469800971041293); + DK(KP098017140, +0.098017140329560601994195563888641845861136673); + DK(KP995184726, +0.995184726672196886244836953109479921575474869); + DK(KP471396736, +0.471396736825997648556387625905254377657460319); + DK(KP881921264, +0.881921264348355029712756863660388349508442621); + DK(KP290284677, +0.290284677254462367636192375817395274691476278); + DK(KP956940335, +0.956940335732208864935797886980269969482849206); + DK(KP195090322, +0.195090322016128267848284868477022240927691618); + DK(KP980785280, +0.980785280403230449126182236134239036973933731); + DK(KP555570233, +0.555570233019602224742830813948532874374937191); + DK(KP831469612, +0.831469612302545237078788377617905756738560812); + DK(KP382683432, +0.382683432365089771728459984030398866761344562); + DK(KP923879532, +0.923879532511286756128183189396788286822416626); + DK(KP707106781, +0.707106781186547524400844362104849039284835938); + int i; + for (i = m - 2; i > 0; i = i - 2, rio = rio + dist, iio = iio - dist, W = W + 126) { + E Tf, T7i, Tfa, ThM, Tgp, ThH, T2c, T5O, T4T, T6n, Tcp, Ted, TcA, TdE, T87; + E T9o, TK, T93, T2P, T4F, Tfo, Thz, T5T, T6j, Tbx, TdI, Tfl, ThA, T7r, T81; + E TbE, TdH, TZ, T94, T38, T4G, Tfv, ThC, T5W, T6k, TbQ, TdK, Tfs, ThD, T7w; + E T82, TbX, TdL, Tu, T84, Tfh, ThG, Tgm, ThN, T2v, T6m, T4K, T5P, Tce, TdF; + E TcD, Tec, T7l, T9p, T1L, T20, T9c, T9d, T9e, T9f, T40, T66, Tg1, Thu, Tg8; + E Thv, Tg5, Thr, T4n, T67, T4j, T69, T4w, T6a, TaT, TdW, Tb8, TdZ, TfU, Ths; + E T7O, T8y, T7T, T8z, Tbc, TdX, Tbj, Te0, T1g, T1v, T97, T98, T99, T9a, T3j; + E T5Z, TfI, Thk, TfP, Thl, TfM, Tho, T3G, T60, T3C, T62, T3P, T63, Tak, TdQ; + E Tav, TdT, TfB, Thn, T7D, T8v, T7I, T8w, TaD, TdP, TaG, TdS; + { + E T3, Tcm, T4O, Tcv, T6, Tcu, T4R, Tcn, Td, Tcy, T2a, Tch, Ta, Tcx, T27; + E Tck; + { + E T1, T2, T4P, T4Q; + T1 = rio[0]; + T2 = iio[-WS(ios, 32)]; + T3 = T1 + T2; + Tcm = T1 - T2; + { + E T4M, T4N, T4, T5; + T4M = iio[0]; + T4N = rio[WS(ios, 32)]; + T4O = T4M - T4N; + Tcv = T4M + T4N; + T4 = rio[WS(ios, 16)]; + T5 = iio[-WS(ios, 48)]; + T6 = T4 + T5; + Tcu = T4 - T5; + } + T4P = iio[-WS(ios, 16)]; + T4Q = rio[WS(ios, 48)]; + T4R = T4P - T4Q; + Tcn = T4P + T4Q; + { + E Tb, Tc, Tcf, T28, T29, Tcg; + Tb = iio[-WS(ios, 56)]; + Tc = rio[WS(ios, 24)]; + Tcf = Tb - Tc; + T28 = iio[-WS(ios, 24)]; + T29 = rio[WS(ios, 56)]; + Tcg = T29 + T28; + Td = Tb + Tc; + Tcy = Tcf + Tcg; + T2a = T28 - T29; + Tch = Tcf - Tcg; + } + { + E T8, T9, Tcj, T25, T26, Tci; + T8 = rio[WS(ios, 8)]; + T9 = iio[-WS(ios, 40)]; + Tcj = T8 - T9; + T25 = iio[-WS(ios, 8)]; + T26 = rio[WS(ios, 40)]; + Tci = T25 + T26; + Ta = T8 + T9; + Tcx = Tcj + Tci; + T27 = T25 - T26; + Tck = Tci - Tcj; + } + } + { + E T7, Te, Tf8, Tf9; + T7 = T3 + T6; + Te = Ta + Td; + Tf = T7 + Te; + T7i = T7 - Te; + Tf8 = Tcv - Tcu; + Tf9 = KP707106781 * (Tck + Tch); + Tfa = Tf8 + Tf9; + ThM = Tf8 - Tf9; + } + { + E Tgn, Tgo, T24, T2b; + Tgn = KP707106781 * (Tcx + Tcy); + Tgo = Tcm + Tcn; + Tgp = Tgn + Tgo; + ThH = Tgo - Tgn; + T24 = T3 - T6; + T2b = T27 - T2a; + T2c = T24 + T2b; + T5O = T24 - T2b; + } + { + E T4L, T4S, Tcl, Tco; + T4L = Td - Ta; + T4S = T4O - T4R; + T4T = T4L + T4S; + T6n = T4S - T4L; + Tcl = KP707106781 * (Tch - Tck); + Tco = Tcm - Tcn; + Tcp = Tcl + Tco; + Ted = Tco - Tcl; + } + { + E Tcw, Tcz, T85, T86; + Tcw = Tcu + Tcv; + Tcz = KP707106781 * (Tcx - Tcy); + TcA = Tcw + Tcz; + TdE = Tcw - Tcz; + T85 = T4O + T4R; + T86 = T27 + T2a; + T87 = T85 - T86; + T9o = T86 + T85; + } + } + { + E TC, Tby, T2x, Tbu, T2N, Tbz, T7o, Tbv, TJ, TbB, TbC, T2E, T2G, Tbp, Tbs; + E T7p, Tfj, Tfk; + { + E Tw, Tx, Ty, Tz, TA, TB; + Tw = rio[WS(ios, 2)]; + Tx = iio[-WS(ios, 34)]; + Ty = Tw + Tx; + Tz = rio[WS(ios, 18)]; + TA = iio[-WS(ios, 50)]; + TB = Tz + TA; + TC = Ty + TB; + Tby = Tz - TA; + T2x = Ty - TB; + Tbu = Tw - Tx; + } + { + E T2H, T2I, T2J, T2K, T2L, T2M; + T2H = iio[-WS(ios, 2)]; + T2I = rio[WS(ios, 34)]; + T2J = T2H - T2I; + T2K = iio[-WS(ios, 18)]; + T2L = rio[WS(ios, 50)]; + T2M = T2K - T2L; + T2N = T2J - T2M; + Tbz = T2H + T2I; + T7o = T2J + T2M; + Tbv = T2K + T2L; + } + { + E TF, Tbr, T2A, Tbq, TI, Tbn, T2D, Tbo; + { + E TD, TE, T2y, T2z; + TD = rio[WS(ios, 10)]; + TE = iio[-WS(ios, 42)]; + TF = TD + TE; + Tbr = TD - TE; + T2y = iio[-WS(ios, 10)]; + T2z = rio[WS(ios, 42)]; + T2A = T2y - T2z; + Tbq = T2y + T2z; + } + { + E TG, TH, T2B, T2C; + TG = iio[-WS(ios, 58)]; + TH = rio[WS(ios, 26)]; + TI = TG + TH; + Tbn = TG - TH; + T2B = iio[-WS(ios, 26)]; + T2C = rio[WS(ios, 58)]; + T2D = T2B - T2C; + Tbo = T2C + T2B; + } + TJ = TF + TI; + TbB = Tbr + Tbq; + TbC = Tbn + Tbo; + T2E = T2A - T2D; + T2G = TI - TF; + Tbp = Tbn - Tbo; + Tbs = Tbq - Tbr; + T7p = T2A + T2D; + } + TK = TC + TJ; + T93 = T7p + T7o; + { + E T2F, T2O, Tfm, Tfn; + T2F = T2x + T2E; + T2O = T2G + T2N; + T2P = FMA(KP923879532, T2F, KP382683432 * T2O); + T4F = FNMS(KP382683432, T2F, KP923879532 * T2O); + Tfm = KP707106781 * (TbB + TbC); + Tfn = Tbu + Tbv; + Tfo = Tfm + Tfn; + Thz = Tfn - Tfm; + } + { + E T5R, T5S, Tbt, Tbw; + T5R = T2x - T2E; + T5S = T2N - T2G; + T5T = FNMS(KP382683432, T5S, KP923879532 * T5R); + T6j = FMA(KP382683432, T5R, KP923879532 * T5S); + Tbt = KP707106781 * (Tbp - Tbs); + Tbw = Tbu - Tbv; + Tbx = Tbt + Tbw; + TdI = Tbw - Tbt; + } + Tfj = Tbz - Tby; + Tfk = KP707106781 * (Tbs + Tbp); + Tfl = Tfj + Tfk; + ThA = Tfj - Tfk; + { + E T7n, T7q, TbA, TbD; + T7n = TC - TJ; + T7q = T7o - T7p; + T7r = T7n + T7q; + T81 = T7q - T7n; + TbA = Tby + Tbz; + TbD = KP707106781 * (TbB - TbC); + TbE = TbA + TbD; + TdH = TbA - TbD; + } + } + { + E TR, TbU, T2Q, TbN, T36, TbV, T7t, TbO, TY, TbR, TbS, T2X, T2Z, TbI, TbL; + E T7u, Tfq, Tfr; + { + E TL, TM, TN, TO, TP, TQ; + TL = iio[-WS(ios, 62)]; + TM = rio[WS(ios, 30)]; + TN = TL + TM; + TO = rio[WS(ios, 14)]; + TP = iio[-WS(ios, 46)]; + TQ = TO + TP; + TR = TN + TQ; + TbU = TL - TM; + T2Q = TN - TQ; + TbN = TO - TP; + } + { + E T30, T31, T32, T33, T34, T35; + T30 = iio[-WS(ios, 30)]; + T31 = rio[WS(ios, 62)]; + T32 = T30 - T31; + T33 = iio[-WS(ios, 14)]; + T34 = rio[WS(ios, 46)]; + T35 = T33 - T34; + T36 = T32 - T35; + TbV = T33 + T34; + T7t = T32 + T35; + TbO = T31 + T30; + } + { + E TU, TbG, T2T, TbH, TX, TbJ, T2W, TbK; + { + E TS, TT, T2R, T2S; + TS = rio[WS(ios, 6)]; + TT = iio[-WS(ios, 38)]; + TU = TS + TT; + TbG = TS - TT; + T2R = iio[-WS(ios, 6)]; + T2S = rio[WS(ios, 38)]; + T2T = T2R - T2S; + TbH = T2R + T2S; + } + { + E TV, TW, T2U, T2V; + TV = iio[-WS(ios, 54)]; + TW = rio[WS(ios, 22)]; + TX = TV + TW; + TbJ = TV - TW; + T2U = iio[-WS(ios, 22)]; + T2V = rio[WS(ios, 54)]; + T2W = T2U - T2V; + TbK = T2V + T2U; + } + TY = TU + TX; + TbR = TbJ - TbK; + TbS = TbH - TbG; + T2X = T2T - T2W; + T2Z = TX - TU; + TbI = TbG + TbH; + TbL = TbJ + TbK; + T7u = T2T + T2W; + } + TZ = TR + TY; + T94 = T7u + T7t; + { + E T2Y, T37, Tft, Tfu; + T2Y = T2Q + T2X; + T37 = T2Z + T36; + T38 = FNMS(KP382683432, T37, KP923879532 * T2Y); + T4G = FMA(KP382683432, T2Y, KP923879532 * T37); + Tft = KP707106781 * (TbI + TbL); + Tfu = TbU + TbV; + Tfv = Tft + Tfu; + ThC = Tfu - Tft; + } + { + E T5U, T5V, TbM, TbP; + T5U = T2Q - T2X; + T5V = T36 - T2Z; + T5W = FMA(KP923879532, T5U, KP382683432 * T5V); + T6k = FNMS(KP382683432, T5U, KP923879532 * T5V); + TbM = KP707106781 * (TbI - TbL); + TbP = TbN - TbO; + TbQ = TbM + TbP; + TdK = TbP - TbM; + } + Tfq = KP707106781 * (TbS + TbR); + Tfr = TbN + TbO; + Tfs = Tfq - Tfr; + ThD = Tfq + Tfr; + { + E T7s, T7v, TbT, TbW; + T7s = TR - TY; + T7v = T7t - T7u; + T7w = T7s - T7v; + T82 = T7s + T7v; + TbT = KP707106781 * (TbR - TbS); + TbW = TbU - TbV; + TbX = TbT + TbW; + TdL = TbW - TbT; + } + } + { + E Ti, T2g, Tl, T2j, T2d, T2k, Tfc, Tfb, Tc5, Tc2, Tp, T2p, Ts, T2s, T2m; + E T2t, Tff, Tfe, Tcc, Tc9; + { + E Tc0, Tc4, Tc3, Tc1; + { + E Tg, Th, T2e, T2f; + Tg = rio[WS(ios, 4)]; + Th = iio[-WS(ios, 36)]; + Ti = Tg + Th; + Tc0 = Tg - Th; + T2e = iio[-WS(ios, 4)]; + T2f = rio[WS(ios, 36)]; + T2g = T2e - T2f; + Tc4 = T2e + T2f; + } + { + E Tj, Tk, T2h, T2i; + Tj = rio[WS(ios, 20)]; + Tk = iio[-WS(ios, 52)]; + Tl = Tj + Tk; + Tc3 = Tj - Tk; + T2h = iio[-WS(ios, 20)]; + T2i = rio[WS(ios, 52)]; + T2j = T2h - T2i; + Tc1 = T2h + T2i; + } + T2d = Ti - Tl; + T2k = T2g - T2j; + Tfc = Tc0 + Tc1; + Tfb = Tc4 - Tc3; + Tc5 = Tc3 + Tc4; + Tc2 = Tc0 - Tc1; + } + { + E Tc7, Tcb, Tca, Tc8; + { + E Tn, To, T2n, T2o; + Tn = iio[-WS(ios, 60)]; + To = rio[WS(ios, 28)]; + Tp = Tn + To; + Tc7 = Tn - To; + T2n = iio[-WS(ios, 28)]; + T2o = rio[WS(ios, 60)]; + T2p = T2n - T2o; + Tcb = T2o + T2n; + } + { + E Tq, Tr, T2q, T2r; + Tq = rio[WS(ios, 12)]; + Tr = iio[-WS(ios, 44)]; + Ts = Tq + Tr; + Tca = Tq - Tr; + T2q = iio[-WS(ios, 12)]; + T2r = rio[WS(ios, 44)]; + T2s = T2q - T2r; + Tc8 = T2q + T2r; + } + T2m = Tp - Ts; + T2t = T2p - T2s; + Tff = Tca + Tcb; + Tfe = Tc7 + Tc8; + Tcc = Tca - Tcb; + Tc9 = Tc7 - Tc8; + } + { + E Tm, Tt, Tfd, Tfg; + Tm = Ti + Tl; + Tt = Tp + Ts; + Tu = Tm + Tt; + T84 = Tt - Tm; + Tfd = FNMS(KP382683432, Tfc, KP923879532 * Tfb); + Tfg = FNMS(KP923879532, Tff, KP382683432 * Tfe); + Tfh = Tfd + Tfg; + ThG = Tfg - Tfd; + } + { + E Tgk, Tgl, T2l, T2u; + Tgk = FMA(KP382683432, Tfb, KP923879532 * Tfc); + Tgl = FMA(KP923879532, Tfe, KP382683432 * Tff); + Tgm = Tgk + Tgl; + ThN = Tgk - Tgl; + T2l = T2d + T2k; + T2u = T2m - T2t; + T2v = KP707106781 * (T2l + T2u); + T6m = KP707106781 * (T2l - T2u); + } + { + E T4I, T4J, Tc6, Tcd; + T4I = T2k - T2d; + T4J = T2m + T2t; + T4K = KP707106781 * (T4I + T4J); + T5P = KP707106781 * (T4J - T4I); + Tc6 = FNMS(KP382683432, Tc5, KP923879532 * Tc2); + Tcd = FMA(KP923879532, Tc9, KP382683432 * Tcc); + Tce = Tc6 + Tcd; + TdF = Tcd - Tc6; + } + { + E TcB, TcC, T7j, T7k; + TcB = FMA(KP923879532, Tc5, KP382683432 * Tc2); + TcC = FNMS(KP382683432, Tc9, KP923879532 * Tcc); + TcD = TcB + TcC; + Tec = TcB - TcC; + T7j = T2g + T2j; + T7k = T2s + T2p; + T7l = T7j - T7k; + T9p = T7j + T7k; + } + } + { + E T1z, T1C, T1D, Tbg, TaQ, T4r, T4u, T7Q, Tbh, TaR, T1G, T3V, T1J, T3Y, T1K; + E T7R, Tbe, Tbd, TaO, TaL, T1S, TfV, TfW, T41, T48, TaW, TaZ, T7L, T1Z, TfY; + E TfZ, T4a, T4h, Tb3, Tb6, T7M; + { + E T1x, T1y, T1A, T1B; + T1x = iio[-WS(ios, 63)]; + T1y = rio[WS(ios, 31)]; + T1z = T1x + T1y; + T1A = rio[WS(ios, 15)]; + T1B = iio[-WS(ios, 47)]; + T1C = T1A + T1B; + T1D = T1z + T1C; + Tbg = T1x - T1y; + TaQ = T1A - T1B; + } + { + E T4p, T4q, T4s, T4t; + T4p = iio[-WS(ios, 31)]; + T4q = rio[WS(ios, 63)]; + T4r = T4p - T4q; + T4s = iio[-WS(ios, 15)]; + T4t = rio[WS(ios, 47)]; + T4u = T4s - T4t; + T7Q = T4r + T4u; + Tbh = T4s + T4t; + TaR = T4q + T4p; + } + { + E TaJ, TaK, TaM, TaN; + { + E T1E, T1F, T3T, T3U; + T1E = rio[WS(ios, 7)]; + T1F = iio[-WS(ios, 39)]; + T1G = T1E + T1F; + TaJ = T1E - T1F; + T3T = iio[-WS(ios, 7)]; + T3U = rio[WS(ios, 39)]; + T3V = T3T - T3U; + TaK = T3T + T3U; + } + { + E T1H, T1I, T3W, T3X; + T1H = iio[-WS(ios, 55)]; + T1I = rio[WS(ios, 23)]; + T1J = T1H + T1I; + TaM = T1H - T1I; + T3W = iio[-WS(ios, 23)]; + T3X = rio[WS(ios, 55)]; + T3Y = T3W - T3X; + TaN = T3X + T3W; + } + T1K = T1G + T1J; + T7R = T3V + T3Y; + Tbe = TaK - TaJ; + Tbd = TaM - TaN; + TaO = TaM + TaN; + TaL = TaJ + TaK; + } + { + E T1O, TaX, T44, TaV, T1R, TaU, T47, TaY; + { + E T1M, T1N, T42, T43; + T1M = rio[WS(ios, 3)]; + T1N = iio[-WS(ios, 35)]; + T1O = T1M + T1N; + TaX = T1M - T1N; + T42 = iio[-WS(ios, 3)]; + T43 = rio[WS(ios, 35)]; + T44 = T42 - T43; + TaV = T42 + T43; + } + { + E T1P, T1Q, T45, T46; + T1P = rio[WS(ios, 19)]; + T1Q = iio[-WS(ios, 51)]; + T1R = T1P + T1Q; + TaU = T1P - T1Q; + T45 = iio[-WS(ios, 19)]; + T46 = rio[WS(ios, 51)]; + T47 = T45 - T46; + TaY = T45 + T46; + } + T1S = T1O + T1R; + TfV = TaV - TaU; + TfW = TaX + TaY; + T41 = T1O - T1R; + T48 = T44 - T47; + TaW = TaU + TaV; + TaZ = TaX - TaY; + T7L = T44 + T47; + } + { + E T1V, Tb4, T4d, Tb2, T1Y, Tb1, T4g, Tb5; + { + E T1T, T1U, T4b, T4c; + T1T = iio[-WS(ios, 59)]; + T1U = rio[WS(ios, 27)]; + T1V = T1T + T1U; + Tb4 = T1T - T1U; + T4b = iio[-WS(ios, 27)]; + T4c = rio[WS(ios, 59)]; + T4d = T4b - T4c; + Tb2 = T4c + T4b; + } + { + E T1W, T1X, T4e, T4f; + T1W = rio[WS(ios, 11)]; + T1X = iio[-WS(ios, 43)]; + T1Y = T1W + T1X; + Tb1 = T1W - T1X; + T4e = iio[-WS(ios, 11)]; + T4f = rio[WS(ios, 43)]; + T4g = T4e - T4f; + Tb5 = T4e + T4f; + } + T1Z = T1V + T1Y; + TfY = Tb4 + Tb5; + TfZ = Tb1 + Tb2; + T4a = T1V - T1Y; + T4h = T4d - T4g; + Tb3 = Tb1 - Tb2; + Tb6 = Tb4 - Tb5; + T7M = T4g + T4d; + } + T1L = T1D + T1K; + T20 = T1S + T1Z; + T9c = T1L - T20; + T9d = T7R + T7Q; + T9e = T7L + T7M; + T9f = T9d - T9e; + { + E T3S, T3Z, TfX, Tg0; + T3S = T1z - T1C; + T3Z = T3V - T3Y; + T40 = T3S + T3Z; + T66 = T3S - T3Z; + TfX = FNMS(KP382683432, TfW, KP923879532 * TfV); + Tg0 = FNMS(KP923879532, TfZ, KP382683432 * TfY); + Tg1 = TfX + Tg0; + Thu = Tg0 - TfX; + } + { + E Tg6, Tg7, Tg3, Tg4; + Tg6 = KP707106781 * (TaL + TaO); + Tg7 = Tbg + Tbh; + Tg8 = Tg6 + Tg7; + Thv = Tg7 - Tg6; + Tg3 = FMA(KP382683432, TfV, KP923879532 * TfW); + Tg4 = FMA(KP923879532, TfY, KP382683432 * TfZ); + Tg5 = Tg3 + Tg4; + Thr = Tg3 - Tg4; + } + { + E T4l, T4m, T49, T4i; + T4l = T48 - T41; + T4m = T4a + T4h; + T4n = KP707106781 * (T4l + T4m); + T67 = KP707106781 * (T4m - T4l); + T49 = T41 + T48; + T4i = T4a - T4h; + T4j = KP707106781 * (T49 + T4i); + T69 = KP707106781 * (T49 - T4i); + } + { + E T4o, T4v, TaP, TaS; + T4o = T1J - T1G; + T4v = T4r - T4u; + T4w = T4o + T4v; + T6a = T4v - T4o; + TaP = KP707106781 * (TaL - TaO); + TaS = TaQ - TaR; + TaT = TaP + TaS; + TdW = TaS - TaP; + } + { + E Tb0, Tb7, TfS, TfT; + Tb0 = FMA(KP923879532, TaW, KP382683432 * TaZ); + Tb7 = FNMS(KP382683432, Tb6, KP923879532 * Tb3); + Tb8 = Tb0 + Tb7; + TdZ = Tb0 - Tb7; + TfS = KP707106781 * (Tbe + Tbd); + TfT = TaQ + TaR; + TfU = TfS - TfT; + Ths = TfS + TfT; + } + { + E T7K, T7N, T7P, T7S; + T7K = T1D - T1K; + T7N = T7L - T7M; + T7O = T7K + T7N; + T8y = T7K - T7N; + T7P = T1Z - T1S; + T7S = T7Q - T7R; + T7T = T7P + T7S; + T8z = T7S - T7P; + } + { + E Tba, Tbb, Tbf, Tbi; + Tba = FNMS(KP382683432, TaW, KP923879532 * TaZ); + Tbb = FMA(KP923879532, Tb6, KP382683432 * Tb3); + Tbc = Tba + Tbb; + TdX = Tbb - Tba; + Tbf = KP707106781 * (Tbd - Tbe); + Tbi = Tbg - Tbh; + Tbj = Tbf + Tbi; + Te0 = Tbi - Tbf; + } + } + { + E T14, T17, T18, Tax, Tas, T3K, T3N, T7F, Tay, Tat, T1b, T3e, T1e, T3h, T1f; + E T7G, TaB, TaA, Taq, Tan, T1n, TfC, TfD, T3k, T3r, Ta8, Tab, T7A, T1u, TfF; + E TfG, T3t, T3A, Taf, Tai, T7B; + { + E T12, T13, T15, T16; + T12 = rio[WS(ios, 1)]; + T13 = iio[-WS(ios, 33)]; + T14 = T12 + T13; + T15 = rio[WS(ios, 17)]; + T16 = iio[-WS(ios, 49)]; + T17 = T15 + T16; + T18 = T14 + T17; + Tax = T15 - T16; + Tas = T12 - T13; + } + { + E T3I, T3J, T3L, T3M; + T3I = iio[-WS(ios, 1)]; + T3J = rio[WS(ios, 33)]; + T3K = T3I - T3J; + T3L = iio[-WS(ios, 17)]; + T3M = rio[WS(ios, 49)]; + T3N = T3L - T3M; + T7F = T3K + T3N; + Tay = T3I + T3J; + Tat = T3L + T3M; + } + { + E Tap, Tao, Tal, Tam; + { + E T19, T1a, T3c, T3d; + T19 = rio[WS(ios, 9)]; + T1a = iio[-WS(ios, 41)]; + T1b = T19 + T1a; + Tap = T19 - T1a; + T3c = iio[-WS(ios, 9)]; + T3d = rio[WS(ios, 41)]; + T3e = T3c - T3d; + Tao = T3c + T3d; + } + { + E T1c, T1d, T3f, T3g; + T1c = iio[-WS(ios, 57)]; + T1d = rio[WS(ios, 25)]; + T1e = T1c + T1d; + Tal = T1c - T1d; + T3f = iio[-WS(ios, 25)]; + T3g = rio[WS(ios, 57)]; + T3h = T3f - T3g; + Tam = T3g + T3f; + } + T1f = T1b + T1e; + T7G = T3e + T3h; + TaB = Tal + Tam; + TaA = Tap + Tao; + Taq = Tao - Tap; + Tan = Tal - Tam; + } + { + E T1j, Ta6, T3n, Taa, T1m, Ta9, T3q, Ta7; + { + E T1h, T1i, T3l, T3m; + T1h = rio[WS(ios, 5)]; + T1i = iio[-WS(ios, 37)]; + T1j = T1h + T1i; + Ta6 = T1h - T1i; + T3l = iio[-WS(ios, 5)]; + T3m = rio[WS(ios, 37)]; + T3n = T3l - T3m; + Taa = T3l + T3m; + } + { + E T1k, T1l, T3o, T3p; + T1k = rio[WS(ios, 21)]; + T1l = iio[-WS(ios, 53)]; + T1m = T1k + T1l; + Ta9 = T1k - T1l; + T3o = iio[-WS(ios, 21)]; + T3p = rio[WS(ios, 53)]; + T3q = T3o - T3p; + Ta7 = T3o + T3p; + } + T1n = T1j + T1m; + TfC = Taa - Ta9; + TfD = Ta6 + Ta7; + T3k = T1j - T1m; + T3r = T3n - T3q; + Ta8 = Ta6 - Ta7; + Tab = Ta9 + Taa; + T7A = T3n + T3q; + } + { + E T1q, Tad, T3w, Tah, T1t, Tag, T3z, Tae; + { + E T1o, T1p, T3u, T3v; + T1o = iio[-WS(ios, 61)]; + T1p = rio[WS(ios, 29)]; + T1q = T1o + T1p; + Tad = T1o - T1p; + T3u = iio[-WS(ios, 29)]; + T3v = rio[WS(ios, 61)]; + T3w = T3u - T3v; + Tah = T3v + T3u; + } + { + E T1r, T1s, T3x, T3y; + T1r = rio[WS(ios, 13)]; + T1s = iio[-WS(ios, 45)]; + T1t = T1r + T1s; + Tag = T1r - T1s; + T3x = iio[-WS(ios, 13)]; + T3y = rio[WS(ios, 45)]; + T3z = T3x - T3y; + Tae = T3x + T3y; + } + T1u = T1q + T1t; + TfF = Tad + Tae; + TfG = Tag + Tah; + T3t = T1q - T1t; + T3A = T3w - T3z; + Taf = Tad - Tae; + Tai = Tag - Tah; + T7B = T3z + T3w; + } + T1g = T18 + T1f; + T1v = T1n + T1u; + T97 = T1g - T1v; + T98 = T7G + T7F; + T99 = T7A + T7B; + T9a = T98 - T99; + { + E T3b, T3i, TfE, TfH; + T3b = T14 - T17; + T3i = T3e - T3h; + T3j = T3b + T3i; + T5Z = T3b - T3i; + TfE = FNMS(KP382683432, TfD, KP923879532 * TfC); + TfH = FNMS(KP923879532, TfG, KP382683432 * TfF); + TfI = TfE + TfH; + Thk = TfH - TfE; + } + { + E TfN, TfO, TfK, TfL; + TfN = KP707106781 * (TaA + TaB); + TfO = Tas + Tat; + TfP = TfN + TfO; + Thl = TfO - TfN; + TfK = FMA(KP382683432, TfC, KP923879532 * TfD); + TfL = FMA(KP923879532, TfF, KP382683432 * TfG); + TfM = TfK + TfL; + Tho = TfK - TfL; + } + { + E T3E, T3F, T3s, T3B; + T3E = T3r - T3k; + T3F = T3t + T3A; + T3G = KP707106781 * (T3E + T3F); + T60 = KP707106781 * (T3F - T3E); + T3s = T3k + T3r; + T3B = T3t - T3A; + T3C = KP707106781 * (T3s + T3B); + T62 = KP707106781 * (T3s - T3B); + } + { + E T3H, T3O, Tac, Taj; + T3H = T1e - T1b; + T3O = T3K - T3N; + T3P = T3H + T3O; + T63 = T3O - T3H; + Tac = FNMS(KP382683432, Tab, KP923879532 * Ta8); + Taj = FMA(KP923879532, Taf, KP382683432 * Tai); + Tak = Tac + Taj; + TdQ = Taj - Tac; + } + { + E Tar, Tau, Tfz, TfA; + Tar = KP707106781 * (Tan - Taq); + Tau = Tas - Tat; + Tav = Tar + Tau; + TdT = Tau - Tar; + Tfz = Tay - Tax; + TfA = KP707106781 * (Taq + Tan); + TfB = Tfz + TfA; + Thn = Tfz - TfA; + } + { + E T7z, T7C, T7E, T7H; + T7z = T18 - T1f; + T7C = T7A - T7B; + T7D = T7z + T7C; + T8v = T7z - T7C; + T7E = T1u - T1n; + T7H = T7F - T7G; + T7I = T7E + T7H; + T8w = T7H - T7E; + } + { + E Taz, TaC, TaE, TaF; + Taz = Tax + Tay; + TaC = KP707106781 * (TaA - TaB); + TaD = Taz + TaC; + TdP = Taz - TaC; + TaE = FMA(KP923879532, Tab, KP382683432 * Ta8); + TaF = FNMS(KP382683432, Taf, KP923879532 * Tai); + TaG = TaE + TaF; + TdS = TaE - TaF; + } + } + { + E T11, T9K, T9T, Ta2, T22, T9Q, T9N, Ta3; + { + E Tv, T10, T9R, T9S; + Tv = Tf + Tu; + T10 = TK + TZ; + T11 = Tv + T10; + T9K = Tv - T10; + T9R = T9p + T9o; + T9S = T93 + T94; + T9T = T9R - T9S; + Ta2 = T9S + T9R; + } + { + E T1w, T21, T9L, T9M; + T1w = T1g + T1v; + T21 = T1L + T20; + T22 = T1w + T21; + T9Q = T21 - T1w; + T9L = T99 + T98; + T9M = T9e + T9d; + T9N = T9L - T9M; + Ta3 = T9L + T9M; + } + rio[0] = T11 + T22; + iio[-WS(ios, 63)] = Ta3 + Ta2; + { + E T9O, T9U, T9J, T9P; + T9O = T9K + T9N; + T9U = T9Q + T9T; + T9J = W[94]; + T9P = W[95]; + rio[WS(ios, 48)] = FNMS(T9P, T9U, T9J * T9O); + iio[-WS(ios, 15)] = FMA(T9P, T9O, T9J * T9U); + } + { + E T9W, T9Y, T9V, T9X; + T9W = T9K - T9N; + T9Y = T9T - T9Q; + T9V = W[30]; + T9X = W[31]; + rio[WS(ios, 16)] = FNMS(T9X, T9Y, T9V * T9W); + iio[-WS(ios, 47)] = FMA(T9X, T9W, T9V * T9Y); + } + { + E Ta0, Ta4, T9Z, Ta1; + Ta0 = T11 - T22; + Ta4 = Ta2 - Ta3; + T9Z = W[62]; + Ta1 = W[63]; + rio[WS(ios, 32)] = FNMS(Ta1, Ta4, T9Z * Ta0); + iio[-WS(ios, 31)] = FMA(Ta1, Ta0, T9Z * Ta4); + } + } + { + E T96, T9y, T9r, T9D, T9h, T9C, T9m, T9z; + { + E T92, T95, T9n, T9q; + T92 = Tf - Tu; + T95 = T93 - T94; + T96 = T92 + T95; + T9y = T92 - T95; + T9n = TZ - TK; + T9q = T9o - T9p; + T9r = T9n + T9q; + T9D = T9q - T9n; + } + { + E T9b, T9g, T9k, T9l; + T9b = T97 + T9a; + T9g = T9c - T9f; + T9h = KP707106781 * (T9b + T9g); + T9C = KP707106781 * (T9b - T9g); + T9k = T9a - T97; + T9l = T9c + T9f; + T9m = KP707106781 * (T9k + T9l); + T9z = KP707106781 * (T9l - T9k); + } + { + E T9i, T9s, T91, T9j; + T9i = T96 + T9h; + T9s = T9m + T9r; + T91 = W[110]; + T9j = W[111]; + rio[WS(ios, 56)] = FNMS(T9j, T9s, T91 * T9i); + iio[-WS(ios, 7)] = FMA(T9j, T9i, T91 * T9s); + } + { + E T9G, T9I, T9F, T9H; + T9G = T9y - T9z; + T9I = T9D - T9C; + T9F = W[78]; + T9H = W[79]; + rio[WS(ios, 40)] = FNMS(T9H, T9I, T9F * T9G); + iio[-WS(ios, 23)] = FMA(T9H, T9G, T9F * T9I); + } + { + E T9u, T9w, T9t, T9v; + T9u = T96 - T9h; + T9w = T9r - T9m; + T9t = W[46]; + T9v = W[47]; + rio[WS(ios, 24)] = FNMS(T9v, T9w, T9t * T9u); + iio[-WS(ios, 39)] = FMA(T9v, T9u, T9t * T9w); + } + { + E T9A, T9E, T9x, T9B; + T9A = T9y + T9z; + T9E = T9C + T9D; + T9x = W[14]; + T9B = W[15]; + rio[WS(ios, 8)] = FNMS(T9B, T9E, T9x * T9A); + iio[-WS(ios, 55)] = FMA(T9B, T9A, T9x * T9E); + } + } + { + E T8u, T8Q, T8J, T8V, T8B, T8U, T8G, T8R; + { + E T8s, T8t, T8H, T8I; + T8s = T7i - T7l; + T8t = KP707106781 * (T82 - T81); + T8u = T8s + T8t; + T8Q = T8s - T8t; + T8H = KP707106781 * (T7r - T7w); + T8I = T87 - T84; + T8J = T8H + T8I; + T8V = T8I - T8H; + } + { + E T8x, T8A, T8E, T8F; + T8x = FNMS(KP382683432, T8w, KP923879532 * T8v); + T8A = FMA(KP923879532, T8y, KP382683432 * T8z); + T8B = T8x + T8A; + T8U = T8A - T8x; + T8E = FMA(KP382683432, T8v, KP923879532 * T8w); + T8F = FNMS(KP382683432, T8y, KP923879532 * T8z); + T8G = T8E + T8F; + T8R = T8E - T8F; + } + { + E T8C, T8K, T8r, T8D; + T8C = T8u + T8B; + T8K = T8G + T8J; + T8r = W[6]; + T8D = W[7]; + rio[WS(ios, 4)] = FNMS(T8D, T8K, T8r * T8C); + iio[-WS(ios, 59)] = FMA(T8D, T8C, T8r * T8K); + } + { + E T8Y, T90, T8X, T8Z; + T8Y = T8Q - T8R; + T90 = T8V - T8U; + T8X = W[38]; + T8Z = W[39]; + rio[WS(ios, 20)] = FNMS(T8Z, T90, T8X * T8Y); + iio[-WS(ios, 43)] = FMA(T8Z, T8Y, T8X * T90); + } + { + E T8M, T8O, T8L, T8N; + T8M = T8u - T8B; + T8O = T8J - T8G; + T8L = W[70]; + T8N = W[71]; + rio[WS(ios, 36)] = FNMS(T8N, T8O, T8L * T8M); + iio[-WS(ios, 27)] = FMA(T8N, T8M, T8L * T8O); + } + { + E T8S, T8W, T8P, T8T; + T8S = T8Q + T8R; + T8W = T8U + T8V; + T8P = W[102]; + T8T = W[103]; + rio[WS(ios, 52)] = FNMS(T8T, T8W, T8P * T8S); + iio[-WS(ios, 11)] = FMA(T8T, T8S, T8P * T8W); + } + } + { + E T7y, T8g, T89, T8l, T7V, T8k, T80, T8h; + { + E T7m, T7x, T83, T88; + T7m = T7i + T7l; + T7x = KP707106781 * (T7r + T7w); + T7y = T7m + T7x; + T8g = T7m - T7x; + T83 = KP707106781 * (T81 + T82); + T88 = T84 + T87; + T89 = T83 + T88; + T8l = T88 - T83; + } + { + E T7J, T7U, T7Y, T7Z; + T7J = FMA(KP923879532, T7D, KP382683432 * T7I); + T7U = FNMS(KP382683432, T7T, KP923879532 * T7O); + T7V = T7J + T7U; + T8k = T7J - T7U; + T7Y = FNMS(KP382683432, T7D, KP923879532 * T7I); + T7Z = FMA(KP382683432, T7O, KP923879532 * T7T); + T80 = T7Y + T7Z; + T8h = T7Z - T7Y; + } + { + E T7W, T8a, T7h, T7X; + T7W = T7y + T7V; + T8a = T80 + T89; + T7h = W[118]; + T7X = W[119]; + rio[WS(ios, 60)] = FNMS(T7X, T8a, T7h * T7W); + iio[-WS(ios, 3)] = FMA(T7X, T7W, T7h * T8a); + } + { + E T8o, T8q, T8n, T8p; + T8o = T8g - T8h; + T8q = T8l - T8k; + T8n = W[86]; + T8p = W[87]; + rio[WS(ios, 44)] = FNMS(T8p, T8q, T8n * T8o); + iio[-WS(ios, 19)] = FMA(T8p, T8o, T8n * T8q); + } + { + E T8c, T8e, T8b, T8d; + T8c = T7y - T7V; + T8e = T89 - T80; + T8b = W[54]; + T8d = W[55]; + rio[WS(ios, 28)] = FNMS(T8d, T8e, T8b * T8c); + iio[-WS(ios, 35)] = FMA(T8d, T8c, T8b * T8e); + } + { + E T8i, T8m, T8f, T8j; + T8i = T8g + T8h; + T8m = T8k + T8l; + T8f = W[22]; + T8j = W[23]; + rio[WS(ios, 12)] = FNMS(T8j, T8m, T8f * T8i); + iio[-WS(ios, 51)] = FMA(T8j, T8i, T8f * T8m); + } + } + { + E T6K, T76, T6Z, T7b, T6R, T7a, T6W, T77; + { + E T6I, T6J, T6X, T6Y; + T6I = T5O - T5P; + T6J = T6j - T6k; + T6K = T6I + T6J; + T76 = T6I - T6J; + T6X = T5W - T5T; + T6Y = T6n - T6m; + T6Z = T6X + T6Y; + T7b = T6Y - T6X; + { + E T6N, T6U, T6Q, T6V; + { + E T6L, T6M, T6O, T6P; + T6L = T5Z - T60; + T6M = T63 - T62; + T6N = FMA(KP831469612, T6L, KP555570233 * T6M); + T6U = FNMS(KP555570233, T6L, KP831469612 * T6M); + T6O = T66 - T67; + T6P = T6a - T69; + T6Q = FNMS(KP555570233, T6P, KP831469612 * T6O); + T6V = FMA(KP555570233, T6O, KP831469612 * T6P); + } + T6R = T6N + T6Q; + T7a = T6N - T6Q; + T6W = T6U + T6V; + T77 = T6V - T6U; + } + } + { + E T6S, T70, T6H, T6T; + T6S = T6K + T6R; + T70 = T6W + T6Z; + T6H = W[114]; + T6T = W[115]; + rio[WS(ios, 58)] = FNMS(T6T, T70, T6H * T6S); + iio[-WS(ios, 5)] = FMA(T6T, T6S, T6H * T70); + } + { + E T7e, T7g, T7d, T7f; + T7e = T76 - T77; + T7g = T7b - T7a; + T7d = W[82]; + T7f = W[83]; + rio[WS(ios, 42)] = FNMS(T7f, T7g, T7d * T7e); + iio[-WS(ios, 21)] = FMA(T7f, T7e, T7d * T7g); + } + { + E T72, T74, T71, T73; + T72 = T6K - T6R; + T74 = T6Z - T6W; + T71 = W[50]; + T73 = W[51]; + rio[WS(ios, 26)] = FNMS(T73, T74, T71 * T72); + iio[-WS(ios, 37)] = FMA(T73, T72, T71 * T74); + } + { + E T78, T7c, T75, T79; + T78 = T76 + T77; + T7c = T7a + T7b; + T75 = W[18]; + T79 = W[19]; + rio[WS(ios, 10)] = FNMS(T79, T7c, T75 * T78); + iio[-WS(ios, 53)] = FMA(T79, T78, T75 * T7c); + } + } + { + E T3a, T52, T4V, T57, T4z, T56, T4E, T53; + { + E T2w, T39, T4H, T4U; + T2w = T2c + T2v; + T39 = T2P + T38; + T3a = T2w + T39; + T52 = T2w - T39; + T4H = T4F + T4G; + T4U = T4K + T4T; + T4V = T4H + T4U; + T57 = T4U - T4H; + { + E T3R, T4C, T4y, T4D; + { + E T3D, T3Q, T4k, T4x; + T3D = T3j + T3C; + T3Q = T3G + T3P; + T3R = FMA(KP980785280, T3D, KP195090322 * T3Q); + T4C = FNMS(KP195090322, T3D, KP980785280 * T3Q); + T4k = T40 + T4j; + T4x = T4n + T4w; + T4y = FNMS(KP195090322, T4x, KP980785280 * T4k); + T4D = FMA(KP195090322, T4k, KP980785280 * T4x); + } + T4z = T3R + T4y; + T56 = T3R - T4y; + T4E = T4C + T4D; + T53 = T4D - T4C; + } + } + { + E T4A, T4W, T23, T4B; + T4A = T3a + T4z; + T4W = T4E + T4V; + T23 = W[122]; + T4B = W[123]; + rio[WS(ios, 62)] = FNMS(T4B, T4W, T23 * T4A); + iio[-WS(ios, 1)] = FMA(T4B, T4A, T23 * T4W); + } + { + E T5a, T5c, T59, T5b; + T5a = T52 - T53; + T5c = T57 - T56; + T59 = W[90]; + T5b = W[91]; + rio[WS(ios, 46)] = FNMS(T5b, T5c, T59 * T5a); + iio[-WS(ios, 17)] = FMA(T5b, T5a, T59 * T5c); + } + { + E T4Y, T50, T4X, T4Z; + T4Y = T3a - T4z; + T50 = T4V - T4E; + T4X = W[58]; + T4Z = W[59]; + rio[WS(ios, 30)] = FNMS(T4Z, T50, T4X * T4Y); + iio[-WS(ios, 33)] = FMA(T4Z, T4Y, T4X * T50); + } + { + E T54, T58, T51, T55; + T54 = T52 + T53; + T58 = T56 + T57; + T51 = W[26]; + T55 = W[27]; + rio[WS(ios, 14)] = FNMS(T55, T58, T51 * T54); + iio[-WS(ios, 49)] = FMA(T55, T54, T51 * T58); + } + } + { + E T5g, T5C, T5v, T5H, T5n, T5G, T5s, T5D; + { + E T5e, T5f, T5t, T5u; + T5e = T2c - T2v; + T5f = T4G - T4F; + T5g = T5e + T5f; + T5C = T5e - T5f; + T5t = T2P - T38; + T5u = T4T - T4K; + T5v = T5t + T5u; + T5H = T5u - T5t; + { + E T5j, T5q, T5m, T5r; + { + E T5h, T5i, T5k, T5l; + T5h = T3j - T3C; + T5i = T3P - T3G; + T5j = FNMS(KP555570233, T5i, KP831469612 * T5h); + T5q = FMA(KP555570233, T5h, KP831469612 * T5i); + T5k = T40 - T4j; + T5l = T4w - T4n; + T5m = FMA(KP831469612, T5k, KP555570233 * T5l); + T5r = FNMS(KP555570233, T5k, KP831469612 * T5l); + } + T5n = T5j + T5m; + T5G = T5m - T5j; + T5s = T5q + T5r; + T5D = T5q - T5r; + } + } + { + E T5o, T5w, T5d, T5p; + T5o = T5g + T5n; + T5w = T5s + T5v; + T5d = W[10]; + T5p = W[11]; + rio[WS(ios, 6)] = FNMS(T5p, T5w, T5d * T5o); + iio[-WS(ios, 57)] = FMA(T5p, T5o, T5d * T5w); + } + { + E T5K, T5M, T5J, T5L; + T5K = T5C - T5D; + T5M = T5H - T5G; + T5J = W[42]; + T5L = W[43]; + rio[WS(ios, 22)] = FNMS(T5L, T5M, T5J * T5K); + iio[-WS(ios, 41)] = FMA(T5L, T5K, T5J * T5M); + } + { + E T5y, T5A, T5x, T5z; + T5y = T5g - T5n; + T5A = T5v - T5s; + T5x = W[74]; + T5z = W[75]; + rio[WS(ios, 38)] = FNMS(T5z, T5A, T5x * T5y); + iio[-WS(ios, 25)] = FMA(T5z, T5y, T5x * T5A); + } + { + E T5E, T5I, T5B, T5F; + T5E = T5C + T5D; + T5I = T5G + T5H; + T5B = W[106]; + T5F = W[107]; + rio[WS(ios, 54)] = FNMS(T5F, T5I, T5B * T5E); + iio[-WS(ios, 9)] = FMA(T5F, T5E, T5B * T5I); + } + } + { + E T5Y, T6w, T6p, T6B, T6d, T6A, T6i, T6x; + { + E T5Q, T5X, T6l, T6o; + T5Q = T5O + T5P; + T5X = T5T + T5W; + T5Y = T5Q + T5X; + T6w = T5Q - T5X; + T6l = T6j + T6k; + T6o = T6m + T6n; + T6p = T6l + T6o; + T6B = T6o - T6l; + { + E T65, T6g, T6c, T6h; + { + E T61, T64, T68, T6b; + T61 = T5Z + T60; + T64 = T62 + T63; + T65 = FNMS(KP195090322, T64, KP980785280 * T61); + T6g = FMA(KP195090322, T61, KP980785280 * T64); + T68 = T66 + T67; + T6b = T69 + T6a; + T6c = FMA(KP980785280, T68, KP195090322 * T6b); + T6h = FNMS(KP195090322, T68, KP980785280 * T6b); + } + T6d = T65 + T6c; + T6A = T6c - T65; + T6i = T6g + T6h; + T6x = T6g - T6h; + } + } + { + E T6e, T6q, T5N, T6f; + T6e = T5Y + T6d; + T6q = T6i + T6p; + T5N = W[2]; + T6f = W[3]; + rio[WS(ios, 2)] = FNMS(T6f, T6q, T5N * T6e); + iio[-WS(ios, 61)] = FMA(T6f, T6e, T5N * T6q); + } + { + E T6E, T6G, T6D, T6F; + T6E = T6w - T6x; + T6G = T6B - T6A; + T6D = W[34]; + T6F = W[35]; + rio[WS(ios, 18)] = FNMS(T6F, T6G, T6D * T6E); + iio[-WS(ios, 45)] = FMA(T6F, T6E, T6D * T6G); + } + { + E T6s, T6u, T6r, T6t; + T6s = T5Y - T6d; + T6u = T6p - T6i; + T6r = W[66]; + T6t = W[67]; + rio[WS(ios, 34)] = FNMS(T6t, T6u, T6r * T6s); + iio[-WS(ios, 29)] = FMA(T6t, T6s, T6r * T6u); + } + { + E T6y, T6C, T6v, T6z; + T6y = T6w + T6x; + T6C = T6A + T6B; + T6v = W[98]; + T6z = W[99]; + rio[WS(ios, 50)] = FNMS(T6z, T6C, T6v * T6y); + iio[-WS(ios, 13)] = FMA(T6z, T6y, T6v * T6C); + } + } + { + E TdO, Tf1, Teq, TeH, Tef, TeW, Ten, TeM, Te3, Ter, Te8, Tem, TeE, Tf0, TeP; + E TeX; + { + E TdG, TeG, TdN, TeF, TdJ, TdM; + TdG = TdE + TdF; + TeG = Ted - Tec; + TdJ = FNMS(KP555570233, TdI, KP831469612 * TdH); + TdM = FMA(KP831469612, TdK, KP555570233 * TdL); + TdN = TdJ + TdM; + TeF = TdM - TdJ; + TdO = TdG + TdN; + Tf1 = TeG - TeF; + Teq = TdG - TdN; + TeH = TeF + TeG; + } + { + E Tee, TeK, Teb, TeL, Te9, Tea; + Tee = Tec + Ted; + TeK = TdE - TdF; + Te9 = FMA(KP555570233, TdH, KP831469612 * TdI); + Tea = FNMS(KP555570233, TdK, KP831469612 * TdL); + Teb = Te9 + Tea; + TeL = Te9 - Tea; + Tef = Teb + Tee; + TeW = TeK - TeL; + Ten = Tee - Teb; + TeM = TeK + TeL; + } + { + E TdV, Te6, Te2, Te7; + { + E TdR, TdU, TdY, Te1; + TdR = TdP + TdQ; + TdU = TdS + TdT; + TdV = FNMS(KP290284677, TdU, KP956940335 * TdR); + Te6 = FMA(KP290284677, TdR, KP956940335 * TdU); + TdY = TdW + TdX; + Te1 = TdZ + Te0; + Te2 = FMA(KP956940335, TdY, KP290284677 * Te1); + Te7 = FNMS(KP290284677, TdY, KP956940335 * Te1); + } + Te3 = TdV + Te2; + Ter = Te6 - Te7; + Te8 = Te6 + Te7; + Tem = Te2 - TdV; + } + { + E TeA, TeN, TeD, TeO; + { + E Tey, Tez, TeB, TeC; + Tey = TdT - TdS; + Tez = TdP - TdQ; + TeA = FNMS(KP471396736, Tez, KP881921264 * Tey); + TeN = FMA(KP881921264, Tez, KP471396736 * Tey); + TeB = TdW - TdX; + TeC = Te0 - TdZ; + TeD = FMA(KP471396736, TeB, KP881921264 * TeC); + TeO = FNMS(KP471396736, TeC, KP881921264 * TeB); + } + TeE = TeA + TeD; + Tf0 = TeN - TeO; + TeP = TeN + TeO; + TeX = TeD - TeA; + } + { + E Te4, Teg, TdD, Te5; + Te4 = TdO + Te3; + Teg = Te8 + Tef; + TdD = W[120]; + Te5 = W[121]; + iio[-WS(ios, 2)] = FMA(TdD, Te4, Te5 * Teg); + rio[WS(ios, 61)] = FNMS(Te5, Te4, TdD * Teg); + } + { + E TeY, Tf2, TeV, TeZ; + TeY = TeW + TeX; + Tf2 = Tf0 + Tf1; + TeV = W[104]; + TeZ = W[105]; + iio[-WS(ios, 10)] = FMA(TeV, TeY, TeZ * Tf2); + rio[WS(ios, 53)] = FNMS(TeZ, TeY, TeV * Tf2); + } + { + E Tf4, Tf6, Tf3, Tf5; + Tf4 = Tf1 - Tf0; + Tf6 = TeW - TeX; + Tf3 = W[40]; + Tf5 = W[41]; + rio[WS(ios, 21)] = FNMS(Tf5, Tf6, Tf3 * Tf4); + iio[-WS(ios, 42)] = FMA(Tf3, Tf6, Tf5 * Tf4); + } + { + E Tei, Tek, Teh, Tej; + Tei = Tef - Te8; + Tek = TdO - Te3; + Teh = W[56]; + Tej = W[57]; + rio[WS(ios, 29)] = FNMS(Tej, Tek, Teh * Tei); + iio[-WS(ios, 34)] = FMA(Teh, Tek, Tej * Tei); + } + { + E Teo, Tes, Tel, Tep; + Teo = Tem + Ten; + Tes = Teq + Ter; + Tel = W[24]; + Tep = W[25]; + rio[WS(ios, 13)] = FNMS(Tep, Tes, Tel * Teo); + iio[-WS(ios, 50)] = FMA(Tel, Tes, Tep * Teo); + } + { + E TeI, TeQ, Tex, TeJ; + TeI = TeE + TeH; + TeQ = TeM + TeP; + Tex = W[8]; + TeJ = W[9]; + rio[WS(ios, 5)] = FNMS(TeJ, TeQ, Tex * TeI); + iio[-WS(ios, 58)] = FMA(Tex, TeQ, TeJ * TeI); + } + { + E TeS, TeU, TeR, TeT; + TeS = TeM - TeP; + TeU = TeH - TeE; + TeR = W[72]; + TeT = W[73]; + iio[-WS(ios, 26)] = FMA(TeR, TeS, TeT * TeU); + rio[WS(ios, 37)] = FNMS(TeT, TeS, TeR * TeU); + } + { + E Teu, Tew, Tet, Tev; + Teu = Teq - Ter; + Tew = Ten - Tem; + Tet = W[88]; + Tev = W[89]; + iio[-WS(ios, 18)] = FMA(Tet, Teu, Tev * Tew); + rio[WS(ios, 45)] = FNMS(Tev, Teu, Tet * Tew); + } + } + { + E Tcr, Tdw, TcX, Td6, TcI, Tdt, TcS, Tdl, Tbm, TcW, TcL, TcT, Tdd, Tdx, Tdi; + E Tds; + { + E Tcq, Td4, TbZ, Td5, TbF, TbY; + Tcq = Tce + Tcp; + Td4 = TcA - TcD; + TbF = FNMS(KP195090322, TbE, KP980785280 * Tbx); + TbY = FMA(KP195090322, TbQ, KP980785280 * TbX); + TbZ = TbF + TbY; + Td5 = TbY - TbF; + Tcr = TbZ + Tcq; + Tdw = Td4 - Td5; + TcX = Tcq - TbZ; + Td6 = Td4 + Td5; + } + { + E TcE, Tdk, TcH, Tdj, TcF, TcG; + TcE = TcA + TcD; + Tdk = Tcp - Tce; + TcF = FMA(KP980785280, TbE, KP195090322 * Tbx); + TcG = FNMS(KP195090322, TbX, KP980785280 * TbQ); + TcH = TcF + TcG; + Tdj = TcF - TcG; + TcI = TcE + TcH; + Tdt = Tdk - Tdj; + TcS = TcE - TcH; + Tdl = Tdj + Tdk; + } + { + E TaI, TcJ, Tbl, TcK; + { + E Taw, TaH, Tb9, Tbk; + Taw = Tak + Tav; + TaH = TaD + TaG; + TaI = FNMS(KP098017140, TaH, KP995184726 * Taw); + TcJ = FMA(KP995184726, TaH, KP098017140 * Taw); + Tb9 = TaT + Tb8; + Tbk = Tbc + Tbj; + Tbl = FMA(KP098017140, Tb9, KP995184726 * Tbk); + TcK = FNMS(KP098017140, Tbk, KP995184726 * Tb9); + } + Tbm = TaI + Tbl; + TcW = TcJ - TcK; + TcL = TcJ + TcK; + TcT = Tbl - TaI; + } + { + E Td9, Tdg, Tdc, Tdh; + { + E Td7, Td8, Tda, Tdb; + Td7 = TaD - TaG; + Td8 = Tav - Tak; + Td9 = FNMS(KP634393284, Td8, KP773010453 * Td7); + Tdg = FMA(KP634393284, Td7, KP773010453 * Td8); + Tda = TaT - Tb8; + Tdb = Tbj - Tbc; + Tdc = FMA(KP773010453, Tda, KP634393284 * Tdb); + Tdh = FNMS(KP634393284, Tda, KP773010453 * Tdb); + } + Tdd = Td9 + Tdc; + Tdx = Tdg - Tdh; + Tdi = Tdg + Tdh; + Tds = Tdc - Td9; + } + { + E Tcs, TcM, Ta5, Tct; + Tcs = Tbm + Tcr; + TcM = TcI + TcL; + Ta5 = W[0]; + Tct = W[1]; + rio[WS(ios, 1)] = FNMS(Tct, TcM, Ta5 * Tcs); + iio[-WS(ios, 62)] = FMA(Ta5, TcM, Tct * Tcs); + } + { + E Tdu, Tdy, Tdr, Tdv; + Tdu = Tds + Tdt; + Tdy = Tdw + Tdx; + Tdr = W[16]; + Tdv = W[17]; + rio[WS(ios, 9)] = FNMS(Tdv, Tdy, Tdr * Tdu); + iio[-WS(ios, 54)] = FMA(Tdr, Tdy, Tdv * Tdu); + } + { + E TdA, TdC, Tdz, TdB; + TdA = Tdw - Tdx; + TdC = Tdt - Tds; + Tdz = W[80]; + TdB = W[81]; + iio[-WS(ios, 22)] = FMA(Tdz, TdA, TdB * TdC); + rio[WS(ios, 41)] = FNMS(TdB, TdA, Tdz * TdC); + } + { + E TcO, TcQ, TcN, TcP; + TcO = TcI - TcL; + TcQ = Tcr - Tbm; + TcN = W[64]; + TcP = W[65]; + iio[-WS(ios, 30)] = FMA(TcN, TcO, TcP * TcQ); + rio[WS(ios, 33)] = FNMS(TcP, TcO, TcN * TcQ); + } + { + E TcU, TcY, TcR, TcV; + TcU = TcS + TcT; + TcY = TcW + TcX; + TcR = W[96]; + TcV = W[97]; + iio[-WS(ios, 14)] = FMA(TcR, TcU, TcV * TcY); + rio[WS(ios, 49)] = FNMS(TcV, TcU, TcR * TcY); + } + { + E Tde, Tdm, Td3, Tdf; + Tde = Td6 + Tdd; + Tdm = Tdi + Tdl; + Td3 = W[112]; + Tdf = W[113]; + iio[-WS(ios, 6)] = FMA(Td3, Tde, Tdf * Tdm); + rio[WS(ios, 57)] = FNMS(Tdf, Tde, Td3 * Tdm); + } + { + E Tdo, Tdq, Tdn, Tdp; + Tdo = Tdl - Tdi; + Tdq = Td6 - Tdd; + Tdn = W[48]; + Tdp = W[49]; + rio[WS(ios, 25)] = FNMS(Tdp, Tdq, Tdn * Tdo); + iio[-WS(ios, 38)] = FMA(Tdn, Tdq, Tdp * Tdo); + } + { + E Td0, Td2, TcZ, Td1; + Td0 = TcX - TcW; + Td2 = TcS - TcT; + TcZ = W[32]; + Td1 = W[33]; + rio[WS(ios, 17)] = FNMS(Td1, Td2, TcZ * Td0); + iio[-WS(ios, 46)] = FMA(TcZ, Td2, Td1 * Td0); + } + } + { + E Tfy, Thd, TgC, TgT, Tgr, Th8, Tgz, TgY, Tgb, TgD, Tgg, Tgy, TgQ, Thc, Th1; + E Th9; + { + E Tfi, TgS, Tfx, TgR, Tfp, Tfw; + Tfi = Tfa + Tfh; + TgS = Tgp - Tgm; + Tfp = FNMS(KP195090322, Tfo, KP980785280 * Tfl); + Tfw = FMA(KP980785280, Tfs, KP195090322 * Tfv); + Tfx = Tfp + Tfw; + TgR = Tfw - Tfp; + Tfy = Tfi + Tfx; + Thd = TgS - TgR; + TgC = Tfi - Tfx; + TgT = TgR + TgS; + } + { + E Tgq, TgW, Tgj, TgX, Tgh, Tgi; + Tgq = Tgm + Tgp; + TgW = Tfa - Tfh; + Tgh = FMA(KP195090322, Tfl, KP980785280 * Tfo); + Tgi = FNMS(KP195090322, Tfs, KP980785280 * Tfv); + Tgj = Tgh + Tgi; + TgX = Tgh - Tgi; + Tgr = Tgj + Tgq; + Th8 = TgW - TgX; + Tgz = Tgq - Tgj; + TgY = TgW + TgX; + } + { + E TfR, Tge, Tga, Tgf; + { + E TfJ, TfQ, Tg2, Tg9; + TfJ = TfB + TfI; + TfQ = TfM + TfP; + TfR = FNMS(KP098017140, TfQ, KP995184726 * TfJ); + Tge = FMA(KP098017140, TfJ, KP995184726 * TfQ); + Tg2 = TfU + Tg1; + Tg9 = Tg5 + Tg8; + Tga = FMA(KP995184726, Tg2, KP098017140 * Tg9); + Tgf = FNMS(KP098017140, Tg2, KP995184726 * Tg9); + } + Tgb = TfR + Tga; + TgD = Tge - Tgf; + Tgg = Tge + Tgf; + Tgy = Tga - TfR; + } + { + E TgM, TgZ, TgP, Th0; + { + E TgK, TgL, TgN, TgO; + TgK = TfP - TfM; + TgL = TfB - TfI; + TgM = FNMS(KP634393284, TgL, KP773010453 * TgK); + TgZ = FMA(KP773010453, TgL, KP634393284 * TgK); + TgN = TfU - Tg1; + TgO = Tg8 - Tg5; + TgP = FMA(KP634393284, TgN, KP773010453 * TgO); + Th0 = FNMS(KP634393284, TgO, KP773010453 * TgN); + } + TgQ = TgM + TgP; + Thc = TgZ - Th0; + Th1 = TgZ + Th0; + Th9 = TgP - TgM; + } + { + E Tgc, Tgs, Tf7, Tgd; + Tgc = Tfy + Tgb; + Tgs = Tgg + Tgr; + Tf7 = W[124]; + Tgd = W[125]; + iio[0] = FMA(Tf7, Tgc, Tgd * Tgs); + rio[WS(ios, 63)] = FNMS(Tgd, Tgc, Tf7 * Tgs); + } + { + E Tha, The, Th7, Thb; + Tha = Th8 + Th9; + The = Thc + Thd; + Th7 = W[108]; + Thb = W[109]; + iio[-WS(ios, 8)] = FMA(Th7, Tha, Thb * The); + rio[WS(ios, 55)] = FNMS(Thb, Tha, Th7 * The); + } + { + E Thg, Thi, Thf, Thh; + Thg = Thd - Thc; + Thi = Th8 - Th9; + Thf = W[44]; + Thh = W[45]; + rio[WS(ios, 23)] = FNMS(Thh, Thi, Thf * Thg); + iio[-WS(ios, 40)] = FMA(Thf, Thi, Thh * Thg); + } + { + E Tgu, Tgw, Tgt, Tgv; + Tgu = Tgr - Tgg; + Tgw = Tfy - Tgb; + Tgt = W[60]; + Tgv = W[61]; + rio[WS(ios, 31)] = FNMS(Tgv, Tgw, Tgt * Tgu); + iio[-WS(ios, 32)] = FMA(Tgt, Tgw, Tgv * Tgu); + } + { + E TgA, TgE, Tgx, TgB; + TgA = Tgy + Tgz; + TgE = TgC + TgD; + Tgx = W[28]; + TgB = W[29]; + rio[WS(ios, 15)] = FNMS(TgB, TgE, Tgx * TgA); + iio[-WS(ios, 48)] = FMA(Tgx, TgE, TgB * TgA); + } + { + E TgU, Th2, TgJ, TgV; + TgU = TgQ + TgT; + Th2 = TgY + Th1; + TgJ = W[12]; + TgV = W[13]; + rio[WS(ios, 7)] = FNMS(TgV, Th2, TgJ * TgU); + iio[-WS(ios, 56)] = FMA(TgJ, Th2, TgV * TgU); + } + { + E Th4, Th6, Th3, Th5; + Th4 = TgY - Th1; + Th6 = TgT - TgQ; + Th3 = W[76]; + Th5 = W[77]; + iio[-WS(ios, 24)] = FMA(Th3, Th4, Th5 * Th6); + rio[WS(ios, 39)] = FNMS(Th5, Th4, Th3 * Th6); + } + { + E TgG, TgI, TgF, TgH; + TgG = TgC - TgD; + TgI = Tgz - Tgy; + TgF = W[92]; + TgH = W[93]; + iio[-WS(ios, 16)] = FMA(TgF, TgG, TgH * TgI); + rio[WS(ios, 47)] = FNMS(TgH, TgG, TgF * TgI); + } + } + { + E ThJ, TiG, Ti7, Tig, ThS, TiD, Ti2, Tiv, Thy, Ti6, ThV, Ti3, Tin, TiH, Tis; + E TiC; + { + E ThI, Tie, ThF, Tif, ThB, ThE; + ThI = ThG + ThH; + Tie = ThM - ThN; + ThB = FNMS(KP555570233, ThA, KP831469612 * Thz); + ThE = FNMS(KP555570233, ThD, KP831469612 * ThC); + ThF = ThB + ThE; + Tif = ThE - ThB; + ThJ = ThF + ThI; + TiG = Tie - Tif; + Ti7 = ThI - ThF; + Tig = Tie + Tif; + } + { + E ThO, Tiu, ThR, Tit, ThP, ThQ; + ThO = ThM + ThN; + Tiu = ThH - ThG; + ThP = FMA(KP831469612, ThA, KP555570233 * Thz); + ThQ = FMA(KP831469612, ThD, KP555570233 * ThC); + ThR = ThP - ThQ; + Tit = ThP + ThQ; + ThS = ThO + ThR; + TiD = Tiu - Tit; + Ti2 = ThO - ThR; + Tiv = Tit + Tiu; + } + { + E Thq, ThT, Thx, ThU; + { + E Thm, Thp, Tht, Thw; + Thm = Thk + Thl; + Thp = Thn + Tho; + Thq = FNMS(KP290284677, Thp, KP956940335 * Thm); + ThT = FMA(KP956940335, Thp, KP290284677 * Thm); + Tht = Thr - Ths; + Thw = Thu + Thv; + Thx = FMA(KP290284677, Tht, KP956940335 * Thw); + ThU = FNMS(KP290284677, Thw, KP956940335 * Tht); + } + Thy = Thq + Thx; + Ti6 = ThT - ThU; + ThV = ThT + ThU; + Ti3 = Thx - Thq; + } + { + E Tij, Tiq, Tim, Tir; + { + E Tih, Tii, Tik, Til; + Tih = Thn - Tho; + Tii = Thl - Thk; + Tij = FNMS(KP471396736, Tii, KP881921264 * Tih); + Tiq = FMA(KP471396736, Tih, KP881921264 * Tii); + Tik = Thv - Thu; + Til = Ths + Thr; + Tim = FNMS(KP881921264, Til, KP471396736 * Tik); + Tir = FMA(KP471396736, Til, KP881921264 * Tik); + } + Tin = Tij + Tim; + TiH = Tiq - Tir; + Tis = Tiq + Tir; + TiC = Tim - Tij; + } + { + E ThK, ThW, Thj, ThL; + ThK = Thy + ThJ; + ThW = ThS + ThV; + Thj = W[4]; + ThL = W[5]; + rio[WS(ios, 3)] = FNMS(ThL, ThW, Thj * ThK); + iio[-WS(ios, 60)] = FMA(Thj, ThW, ThL * ThK); + } + { + E TiE, TiI, TiB, TiF; + TiE = TiC + TiD; + TiI = TiG + TiH; + TiB = W[20]; + TiF = W[21]; + rio[WS(ios, 11)] = FNMS(TiF, TiI, TiB * TiE); + iio[-WS(ios, 52)] = FMA(TiB, TiI, TiF * TiE); + } + { + E TiK, TiM, TiJ, TiL; + TiK = TiG - TiH; + TiM = TiD - TiC; + TiJ = W[84]; + TiL = W[85]; + iio[-WS(ios, 20)] = FMA(TiJ, TiK, TiL * TiM); + rio[WS(ios, 43)] = FNMS(TiL, TiK, TiJ * TiM); + } + { + E ThY, Ti0, ThX, ThZ; + ThY = ThS - ThV; + Ti0 = ThJ - Thy; + ThX = W[68]; + ThZ = W[69]; + iio[-WS(ios, 28)] = FMA(ThX, ThY, ThZ * Ti0); + rio[WS(ios, 35)] = FNMS(ThZ, ThY, ThX * Ti0); + } + { + E Ti4, Ti8, Ti1, Ti5; + Ti4 = Ti2 + Ti3; + Ti8 = Ti6 + Ti7; + Ti1 = W[100]; + Ti5 = W[101]; + iio[-WS(ios, 12)] = FMA(Ti1, Ti4, Ti5 * Ti8); + rio[WS(ios, 51)] = FNMS(Ti5, Ti4, Ti1 * Ti8); + } + { + E Tio, Tiw, Tid, Tip; + Tio = Tig + Tin; + Tiw = Tis + Tiv; + Tid = W[116]; + Tip = W[117]; + iio[-WS(ios, 4)] = FMA(Tid, Tio, Tip * Tiw); + rio[WS(ios, 59)] = FNMS(Tip, Tio, Tid * Tiw); + } + { + E Tiy, TiA, Tix, Tiz; + Tiy = Tiv - Tis; + TiA = Tig - Tin; + Tix = W[52]; + Tiz = W[53]; + rio[WS(ios, 27)] = FNMS(Tiz, TiA, Tix * Tiy); + iio[-WS(ios, 36)] = FMA(Tix, TiA, Tiz * Tiy); + } + { + E Tia, Tic, Ti9, Tib; + Tia = Ti7 - Ti6; + Tic = Ti2 - Ti3; + Ti9 = W[36]; + Tib = W[37]; + rio[WS(ios, 19)] = FNMS(Tib, Tic, Ti9 * Tia); + iio[-WS(ios, 44)] = FMA(Ti9, Tic, Tib * Tia); + } + } + } + return W; +} + +static const tw_instr twinstr[] = { + {TW_FULL, 0, 64}, + {TW_NEXT, 1, 0} +}; + +static const hc2hc_desc desc = { 64, "hb_64", twinstr, {808, 270, 230, 0}, &GENUS, 0, 0, 0 }; + +void X(codelet_hb_64) (planner *p) { + X(khc2hc_dif_register) (p, hb_64, &desc); +} |