t2_64.c

最新的FFT程序

						       T98 = T96 + T97;
						       Te1 = T96 - T97;
						       T2G = T2o + T2F;
						       Tge = T2o - T2F;
						       Tgg = T9b + T9d;
						       T9e = T9b - T9d;
						       T2K = Tr * T2J;
						       T9j = Tt * T2J;
						  }
						  T9f = T99 - T9e;
						  Te0 = T99 + T9e;
						  Tgh = Tgf - Tgg;
						  TiK = Tgf + Tgg;
						  T2M = FNMS(Tt, T2L, T2K);
						  T9k = FMA(Tr, T2L, T9j);
					     }
					     T36 = ii[WS(ios, 50)];
					     T32 = ri[WS(ios, 50)];
					     T2O = ri[WS(ios, 34)];
					     T2R = ii[WS(ios, 34)];
					     T2Y = ii[WS(ios, 18)];
					     {
						  E T9J, T33, T2P, T9l;
						  T9J = T35 * T32;
						  T33 = T31 * T32;
						  T2P = T2N * T2O;
						  T9l = T2N * T2R;
						  T9K = FMA(T31, T36, T9J);
						  T37 = FNMS(T35, T36, T33);
						  T2S = FMA(T2Q, T2R, T2P);
						  T9m = FNMS(T2Q, T2O, T9l);
						  T2V = ri[WS(ios, 18)];
					     }
					}
					{
					     E T71, T6X, T6K, T6O;
					     {
						  E T6E, T9G, T2T, T9n, Tgl, T6F, T6G, T9I, T2Z, T9H, T2W, TbC;
						  T6E = ri[WS(ios, 63)];
						  T9G = T2M - T2S;
						  T2T = T2M + T2S;
						  T9n = T9k - T9m;
						  Tgl = T9k + T9m;
						  T9H = T2X * T2V;
						  T2W = T2U * T2V;
						  T6F = TF * T6E;
						  T6G = ii[WS(ios, 63)];
						  T9I = FMA(T2U, T2Y, T9H);
						  T2Z = FNMS(T2X, T2Y, T2W);
						  T71 = ii[WS(ios, 47)];
						  T6H = FMA(TN, T6G, T6F);
						  TbC = TF * T6G;
						  {
						       E Tgm, T9L, T38, T9o;
						       Tgm = T9I + T9K;
						       T9L = T9I - T9K;
						       T38 = T2Z + T37;
						       T9o = T2Z - T37;
						       TbD = FNMS(TN, T6E, TbC);
						       Te8 = T9G + T9L;
						       T9M = T9G - T9L;
						       TiN = Tgl + Tgm;
						       Tgn = Tgl - Tgm;
						       Te5 = T9n - T9o;
						       T9p = T9n + T9o;
						       Tgq = T2T - T38;
						       T39 = T2T + T38;
						       T6X = ri[WS(ios, 47)];
						  }
					     }
					     T6K = ri[WS(ios, 31)];
					     T6O = ii[WS(ios, 31)];
					     T6T = ii[WS(ios, 15)];
					     {
						  E Tcy, T6Y, T6L, TbE;
						  Tcy = T70 * T6X;
						  T6Y = T6W * T6X;
						  T6L = T6J * T6K;
						  TbE = T6J * T6O;
						  Tcz = FMA(T6W, T71, Tcy);
						  T72 = FNMS(T70, T71, T6Y);
						  T6P = FMA(T6N, T6O, T6L);
						  TbF = FNMS(T6N, T6K, TbE);
						  T6R = ri[WS(ios, 15)];
					     }
					}
				   }
			      }
			      {
				   E T25, T8Y, T1W, T8U, T1P, Tga, T8R, T8V, T1T, T4P, T4R;
				   {
					E T8y, T16, T1q, T8I, T1z, T1g, T8A, T1m;
					{
					     E T1y, T1u, T1b, T1f;
					     {
						  E T15, Tcv, T6Q, TbG, Tha, Tcw, T6S, T11;
						  T15 = ii[WS(ios, 8)];
						  Tcv = T6H - T6P;
						  T6Q = T6H + T6P;
						  TbG = TbD - TbF;
						  Tha = TbD + TbF;
						  Tcw = TP * T6R;
						  T6S = TL * T6R;
						  T11 = ri[WS(ios, 8)];
						  T1y = ii[WS(ios, 24)];
						  {
						       E Tcx, T6U, T8x, T12;
						       Tcx = FMA(TL, T6T, Tcw);
						       T6U = FNMS(TP, T6T, T6S);
						       T8x = T14 * T11;
						       T12 = T10 * T11;
						       {
							    E TcA, Thb, TbH, T73;
							    TcA = Tcx - Tcz;
							    Thb = Tcx + Tcz;
							    TbH = T6U - T72;
							    T73 = T6U + T72;
							    T8y = FMA(T10, T15, T8x);
							    T16 = FNMS(T14, T15, T12);
							    TeP = Tcv + TcA;
							    TcB = Tcv - TcA;
							    Tja = Tha + Thb;
							    Thc = Tha - Thb;
							    TeE = TbG - TbH;
							    TbI = TbG + TbH;
							    Thr = T6Q - T73;
							    T74 = T6Q + T73;
							    T1u = ri[WS(ios, 24)];
						       }
						  }
					     }
					     T1b = ri[WS(ios, 40)];
					     T1f = ii[WS(ios, 40)];
					     T1q = ii[WS(ios, 56)];
					     {
						  E T8H, T1v, T1c, T8z;
						  T8H = T1x * T1u;
						  T1v = T1t * T1u;
						  T1c = T1a * T1b;
						  T8z = T1a * T1f;
						  T8I = FMA(T1t, T1y, T8H);
						  T1z = FNMS(T1x, T1y, T1v);
						  T1g = FMA(T1e, T1f, T1c);
						  T8A = FNMS(T1e, T1b, T8z);
						  T1m = ri[WS(ios, 56)];
					     }
					}
					{
					     E T1G, T8O, T1O, T8Q, T1S, T20, T24;
					     {
						  E T1D, T1h, Tg5, T1E, T1F, T8J, Tg6, T8E, T1A, T8N, T8G, T1r;
						  T1D = ri[WS(ios, 4)];
						  {
						       E T8C, T8B, T8F, T1n;
						       T8C = T16 - T1g;
						       T1h = T16 + T1g;
						       T8B = T8y - T8A;
						       Tg5 = T8y + T8A;
						       T8F = T1p * T1m;
						       T1n = T1l * T1m;
						       TdT = T8C + T8B;
						       T8D = T8B - T8C;
						       T8G = FMA(T1l, T1q, T8F);
						       T1r = FNMS(T1p, T1q, T1n);
						       T1E = T9 * T1D;
						  }
						  T1F = ii[WS(ios, 4)];
						  T8J = T8G - T8I;
						  Tg6 = T8G + T8I;
						  T8E = T1r - T1z;
						  T1A = T1r + T1z;
						  T8N = T9 * T1F;
						  Tg7 = Tg5 - Tg6;
						  Tk7 = Tg5 + Tg6;
						  T8K = T8E + T8J;
						  TdU = T8E - T8J;
						  T1B = T1h + T1A;
						  TkD = T1A - T1h;
						  T1G = FMA(Td, T1F, T1E);
						  T8O = FNMS(Td, T1D, T8N);
					     }
					     T20 = ri[WS(ios, 52)];
					     T24 = ii[WS(ios, 52)];
					     {
						  E T1J, T1N, T21, T8X, T1K, T8P;
						  T1J = ri[WS(ios, 36)];
						  T1N = ii[WS(ios, 36)];
						  T21 = T1Z * T20;
						  T8X = T1Z * T24;
						  T1K = T1I * T1J;
						  T8P = T1I * T1N;
						  T25 = FMA(T23, T24, T21);
						  T8Y = FNMS(T23, T20, T8X);
						  T1O = FMA(T1M, T1N, T1K);
						  T8Q = FNMS(T1M, T1J, T8P);
					     }
					     T1W = ii[WS(ios, 20)];
					     T1S = ri[WS(ios, 20)];
					     T8U = T1G - T1O;
					     T1P = T1G + T1O;
					     Tga = T8O + T8Q;
					     T8R = T8O - T8Q;
					     T8V = T1V * T1S;
					     T1T = T1R * T1S;
					}
				   }
				   {
					E T9T, T3T, T45, T4i, Taj, T3Z, T9V, T43, Tag;
					{
					     E T4d, T4e, T4h, T3V, T3Y;
					     {
						  E T3S, T8W, T1X, T3O, T8Z, Tgb;
						  T3S = ii[WS(ios, 62)];
						  T8W = FMA(T1R, T1W, T8V);
						  T1X = FNMS(T1V, T1W, T1T);
						  T3O = ri[WS(ios, 62)];
						  T4d = ri[WS(ios, 46)];
						  T8Z = T8W - T8Y;
						  Tgb = T8W + T8Y;
						  {
						       E T8S, T26, T9S, T3P;
						       T8S = T1X - T25;
						       T26 = T1X + T25;
						       T9S = T3R * T3O;
						       T3P = T3N * T3O;
						       TdX = T8U + T8Z;
						       T90 = T8U - T8Z;
						       TiJ = Tga + Tgb;
						       Tgc = Tga - Tgb;
						       TdY = T8R - T8S;
						       T8T = T8R + T8S;
						       Tg9 = T1P - T26;
						       T27 = T1P + T26;
						       T9T = FMA(T3N, T3S, T9S);
						       T3T = FNMS(T3R, T3S, T3P);
						       T4e = T4c * T4d;
						  }
						  T4h = ii[WS(ios, 46)];
					     }
					     T3V = ri[WS(ios, 30)];
					     T3Y = ii[WS(ios, 30)];
					     T45 = ii[WS(ios, 14)];
					     T4i = FMA(T4g, T4h, T4e);
					     {
						  E Tai, T3W, T9U, T42;
						  Tai = T4c * T4h;
						  T3W = T3U * T3V;
						  T9U = T3U * T3Y;
						  T42 = ri[WS(ios, 14)];
						  Taj = FNMS(T4g, T4d, Tai);
						  T3Z = FMA(T3X, T3Y, T3W);
						  T9V = FNMS(T3X, T3V, T9U);
						  T43 = T41 * T42;
						  Tag = T44 * T42;
					     }
					}
					{
					     E Taf, T40, Tgw, T9W, T46, Tah;
					     Taf = T3T - T3Z;
					     T40 = T3T + T3Z;
					     Tgw = T9T + T9V;
					     T9W = T9T - T9V;
					     T46 = FNMS(T44, T45, T43);
					     Tah = FMA(T41, T45, Tag);
					     {
						  E T4j, T9X, Tak, Tgx;
						  T4j = T46 + T4i;
						  T9X = T46 - T4i;
						  Tak = Tah - Taj;
						  Tgx = Tah + Taj;
						  T9Y = T9W + T9X;
						  Tec = T9W - T9X;
						  T4k = T40 + T4j;
						  TgB = T40 - T4j;
						  Tgy = Tgw - Tgx;
						  TiT = Tgw + Tgx;
						  Tef = Taf + Tak;
						  Tal = Taf - Tak;
					     }
					}
				   }
				   T4P = ri[WS(ios, 1)];
				   T4R = ii[WS(ios, 1)];
				   {
					E T56, T5a, T4Q, Tat, T57, Tbp;
					T56 = ri[WS(ios, 49)];
					T5a = ii[WS(ios, 49)];
					T4Q = T5 * T4P;
					Tat = T5 * T4R;
					T57 = T55 * T56;
					Tbp = T55 * T5a;
					T4S = FMA(T8, T4R, T4Q);
					Tau = FNMS(T8, T4P, Tat);
					T5b = FMA(T59, T5a, T57);
					Tbq = FNMS(T59, T56, Tbp);
				   }
				   T4X = ii[WS(ios, 33)];
				   T4U = ri[WS(ios, 33)];
				   T52 = ii[WS(ios, 17)];
				   T50 = ri[WS(ios, 17)];
			      }
			 }
			 {
			      E T7i, TbU, T7t, T7f, TbJ, TbO, Ths, T7n, TbR;
			      {
				   E Ta6, Tab, Tan, Ta5;
				   {
					E T9x, T9C, T9O, T9w;
					{
					     E T3f, T9s, T9B, T3I, T3x, T3B, T3t, T9u, T3y, T9y, T9q, T3u;
					     {
						  E T3b, T3c, T3e, T3H, T3E, T3s, T3p, T9t;
						  T3b = ri[WS(ios, 10)];
						  {
						       E Tav, T4V, Tbn, T51;
						       Tav = T4W * T4U;
						       T4V = T4T * T4U;
						       Tbn = T4b * T50;
						       T51 = T48 * T50;
						       {
							    E Taw, T4Y, Tbo, T53;
							    Taw = FMA(T4T, T4X, Tav);
							    T4Y = FNMS(T4W, T4X, T4V);
							    Tbo = FMA(T48, T52, Tbn);
							    T53 = FNMS(T4b, T52, T51);
							    {
								 E Tax, TgJ, Tbm, T4Z;
								 Tax = Tau - Taw;
								 TgJ = Tau + Taw;
								 Tbm = T4S - T4Y;
								 T4Z = T4S + T4Y;
								 {
								      E Tbr, TgK, Tay, T5c;
								      Tbr = Tbo - Tbq;
								      TgK = Tbo + Tbq;
								      Tay = T53 - T5b;
								      T5c = T53 + T5b;
								      Tew = Tbm + Tbr;
								      Tbs = Tbm - Tbr;
								      TiZ = TgJ + TgK;
								      TgL = TgJ - TgK;
								      Tel = Tax - Tay;
								      Taz = Tax + Tay;
								      Th0 = T4Z - T5c;
								      T5d = T4Z + T5c;
								      T3c = T3a * T3b;
								 }
							    }
						       }
						  }
						  T3e = ii[WS(ios, 10)];
						  T3H = ii[WS(ios, 26)];
						  T3E = ri[WS(ios, 26)];
						  T3s = ii[WS(ios, 42)];
						  T3f = FMA(T3d, T3e, T3c);
						  {
						       E T9r, T9A, T3F, T3o;
						       T9r = T3a * T3e;
						       T9A = T3G * T3E;
						       T3F = T3D * T3E;
						       T3o = ri[WS(ios, 42)];
						       T9s = FNMS(T3d, T3b, T9r);
						       T9B = FMA(T3D, T3H, T9A);
						       T3I = FNMS(T3G, T3H, T3F);
						       T3p = T3n * T3o;
						       T9t = T3r * T3o;
						  }
						  T3x = ri[WS(ios, 58)];
						  T3B = ii[WS(ios, 58)];
						  T3t = FNMS(T3r, T3s, T3p);
						  T9u = FMA(T3n, T3s, T9t);
						  T3y = T3w * T3x;
						  T9y = T3w * T3B;
					     }
					     T9q = T3f - T3t;
					     T3u = T3f + T3t;
					     {
						  E Tgr, T9v, T3C, T9z, T3J, Tgs;
						  Tgr = T9s + T9u;
						  T9v = T9s - T9u;
						  T3C = FMA(T3A, T3B, T3y);
						  T9z = FNMS(T3A, T3x, T9y);
						  T9x = T3C - T3I;
						  T3J = T3C + T3I;
						  Tgs = T9z + T9B;
						  T9C = T9z - T9B;
						  T9O = T9v - T9q;
						  T9w = T9q + T9v;
						  Tgo = T3J - T3u;
						  T3K = T3u + T3J;
						  TiO = Tgr + Tgs;
						  Tgt = Tgr - Tgs;
					     }
					}
					{
					     E Ta1, T4o, T4C, T4J, Taa, Ta3, T4u, T4z, Ta7;
					     {
						  E T4F, T4G, T4I, T4n, T4l, T4t, T4q;
						  T4n = ii[WS(ios, 6)];
						  T4l = ri[WS(ios, 6)];
						  T4F = ri[WS(ios, 22)];
						  {
						       E T9N, T9D, Ta0, T4m;
						       T9N = T9x + T9C;
						       T9D = T9x - T9C;
						       Ta0 = T3i * T4l;
						       T4m = T3g * T4l;
						       Te6 = T9O + T9N;
						       T9P = T9N - T9O;
						       Te9 = T9w + T9D;
						       T9E = T9w - T9D;
						       Ta1 = FMA(T3g, T4n, Ta0);
						       T4o = FNMS(T3i, T4n, T4m);
						       T4G = T4E * T4F;
						  }
						  T4I = ii[WS(ios, 22)];
						  T4t = ii[WS(ios, 38)];
						  T4q = ri[WS(ios, 38)];
						  T4C = ii[WS(ios, 54)];
						  T4J = FMA(T4H, T4I, T4G);