t2_64.c

最新的FFT程序

						  {
						       E Ta9, Ta2, T4r, T4y;
						       Ta9 = T4E * T4I;
						       Ta2 = T4s * T4q;
						       T4r = T4p * T4q;
						       T4y = ri[WS(ios, 54)];
						       Taa = FNMS(T4H, T4F, Ta9);
						       Ta3 = FMA(T4p, T4t, Ta2);
						       T4u = FNMS(T4s, T4t, T4r);
						       T4z = T4x * T4y;
						       Ta7 = T4B * T4y;
						  }
					     }
					     {
						  E Ta4, TgC, T4v, T9Z, T4D, Ta8, T4K, TgD;
						  Ta4 = Ta1 - Ta3;
						  TgC = Ta1 + Ta3;
						  T4v = T4o + T4u;
						  T9Z = T4o - T4u;
						  T4D = FNMS(T4B, T4C, T4z);
						  Ta8 = FMA(T4x, T4C, Ta7);
						  Ta6 = T4D - T4J;
						  T4K = T4D + T4J;
						  TgD = Ta8 + Taa;
						  Tab = Ta8 - Taa;
						  Tan = Ta4 - T9Z;
						  Ta5 = T9Z + Ta4;
						  Tgz = T4K - T4v;
						  T4L = T4v + T4K;
						  TiU = TgC + TgD;
						  TgE = TgC - TgD;
					     }
					}
				   }
				   {
					E T5h, TaC, T5t, T5G, T5u, TaL, TaE, T5p, T5x;
					{
					     E T5B, T5C, T5F, T5e, T5g;
					     T5e = ri[WS(ios, 9)];
					     T5g = ii[WS(ios, 9)];
					     T5B = ri[WS(ios, 25)];
					     {
						  E Tam, Tac, T5f, TaB;
						  Tam = Ta6 + Tab;
						  Tac = Ta6 - Tab;
						  T5f = T3 * T5e;
						  TaB = T3 * T5g;
						  Ted = Tan + Tam;
						  Tao = Tam - Tan;
						  Teg = Ta5 + Tac;
						  Tad = Ta5 - Tac;
						  T5h = FMA(Tb, T5g, T5f);
						  TaC = FNMS(Tb, T5e, TaB);
						  T5C = T5A * T5B;
					     }
					     T5F = ii[WS(ios, 25)];
					     {
						  E T5o, T5k, TaK, TaD, T5l;
						  T5o = ii[WS(ios, 41)];
						  T5k = ri[WS(ios, 41)];
						  T5t = ri[WS(ios, 57)];
						  T5G = FMA(T5E, T5F, T5C);
						  TaK = T5A * T5F;
						  TaD = T5n * T5k;
						  T5l = T5j * T5k;
						  T5u = T5s * T5t;
						  TaL = FNMS(T5E, T5B, TaK);
						  TaE = FMA(T5j, T5o, TaD);
						  T5p = FNMS(T5n, T5o, T5l);
						  T5x = ii[WS(ios, 57)];
					     }
					}
					{
					     E T78, TbL, T7e, TbN, T7j, T7m;
					     {
						  E T7s, T7p, T7a, T7d;
						  {
						       E T75, Th1, TaG, Tbu, TaH, T76, T77, TaM, Th2, TbK, T5q, T5H, TaJ, Tbt, TaN;
						       T75 = ri[WS(ios, 7)];
						       {
							    E TaF, TaA, T5y, TaI;
							    Th1 = TaC + TaE;
							    TaF = TaC - TaE;
							    T5q = T5h + T5p;
							    TaA = T5h - T5p;
							    T5y = FMA(T5w, T5x, T5u);
							    TaI = T5s * T5x;
							    TaG = TaA + TaF;
							    Tbu = TaF - TaA;
							    T5H = T5y + T5G;
							    TaH = T5y - T5G;
							    TaJ = FNMS(T5w, T5t, TaI);
							    T76 = T1i * T75;
						       }
						       T77 = ii[WS(ios, 7)];
						       TgM = T5H - T5q;
						       T5I = T5q + T5H;
						       TaM = TaJ - TaL;
						       Th2 = TaJ + TaL;
						       TbK = T1i * T77;
						       T78 = FMA(T1k, T77, T76);
						       Tbt = TaH + TaM;
						       TaN = TaH - TaM;
						       Tj0 = Th1 + Th2;
						       Th3 = Th1 - Th2;
						       Tem = Tbu + Tbt;
						       Tbv = Tbt - Tbu;
						       Tex = TaG + TaN;
						       TaO = TaG - TaN;
						       TbL = FNMS(T1k, T75, TbK);
						  }
						  T7s = ii[WS(ios, 23)];
						  T7p = ri[WS(ios, 23)];
						  T7a = ri[WS(ios, 39)];
						  T7d = ii[WS(ios, 39)];
						  T7i = ri[WS(ios, 55)];
						  {
						       E TbT, T7q, T7b, TbM;
						       TbT = T7r * T7p;
						       T7q = T7o * T7p;
						       T7b = T79 * T7a;
						       TbM = T79 * T7d;
						       TbU = FMA(T7o, T7s, TbT);
						       T7t = FNMS(T7r, T7s, T7q);
						       T7e = FMA(T7c, T7d, T7b);
						       TbN = FNMS(T7c, T7a, TbM);
						       T7j = T7h * T7i;
						  }
						  T7m = ii[WS(ios, 55)];
					     }
					     T7f = T78 + T7e;
					     TbJ = T78 - T7e;
					     TbO = TbL - TbN;
					     Ths = TbL + TbN;
					     T7n = FMA(T7l, T7m, T7j);
					     TbR = T7h * T7m;
					}
				   }
			      }
			      {
				   E T6g, TaY, T6o, T6z, T6p, TaU, Tb0, T6m, T6q;
				   {
					E T5V, Tb9, T66, T5U, Tb5, Tbg, TgO, T5Y, Tb6, T6u, T6v, T6y;
					{
					     E Tbd, T5N, T5T, Tbf, T5W, T5X;
					     {
						  E T65, T61, T5P, T5S;
						  {
						       E T5M, TbQ, TcD, TbP, TbV, Tht, T5K, T7u, TbS;
						       T5M = ii[WS(ios, 5)];
						       T7u = T7n + T7t;
						       TbQ = T7n - T7t;
						       TbS = FNMS(T7l, T7i, TbR);
						       TcD = TbO - TbJ;
						       TbP = TbJ + TbO;
						       Thd = T7u - T7f;
						       T7v = T7f + T7u;
						       TbV = TbS - TbU;
						       Tht = TbS + TbU;
						       T5K = ri[WS(ios, 5)];
						       T65 = ii[WS(ios, 53)];
						       {
							    E TcC, TbW, Tbc, T5L;
							    TcC = TbQ + TbV;
							    TbW = TbQ - TbV;
							    Tjb = Ths + Tht;
							    Thu = Ths - Tht;
							    Tbc = Ti * T5K;
							    T5L = Te * T5K;
							    TeF = TcD + TcC;
							    TcE = TcC - TcD;
							    TeQ = TbP + TbW;
							    TbX = TbP - TbW;
							    Tbd = FMA(Te, T5M, Tbc);
							    T5N = FNMS(Ti, T5M, T5L);
							    T61 = ri[WS(ios, 53)];
						       }
						  }
						  T5P = ri[WS(ios, 37)];
						  T5S = ii[WS(ios, 37)];
						  T5V = ri[WS(ios, 21)];
						  {
						       E Tb8, T62, T5Q, Tbe;
						       Tb8 = T64 * T61;
						       T62 = T60 * T61;
						       T5Q = T5O * T5P;
						       Tbe = T5O * T5S;
						       Tb9 = FMA(T60, T65, Tb8);
						       T66 = FNMS(T64, T65, T62);
						       T5T = FMA(T5R, T5S, T5Q);
						       Tbf = FNMS(T5R, T5P, Tbe);
						       T5W = T3j * T5V;
						  }
						  T5X = ii[WS(ios, 21)];
					     }
					     T5U = T5N + T5T;
					     Tb5 = T5N - T5T;
					     Tbg = Tbd - Tbf;
					     TgO = Tbd + Tbf;
					     T5Y = FMA(T3m, T5X, T5W);
					     Tb6 = T3j * T5X;
					}
					{
					     E T6b, T6c, Tbi, Teo, TgR, TgP, Tba, T6f, T67, Tbh, Tb7;
					     T6b = ri[WS(ios, 61)];
					     T67 = T5Y + T66;
					     Tbh = T5Y - T66;
					     Tb7 = FNMS(T3m, T5V, Tb6);
					     T6c = T6a * T6b;
					     Tbi = Tbg + Tbh;
					     Teo = Tbg - Tbh;
					     TgR = T5U - T67;
					     T68 = T5U + T67;
					     TgP = Tb7 + Tb9;
					     Tba = Tb7 - Tb9;
					     T6f = ii[WS(ios, 61)];
					     T6u = ri[WS(ios, 45)];
					     {
						  E Tbb, Tep, TgQ, TaX;
						  Tbb = Tb5 - Tba;
						  Tep = Tb5 + Tba;
						  TgQ = TgO - TgP;
						  Tj5 = TgO + TgP;
						  T6g = FMA(T6e, T6f, T6c);
						  TaX = T6a * T6f;
						  Tez = FMA(KP414213562, Teo, Tep);
						  Teq = FNMS(KP414213562, Tep, Teo);
						  Tbj = FNMS(KP414213562, Tbi, Tbb);
						  Tbx = FMA(KP414213562, Tbb, Tbi);
						  TgS = TgQ - TgR;
						  Th5 = TgR + TgQ;
						  TaY = FNMS(T6e, T6b, TaX);
						  T6v = T6t * T6u;
					     }
					     T6y = ii[WS(ios, 45)];
					}
					{
					     E T6l, T6i, TaT, TaZ, T6j;
					     T6l = ii[WS(ios, 29)];
					     T6i = ri[WS(ios, 29)];
					     T6o = ri[WS(ios, 13)];
					     T6z = FMA(T6x, T6y, T6v);
					     TaT = T6t * T6y;
					     TaZ = T6k * T6i;
					     T6j = T6h * T6i;
					     T6p = T17 * T6o;
					     TaU = FNMS(T6x, T6u, TaT);
					     Tb0 = FMA(T6h, T6l, TaZ);
					     T6m = FNMS(T6k, T6l, T6j);
					     T6q = ii[WS(ios, 13)];
					}
				   }
				   {
					E T7I, Tci, T7T, Thf, Tcp, Tce, T7H, T7L, Tcf;
					{
					     E T7A, Tcm, Tco, T7G, T7J, T7K;
					     {
						  E T7S, T7O, T7F, T7C;
						  {
						       E TgU, TaQ, Tb3, Ter, TgT, TgV, TaV, TaW, Tes, TgW;
						       {
							    E T7x, Tb1, T6n, T7y, T7z, T6A, Tb2, TaS, Tcl, T6r, TaR;
							    T7x = ri[WS(ios, 3)];
							    TgU = TaY + Tb0;
							    Tb1 = TaY - Tb0;
							    T6n = T6g + T6m;
							    TaQ = T6g - T6m;
							    T6r = FMA(T19, T6q, T6p);
							    TaR = T17 * T6q;
							    T7y = T4 * T7x;
							    T7z = ii[WS(ios, 3)];
							    T6A = T6r + T6z;
							    Tb2 = T6r - T6z;
							    TaS = FNMS(T19, T6o, TaR);
							    Tcl = T4 * T7z;
							    Tb3 = Tb1 + Tb2;
							    Ter = Tb1 - Tb2;
							    TgT = T6n - T6A;
							    T6B = T6n + T6A;
							    TgV = TaS + TaU;
							    TaV = TaS - TaU;
							    T7A = FMA(T7, T7z, T7y);
							    Tcm = FNMS(T7, T7x, Tcl);
						       }
						       TaW = TaQ - TaV;
						       Tes = TaQ + TaV;
						       TgW = TgU - TgV;
						       Tj6 = TgU + TgV;
						       TeA = FNMS(KP414213562, Ter, Tes);
						       Tet = FMA(KP414213562, Tes, Ter);
						       Tb4 = FMA(KP414213562, Tb3, TaW);
						       Tby = FNMS(KP414213562, TaW, Tb3);
						       TgX = TgT + TgW;
						       Th6 = TgT - TgW;
						       T7S = ii[WS(ios, 51)];
						       T7O = ri[WS(ios, 51)];
						  }
						  T7F = ii[WS(ios, 35)];
						  T7C = ri[WS(ios, 35)];
						  T7I = ri[WS(ios, 19)];
						  {
						       E Tch, T7P, Tcn, T7D;
						       Tch = T7R * T7O;
						       T7P = T7N * T7O;
						       Tcn = T7E * T7C;
						       T7D = T7B * T7C;
						       Tci = FMA(T7N, T7S, Tch);
						       T7T = FNMS(T7R, T7S, T7P);
						       Tco = FMA(T7B, T7F, Tcn);
						       T7G = FNMS(T7E, T7F, T7D);
						       T7J = T2u * T7I;
						  }
						  T7K = ii[WS(ios, 19)];
					     }
					     Thf = Tcm + Tco;
					     Tcp = Tcm - Tco;
					     Tce = T7A - T7G;
					     T7H = T7A + T7G;
					     T7L = FMA(T2x, T7K, T7J);
					     Tcf = T2u * T7K;
					}
					{
					     E T83, Tc7, T87, Tc9, T89, T8f, T8j;
					     {
						  E T82, Tcr, TeH, Thi, Thg, Tcj, T7Z, Tc6, Tck, TeI, Thh;
						  T82 = ii[WS(ios, 59)];
						  {
						       E T7U, Tcq, Tcg, T7Y;
						       T7U = T7L + T7T;
						       Tcq = T7L - T7T;
						       Tcg = FNMS(T2x, T7I, Tcf);
						       T7Y = ri[WS(ios, 59)];
						       Tcr = Tcp + Tcq;
						       TeH = Tcp - Tcq;
						       Thi = T7H - T7U;
						       T7V = T7H + T7U;
						       Thg = Tcg + Tci;
						       Tcj = Tcg - Tci;
						       T7Z = T7X * T7Y;
						       Tc6 = T81 * T7Y;
						  }
						  Tck = Tce - Tcj;
						  TeI = Tce + Tcj;
						  Thh = Thf - Thg;
						  Tjg = Thf + Thg;
						  TeS = FMA(KP414213562, TeH, TeI);
						  TeJ = FNMS(KP414213562, TeI, TeH);
						  Tcs = FNMS(KP414213562, Tcr, Tck);
						  TcG = FMA(KP414213562, Tck, Tcr);
						  Thj = Thh - Thi;
						  Thw = Thi + Thh;
						  T83 = FNMS(T81, T82, T7Z);
						  Tc7 = FMA(T7X, T82, Tc6);
					     }
					     T8f = ri[WS(ios, 43)];
					     T8j = ii[WS(ios, 43)];
					     {
						  E T84, T86, T8g, Tc2, T85, Tc8;
						  T84 = ri[WS(ios, 27)];
						  T86 = ii[WS(ios, 27)];
						  T8g = T8e * T8f;
						  Tc2 = T8e * T8j;
						  T85 = T2 * T84;
						  Tc8 = T2 * T86;
						  T8k = FMA(T8i, T8j, T8g);
						  Tc3 = FNMS(T8i, T8f, Tc2);
						  T87 = FMA(Tg, T86, T85);
						  Tc9 = FNMS(Tg, T84, Tc8);
					     }
					     T8b = ii[WS(ios, 11)];
					     T89 = ri[WS(ios, 11)];
					     TbZ = T83 - T87;
					     T88 = T83 + T87;
					     Thl = Tc7 + Tc9;
					     Tca = Tc7 - Tc9;
					     Tc0 = Tx * T89;
					     T8a = Tu * T89;
					}
				   }
			      }
			 }
		    }
		    {
			 E TeM, TeT, TcH, Tcd, Thx, Tho, Tkw, Tkv, Tl6, Tl5;
			 {
			      E TiI, Tkp, TiQ, TiS, TiL, Tkq, TiP, TiV, Tjf, Tjd, Tjc, Tji, Tj4, Tj2, Tj1;
			      E Tj7, Tkh, Tki;
			      {
				   E TjG, T2I, Tkj, T4N, Tkk, Tkf, Tk5, TjJ, T8o, Tk2, TjL, T6D, TjY, TjU, Tk1;
				   E TjO;
				   {
					E T3L, Tjh, T8m, T4M, Tk6, Tke, TjH, TjI;
					{
					     E T1C, Tc1, T8c, T2H, Tc4, Thm;
					     TiI = TY - T1B;
					     T1C = TY + T1B;
					     Tc1 = FMA(Tu, T8b, Tc0);
					     T8c = FNMS(Tx, T8b, T8a);
					     T2H = T27 + T2G;
					     Tkp = T2G - T27;
					     TiQ = T39 - T3K;
					     T3L = T39 + T3K;
					     Tc4 = Tc1 - Tc3;
					     Thm = Tc1 + Tc3;
					     {
						  E Tcb, T8l, TeL, Tc5;
						  Tcb = T8c - T8k;
						  T8l = T8c + T8k;
						  TjG = T1C - T2H;
						  T2I = T1C + T2H;
						  TeL = TbZ + Tc4;
						  Tc5 = TbZ - Tc4;
						  {
						       E Thn, TeK, Tcc, Thk;
						       Tjh = Thl + Thm;
						       Thn = Thl - Thm;
						       TeK = Tca - Tcb;
						       Tcc = Tca + Tcb;
						       T8m = T88 + T8l;
						       Thk = T88 - T8l;