m1_64.c

fftw-3.0.1

		    T7W = T3Q - T3P;
		    T7X = KP707106781 * (T43 + T44);
		    T7Y = T7W - T7X;
		    T9o = T7W + T7X;
		    T42 = T40 - T41;
		    T45 = KP707106781 * (T43 - T44);
		    T46 = T42 - T45;
		    T6K = T42 + T45;
	       }
	  }
	  {
	       E T14, T4P, T4d, TaG, T17, T4a, T4S, TaH, T1e, TaZ, T4j, T4V, T1b, TaY, T4o;
	       E T4U;
	       {
		    E T12, T13, T4Q, T4R;
		    T12 = ri[WS(is, 1)];
		    T13 = ri[WS(is, 33)];
		    T14 = T12 + T13;
		    T4P = T12 - T13;
		    {
			 E T4b, T4c, T15, T16;
			 T4b = ii[WS(is, 1)];
			 T4c = ii[WS(is, 33)];
			 T4d = T4b - T4c;
			 TaG = T4b + T4c;
			 T15 = ri[WS(is, 17)];
			 T16 = ri[WS(is, 49)];
			 T17 = T15 + T16;
			 T4a = T15 - T16;
		    }
		    T4Q = ii[WS(is, 17)];
		    T4R = ii[WS(is, 49)];
		    T4S = T4Q - T4R;
		    TaH = T4Q + T4R;
		    {
			 E T1c, T1d, T4f, T4g, T4h, T4i;
			 T1c = ri[WS(is, 57)];
			 T1d = ri[WS(is, 25)];
			 T4f = T1c - T1d;
			 T4g = ii[WS(is, 57)];
			 T4h = ii[WS(is, 25)];
			 T4i = T4g - T4h;
			 T1e = T1c + T1d;
			 TaZ = T4g + T4h;
			 T4j = T4f - T4i;
			 T4V = T4f + T4i;
		    }
		    {
			 E T19, T1a, T4k, T4l, T4m, T4n;
			 T19 = ri[WS(is, 9)];
			 T1a = ri[WS(is, 41)];
			 T4k = T19 - T1a;
			 T4l = ii[WS(is, 9)];
			 T4m = ii[WS(is, 41)];
			 T4n = T4l - T4m;
			 T1b = T19 + T1a;
			 TaY = T4l + T4m;
			 T4o = T4k + T4n;
			 T4U = T4n - T4k;
		    }
	       }
	       {
		    E T18, T1f, TaX, Tb0;
		    T18 = T14 + T17;
		    T1f = T1b + T1e;
		    T1g = T18 + T1f;
		    Tdp = T18 - T1f;
		    TaX = T14 - T17;
		    Tb0 = TaY - TaZ;
		    Tb1 = TaX - Tb0;
		    Tcm = TaX + Tb0;
	       }
	       {
		    E Tdk, Tdl, T4e, T4p;
		    Tdk = TaG + TaH;
		    Tdl = TaY + TaZ;
		    Tdm = Tdk - Tdl;
		    Tej = Tdk + Tdl;
		    T4e = T4a + T4d;
		    T4p = KP707106781 * (T4j - T4o);
		    T4q = T4e - T4p;
		    T6R = T4e + T4p;
	       }
	       {
		    E T4T, T4W, T8d, T8e;
		    T4T = T4P - T4S;
		    T4W = KP707106781 * (T4U - T4V);
		    T4X = T4T - T4W;
		    T6O = T4T + T4W;
		    T8d = T4P + T4S;
		    T8e = KP707106781 * (T4o + T4j);
		    T8f = T8d - T8e;
		    T9s = T8d + T8e;
	       }
	       {
		    E TaI, TaJ, T82, T83;
		    TaI = TaG - TaH;
		    TaJ = T1e - T1b;
		    TaK = TaI - TaJ;
		    Tcp = TaJ + TaI;
		    T82 = T4d - T4a;
		    T83 = KP707106781 * (T4U + T4V);
		    T84 = T82 - T83;
		    T9v = T82 + T83;
	       }
	  }
	  {
	       E T1j, TaR, T1m, TaS, T4G, T4L, TaT, TaQ, T89, T88, T1q, TaM, T1t, TaN, T4v;
	       E T4A, TaO, TaL, T86, T85;
	       {
		    E T4H, T4F, T4C, T4K;
		    {
			 E T1h, T1i, T4D, T4E;
			 T1h = ri[WS(is, 5)];
			 T1i = ri[WS(is, 37)];
			 T1j = T1h + T1i;
			 T4H = T1h - T1i;
			 T4D = ii[WS(is, 5)];
			 T4E = ii[WS(is, 37)];
			 T4F = T4D - T4E;
			 TaR = T4D + T4E;
		    }
		    {
			 E T1k, T1l, T4I, T4J;
			 T1k = ri[WS(is, 21)];
			 T1l = ri[WS(is, 53)];
			 T1m = T1k + T1l;
			 T4C = T1k - T1l;
			 T4I = ii[WS(is, 21)];
			 T4J = ii[WS(is, 53)];
			 T4K = T4I - T4J;
			 TaS = T4I + T4J;
		    }
		    T4G = T4C + T4F;
		    T4L = T4H - T4K;
		    TaT = TaR - TaS;
		    TaQ = T1j - T1m;
		    T89 = T4H + T4K;
		    T88 = T4F - T4C;
	       }
	       {
		    E T4r, T4z, T4w, T4u;
		    {
			 E T1o, T1p, T4x, T4y;
			 T1o = ri[WS(is, 61)];
			 T1p = ri[WS(is, 29)];
			 T1q = T1o + T1p;
			 T4r = T1o - T1p;
			 T4x = ii[WS(is, 61)];
			 T4y = ii[WS(is, 29)];
			 T4z = T4x - T4y;
			 TaM = T4x + T4y;
		    }
		    {
			 E T1r, T1s, T4s, T4t;
			 T1r = ri[WS(is, 13)];
			 T1s = ri[WS(is, 45)];
			 T1t = T1r + T1s;
			 T4w = T1r - T1s;
			 T4s = ii[WS(is, 13)];
			 T4t = ii[WS(is, 45)];
			 T4u = T4s - T4t;
			 TaN = T4s + T4t;
		    }
		    T4v = T4r - T4u;
		    T4A = T4w + T4z;
		    TaO = TaM - TaN;
		    TaL = T1q - T1t;
		    T86 = T4z - T4w;
		    T85 = T4r + T4u;
	       }
	       {
		    E T1n, T1u, Tb2, Tb3;
		    T1n = T1j + T1m;
		    T1u = T1q + T1t;
		    T1v = T1n + T1u;
		    Tdn = T1u - T1n;
		    Tb2 = TaT - TaQ;
		    Tb3 = TaL + TaO;
		    Tb4 = KP707106781 * (Tb2 - Tb3);
		    Tcq = KP707106781 * (Tb2 + Tb3);
	       }
	       {
		    E Tdq, Tdr, T4B, T4M;
		    Tdq = TaR + TaS;
		    Tdr = TaM + TaN;
		    Tds = Tdq - Tdr;
		    Tek = Tdq + Tdr;
		    T4B = FNMS(KP923879532, T4A, KP382683432 * T4v);
		    T4M = FMA(KP923879532, T4G, KP382683432 * T4L);
		    T4N = T4B - T4M;
		    T6P = T4M + T4B;
	       }
	       {
		    E T4Y, T4Z, T8g, T8h;
		    T4Y = FNMS(KP923879532, T4L, KP382683432 * T4G);
		    T4Z = FMA(KP382683432, T4A, KP923879532 * T4v);
		    T50 = T4Y - T4Z;
		    T6S = T4Y + T4Z;
		    T8g = FNMS(KP382683432, T89, KP923879532 * T88);
		    T8h = FMA(KP923879532, T86, KP382683432 * T85);
		    T8i = T8g - T8h;
		    T9w = T8g + T8h;
	       }
	       {
		    E TaP, TaU, T87, T8a;
		    TaP = TaL - TaO;
		    TaU = TaQ + TaT;
		    TaV = KP707106781 * (TaP - TaU);
		    Tcn = KP707106781 * (TaU + TaP);
		    T87 = FNMS(KP382683432, T86, KP923879532 * T85);
		    T8a = FMA(KP382683432, T88, KP923879532 * T89);
		    T8b = T87 - T8a;
		    T9t = T8a + T87;
	       }
	  }
	  {
	       E T1O, Tbc, T1R, Tbd, T5o, T5t, Tbf, Tbe, T8p, T8o, T1V, Tbi, T1Y, Tbj, T5z;
	       E T5E, Tbk, Tbh, T8s, T8r;
	       {
		    E T5p, T5n, T5k, T5s;
		    {
			 E T1M, T1N, T5l, T5m;
			 T1M = ri[WS(is, 3)];
			 T1N = ri[WS(is, 35)];
			 T1O = T1M + T1N;
			 T5p = T1M - T1N;
			 T5l = ii[WS(is, 3)];
			 T5m = ii[WS(is, 35)];
			 T5n = T5l - T5m;
			 Tbc = T5l + T5m;
		    }
		    {
			 E T1P, T1Q, T5q, T5r;
			 T1P = ri[WS(is, 19)];
			 T1Q = ri[WS(is, 51)];
			 T1R = T1P + T1Q;
			 T5k = T1P - T1Q;
			 T5q = ii[WS(is, 19)];
			 T5r = ii[WS(is, 51)];
			 T5s = T5q - T5r;
			 Tbd = T5q + T5r;
		    }
		    T5o = T5k + T5n;
		    T5t = T5p - T5s;
		    Tbf = T1O - T1R;
		    Tbe = Tbc - Tbd;
		    T8p = T5p + T5s;
		    T8o = T5n - T5k;
	       }
	       {
		    E T5A, T5y, T5v, T5D;
		    {
			 E T1T, T1U, T5w, T5x;
			 T1T = ri[WS(is, 59)];
			 T1U = ri[WS(is, 27)];
			 T1V = T1T + T1U;
			 T5A = T1T - T1U;
			 T5w = ii[WS(is, 59)];
			 T5x = ii[WS(is, 27)];
			 T5y = T5w - T5x;
			 Tbi = T5w + T5x;
		    }
		    {
			 E T1W, T1X, T5B, T5C;
			 T1W = ri[WS(is, 11)];
			 T1X = ri[WS(is, 43)];
			 T1Y = T1W + T1X;
			 T5v = T1W - T1X;
			 T5B = ii[WS(is, 11)];
			 T5C = ii[WS(is, 43)];
			 T5D = T5B - T5C;
			 Tbj = T5B + T5C;
		    }
		    T5z = T5v + T5y;
		    T5E = T5A - T5D;
		    Tbk = Tbi - Tbj;
		    Tbh = T1V - T1Y;
		    T8s = T5A + T5D;
		    T8r = T5y - T5v;
	       }
	       {
		    E T1S, T1Z, Tbt, Tbu;
		    T1S = T1O + T1R;
		    T1Z = T1V + T1Y;
		    T20 = T1S + T1Z;
		    TdD = T1Z - T1S;
		    Tbt = Tbh - Tbk;
		    Tbu = Tbf + Tbe;
		    Tbv = KP707106781 * (Tbt - Tbu);
		    Tcu = KP707106781 * (Tbu + Tbt);
	       }
	       {
		    E Tdw, Tdx, T5u, T5F;
		    Tdw = Tbc + Tbd;
		    Tdx = Tbi + Tbj;
		    Tdy = Tdw - Tdx;
		    Tep = Tdw + Tdx;
		    T5u = FNMS(KP923879532, T5t, KP382683432 * T5o);
		    T5F = FMA(KP382683432, T5z, KP923879532 * T5E);
		    T5G = T5u - T5F;
		    T6Z = T5u + T5F;
	       }
	       {
		    E T5R, T5S, T8z, T8A;
		    T5R = FNMS(KP923879532, T5z, KP382683432 * T5E);
		    T5S = FMA(KP923879532, T5o, KP382683432 * T5t);
		    T5T = T5R - T5S;
		    T6W = T5S + T5R;
		    T8z = FNMS(KP382683432, T8r, KP923879532 * T8s);
		    T8A = FMA(KP382683432, T8o, KP923879532 * T8p);
		    T8B = T8z - T8A;
		    T9A = T8A + T8z;
	       }
	       {
		    E Tbg, Tbl, T8q, T8t;
		    Tbg = Tbe - Tbf;
		    Tbl = Tbh + Tbk;
		    Tbm = KP707106781 * (Tbg - Tbl);
		    Tcx = KP707106781 * (Tbg + Tbl);
		    T8q = FNMS(KP382683432, T8p, KP923879532 * T8o);
		    T8t = FMA(KP923879532, T8r, KP382683432 * T8s);
		    T8u = T8q - T8t;
		    T9D = T8q + T8t;
	       }
	  }
	  {
	       E T11, TeD, TeG, TeI, T22, T23, T34, TeH;
	       {
		    E Tv, T10, TeE, TeF;
		    Tv = Tf + Tu;
		    T10 = TK + TZ;
		    T11 = Tv + T10;
		    TeD = Tv - T10;
		    TeE = Tej + Tek;
		    TeF = Teo + Tep;
		    TeG = TeE - TeF;
		    TeI = TeE + TeF;
	       }
	       {
		    E T1w, T21, T2y, T33;
		    T1w = T1g + T1v;
		    T21 = T1L + T20;
		    T22 = T1w + T21;
		    T23 = T21 - T1w;
		    T2y = T2i + T2x;
		    T33 = T2N + T32;
		    T34 = T2y - T33;
		    TeH = T2y + T33;
	       }
	       ro[WS(os, 32)] = T11 - T22;
	       io[WS(os, 32)] = TeH - TeI;
	       ro[0] = T11 + T22;
	       io[0] = TeH + TeI;
	       io[WS(os, 16)] = T23 + T34;
	       ro[WS(os, 16)] = TeD + TeG;
	       io[WS(os, 48)] = T34 - T23;
	       ro[WS(os, 48)] = TeD - TeG;
	  }
	  {
	       E Teh, Tex, Tev, TeB, Tem, Tey, Ter, Tez;
	       {
		    E Tef, Teg, Tet, Teu;
		    Tef = Tf - Tu;
		    Teg = T2N - T32;
		    Teh = Tef + Teg;
		    Tex = Tef - Teg;
		    Tet = T2i - T2x;
		    Teu = TZ - TK;
		    Tev = Tet - Teu;
		    TeB = Teu + Tet;
	       }
	       {
		    E Tei, Tel, Ten, Teq;
		    Tei = T1g - T1v;
		    Tel = Tej - Tek;
		    Tem = Tei + Tel;
		    Tey = Tel - Tei;
		    Ten = T1L - T20;
		    Teq = Teo - Tep;
		    Ter = Ten - Teq;
		    Tez = Ten + Teq;
	       }
	       {
		    E Tes, TeC, Tew, TeA;
		    Tes = KP707106781 * (Tem + Ter);
		    ro[WS(os, 40)] = Teh - Tes;
		    ro[WS(os, 8)] = Teh + Tes;
		    TeC = KP707106781 * (Tey + Tez);
		    io[WS(os, 40)] = TeB - TeC;
		    io[WS(os, 8)] = TeB + TeC;
		    Tew = KP707106781 * (Ter - Tem);
		    io[WS(os, 56)] = Tev - Tew;
		    io[WS(os, 24)] = Tev + Tew;
		    TeA = KP707106781 * (Tey - Tez);
		    ro[WS(os, 56)] = Tex - TeA;
		    ro[WS(os, 24)] = Tex + TeA;
	       }
	  }
	  {
	       E Tdb, TdV, Te5, TdJ, Tdi, Te6, Te3, Teb, TdM, TdW, Tdu, TdQ, Te0, Tea, TdF;
	       E TdR;
	       {
		    E Tde, Tdh, Tdo, Tdt;
		    Tdb = Td9 - Tda;
		    TdV = Td9 + Tda;
		    Te5 = TdI + TdH;
		    TdJ = TdH - TdI;
		    Tde = Tdc - Tdd;
		    Tdh = Tdf + Tdg;
		    Tdi = KP707106781 * (Tde - Tdh);
		    Te6 = KP707106781 * (Tde + Tdh);
		    {
			 E Te1, Te2, TdK, TdL;
			 Te1 = Tdv + Tdy;
			 Te2 = TdD + TdC;
			 Te3 = FNMS(KP382683432, Te2, KP923879532 * Te1);
			 Teb = FMA(KP923879532, Te2, KP382683432 * Te1);
			 TdK = Tdf - Tdg;
			 TdL = Tdd + Tdc;
			 TdM = KP707106781 * (TdK - TdL);
			 TdW = KP707106781 * (TdL + TdK);
		    }
		    Tdo = Tdm - Tdn;
		    Tdt = Tdp - Tds;
		    Tdu = FMA(KP923879532, Tdo, KP382683432 * Tdt);
		    TdQ = FNMS(KP923879532, Tdt, KP382683432 * Tdo);
		    {
			 E TdY, TdZ, Tdz, TdE;
			 TdY = Tdn + Tdm;
			 TdZ = Tdp + Tds;
			 Te0 = FMA(KP382683432, TdY, KP923879532 * TdZ);
			 Tea = FNMS(KP382683432, TdZ, KP923879532 * TdY);
			 Tdz = Tdv - Tdy;
			 TdE = TdC - TdD;
			 TdF = FNMS(KP923879532, TdE, KP382683432 * Tdz);
			 TdR = FMA(KP382683432, TdE, KP923879532 * Tdz);
		    }
	       }
	       {
		    E Tdj, TdG, TdT, TdU;
		    Tdj = Tdb + Tdi;
		    TdG = Tdu + TdF;
		    ro[WS(os, 44)] = Tdj - TdG;
		    ro[WS(os, 12)] = Tdj + TdG;
		    TdT = TdJ + TdM;
		    TdU = TdQ + TdR;
		    io[WS(os, 44)] = TdT - TdU;
		    io[WS(os, 12)] = TdT + TdU;
	       }
	       {
		    E TdN, TdO, TdP, TdS;
		    TdN = TdJ - TdM;
		    TdO = TdF - Tdu;
		    io[WS(os, 60)] = TdN - TdO;
		    io[WS(os, 28)] = TdN + TdO;
		    TdP = Tdb - Tdi;
		    TdS = TdQ - TdR;
		    ro[WS(os, 60)] = TdP - TdS;
		    ro[WS(os, 28)] = TdP + TdS;
	       }
	       {
		    E TdX, Te4, Ted, Tee;
		    TdX = TdV + TdW;
		    Te4 = Te0 + Te3;
		    ro[WS(os, 36)] = TdX - Te4;
		    ro[WS(os, 4)] = TdX + Te4;
		    Ted = Te5 + Te6;
		    Tee = Tea + Teb;
		    io[WS(os, 36)] = Ted - Tee;
		    io[WS(os, 4)] = Ted + Tee;
	       }
	       {
		    E Te7, Te8, Te9, Tec;
		    Te7 = Te5 - Te6;
		    Te8 = Te3 - Te0;
		    io[WS(os, 52)] = Te7 - Te8;
		    io[WS(os, 20)] = Te7 + Te8;
		    Te9 = TdV - TdW;
		    Tec = Tea - Teb;
		    ro[WS(os, 52)] = Te9 - Tec;
		    ro[WS(os, 20)] = Te9 + Tec;
	       }
	  }
	  {
	       E Tcd, TcP, TcD, TcZ, Tck, Td0, TcX, Td5, Tcs, TcK, TcG, TcQ, TcU, Td4, Tcz;
	       E TcL, Tcc, TcC;
	       Tcc = KP707106781 * (TbD + TbC);
	       Tcd = Tcb - Tcc;
	       TcP = Tcb + Tcc;
	       TcC = KP707106781 * (Tak + Tan);
	       TcD = TcB - TcC;
	       TcZ = TcB + TcC;
	       {
		    E Tcg, Tcj, TcV, TcW;
		    Tcg = FNMS(KP382683432, Tcf, KP923879532 * Tce);
		    Tcj = FMA(KP923879532, Tch, KP382683432 * Tci);
		    Tck = Tcg - Tcj;
		    Td0 = Tcg + Tcj;
		    TcV = Tct + Tcu;
		    TcW = Tcw + Tcx;
		    TcX = FNMS(KP195090322, TcW, KP980785280 * TcV);
		    Td5 = FMA(KP195090322, TcV, KP980785280 * TcW);
	       }
	       {
		    E Tco, Tcr, TcE, TcF;
		    Tco = Tcm - Tcn;
		    Tcr = Tcp - Tcq;
		    Tcs = FMA(KP555570233, Tco, KP831469612 * Tcr);
		    TcK = FNMS(KP831469612, Tco, KP555570233 * Tcr);
		    TcE = FNMS(KP382683432, Tch, KP923879532 * Tci);
		    TcF = FMA(KP382683432, Tce, KP923879532 * Tcf);
		    TcG = TcE - TcF;
		    TcQ = TcF + TcE;
	       }
	       {
		    E TcS, TcT, Tcv, Tcy;
		    TcS = Tcm + Tcn;
		    TcT = Tcp + Tcq;
		    TcU = FMA(KP980785280, TcS, KP195090322 * TcT);
		    Td4 = FNMS(KP195090322, TcS, KP980785280 * TcT);
		    Tcv = Tct - Tcu;
		    Tcy = Tcw - Tcx;