summaryrefslogtreecommitdiff
path: root/src/fftw3/rdft/codelets/hc2r/hb_64.c
diff options
context:
space:
mode:
authorscuri <scuri>2009-08-20 12:35:06 +0000
committerscuri <scuri>2009-08-20 12:35:06 +0000
commit5d735255ddd3cb2f547abd3d03969af3fb7eb04e (patch)
tree8fb66510bc625bb1b08ccb41f1b83fb0f7cb8f19 /src/fftw3/rdft/codelets/hc2r/hb_64.c
parent35733b87eed86e5228f12fa10c98a3d9d22a6073 (diff)
*** empty log message ***
Diffstat (limited to 'src/fftw3/rdft/codelets/hc2r/hb_64.c')
-rw-r--r--src/fftw3/rdft/codelets/hc2r/hb_64.c1972
1 files changed, 0 insertions, 1972 deletions
diff --git a/src/fftw3/rdft/codelets/hc2r/hb_64.c b/src/fftw3/rdft/codelets/hc2r/hb_64.c
deleted file mode 100644
index d3f9afc..0000000
--- a/src/fftw3/rdft/codelets/hc2r/hb_64.c
+++ /dev/null
@@ -1,1972 +0,0 @@
-/*
- * 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);
-}