mhc2r_128.c

fftw-3.0.1

		    T5e = T5c - T5d;
		    T5h = FNMS(KP1_913880671, T5g, KP580569354 * T5f);
		    O[WS(os, 90)] = T5e - T5h;
		    O[WS(os, 26)] = T5e + T5h;
	       }
	  }
	  {
	       E T7o, T7C, T7y, T7G, T7r, T7D, T7v, T7F;
	       {
		    E T7m, T7n, T7w, T7x;
		    T7m = T6q + T6r;
		    T7n = KP1_414213562 * (T6t + T6w);
		    T7o = T7m - T7n;
		    T7C = T7m + T7n;
		    T7w = KP707106781 * (T6U - T6Z);
		    T7x = T78 - T75;
		    T7y = T7w + T7x;
		    T7G = T7x - T7w;
	       }
	       {
		    E T7p, T7q, T7t, T7u;
		    T7p = T6z - T6C;
		    T7q = T6H - T6E;
		    T7r = FNMS(KP1_847759065, T7q, KP765366864 * T7p);
		    T7D = FMA(KP1_847759065, T7p, KP765366864 * T7q);
		    T7t = T6L - T6O;
		    T7u = KP707106781 * (T73 - T72);
		    T7v = T7t + T7u;
		    T7F = T7t - T7u;
	       }
	       {
		    E T7s, T7z, T7I, T7J;
		    T7s = T7o + T7r;
		    T7z = FNMS(KP1_111140466, T7y, KP1_662939224 * T7v);
		    O[WS(os, 76)] = T7s - T7z;
		    O[WS(os, 12)] = T7s + T7z;
		    T7I = T7C + T7D;
		    T7J = FMA(KP1_961570560, T7F, KP390180644 * T7G);
		    O[WS(os, 60)] = T7I - T7J;
		    O[WS(os, 124)] = T7I + T7J;
	       }
	       {
		    E T7A, T7B, T7E, T7H;
		    T7A = T7o - T7r;
		    T7B = FMA(KP1_111140466, T7v, KP1_662939224 * T7y);
		    O[WS(os, 44)] = T7A - T7B;
		    O[WS(os, 108)] = T7A + T7B;
		    T7E = T7C - T7D;
		    T7H = FNMS(KP1_961570560, T7G, KP390180644 * T7F);
		    O[WS(os, 92)] = T7E - T7H;
		    O[WS(os, 28)] = T7E + T7H;
	       }
	  }
	  {
	       E T5q, T5U, T5Q, T5Y, T5x, T5V, T5J, T5X;
	       {
		    E T5m, T5p, T5M, T5P;
		    T5m = T5k - T5l;
		    T5p = FNMS(KP1_847759065, T5o, KP765366864 * T5n);
		    T5q = T5m + T5p;
		    T5U = T5m - T5p;
		    T5M = T5K + T5L;
		    T5P = T5N + T5O;
		    T5Q = T5M + T5P;
		    T5Y = T5P - T5M;
	       }
	       {
		    E T5t, T5w, T5B, T5I;
		    T5t = T5r + T5s;
		    T5w = T5u + T5v;
		    T5x = FNMS(KP1_111140466, T5w, KP1_662939224 * T5t);
		    T5V = FMA(KP1_111140466, T5t, KP1_662939224 * T5w);
		    T5B = T5z + T5A;
		    T5I = T5E + T5H;
		    T5J = T5B + T5I;
		    T5X = T5B - T5I;
	       }
	       {
		    E T5y, T5R, T60, T61;
		    T5y = T5q + T5x;
		    T5R = FNMS(KP580569354, T5Q, KP1_913880671 * T5J);
		    O[WS(os, 70)] = T5y - T5R;
		    O[WS(os, 6)] = T5y + T5R;
		    T60 = T5U + T5V;
		    T61 = FMA(KP1_763842528, T5X, KP942793473 * T5Y);
		    O[WS(os, 54)] = T60 - T61;
		    O[WS(os, 118)] = T60 + T61;
	       }
	       {
		    E T5S, T5T, T5W, T5Z;
		    T5S = T5q - T5x;
		    T5T = FMA(KP580569354, T5J, KP1_913880671 * T5Q);
		    O[WS(os, 38)] = T5S - T5T;
		    O[WS(os, 102)] = T5S + T5T;
		    T5W = T5U - T5V;
		    T5Z = FNMS(KP1_763842528, T5Y, KP942793473 * T5X);
		    O[WS(os, 86)] = T5W - T5Z;
		    O[WS(os, 22)] = T5W + T5Z;
	       }
	  }
	  {
	       E T64, T6i, T6e, T6m, T67, T6j, T6b, T6l;
	       {
		    E T62, T63, T6c, T6d;
		    T62 = T5k + T5l;
		    T63 = FMA(KP1_847759065, T5n, KP765366864 * T5o);
		    T64 = T62 - T63;
		    T6i = T62 + T63;
		    T6c = T5E - T5H;
		    T6d = T5O - T5N;
		    T6e = T6c + T6d;
		    T6m = T6d - T6c;
	       }
	       {
		    E T65, T66, T69, T6a;
		    T65 = T5r - T5s;
		    T66 = T5v - T5u;
		    T67 = FNMS(KP1_961570560, T66, KP390180644 * T65);
		    T6j = FMA(KP1_961570560, T65, KP390180644 * T66);
		    T69 = T5z - T5A;
		    T6a = T5L - T5K;
		    T6b = T69 + T6a;
		    T6l = T69 - T6a;
	       }
	       {
		    E T68, T6f, T6o, T6p;
		    T68 = T64 + T67;
		    T6f = FNMS(KP1_268786568, T6e, KP1_546020906 * T6b);
		    O[WS(os, 78)] = T68 - T6f;
		    O[WS(os, 14)] = T68 + T6f;
		    T6o = T6i + T6j;
		    T6p = FMA(KP1_990369453, T6l, KP196034280 * T6m);
		    O[WS(os, 62)] = T6o - T6p;
		    O[WS(os, 126)] = T6o + T6p;
	       }
	       {
		    E T6g, T6h, T6k, T6n;
		    T6g = T64 - T67;
		    T6h = FMA(KP1_268786568, T6b, KP1_546020906 * T6e);
		    O[WS(os, 46)] = T6g - T6h;
		    O[WS(os, 110)] = T6g + T6h;
		    T6k = T6i - T6j;
		    T6n = FNMS(KP1_990369453, T6m, KP196034280 * T6l);
		    O[WS(os, 94)] = T6k - T6n;
		    O[WS(os, 30)] = T6k + T6n;
	       }
	  }
	  {
	       E T8Z, Tb1, T9C, Tb2, Tbe, Tbq, Tbb, Tbp, TaX, Tbs, Tb5, Tbi, TaI, Tbt, Tb4;
	       E Tbl;
	       {
		    E T8F, T8Y, Tb9, Tba;
		    T8F = T8x + T8E;
		    T8Y = FNMS(KP390180644, T8X, KP1_961570560 * T8Q);
		    T8Z = T8F + T8Y;
		    Tb1 = T8F - T8Y;
		    {
			 E T9q, T9B, Tbc, Tbd;
			 T9q = T9a + T9p;
			 T9B = T9t + T9A;
			 T9C = FNMS(KP196034280, T9B, KP1_990369453 * T9q);
			 Tb2 = FMA(KP196034280, T9q, KP1_990369453 * T9B);
			 Tbc = T9a - T9p;
			 Tbd = T9A - T9t;
			 Tbe = FNMS(KP1_546020906, Tbd, KP1_268786568 * Tbc);
			 Tbq = FMA(KP1_546020906, Tbc, KP1_268786568 * Tbd);
		    }
		    Tb9 = T8x - T8E;
		    Tba = FMA(KP390180644, T8Q, KP1_961570560 * T8X);
		    Tbb = Tb9 - Tba;
		    Tbp = Tb9 + Tba;
		    {
			 E TaW, Tbg, TaL, Tbh, TaJ, TaK;
			 TaW = TaO + TaV;
			 Tbg = T9O - Ta3;
			 TaJ = FMA(KP195090322, Taf, KP980785280 * Tam);
			 TaK = FNMS(KP195090322, Tay, KP980785280 * TaF);
			 TaL = TaJ + TaK;
			 Tbh = TaK - TaJ;
			 TaX = TaL + TaW;
			 Tbs = Tbg - Tbh;
			 Tb5 = TaW - TaL;
			 Tbi = Tbg + Tbh;
		    }
		    {
			 E Ta4, Tbk, TaH, Tbj, Tan, TaG;
			 Ta4 = T9O + Ta3;
			 Tbk = TaV - TaO;
			 Tan = FNMS(KP195090322, Tam, KP980785280 * Taf);
			 TaG = FMA(KP980785280, Tay, KP195090322 * TaF);
			 TaH = Tan + TaG;
			 Tbj = Tan - TaG;
			 TaI = Ta4 + TaH;
			 Tbt = Tbk - Tbj;
			 Tb4 = Ta4 - TaH;
			 Tbl = Tbj + Tbk;
		    }
	       }
	       {
		    E T9D, TaY, Tbr, Tbu;
		    T9D = T8Z + T9C;
		    TaY = FNMS(KP098135348, TaX, KP1_997590912 * TaI);
		    O[WS(os, 65)] = T9D - TaY;
		    O[WS(os, 1)] = T9D + TaY;
		    Tbr = Tbp - Tbq;
		    Tbu = FNMS(KP1_883088130, Tbt, KP673779706 * Tbs);
		    O[WS(os, 89)] = Tbr - Tbu;
		    O[WS(os, 25)] = Tbr + Tbu;
	       }
	       {
		    E Tbv, Tbw, TaZ, Tb0;
		    Tbv = Tbp + Tbq;
		    Tbw = FMA(KP1_883088130, Tbs, KP673779706 * Tbt);
		    O[WS(os, 57)] = Tbv - Tbw;
		    O[WS(os, 121)] = Tbv + Tbw;
		    TaZ = T8Z - T9C;
		    Tb0 = FMA(KP098135348, TaI, KP1_997590912 * TaX);
		    O[WS(os, 33)] = TaZ - Tb0;
		    O[WS(os, 97)] = TaZ + Tb0;
	       }
	       {
		    E Tb3, Tb6, Tbf, Tbm;
		    Tb3 = Tb1 - Tb2;
		    Tb6 = FNMS(KP1_481902250, Tb5, KP1_343117909 * Tb4);
		    O[WS(os, 81)] = Tb3 - Tb6;
		    O[WS(os, 17)] = Tb3 + Tb6;
		    Tbf = Tbb + Tbe;
		    Tbm = FNMS(KP855110186, Tbl, KP1_807978586 * Tbi);
		    O[WS(os, 73)] = Tbf - Tbm;
		    O[WS(os, 9)] = Tbf + Tbm;
	       }
	       {
		    E Tbn, Tbo, Tb7, Tb8;
		    Tbn = Tbb - Tbe;
		    Tbo = FMA(KP855110186, Tbi, KP1_807978586 * Tbl);
		    O[WS(os, 41)] = Tbn - Tbo;
		    O[WS(os, 105)] = Tbn + Tbo;
		    Tb7 = Tb1 + Tb2;
		    Tb8 = FMA(KP1_481902250, Tb4, KP1_343117909 * Tb5);
		    O[WS(os, 49)] = Tb7 - Tb8;
		    O[WS(os, 113)] = Tb7 + Tb8;
	       }
	  }
	  {
	       E TcR, TdR, Tda, TdS, Te4, Teg, Te1, Tef, TdN, Tei, TdV, Te8, TdC, Tej, TdU;
	       E Teb;
	       {
		    E TcJ, TcQ, TdZ, Te0;
		    TcJ = TcF + TcI;
		    TcQ = FNMS(KP1_111140466, TcP, KP1_662939224 * TcM);
		    TcR = TcJ + TcQ;
		    TdR = TcJ - TcQ;
		    {
			 E Td2, Td9, Te2, Te3;
			 Td2 = TcU + Td1;
			 Td9 = Td5 + Td8;
			 Tda = FNMS(KP580569354, Td9, KP1_913880671 * Td2);
			 TdS = FMA(KP580569354, Td2, KP1_913880671 * Td9);
			 Te2 = TcU - Td1;
			 Te3 = Td8 - Td5;
			 Te4 = FNMS(KP1_763842528, Te3, KP942793473 * Te2);
			 Teg = FMA(KP1_763842528, Te2, KP942793473 * Te3);
		    }
		    TdZ = TcF - TcI;
		    Te0 = FMA(KP1_111140466, TcM, KP1_662939224 * TcP);
		    Te1 = TdZ - Te0;
		    Tef = TdZ + Te0;
		    {
			 E TdM, Te6, TdF, Te7, TdD, TdE;
			 TdM = TdI + TdL;
			 Te6 = Tde - Tdl;
			 TdD = FMA(KP555570233, Tdp, KP831469612 * Tds);
			 TdE = FNMS(KP555570233, Tdw, KP831469612 * Tdz);
			 TdF = TdD + TdE;
			 Te7 = TdE - TdD;
			 TdN = TdF + TdM;
			 Tei = Te6 - Te7;
			 TdV = TdM - TdF;
			 Te8 = Te6 + Te7;
		    }
		    {
			 E Tdm, Tea, TdB, Te9, Tdt, TdA;
			 Tdm = Tde + Tdl;
			 Tea = TdL - TdI;
			 Tdt = FNMS(KP555570233, Tds, KP831469612 * Tdp);
			 TdA = FMA(KP831469612, Tdw, KP555570233 * Tdz);
			 TdB = Tdt + TdA;
			 Te9 = Tdt - TdA;
			 TdC = Tdm + TdB;
			 Tej = Tea - Te9;
			 TdU = Tdm - TdB;
			 Teb = Te9 + Tea;
		    }
	       }
	       {
		    E Tdb, TdO, Teh, Tek;
		    Tdb = TcR + Tda;
		    TdO = FNMS(KP293460948, TdN, KP1_978353019 * TdC);
		    O[WS(os, 67)] = Tdb - TdO;
		    O[WS(os, 3)] = Tdb + TdO;
		    Teh = Tef - Teg;
		    Tek = FNMS(KP1_940062506, Tej, KP485960359 * Tei);
		    O[WS(os, 91)] = Teh - Tek;
		    O[WS(os, 27)] = Teh + Tek;
	       }
	       {
		    E Tel, Tem, TdP, TdQ;
		    Tel = Tef + Teg;
		    Tem = FMA(KP1_940062506, Tei, KP485960359 * Tej);
		    O[WS(os, 59)] = Tel - Tem;
		    O[WS(os, 123)] = Tel + Tem;
		    TdP = TcR - Tda;
		    TdQ = FMA(KP293460948, TdC, KP1_978353019 * TdN);
		    O[WS(os, 35)] = TdP - TdQ;
		    O[WS(os, 99)] = TdP + TdQ;
	       }
	       {
		    E TdT, TdW, Te5, Tec;
		    TdT = TdR - TdS;
		    TdW = FNMS(KP1_606415062, TdV, KP1_191398608 * TdU);
		    O[WS(os, 83)] = TdT - TdW;
		    O[WS(os, 19)] = TdT + TdW;
		    Te5 = Te1 + Te4;
		    Tec = FNMS(KP1_028205488, Teb, KP1_715457220 * Te8);
		    O[WS(os, 75)] = Te5 - Tec;
		    O[WS(os, 11)] = Te5 + Tec;
	       }
	       {
		    E Ted, Tee, TdX, TdY;
		    Ted = Te1 - Te4;
		    Tee = FMA(KP1_028205488, Te8, KP1_715457220 * Teb);
		    O[WS(os, 43)] = Ted - Tee;
		    O[WS(os, 107)] = Ted + Tee;
		    TdX = TdR + TdS;
		    TdY = FMA(KP1_606415062, TdU, KP1_191398608 * TdV);
		    O[WS(os, 51)] = TdX - TdY;
		    O[WS(os, 115)] = TdX + TdY;
	       }
	  }
	  {
	       E TbD, Tc7, TbK, Tc8, Tck, Tcw, Tch, Tcv, Tc3, Tcy, Tcb, Tco, TbW, Tcz, Tca;
	       E Tcr;
	       {
		    E Tbz, TbC, Tcf, Tcg;
		    Tbz = Tbx - Tby;
		    TbC = FNMS(KP1_662939224, TbB, KP1_111140466 * TbA);
		    TbD = Tbz + TbC;
		    Tc7 = Tbz - TbC;
		    {
			 E TbG, TbJ, Tci, Tcj;
			 TbG = TbE + TbF;
			 TbJ = TbH + TbI;
			 TbK = FNMS(KP942793473, TbJ, KP1_763842528 * TbG);
			 Tc8 = FMA(KP942793473, TbG, KP1_763842528 * TbJ);
			 Tci = TbE - TbF;
			 Tcj = TbI - TbH;
			 Tck = FNMS(KP1_913880671, Tcj, KP580569354 * Tci);
			 Tcw = FMA(KP1_913880671, Tci, KP580569354 * Tcj);
		    }
		    Tcf = Tbx + Tby;
		    Tcg = FMA(KP1_662939224, TbA, KP1_111140466 * TbB);
		    Tch = Tcf - Tcg;
		    Tcv = Tcf + Tcg;
		    {
			 E Tc2, Tcm, TbZ, Tcn, TbX, TbY;
			 Tc2 = Tc0 + Tc1;
			 Tcm = TbM - TbN;
			 TbX = FMA(KP831469612, TbP, KP555570233 * TbQ);
			 TbY = FNMS(KP831469612, TbS, KP555570233 * TbT);
			 TbZ = TbX + TbY;
			 Tcn = TbY - TbX;
			 Tc3 = TbZ + Tc2;
			 Tcy = Tcm - Tcn;
			 Tcb = Tc2 - TbZ;
			 Tco = Tcm + Tcn;
		    }
		    {
			 E TbO, Tcq, TbV, Tcp, TbR, TbU;
			 TbO = TbM + TbN;
			 Tcq = Tc1 - Tc0;
			 TbR = FNMS(KP831469612, TbQ, KP555570233 * TbP);
			 TbU = FMA(KP555570233, TbS, KP831469612 * TbT);
			 TbV = TbR + TbU;
			 Tcp = TbR - TbU;
			 TbW = TbO + TbV;
			 Tcz = Tcq - Tcp;
			 Tca = TbO - TbV;
			 Tcr = Tcp + Tcq;
		    }
	       }
	       {
		    E TbL, Tc4, Tcx, TcA;
		    TbL = TbD + TbK;
		    Tc4 = FNMS(KP485960359, Tc3, KP1_940062506 * TbW);
		    O[WS(os, 69)] = TbL - Tc4;
		    O[WS(os, 5)] = TbL + Tc4;
		    Tcx = Tcv - Tcw;
		    TcA = FNMS(KP1_978353019, Tcz, KP293460948 * Tcy);
		    O[WS(os, 93)] = Tcx - TcA;
		    O[WS(os, 29)] = Tcx + TcA;
	       }
	       {
		    E TcB, TcC, Tc5, Tc6;
		    TcB = Tcv + Tcw;
		    TcC = FMA(KP1_978353019, Tcy, KP293460948 * Tcz);
		    O[WS(os, 61)] = TcB - TcC;
		    O[WS(os, 125)] = TcB + TcC;
		    Tc5 = TbD - TbK;
		    Tc6 = FMA(KP485960359, TbW, KP1_940062506 * Tc3);
		    O[WS(os, 37)] = Tc5 - Tc6;
		    O[WS(os, 101)] = Tc5 + Tc6;
	       }
	       {
		    E Tc9, Tcc, Tcl, Tcs;
		    Tc9 = Tc7 - Tc8;
		    Tcc = FNMS(KP1_715457220, Tcb, KP1_028205488 * Tca);
		    O[WS(os, 85)] = Tc9 - Tcc;
		    O[WS(os, 21)] = Tc9 + Tcc;
		    Tcl = Tch + Tck;
		    Tcs = FNMS(KP1_191398608, Tcr, KP1_606415062 * Tco);
		    O[WS(os, 77)] = Tcl - Tcs;
		    O[WS(os, 13)] = Tcl + Tcs;
	       }
	       {
		    E Tct, Tcu, Tcd, Tce;
		    Tct = Tch - Tck;
		    Tcu = FMA(KP1_191398608, Tco, KP1_606415062 * Tcr);
		    O[WS(os, 45)] = Tct - Tcu;
		    O[WS(os, 109)] = Tct + Tcu;
		    Tcd = Tc7 + Tc8;
		    Tce = FMA(KP1_715457220, Tca, KP1_028205488 * Tcb);
		    O[WS(os, 53)] = Tcd - Tce;
		    O[WS(os, 117)] = Tcd + Tce;
	       }
	  }
	  {
	       E Tet, TeX, TeA, TeY, Tfa, Tfm, Tf7, Tfl, TeT, Tfo, Tf1, Tfe, TeM, Tfp, Tf0;
	       E Tfh;
	       {
		    E Tep, Tes, Tf5, Tf6;
		    Tep = Ten - Teo;
		    Tes = FNMS(KP1_961570560, Ter, KP390180644 * Teq);
		    Tet = Tep + Tes;
		    TeX = Tep - Tes;
		    {
			 E Tew, Tez, Tf8, Tf9;
			 Tew = Teu - Tev;
			 Tez = Tex + Tey;
			 TeA = FNMS(KP1_268786568, Tez, KP1_546020906 * Tew);
			 TeY = FMA(KP1_268786568, Tew, KP1_546020906 * Tez);
			 Tf8 = Teu + Tev;
			 Tf9 = Tey - Tex;
			 Tfa = FNMS(KP1_990369453, Tf9, KP196034280 * Tf8);
			 Tfm = FMA(KP1_990369453, Tf8, KP196034280 * Tf9);
		    }
		    Tf5 = Ten + Teo;
		    Tf6 = FMA(KP1_961570560, Teq, KP390180644 * Ter);
		    Tf7 = Tf5 - Tf6;
		    Tfl = Tf5 + Tf6;
		    {
			 E TeS, Tfc, TeP, Tfd, TeN, TeO;
			 TeS = TeQ + TeR;
			 Tfc = TeC + TeD;
			 TeN = FMA(KP980785280, TeF, KP195090322 * TeG);
			 TeO = FMA(KP980785280, TeI, KP195090322 * TeJ);
			 TeP = TeN - TeO;
			 Tfd = TeN + TeO;
			 TeT = TeP + TeS;
			 Tfo = Tfc + Tfd;
			 Tf1 = TeS - TeP;
			 Tfe = Tfc - Tfd;
		    }
		    {
			 E TeE, Tfg, TeL, Tff, TeH, TeK;
			 TeE = TeC - TeD;
			 Tfg = TeR - TeQ;
			 TeH = FNMS(KP980785280, TeG, KP195090322 * TeF);
			 TeK = FNMS(KP980785280, TeJ, KP195090322 * TeI);
			 TeL = TeH + TeK;
			 Tff = TeH - TeK;
			 TeM = TeE + TeL;
			 Tfp = Tfg - Tff;
			 Tf0 = TeE - TeL;
			 Tfh = Tff + Tfg;
		    }
	       }
	       {
		    E TeB, TeU, Tfn, Tfq;
		    TeB = Tet + TeA;
		    TeU = FNMS(KP673779706, TeT, KP1_883088130 * TeM);
		    O[WS(os, 71)] = TeB - TeU;
		    O[WS(os, 7)] = TeB + TeU;
		    Tfn = Tfl - Tfm;
		    Tfq = FNMS(KP1_997590912, Tfp, KP098135348 * Tfo);
		    O[WS(os, 95)] = Tfn - Tfq;
		    O[WS(os, 31)] = Tfn + Tfq;
	       }
	       {
		    E Tfr, Tfs, TeV, TeW;
		    Tfr = Tfl + Tfm;
		    Tfs = FMA(KP1_997590912, Tfo, KP098135348 * Tfp);
		    O[WS(os, 63)] = Tfr - Tfs;
		    O[WS(os, 127)] = Tfr + Tfs;
		    TeV = Tet - TeA;
		    TeW = FMA(KP673779706, TeM, KP1_883088130 * TeT);
		    O[WS(os, 39)] = TeV - TeW;
		    O[WS(os, 103)] = TeV + TeW;
	       }
	       {
		    E TeZ, Tf2, Tfb, Tfi;
		    TeZ = TeX - TeY;
		    Tf2 = FNMS(KP1_807978586, Tf1, KP855110186 * Tf0);
		    O[WS(os, 87)] = TeZ - Tf2;
		    O[WS(os, 23)] = TeZ + Tf2;
		    Tfb = Tf7 + Tfa;
		    Tfi = FNMS(KP1_343117909, Tfh, KP1_481902250 * Tfe);
		    O[WS(os, 79)] = Tfb - Tfi;
		    O[WS(os, 15)] = Tfb + Tfi;
	       }
	       {
		    E Tfj, Tfk, Tf3, Tf4;
		    Tfj = Tf7 - Tfa;
		    Tfk = FMA(KP1_343117909, Tfe, KP1_481902250 * Tfh);
		    O[WS(os, 47)] = Tfj - Tfk;
		    O[WS(os, 111)] = Tfj + Tfk;
		    Tf3 = TeX + TeY;
		    Tf4 = FMA(KP1_807978586, Tf0, KP855110186 * Tf1);
		    O[WS(os, 55)] = Tf3 - Tf4;
		    O[WS(os, 119)] = Tf3 + Tf4;
	       }
	  }
     }
}

static void mhc2r_128(const R *ri, const R *ii, R *O, stride ris, stride iis, stride os, int v, int ivs, int ovs)
{
     int i;
     for (i = v; i > 0; --i) {
	  mhc2r_128_0(ri, ii, O, ris, iis, os);
	  ri += ivs;
	  ii += ivs;
	  O += ovs;
     }
}

static const khc2r_desc desc = { 128, "mhc2r_128", {812, 198, 144, 0}, &GENUS, 0, 0, 0, 0, 0 };

void X(codelet_mhc2r_128) (planner *p) {
     X(khc2r_register) (p, mhc2r_128, &desc);
}