mr2hc_128.c

fftw-3.0.1

		    T4n = T4j - T4m;
		    T4q = T4o - T4p;
		    T4r = T4n + T4q;
		    T6F = T4q - T4n;
	       }
	       {
		    E T7S, T7T, Tak, Tan;
		    T7S = T3k - T3r;
		    T7T = T4p + T4o;
		    T7U = T7S + T7T;
		    T9l = T7S - T7T;
		    Tak = Tai + Taj;
		    Tan = Tal + Tam;
		    Tao = Tak - Tan;
		    TbW = Tak + Tan;
	       }
	       {
		    E TaH, TaK, T83, T84;
		    TaH = TaF + TaG;
		    TaK = TaI + TaJ;
		    TaL = TaH - TaK;
		    TbX = TaK + TaH;
		    T83 = T4m + T4j;
		    T84 = T3G - T3z;
		    T85 = T83 + T84;
		    T9j = T84 - T83;
	       }
	  }
	  {
	       E T1z, T2V, T1C, T39, T1G, T38, T1J, T2Y, T1O, T1R, T32, Td0, T3c, T1V, T1Y;
	       E T35, Td1, T3b;
	       {
		    E T1x, T1y, T1A, T1B;
		    T1x = I[WS(is, 126)];
		    T1y = I[WS(is, 62)];
		    T1z = T1x + T1y;
		    T2V = T1x - T1y;
		    T1A = I[WS(is, 30)];
		    T1B = I[WS(is, 94)];
		    T1C = T1A + T1B;
		    T39 = T1A - T1B;
	       }
	       {
		    E T1E, T1F, T2W, T1H, T1I, T2X;
		    T1E = I[WS(is, 14)];
		    T1F = I[WS(is, 78)];
		    T2W = T1E - T1F;
		    T1H = I[WS(is, 110)];
		    T1I = I[WS(is, 46)];
		    T2X = T1H - T1I;
		    T1G = T1E + T1F;
		    T38 = KP707106781 * (T2X - T2W);
		    T1J = T1H + T1I;
		    T2Y = KP707106781 * (T2W + T2X);
	       }
	       {
		    E T30, T31, T33, T34;
		    {
			 E T1M, T1N, T1P, T1Q;
			 T1M = I[WS(is, 6)];
			 T1N = I[WS(is, 70)];
			 T1O = T1M + T1N;
			 T30 = T1M - T1N;
			 T1P = I[WS(is, 38)];
			 T1Q = I[WS(is, 102)];
			 T1R = T1P + T1Q;
			 T31 = T1P - T1Q;
		    }
		    T32 = FNMS(KP382683432, T31, KP923879532 * T30);
		    Td0 = T1O - T1R;
		    T3c = FMA(KP382683432, T30, KP923879532 * T31);
		    {
			 E T1T, T1U, T1W, T1X;
			 T1T = I[WS(is, 118)];
			 T1U = I[WS(is, 54)];
			 T1V = T1T + T1U;
			 T33 = T1T - T1U;
			 T1W = I[WS(is, 22)];
			 T1X = I[WS(is, 86)];
			 T1Y = T1W + T1X;
			 T34 = T1W - T1X;
		    }
		    T35 = FMA(KP923879532, T33, KP382683432 * T34);
		    Td1 = T1V - T1Y;
		    T3b = FNMS(KP923879532, T34, KP382683432 * T33);
	       }
	       {
		    E T1D, T1K, TcZ, Td2;
		    T1D = T1z + T1C;
		    T1K = T1G + T1J;
		    T1L = T1D + T1K;
		    Tad = T1D - T1K;
		    TcZ = T1z - T1C;
		    Td2 = KP707106781 * (Td0 + Td1);
		    Td3 = TcZ + Td2;
		    Tew = TcZ - Td2;
	       }
	       {
		    E Td4, Td5, T1S, T1Z;
		    Td4 = T1J - T1G;
		    Td5 = KP707106781 * (Td1 - Td0);
		    Td6 = Td4 + Td5;
		    Tex = Td5 - Td4;
		    T1S = T1O + T1R;
		    T1Z = T1V + T1Y;
		    T20 = T1S + T1Z;
		    Tae = T1Z - T1S;
	       }
	       {
		    E T2Z, T36, T7J, T7K;
		    T2Z = T2V + T2Y;
		    T36 = T32 + T35;
		    T37 = T2Z + T36;
		    T6x = T2Z - T36;
		    T7J = T2V - T2Y;
		    T7K = T3c + T3b;
		    T7L = T7J + T7K;
		    T9a = T7J - T7K;
	       }
	       {
		    E T7M, T7N, T3a, T3d;
		    T7M = T39 + T38;
		    T7N = T35 - T32;
		    T7O = T7M + T7N;
		    T9b = T7N - T7M;
		    T3a = T38 - T39;
		    T3d = T3b - T3c;
		    T3e = T3a + T3d;
		    T6w = T3d - T3a;
	       }
	  }
	  {
	       E T3L, Tdf, T3X, Tar, T42, Tdi, T4e, Tay, T3S, Tdg, T3U, Tau, T49, Tdj, T4b;
	       E TaB, Tdh, Tdk;
	       {
		    E T3J, T3K, Tap, T3V, T3W, Taq;
		    T3J = I[WS(is, 5)];
		    T3K = I[WS(is, 69)];
		    Tap = T3J + T3K;
		    T3V = I[WS(is, 37)];
		    T3W = I[WS(is, 101)];
		    Taq = T3V + T3W;
		    T3L = T3J - T3K;
		    Tdf = Tap - Taq;
		    T3X = T3V - T3W;
		    Tar = Tap + Taq;
	       }
	       {
		    E T40, T41, Taw, T4c, T4d, Tax;
		    T40 = I[WS(is, 125)];
		    T41 = I[WS(is, 61)];
		    Taw = T40 + T41;
		    T4c = I[WS(is, 29)];
		    T4d = I[WS(is, 93)];
		    Tax = T4c + T4d;
		    T42 = T40 - T41;
		    Tdi = Taw - Tax;
		    T4e = T4c - T4d;
		    Tay = Taw + Tax;
	       }
	       {
		    E T3O, Tas, T3R, Tat;
		    {
			 E T3M, T3N, T3P, T3Q;
			 T3M = I[WS(is, 21)];
			 T3N = I[WS(is, 85)];
			 T3O = T3M - T3N;
			 Tas = T3M + T3N;
			 T3P = I[WS(is, 117)];
			 T3Q = I[WS(is, 53)];
			 T3R = T3P - T3Q;
			 Tat = T3P + T3Q;
		    }
		    T3S = KP707106781 * (T3O + T3R);
		    Tdg = Tat - Tas;
		    T3U = KP707106781 * (T3R - T3O);
		    Tau = Tas + Tat;
	       }
	       {
		    E T45, Taz, T48, TaA;
		    {
			 E T43, T44, T46, T47;
			 T43 = I[WS(is, 13)];
			 T44 = I[WS(is, 77)];
			 T45 = T43 - T44;
			 Taz = T43 + T44;
			 T46 = I[WS(is, 109)];
			 T47 = I[WS(is, 45)];
			 T48 = T46 - T47;
			 TaA = T46 + T47;
		    }
		    T49 = KP707106781 * (T45 + T48);
		    Tdj = TaA - Taz;
		    T4b = KP707106781 * (T48 - T45);
		    TaB = Taz + TaA;
	       }
	       TbZ = Tar + Tau;
	       Tc0 = Tay + TaB;
	       {
		    E T3T, T3Y, Tdq, Tdr;
		    T3T = T3L + T3S;
		    T3Y = T3U - T3X;
		    T3Z = FMA(KP980785280, T3T, KP195090322 * T3Y);
		    T4s = FNMS(KP195090322, T3T, KP980785280 * T3Y);
		    Tdq = FNMS(KP382683432, Tdf, KP923879532 * Tdg);
		    Tdr = FMA(KP382683432, Tdi, KP923879532 * Tdj);
		    Tds = Tdq + Tdr;
		    TeI = Tdr - Tdq;
	       }
	       {
		    E T4a, T4f, T7Y, T7Z;
		    T4a = T42 + T49;
		    T4f = T4b - T4e;
		    T4g = FNMS(KP195090322, T4f, KP980785280 * T4a);
		    T4t = FMA(KP195090322, T4a, KP980785280 * T4f);
		    T7Y = T42 - T49;
		    T7Z = T4e + T4b;
		    T80 = FNMS(KP555570233, T7Z, KP831469612 * T7Y);
		    T87 = FMA(KP555570233, T7Y, KP831469612 * T7Z);
	       }
	       Tdh = FMA(KP923879532, Tdf, KP382683432 * Tdg);
	       Tdk = FNMS(KP382683432, Tdj, KP923879532 * Tdi);
	       Tdl = Tdh + Tdk;
	       TeE = Tdk - Tdh;
	       {
		    E T7V, T7W, Tav, TaC;
		    T7V = T3L - T3S;
		    T7W = T3X + T3U;
		    T7X = FMA(KP831469612, T7V, KP555570233 * T7W);
		    T86 = FNMS(KP555570233, T7V, KP831469612 * T7W);
		    Tav = Tar - Tau;
		    TaC = Tay - TaB;
		    TaD = KP707106781 * (Tav + TaC);
		    TaM = KP707106781 * (TaC - Tav);
	       }
	  }
	  {
	       E T50, TdA, T5c, TaY, T5h, TdD, T5t, Tb5, T57, TdB, T59, Tb1, T5o, TdE, T5q;
	       E Tb8, TdC, TdF;
	       {
		    E T4Y, T4Z, TaW, T5a, T5b, TaX;
		    T4Y = I[WS(is, 3)];
		    T4Z = I[WS(is, 67)];
		    TaW = T4Y + T4Z;
		    T5a = I[WS(is, 35)];
		    T5b = I[WS(is, 99)];
		    TaX = T5a + T5b;
		    T50 = T4Y - T4Z;
		    TdA = TaW - TaX;
		    T5c = T5a - T5b;
		    TaY = TaW + TaX;
	       }
	       {
		    E T5f, T5g, Tb3, T5r, T5s, Tb4;
		    T5f = I[WS(is, 123)];
		    T5g = I[WS(is, 59)];
		    Tb3 = T5f + T5g;
		    T5r = I[WS(is, 27)];
		    T5s = I[WS(is, 91)];
		    Tb4 = T5r + T5s;
		    T5h = T5f - T5g;
		    TdD = Tb3 - Tb4;
		    T5t = T5r - T5s;
		    Tb5 = Tb3 + Tb4;
	       }
	       {
		    E T53, TaZ, T56, Tb0;
		    {
			 E T51, T52, T54, T55;
			 T51 = I[WS(is, 19)];
			 T52 = I[WS(is, 83)];
			 T53 = T51 - T52;
			 TaZ = T51 + T52;
			 T54 = I[WS(is, 115)];
			 T55 = I[WS(is, 51)];
			 T56 = T54 - T55;
			 Tb0 = T54 + T55;
		    }
		    T57 = KP707106781 * (T53 + T56);
		    TdB = Tb0 - TaZ;
		    T59 = KP707106781 * (T56 - T53);
		    Tb1 = TaZ + Tb0;
	       }
	       {
		    E T5k, Tb6, T5n, Tb7;
		    {
			 E T5i, T5j, T5l, T5m;
			 T5i = I[WS(is, 11)];
			 T5j = I[WS(is, 75)];
			 T5k = T5i - T5j;
			 Tb6 = T5i + T5j;
			 T5l = I[WS(is, 107)];
			 T5m = I[WS(is, 43)];
			 T5n = T5l - T5m;
			 Tb7 = T5l + T5m;
		    }
		    T5o = KP707106781 * (T5k + T5n);
		    TdE = Tb7 - Tb6;
		    T5q = KP707106781 * (T5n - T5k);
		    Tb8 = Tb6 + Tb7;
	       }
	       Tc6 = TaY + Tb1;
	       Tc7 = Tb5 + Tb8;
	       {
		    E T58, T5d, TdL, TdM;
		    T58 = T50 + T57;
		    T5d = T59 - T5c;
		    T5e = FMA(KP980785280, T58, KP195090322 * T5d);
		    T5H = FNMS(KP195090322, T58, KP980785280 * T5d);
		    TdL = FNMS(KP382683432, TdA, KP923879532 * TdB);
		    TdM = FMA(KP382683432, TdD, KP923879532 * TdE);
		    TdN = TdL + TdM;
		    TeM = TdM - TdL;
	       }
	       {
		    E T5p, T5u, T8h, T8i;
		    T5p = T5h + T5o;
		    T5u = T5q - T5t;
		    T5v = FNMS(KP195090322, T5u, KP980785280 * T5p);
		    T5I = FMA(KP195090322, T5p, KP980785280 * T5u);
		    T8h = T5h - T5o;
		    T8i = T5t + T5q;
		    T8j = FNMS(KP555570233, T8i, KP831469612 * T8h);
		    T8q = FMA(KP555570233, T8h, KP831469612 * T8i);
	       }
	       TdC = FMA(KP923879532, TdA, KP382683432 * TdB);
	       TdF = FNMS(KP382683432, TdE, KP923879532 * TdD);
	       TdG = TdC + TdF;
	       TeO = TdF - TdC;
	       {
		    E T8e, T8f, Tb2, Tb9;
		    T8e = T50 - T57;
		    T8f = T5c + T59;
		    T8g = FMA(KP831469612, T8e, KP555570233 * T8f);
		    T8p = FNMS(KP555570233, T8e, KP831469612 * T8f);
		    Tb2 = TaY - Tb1;
		    Tb9 = Tb5 - Tb8;
		    Tba = KP707106781 * (Tb2 + Tb9);
		    Tbj = KP707106781 * (Tb9 - Tb2);
	       }
	  }
	  {
	       E T11, TbV, Tc9, Tcf, T22, Tcb, Tc2, Tce;
	       {
		    E Tv, T10, Tc5, Tc8;
		    Tv = Tf + Tu;
		    T10 = TK + TZ;
		    T11 = Tv + T10;
		    TbV = Tv - T10;
		    Tc5 = Tc3 + Tc4;
		    Tc8 = Tc6 + Tc7;
		    Tc9 = Tc5 - Tc8;
		    Tcf = Tc5 + Tc8;
	       }
	       {
		    E T1w, T21, TbY, Tc1;
		    T1w = T1g + T1v;
		    T21 = T1L + T20;
		    T22 = T1w + T21;
		    Tcb = T21 - T1w;
		    TbY = TbW + TbX;
		    Tc1 = TbZ + Tc0;
		    Tc2 = TbY - Tc1;
		    Tce = TbY + Tc1;
	       }
	       ro[WS(ros, 32)] = T11 - T22;
	       io[WS(ios, 32)] = Tcf - Tce;
	       {
		    E Tca, Tcc, Tcd, Tcg;
		    Tca = KP707106781 * (Tc2 + Tc9);
		    ro[WS(ros, 48)] = TbV - Tca;
		    ro[WS(ros, 16)] = TbV + Tca;
		    Tcc = KP707106781 * (Tc9 - Tc2);
		    io[WS(ios, 16)] = Tcb + Tcc;
		    io[WS(ios, 48)] = Tcc - Tcb;
		    Tcd = T11 + T22;
		    Tcg = Tce + Tcf;
		    ro[WS(ros, 64)] = Tcd - Tcg;
		    ro[0] = Tcd + Tcg;
	       }
	  }
	  {
	       E Tch, Tcu, Tck, Tct, Tco, Tcy, Tcr, Tcz, Tci, Tcj;
	       Tch = Tf - Tu;
	       Tcu = TZ - TK;
	       Tci = T1g - T1v;
	       Tcj = T1L - T20;
	       Tck = KP707106781 * (Tci + Tcj);
	       Tct = KP707106781 * (Tcj - Tci);
	       {
		    E Tcm, Tcn, Tcp, Tcq;
		    Tcm = TbW - TbX;
		    Tcn = Tc0 - TbZ;
		    Tco = FMA(KP923879532, Tcm, KP382683432 * Tcn);
		    Tcy = FNMS(KP382683432, Tcm, KP923879532 * Tcn);
		    Tcp = Tc3 - Tc4;
		    Tcq = Tc7 - Tc6;
		    Tcr = FNMS(KP382683432, Tcq, KP923879532 * Tcp);
		    Tcz = FMA(KP382683432, Tcp, KP923879532 * Tcq);
	       }
	       {
		    E Tcl, Tcs, Tcx, TcA;
		    Tcl = Tch + Tck;
		    Tcs = Tco + Tcr;
		    ro[WS(ros, 56)] = Tcl - Tcs;
		    ro[WS(ros, 8)] = Tcl + Tcs;
		    Tcx = Tcu + Tct;
		    TcA = Tcy + Tcz;
		    io[WS(ios, 8)] = Tcx + TcA;
		    io[WS(ios, 56)] = TcA - Tcx;
	       }
	       {
		    E Tcv, Tcw, TcB, TcC;
		    Tcv = Tct - Tcu;
		    Tcw = Tcr - Tco;
		    io[WS(ios, 24)] = Tcv + Tcw;
		    io[WS(ios, 40)] = Tcw - Tcv;
		    TcB = Tch - Tck;
		    TcC = Tcz - Tcy;
		    ro[WS(ros, 40)] = TcB - TcC;
		    ro[WS(ros, 24)] = TcB + TcC;
	       }
	  }
	  {
	       E Ta9, TbB, Tbs, TbM, Tag, TbL, TbJ, TbR, TaO, Tbw, Tbp, TbC, TbG, TbQ, Tbl;
	       E Tbx, Ta8, Tbr;
	       Ta8 = KP707106781 * (Ta6 + Ta7);
	       Ta9 = Ta5 + Ta8;
	       TbB = Ta5 - Ta8;
	       Tbr = KP707106781 * (Ta7 - Ta6);
	       Tbs = Tbq + Tbr;
	       TbM = Tbr - Tbq;
	       {
		    E Tac, Taf, TbH, TbI;
		    Tac = FMA(KP923879532, Taa, KP382683432 * Tab);
		    Taf = FNMS(KP382683432, Tae, KP923879532 * Tad);
		    Tag = Tac + Taf;
		    TbL = Taf - Tac;
		    TbH = TaV - Tba;
		    TbI = Tbj - Tbi;
		    TbJ = FNMS(KP555570233, TbI, KP831469612 * TbH);
		    TbR = FMA(KP555570233, TbH, KP831469612 * TbI);
	       }
	       {
		    E TaE, TaN, Tbn, Tbo;
		    TaE = Tao + TaD;
		    TaN = TaL + TaM;
		    TaO = FMA(KP980785280, TaE, KP195090322 * TaN);
		    Tbw = FNMS(KP195090322, TaE, KP980785280 * TaN);
		    Tbn = FNMS(KP382683432, Taa, KP923879532 * Tab);
		    Tbo = FMA(KP382683432, Tad, KP923879532 * Tae);
		    Tbp = Tbn + Tbo;
		    TbC = Tbo - Tbn;
	       }
	       {
		    E TbE, TbF, Tbb, Tbk;
		    TbE = Tao - TaD;
		    TbF = TaM - TaL;
		    TbG = FMA(KP831469612, TbE, KP555570233 * TbF);
		    TbQ = FNMS(KP555570233, TbE, KP831469612 * TbF);
		    Tbb = TaV + Tba;
		    Tbk = Tbi + Tbj;
		    Tbl = FNMS(KP195090322, Tbk, KP980785280 * Tbb);
		    Tbx = FMA(KP195090322, Tbb, KP980785280 * Tbk);
	       }
	       {
		    E Tah, Tbm, Tbv, Tby;
		    Tah = Ta9 + Tag;
		    Tbm = TaO + Tbl;
		    ro[WS(ros, 60)] = Tah - Tbm;
		    ro[WS(ros, 4)] = Tah + Tbm;
		    Tbv = Tbs + Tbp;
		    Tby = Tbw + Tbx;
		    io[WS(ios, 4)] = Tbv + Tby;
		    io[WS(ios, 60)] = Tby - Tbv;
	       }
	       {
		    E Tbt, Tbu, Tbz, TbA;
		    Tbt = Tbp - Tbs;
		    Tbu = Tbl - TaO;
		    io[WS(ios, 28)] = Tbt + Tbu;
		    io[WS(ios, 36)] = Tbu - Tbt;
		    Tbz = Ta9 - Tag;
		    TbA = Tbx - Tbw;
		    ro[WS(ros, 36)] = Tbz - TbA;
		    ro[WS(ros, 28)] = Tbz + TbA;
	       }
	       {
		    E TbD, TbK, TbP, TbS;
		    TbD = TbB + TbC;
		    TbK = TbG + TbJ;
		    ro[WS(ros, 52)] = TbD - TbK;
		    ro[WS(ros, 12)] = TbD + TbK;
		    TbP = TbM + TbL;
		    TbS = TbQ + TbR;
		    io[WS(ios, 12)] = TbP + TbS;
		    io[WS(ios, 52)] = TbS - TbP;
	       }
	       {
		    E TbN, TbO, TbT, TbU;
		    TbN = TbL - TbM;
		    TbO = TbJ - TbG;
		    io[WS(ios, 20)] = TbN + TbO;
		    io[WS(ios, 44)] = TbO - TbN;
		    TbT = TbB - TbC;
		    TbU = TbR - TbQ;
		    ro[WS(ros, 44)] = TbT - TbU;
		    ro[WS(ros, 20)] = TbT + TbU;
	       }
	  }
	  {
	       E Tev, Tf7, Tfc, Tfm, Tff, Tfn, TeC, Tfh, TeK, Tf2, TeV, Tf8, TeY, Tfi, TeR;
	       E Tf3;
	       {
		    E Tet, Teu, Tfa, Tfb;
		    Tet = TcD - TcG;
		    Teu = TdY - TdX;
		    Tev = Tet - Teu;
		    Tf7 = Tet + Teu;
		    Tfa = TeF + TeE;
		    Tfb = TeH + TeI;
		    Tfc = FMA(KP290284677, Tfa, KP956940335 * Tfb);
		    Tfm = FNMS(KP290284677, Tfb, KP956940335 * Tfa);
	       }
	       {
		    E Tfd, Tfe, Tey, TeB;
		    Tfd = TeL + TeM;
		    Tfe = TeP + TeO;
		    Tff = FNMS(KP290284677, Tfe, KP956940335 * Tfd);
		    Tfn = FMA(KP956940335, Tfe, KP290284677 * Tfd);
		    Tey = FMA(KP555570233, Tew, KP831469612 * Tex);
		    TeB = FNMS(KP555570233, TeA, KP831469612 * Tez);
		    TeC = Tey - TeB;
		    Tfh = TeB + Tey;
	       }
	       {
		    E TeG, TeJ, TeT, TeU;
		    TeG = TeE - TeF;
		    TeJ = TeH - TeI;
		    TeK = FMA(KP471396736, TeG, KP881921264 * TeJ);
		    Tf2 = FNMS(KP471396736, TeJ, KP881921264 * TeG);
		    TeT = FNMS(KP555570233, Tex, KP831469612 * Tew);
		    TeU = FMA(KP831469612, TeA, KP555570233 * Tez);
		    TeV = TeT - TeU;
		    Tf8 = TeU + TeT;
	       }
	       {
		    E TeW, TeX, TeN, TeQ;
		    TeW = TcN - TcK;
		    TeX = TdV - TdU;