hf2_64.c

快速傅立叶变换库函数

					     TbM = T7o * T7s;
					     T7b = T79 * T7a;
					     TbR = FNMS(T1k, T75, TbQ);
					     T7t = FMA(T7r, T7s, T7q);
					     TbN = FNMS(T7r, T7p, TbM);
					     T7d = iio[-WS(ios, 24)];
					     T7i = rio[WS(ios, 55)];
					     T7m = iio[-WS(ios, 8)];
					}
				   }
				   {
					E T5k, T5h, T5l, TaI, T5G, TaE, T5o, T5t, T5x;
					{
					     E T5e, T5f, T5g, T5B, T5F, TaH, T5C, TaD;
					     {
						  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 = T8 * 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)];
					     T5k = rio[WS(ios, 41)];
					     T5h = FMA(Tc, T5g, T5f);
					     TaH = T8 * T5g;
					     T5C = T5A * T5B;
					     TaD = T5A * T5F;
					     T5l = T5j * T5k;
					     TaI = FNMS(Tc, T5e, TaH);
					     T5G = FMA(T5E, T5F, T5C);
					     TaE = FNMS(T5E, T5B, TaD);
					     T5o = iio[-WS(ios, 22)];
					     T5t = rio[WS(ios, 57)];
					     T5x = iio[-WS(ios, 6)];
					}
					{
					     E T4l, T4m, T4n, T4F, T4I, Ta6, T4G, Ta2;
					     {
						  E T5q, TaM, TaK, T5y, TaC, TgK, TaL;
						  T4l = rio[WS(ios, 6)];
						  {
						       E T5p, TaJ, T5u, TaB;
						       T5p = FMA(T5n, T5o, T5l);
						       TaJ = T5j * T5o;
						       T5u = T5s * T5t;
						       TaB = T5s * T5x;
						       T5q = T5h + T5p;
						       TaM = T5h - T5p;
						       TaK = FNMS(T5n, T5k, TaJ);
						       T5y = FMA(T5w, T5x, T5u);
						       TaC = FNMS(T5w, T5t, TaB);
						       T4m = T3g * T4l;
						  }
						  TgK = TaI + TaK;
						  TaL = TaI - TaK;
						  {
						       E T5H, TaA, TgL, TaF;
						       T5H = T5y + T5G;
						       TaA = T5y - T5G;
						       TgL = TaC + TaE;
						       TaF = TaC - TaE;
						       {
							    E TaN, Tbt, TaG, Tbu;
							    TaN = TaL - TaM;
							    Tbt = TaM + TaL;
							    T5I = T5q + T5H;
							    Th3 = T5H - T5q;
							    TaG = TaA + TaF;
							    Tbu = TaA - TaF;
							    TgM = TgK - TgL;
							    Tj0 = TgK + TgL;
							    Tbv = Tbt - Tbu;
							    Tem = Tbt + Tbu;
							    TaO = TaG - TaN;
							    Tex = TaN + TaG;
							    T4n = iio[-WS(ios, 57)];
						       }
						  }
					     }
					     T4F = rio[WS(ios, 22)];
					     T4I = iio[-WS(ios, 41)];
					     T4q = rio[WS(ios, 38)];
					     T4o = FMA(T3i, T4n, T4m);
					     Ta6 = T3g * T4n;
					     T4G = T4E * T4F;
					     Ta2 = T4E * T4I;
					     T4r = T4p * T4q;
					     Ta7 = FNMS(T3i, T4l, Ta6);
					     T4J = FMA(T4H, T4I, T4G);
					     Ta3 = FNMS(T4H, T4F, Ta2);
					     T4t = iio[-WS(ios, 25)];
					     T4y = rio[WS(ios, 54)];
					     T4C = iio[-WS(ios, 9)];
					}
				   }
			      }
			      {
				   E T84, T83, T85, Tc1, T8k, Tca, T86, T89, T8b;
				   {
					E T7C, T7A, T7D, Tcg, T7T, Tcp, T7F, T7I, T7K;
					{
					     E T7x, T7y, T7z, T7O, T7S, Tcf, T7P, Tco;
					     {
						  E T4v, Tab, Ta9, T4D, Ta1, Tgx, Taa;
						  T7x = rio[WS(ios, 3)];
						  {
						       E T4u, Ta8, T4z, Ta0;
						       T4u = FMA(T4s, T4t, T4r);
						       Ta8 = T4p * T4t;
						       T4z = T4x * T4y;
						       Ta0 = T4x * T4C;
						       T4v = T4o + T4u;
						       Tab = T4o - T4u;
						       Ta9 = FNMS(T4s, T4q, Ta8);
						       T4D = FMA(T4B, T4C, T4z);
						       Ta1 = FNMS(T4B, T4y, Ta0);
						       T7y = T3 * T7x;
						  }
						  Tgx = Ta7 + Ta9;
						  Taa = Ta7 - Ta9;
						  {
						       E T4K, T9Z, Tgy, Ta4;
						       T4K = T4D + T4J;
						       T9Z = T4D - T4J;
						       Tgy = Ta1 + Ta3;
						       Ta4 = Ta1 - Ta3;
						       {
							    E Tac, Tam, Ta5, Tan;
							    Tac = Taa - Tab;
							    Tam = Tab + Taa;
							    T4L = T4v + T4K;
							    TgE = T4K - T4v;
							    Ta5 = T9Z + Ta4;
							    Tan = T9Z - Ta4;
							    Tgz = Tgx - Tgy;
							    TiU = Tgx + Tgy;
							    Tao = Tam - Tan;
							    Ted = Tam + Tan;
							    Tad = Ta5 - Tac;
							    Teg = Tac + Ta5;
							    T7z = iio[-WS(ios, 60)];
						       }
						  }
					     }
					     T7O = rio[WS(ios, 51)];
					     T7S = iio[-WS(ios, 12)];
					     T7C = rio[WS(ios, 35)];
					     T7A = FMA(T6, T7z, T7y);
					     Tcf = T3 * T7z;
					     T7P = T7N * T7O;
					     Tco = T7N * T7S;
					     T7D = T7B * T7C;
					     Tcg = FNMS(T6, T7x, Tcf);
					     T7T = FMA(T7R, T7S, T7P);
					     Tcp = FNMS(T7R, T7O, Tco);
					     T7F = iio[-WS(ios, 28)];
					     T7I = rio[WS(ios, 19)];
					     T7K = iio[-WS(ios, 44)];
					}
					{
					     E T7Y, T7Z, T82, T8f, T8j, Tc0, T8g, Tc9;
					     {
						  E T7H, Tcl, Tci, T7L, Tcn;
						  T7Y = rio[WS(ios, 59)];
						  {
						       E T7G, Tch, T7J, Tcm;
						       T7G = FMA(T7E, T7F, T7D);
						       Tch = T7B * T7F;
						       T7J = T2u * T7I;
						       Tcm = T2u * T7K;
						       T7H = T7A + T7G;
						       Tcl = T7A - T7G;
						       Tci = FNMS(T7E, T7C, Tch);
						       T7L = FMA(T2x, T7K, T7J);
						       Tcn = FNMS(T2x, T7I, Tcm);
						       T7Z = T7X * T7Y;
						  }
						  {
						       E Thg, Tcj, T7U, Tce;
						       Thg = Tcg + Tci;
						       Tcj = Tcg - Tci;
						       T7U = T7L + T7T;
						       Tce = T7T - T7L;
						       {
							    E Thh, Tcq, Tck, TeI;
							    Thh = Tcn + Tcp;
							    Tcq = Tcn - Tcp;
							    Tck = Tce - Tcj;
							    TeI = Tcj + Tce;
							    {
								 E Thf, Tcr, TeH, Thi;
								 Thf = T7H - T7U;
								 T7V = T7H + T7U;
								 Tcr = Tcl - Tcq;
								 TeH = Tcl + Tcq;
								 Thi = Thg - Thh;
								 Tjg = Thg + Thh;
								 TeS = FNMS(KP414213562, TeH, TeI);
								 TeJ = FMA(KP414213562, TeI, TeH);
								 Tcs = FNMS(KP414213562, Tcr, Tck);
								 TcH = FMA(KP414213562, Tck, Tcr);
								 Thj = Thf + Thi;
								 Thw = Thi - Thf;
								 T82 = iio[-WS(ios, 4)];
							    }
						       }
						  }
					     }
					     T8f = rio[WS(ios, 43)];
					     T8j = iio[-WS(ios, 20)];
					     T84 = rio[WS(ios, 27)];
					     T83 = FMA(T81, T82, T7Z);
					     Tc0 = T7X * T82;
					     T8g = T8e * T8f;
					     Tc9 = T8e * T8j;
					     T85 = Te * T84;
					     Tc1 = FNMS(T81, T7Y, Tc0);
					     T8k = FMA(T8i, T8j, T8g);
					     Tca = FNMS(T8i, T8f, Tc9);
					     T86 = iio[-WS(ios, 36)];
					     T89 = rio[WS(ios, 11)];
					     T8b = iio[-WS(ios, 52)];
					}
				   }
				   {
					E T5P, T5N, T5Q, TaR, T66, Tb1, T5S, T5V, T5X;
					{
					     E T5K, T5L, T5M, T61, T65, TaQ, T62, Tb0;
					     {
						  E T88, Tc6, Tc3, T8c, Tc8;
						  T5K = rio[WS(ios, 5)];
						  {
						       E T87, Tc2, T8a, Tc7;
						       T87 = FMA(Ti, T86, T85);
						       Tc2 = Te * T86;
						       T8a = Tu * T89;
						       Tc7 = Tu * T8b;
						       T88 = T83 + T87;
						       Tc6 = T83 - T87;
						       Tc3 = FNMS(Ti, T84, Tc2);
						       T8c = FMA(Tx, T8b, T8a);
						       Tc8 = FNMS(Tx, T89, Tc7);
						       T5L = Td * T5K;
						  }
						  {
						       E Thl, Tc4, T8l, TbZ;
						       Thl = Tc1 + Tc3;
						       Tc4 = Tc1 - Tc3;
						       T8l = T8c + T8k;
						       TbZ = T8k - T8c;
						       {
							    E Thm, Tcb, Tc5, TeL;
							    Thm = Tc8 + Tca;
							    Tcb = Tc8 - Tca;
							    Tc5 = TbZ - Tc4;
							    TeL = Tc4 + TbZ;
							    {
								 E Thk, Tcc, TeK, Thn;
								 Thk = T88 - T8l;
								 T8m = T88 + T8l;
								 Tcc = Tc6 - Tcb;
								 TeK = Tc6 + Tcb;
								 Thn = Thl - Thm;
								 Tjh = Thl + Thm;
								 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;
								 T5M = iio[-WS(ios, 58)];
							    }
						       }
						  }
					     }
					     T61 = rio[WS(ios, 53)];
					     T65 = iio[-WS(ios, 10)];
					     T5P = rio[WS(ios, 37)];
					     T5N = FMA(Th, T5M, T5L);
					     TaQ = Td * T5M;
					     T62 = T60 * T61;
					     Tb0 = T60 * T65;
					     T5Q = T5O * T5P;
					     TaR = FNMS(Th, T5K, TaQ);
					     T66 = FMA(T64, T65, T62);
					     Tb1 = FNMS(T64, T61, Tb0);
					     T5S = iio[-WS(ios, 26)];
					     T5V = rio[WS(ios, 21)];
					     T5X = iio[-WS(ios, 42)];
					}
					{
					     E T6b, T6c, T6f, T6u, T6y, Tb5, T6v, Tbf;
					     {
						  E T5U, TaX, TaT, T5Y, TaZ;
						  T6b = rio[WS(ios, 61)];
						  {
						       E T5T, TaS, T5W, TaY;
						       T5T = FMA(T5R, T5S, T5Q);
						       TaS = T5O * T5S;
						       T5W = T3j * T5V;
						       TaY = T3j * T5X;
						       T5U = T5N + T5T;
						       TaX = T5N - T5T;
						       TaT = FNMS(T5R, T5P, TaS);
						       T5Y = FMA(T3m, T5X, T5W);
						       TaZ = FNMS(T3m, T5V, TaY);
						       T6c = T6a * T6b;
						  }
						  {
						       E TgP, TaU, T67, TaV;
						       TgP = TaR + TaT;
						       TaU = TaR - TaT;
						       T67 = T5Y + T66;
						       TaV = T5Y - T66;
						       {
							    E TgQ, Tb2, TaW, Tep;
							    TgQ = TaZ + Tb1;
							    Tb2 = TaZ - Tb1;
							    TaW = TaU + TaV;
							    Tep = TaU - TaV;
							    {
								 E TgO, Tb3, Teo, TgR;
								 TgO = T5U - T67;
								 T68 = T5U + T67;
								 Tb3 = TaX - Tb2;
								 Teo = TaX + Tb2;
								 TgR = TgP - TgQ;
								 Tj5 = TgP + TgQ;
								 Tez = FNMS(KP414213562, Teo, Tep);
								 Teq = FMA(KP414213562, Tep, Teo);
								 Tb4 = FMA(KP414213562, Tb3, TaW);
								 Tby = FNMS(KP414213562, TaW, Tb3);
								 TgS = TgO + TgR;
								 Th5 = TgR - TgO;
								 T6f = iio[-WS(ios, 2)];
							    }
						       }
						  }
					     }
					     T6u = rio[WS(ios, 45)];
					     T6y = iio[-WS(ios, 18)];
					     T6i = rio[WS(ios, 29)];
					     T6g = FMA(T6e, T6f, T6c);
					     Tb5 = T6a * T6f;
					     T6v = T6t * T6u;
					     Tbf = T6t * T6y;
					     T6j = T6h * T6i;
					     Tb6 = FNMS(T6e, T6b, Tb5);
					     T6z = FMA(T6x, T6y, T6v);
					     Tbg = FNMS(T6x, T6u, Tbf);
					     T6l = iio[-WS(ios, 34)];
					     T6o = rio[WS(ios, 13)];
					     T6q = iio[-WS(ios, 50)];
					}
				   }
			      }
			 }
		    }
		    {
			 E TeA, Tet, Tbj, Tbx, TgX, Th6, Tkw, Tkv, Tl6, Tl5;
			 {
			      E TiI, Tkp, TiQ, TiS, TiL, Tkq, TiP, TiV, Tjf, Tja, Tjd, 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 T6B, Tj6, T3L, T4M, Tk6, Tke, TjH, TjI;
					{
					     E T1C, T6n, Tbc, Tb8, T6r, Tbe, T2H;
					     T1C = TY + T1B;
					     TiI = TY - T1B;
					     {
						  E T6m, Tb7, T6p, Tbd;
						  T6m = FMA(T6k, T6l, T6j);
						  Tb7 = T6h * T6l;
						  T6p = T17 * T6o;
						  Tbd = T17 * T6q;
						  T6n = T6g + T6m;
						  Tbc = T6g - T6m;
						  Tb8 = FNMS(T6k, T6i, Tb7);
						  T6r = FMA(T19, T6q, T6p);
						  Tbe = FNMS(T19, T6o, Tbd);
						  T2H = T27 + T2G;