hf2_64.c

最新的FFT程序

					     {
						  E T7f, TbV, TbT, T7n, TbL, Thb, TbU;
						  T5e = rio[WS(ios, 9)];
						  {
						       E T7e, TbS, T7j, TbK;
						       T7e = FMA(T7c, T7d, T7b);
						       TbS = T79 * T7d;
						       T7j = T7h * T7i;
						       TbK = T7h * T7m;
						       T7f = T78 + T7e;
						       TbV = T78 - T7e;
						       TbT = FNMS(T7c, T7a, TbS);
						       T7n = FMA(T7l, T7m, T7j);
						       TbL = FNMS(T7l, T7i, TbK);
						       T5f = T3 * T5e;
						  }
						  Thb = TbR + TbT;
						  TbU = TbR - TbT;
						  {
						       E T7u, TbJ, Thc, TbO;
						       T7u = T7n + T7t;
						       TbJ = T7n - T7t;
						       Thc = TbL + TbN;
						       TbO = TbL - TbN;
						       {
							    E TbW, TcD, TbP, TcC;
							    TbW = TbU - TbV;
							    TcD = TbV + TbU;
							    T7v = T7f + T7u;
							    Thu = T7u - T7f;
							    TbP = TbJ + TbO;
							    TcC = TbJ - TbO;
							    Thd = Thb - Thc;
							    Tjc = Thb + Thc;
							    TcE = TcC - TcD;
							    TeF = TcD + TcC;
							    TbX = TbP - TbW;
							    TeQ = TbW + TbP;
							    T5g = iio[-WS(ios, 54)];
						       }
						  }
					     }
					     T5B = rio[WS(ios, 25)];
					     T5F = iio[-WS(ios, 38)];
					     T5o = iio[-WS(ios, 22)];
					     T5h = FMA(Tb, T5g, T5f);
					     {
						  E TaH, T5C, TaD, T5k;
						  TaH = T3 * T5g;
						  T5C = T5A * T5B;
						  TaD = T5A * T5F;
						  T5k = rio[WS(ios, 41)];
						  TaI = FNMS(Tb, T5e, TaH);
						  T5G = FMA(T5E, T5F, T5C);
						  TaE = FNMS(T5E, T5B, TaD);
						  T5l = T5j * T5k;
						  TaJ = T5n * T5k;
					     }
					     T5t = rio[WS(ios, 57)];
					     T5x = iio[-WS(ios, 6)];
					     T5p = FNMS(T5n, T5o, T5l);
					     TaK = FMA(T5j, T5o, TaJ);
					     T5u = T5s * T5t;
					     TaB = T5s * T5x;
					}
					TaM = T5h - T5p;
					T5q = T5h + T5p;
					{
					     E TgK, TaL, T5y, TaC, TaN, TaG;
					     TgK = TaI + TaK;
					     TaL = TaI - TaK;
					     T5y = FMA(T5w, T5x, T5u);
					     TaC = FNMS(T5w, T5t, TaB);
					     {
						  E T5H, TaA, TgL, TaF;
						  T5H = T5y + T5G;
						  TaA = T5y - T5G;
						  TgL = TaC + TaE;
						  TaF = TaC - TaE;
						  Tbt = TaM + TaL;
						  TaN = TaL - TaM;
						  Th3 = T5H - T5q;
						  T5I = T5q + T5H;
						  Tbu = TaA - TaF;
						  TaG = TaA + TaF;
						  Tj0 = TgK + TgL;
						  TgM = TgK - TgL;
					     }
					     Tex = TaN + TaG;
					     TaO = TaG - TaN;
					}
				   }
				   {
					E T4n, T4l, T4F, T4I;
					T4n = iio[-WS(ios, 57)];
					Tem = Tbt + Tbu;
					Tbv = Tbt - Tbu;
					T4l = rio[WS(ios, 6)];
					T4F = rio[WS(ios, 22)];
					T4I = iio[-WS(ios, 41)];
					T4t = iio[-WS(ios, 25)];
					{
					     E Ta6, T4m, T4G, Ta2;
					     Ta6 = T3i * T4l;
					     T4m = T3g * T4l;
					     T4G = T4E * T4F;
					     Ta2 = T4E * T4I;
					     Ta7 = FMA(T3g, T4n, Ta6);
					     T4o = FNMS(T3i, T4n, T4m);
					     T4J = FMA(T4H, T4I, T4G);
					     Ta3 = FNMS(T4H, T4F, Ta2);
					     T4q = rio[WS(ios, 38)];
					}
					T4C = iio[-WS(ios, 9)];
					T4y = rio[WS(ios, 54)];
				   }
			      }
			      {
				   E T84, Tc1, T83, T8k, Tca, T85, T86, T8b, T89;
				   {
					E Tck, TeI, Thf, Thi, TeH, Tcr;
					{
					     E T7A, Tcg, Tcp, T7T, T7I, T7K, T7G, Tci, T7J, Tcm;
					     {
						  E T7x, T7y, T7z, T7S, T7O, T7F, T7D, Tch;
						  T7x = rio[WS(ios, 3)];
						  {
						       E Ta8, T4r, Ta0, T4z;
						       Ta8 = T4s * T4q;
						       T4r = T4p * T4q;
						       Ta0 = T4B * T4y;
						       T4z = T4x * T4y;
						       {
							    E Ta9, T4u, Ta1, T4D;
							    Ta9 = FMA(T4p, T4t, Ta8);
							    T4u = FNMS(T4s, T4t, T4r);
							    Ta1 = FMA(T4x, T4C, Ta0);
							    T4D = FNMS(T4B, T4C, T4z);
							    {
								 E Taa, Tgx, Tab, T4v;
								 Taa = Ta7 - Ta9;
								 Tgx = Ta7 + Ta9;
								 Tab = T4o - T4u;
								 T4v = T4o + T4u;
								 {
								      E Ta4, Tgy, T9Z, T4K;
								      Ta4 = Ta1 - Ta3;
								      Tgy = Ta1 + Ta3;
								      T9Z = T4D - T4J;
								      T4K = T4D + T4J;
								      {
									   E Tam, Tac, Tan, Ta5;
									   Tam = Tab + Taa;
									   Tac = Taa - Tab;
									   TiU = Tgx + Tgy;
									   Tgz = Tgx - Tgy;
									   Tan = T9Z - Ta4;
									   Ta5 = T9Z + Ta4;
									   TgE = T4K - T4v;
									   T4L = T4v + T4K;
									   Ted = Tam + Tan;
									   Tao = Tam - Tan;
									   Teg = Tac + Ta5;
									   Tad = Ta5 - Tac;
									   T7y = T4 * T7x;
								      }
								 }
							    }
						       }
						  }
						  T7z = iio[-WS(ios, 60)];
						  T7S = iio[-WS(ios, 12)];
						  T7O = rio[WS(ios, 51)];
						  T7F = iio[-WS(ios, 28)];
						  T7A = FMA(T7, T7z, T7y);
						  {
						       E Tcf, Tco, T7P, T7C;
						       Tcf = T4 * T7z;
						       Tco = T7R * T7O;
						       T7P = T7N * T7O;
						       T7C = rio[WS(ios, 35)];
						       Tcg = FNMS(T7, T7x, Tcf);
						       Tcp = FMA(T7N, T7S, Tco);
						       T7T = FNMS(T7R, T7S, T7P);
						       T7D = T7B * T7C;
						       Tch = T7E * T7C;
						  }
						  T7I = rio[WS(ios, 19)];
						  T7K = iio[-WS(ios, 44)];
						  T7G = FNMS(T7E, T7F, T7D);
						  Tci = FMA(T7B, T7F, Tch);
						  T7J = T2u * T7I;
						  Tcm = T2u * T7K;
					     }
					     {
						  E Tcl, T7H, Thg, Tcj, T7L, Tcn;
						  Tcl = T7A - T7G;
						  T7H = T7A + T7G;
						  Thg = Tcg + Tci;
						  Tcj = Tcg - Tci;
						  T7L = FMA(T2x, T7K, T7J);
						  Tcn = FNMS(T2x, T7I, Tcm);
						  {
						       E T7U, Tce, Tcq, Thh;
						       T7U = T7L + T7T;
						       Tce = T7T - T7L;
						       Tcq = Tcn - Tcp;
						       Thh = Tcn + Tcp;
						       Tck = Tce - Tcj;
						       TeI = Tcj + Tce;
						       Thf = T7H - T7U;
						       T7V = T7H + T7U;
						       Thi = Thg - Thh;
						       Tjg = Thg + Thh;
						       TeH = Tcl + Tcq;
						       Tcr = Tcl - Tcq;
						  }
					     }
					}
					{
					     E T82, T7Y, T8f, T8j;
					     T82 = iio[-WS(ios, 4)];
					     Thw = Thi - Thf;
					     Thj = Thf + Thi;
					     TcH = FMA(KP414213562, Tck, Tcr);
					     Tcs = FNMS(KP414213562, Tcr, Tck);
					     TeJ = FMA(KP414213562, TeI, TeH);
					     TeS = FNMS(KP414213562, TeH, TeI);
					     T7Y = rio[WS(ios, 59)];
					     T8f = rio[WS(ios, 43)];
					     T8j = iio[-WS(ios, 20)];
					     T84 = rio[WS(ios, 27)];
					     {
						  E Tc0, T7Z, T8g, Tc9;
						  Tc0 = T81 * T7Y;
						  T7Z = T7X * T7Y;
						  T8g = T8e * T8f;
						  Tc9 = T8e * T8j;
						  Tc1 = FMA(T7X, T82, Tc0);
						  T83 = FNMS(T81, T82, T7Z);
						  T8k = FMA(T8i, T8j, T8g);
						  Tca = FNMS(T8i, T8f, Tc9);
						  T85 = T2 * T84;
					     }
					     T86 = iio[-WS(ios, 36)];
					     T8b = iio[-WS(ios, 52)];
					     T89 = rio[WS(ios, 11)];
					}
				   }
				   {
					E Tb1, T66, TaX, T5U, TgP, TaU, T5Y, TaZ;
					{
					     E T5N, TaR, T5T, T5V, T5X, TaT, T5W, TaY;
					     {
						  E T65, T61, T5P, TaS;
						  {
						       E T5M, T5L, TaQ, T87, Tc2;
						       T5M = iio[-WS(ios, 58)];
						       T87 = FMA(Tg, T86, T85);
						       Tc2 = T2 * T86;
						       {
							    E Tc7, T8a, T88, Tc6;
							    Tc7 = Tx * T89;
							    T8a = Tu * T89;
							    T88 = T83 + T87;
							    Tc6 = T83 - T87;
							    {
								 E Tc3, Tc8, T8c, T5K;
								 Tc3 = FNMS(Tg, T84, Tc2);
								 Tc8 = FMA(Tu, T8b, Tc7);
								 T8c = FNMS(Tx, T8b, T8a);
								 T5K = rio[WS(ios, 5)];
								 {
								      E Thl, Tc4, Thm, Tcb;
								      Thl = Tc1 + Tc3;
								      Tc4 = Tc1 - Tc3;
								      Thm = Tc8 + Tca;
								      Tcb = Tc8 - Tca;
								      {
									   E T8l, TbZ, Tcc, TeK;
									   T8l = T8c + T8k;
									   TbZ = T8k - T8c;
									   Tcc = Tc6 - Tcb;
									   TeK = Tc6 + Tcb;
									   {
										E Thn, Tc5, TeL, Thk;
										Thn = Thl - Thm;
										Tjh = Thl + Thm;
										Tc5 = TbZ - Tc4;
										TeL = Tc4 + TbZ;
										Thk = T88 - T8l;
										T8m = T88 + T8l;
										TeT = FMA(KP414213562, TeK, TeL);
										TeM = FNMS(KP414213562, TeL, TeK);
										Tcd = FMA(KP414213562, Tcc, Tc5);
										TcG = FNMS(KP414213562, Tc5, Tcc);
										Tho = Thk - Thn;
										Thx = Thk + Thn;
										T5L = Te * T5K;
										TaQ = Ti * T5K;
									   }
								      }
								 }
							    }
						       }
						       T65 = iio[-WS(ios, 10)];
						       T61 = rio[WS(ios, 53)];
						       T5N = FNMS(Ti, T5M, T5L);
						       TaR = FMA(Te, T5M, TaQ);
						  }
						  {
						       E T5S, T5Q, Tb0, T62;
						       T5P = rio[WS(ios, 37)];
						       Tb0 = T64 * T61;
						       T62 = T60 * T61;
						       T5S = iio[-WS(ios, 26)];
						       T5Q = T5O * T5P;
						       Tb1 = FMA(T60, T65, Tb0);
						       T66 = FNMS(T64, T65, T62);
						       TaS = T5O * T5S;
						       T5T = FMA(T5R, T5S, T5Q);
						  }
						  T5V = rio[WS(ios, 21)];
						  T5X = iio[-WS(ios, 42)];
						  TaT = FNMS(T5R, T5P, TaS);
					     }
					     TaX = T5N - T5T;
					     T5U = T5N + T5T;
					     T5W = T3j * T5V;
					     TaY = T3j * T5X;
					     TgP = TaR + TaT;
					     TaU = TaR - TaT;
					     T5Y = FMA(T3m, T5X, T5W);
					     TaZ = FNMS(T3m, T5V, TaY);
					}
					{
					     E T6g, Tb6, T6z, Tbg, T6o, T6q, T6m, Tb8, T6p, Tbd;
					     {
						  E T6b, T6c, T6f, T6u, T6y, TaV, T67, T6l, T6j, Tb7;
						  T6b = rio[WS(ios, 61)];
						  TaV = T5Y - T66;
						  T67 = T5Y + T66;
						  {
						       E Tb2, TgQ, Tep, TaW;
						       Tb2 = TaZ - Tb1;
						       TgQ = TaZ + Tb1;
						       Tep = TaU - TaV;
						       TaW = TaU + TaV;
						       {
							    E TgO, Teo, Tb3, TgR;
							    T68 = T5U + T67;
							    TgO = T5U - T67;
							    Teo = TaX + Tb2;
							    Tb3 = TaX - Tb2;
							    Tj5 = TgP + TgQ;
							    TgR = TgP - TgQ;
							    Teq = FMA(KP414213562, Tep, Teo);
							    Tez = FNMS(KP414213562, Teo, Tep);
							    Tby = FNMS(KP414213562, TaW, Tb3);
							    Tb4 = FMA(KP414213562, Tb3, TaW);
							    Th5 = TgR - TgO;
							    TgS = TgO + TgR;
							    T6c = T6a * T6b;
						       }
						  }
						  T6f = iio[-WS(ios, 2)];
						  T6u = rio[WS(ios, 45)];
						  T6y = iio[-WS(ios, 18)];
						  T6l = iio[-WS(ios, 34)];
						  T6g = FMA(T6e, T6f, T6c);
						  {
						       E Tb5, T6v, Tbf, T6i;
						       Tb5 = T6a * T6f;
						       T6v = T6t * T6u;
						       Tbf = T6t * T6y;
						       T6i = rio[WS(ios, 29)];
						       Tb6 = FNMS(T6e, T6b, Tb5);
						       T6z = FMA(T6x, T6y, T6v);
						       Tbg = FNMS(T6x, T6u, Tbf);
						       T6j = T6h * T6i;
						       Tb7 = T6k * T6i;
						  }
						  T6o = rio[WS(ios, 13)];
						  T6q = iio[-WS(ios, 50)];
						  T6m = FNMS(T6k, T6l, T6j);
						  Tb8 = FMA(T6h, T6l, Tb7);
						  T6p = T17 * T6o;
						  Tbd = T17 * T6q;
					     }
					     {
						  E Tbc, T6n, TgU, Tb9, T6r, Tbe;
						  Tbc = T6g - T6m;
						  T6n = T6g + T6m;
						  TgU = Tb6 + Tb8;
						  Tb9 = Tb6 - Tb8;
						  T6r = FMA(T19, T6q, T6p);
						  Tbe = FNMS(T19, T6o, Tbd);
						  {
						       E T6A, Tba, Tbh, TgV;
						       T6A = T6r + T6z;
						       Tba = T6r - T6z;
						       Tbh = Tbe - Tbg;
						       TgV = Tbe + Tbg;
						       Tbb = Tb9 + Tba;
						       Tes = Tb9 - Tba;
						       TgT = T6n - T6A;
						       T6B = T6n + T6A;
						       TgW = TgU - TgV;
						       Tj6 = TgU + TgV;
						       Ter = Tbc + Tbh;
						       Tbi = Tbc - Tbh;
						  }
					     }
					}
				   }
			      }
			 }
		    }
		    {
			 E Th6, TgX, Tbx, Tbj, Tet, TeA, Tkw, Tkv, Tl6, Tl5;
			 {