t2_64.c

快速fft变换

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