Logo Search packages:      
Sourcecode: fftw3 version File versions  Download package

hb_64.c

/*
 * 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: algsimp.ml,v 1.7 2003/03/15 20:29:42 stevenj Exp $
 * $Id: fft.ml,v 1.2 2003/03/15 20:29:42 stevenj Exp $
 * $Id: gen_hc2hc.ml,v 1.9 2003/04/17 19:25:50 athena 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);
}

Generated by  Doxygen 1.6.0   Back to index