t2_64.c

快速fft变换

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