hf2_64.c

最新的FFT程序

						       T12 = T10 * T11;
						       TeE = TbC + TbH;
						       TbI = TbC - TbH;
						       Tjb = Thr + Ths;
						       Tht = Thr - Ths;
						       TeP = TcA + Tcv;
						       TcB = Tcv - TcA;
						       Tha = T6Q - T73;
						       T74 = T6Q + T73;
						       T8y = FMA(T10, T15, T8x);
						       T16 = FNMS(T14, T15, T12);
						       T1u = rio[WS(ios, 24)];
						  }
					     }
					     T1b = rio[WS(ios, 40)];
					     T1f = iio[-WS(ios, 23)];
					     T1q = iio[-WS(ios, 7)];
					     {
						  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 = rio[WS(ios, 56)];
					     }
					}
					{
					     E T93, T2f, T2n, T95, T2q, T2r;
					     {
						  E T2D, T2z, T2i, T2m;
						  {
						       E T2e, T1h, Tg5, T8G, T1r, T2a, T8J, Tg6;
						       T2e = iio[-WS(ios, 3)];
						       {
							    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);
							    T2a = rio[WS(ios, 60)];
						       }
						       T2D = iio[-WS(ios, 19)];
						       T8J = T8G - T8I;
						       Tg6 = T8G + T8I;
						       {
							    E T8E, T1A, T92, T2b;
							    T8E = T1r - T1z;
							    T1A = T1r + T1z;
							    T92 = T2d * T2a;
							    T2b = T29 * T2a;
							    Tk7 = Tg5 + Tg6;
							    Tg7 = Tg5 - Tg6;
							    TdU = T8E - T8J;
							    T8K = T8E + T8J;
							    TkD = T1A - T1h;
							    T1B = T1h + T1A;
							    T93 = FMA(T29, T2e, T92);
							    T2f = FNMS(T2d, T2e, T2b);
							    T2z = rio[WS(ios, 44)];
						       }
						  }
						  T2i = rio[WS(ios, 28)];
						  T2m = iio[-WS(ios, 35)];
						  T2p = rio[WS(ios, 12)];
						  {
						       E T9c, T2A, T2j, T94;
						       T9c = T2C * T2z;
						       T2A = T2y * T2z;
						       T2j = T2h * T2i;
						       T94 = T2h * T2m;
						       T9d = FMA(T2y, T2D, T9c);
						       T2E = FNMS(T2C, T2D, T2A);
						       T2n = FMA(T2l, T2m, T2j);
						       T95 = FNMS(T2l, T2i, T94);
						       T2q = TH * T2p;
						  }
						  T2r = iio[-WS(ios, 51)];
					     }
					     T2o = T2f + T2n;
					     T99 = T2f - T2n;
					     T96 = T93 - T95;
					     Tgf = T93 + T95;
					     T2s = FMA(TK, T2r, T2q);
					     T9a = TH * T2r;
					}
				   }
			      }
			      {
				   E T4S, Tbn, T5b, T52, Tax, T50, Tbp, T4Y;
				   {
					E T4P, T56, Tbm, T57, T5a;
					{
					     E T4Q, Tgg, T9e, T4R, T2F, T97, T9b;
					     T4P = rio[WS(ios, 1)];
					     T2F = T2s + T2E;
					     T97 = T2s - T2E;
					     T9b = FNMS(TK, T2p, T9a);
					     T4Q = T5 * T4P;
					     T98 = T96 + T97;
					     Te1 = T96 - T97;
					     T2G = T2o + T2F;
					     Tge = T2o - T2F;
					     Tgg = T9b + T9d;
					     T9e = T9b - T9d;
					     T4R = iio[-WS(ios, 62)];
					     T56 = rio[WS(ios, 49)];
					     T9f = T99 - T9e;
					     Te0 = T99 + T9e;
					     Tgh = Tgf - Tgg;
					     TiK = Tgf + Tgg;
					     T4S = FMA(T8, T4R, T4Q);
					     Tbm = T5 * T4R;
					     T57 = T55 * T56;
					     T5a = iio[-WS(ios, 14)];
					}
					{
					     E T4X, T4U, Taw, Tbo, T4V;
					     T4X = iio[-WS(ios, 30)];
					     Tbn = FNMS(T8, T4P, Tbm);
					     T4U = rio[WS(ios, 33)];
					     T5b = FMA(T59, T5a, T57);
					     Taw = T55 * T5a;
					     T52 = iio[-WS(ios, 46)];
					     Tbo = T4W * T4U;
					     T4V = T4T * T4U;
					     Tax = FNMS(T59, T56, Taw);
					     T50 = rio[WS(ios, 17)];
					     Tbp = FMA(T4T, T4X, Tbo);
					     T4Y = FNMS(T4W, T4X, T4V);
					}
				   }
				   {
					E Tag, T3T, T45, T4i, T9W, T3Z, Tai, T43, T9T;
					{
					     E T4d, T4e, T4h, T3V, T3Y;
					     {
						  E T3S, Tbq, Th0, Tat, T4Z, Tav, T53, T3O, Tau, T51, Tay, Th1;
						  T3S = iio[-WS(ios, 1)];
						  Tau = T4b * T50;
						  T51 = T48 * T50;
						  Tbq = Tbn - Tbp;
						  Th0 = Tbn + Tbp;
						  Tat = T4S - T4Y;
						  T4Z = T4S + T4Y;
						  Tav = FMA(T48, T52, Tau);
						  T53 = FNMS(T4b, T52, T51);
						  T3O = rio[WS(ios, 62)];
						  T4d = rio[WS(ios, 46)];
						  Tay = Tav - Tax;
						  Th1 = Tav + Tax;
						  {
						       E Tbr, T5c, Taf, T3P;
						       Tbr = T53 - T5b;
						       T5c = T53 + T5b;
						       Taf = T3R * T3O;
						       T3P = T3N * T3O;
						       Tel = Tat + Tay;
						       Taz = Tat - Tay;
						       TiZ = Th0 + Th1;
						       Th2 = Th0 - Th1;
						       Tew = Tbq - Tbr;
						       Tbs = Tbq + Tbr;
						       TgJ = T4Z - T5c;
						       T5d = T4Z + T5c;
						       Tag = FMA(T3N, T3S, Taf);
						       T3T = FNMS(T3R, T3S, T3P);
						       T4e = T4c * T4d;
						  }
						  T4h = iio[-WS(ios, 17)];
					     }
					     T3V = rio[WS(ios, 30)];
					     T3Y = iio[-WS(ios, 33)];
					     T45 = iio[-WS(ios, 49)];
					     T4i = FMA(T4g, T4h, T4e);
					     {
						  E T9V, T3W, Tah, T42;
						  T9V = T4c * T4h;
						  T3W = T3U * T3V;
						  Tah = T3U * T3Y;
						  T42 = rio[WS(ios, 14)];
						  T9W = FNMS(T4g, T4d, T9V);
						  T3Z = FMA(T3X, T3Y, T3W);
						  Tai = FNMS(T3X, T3V, Tah);
						  T43 = T41 * T42;
						  T9T = T44 * T42;
					     }
					}
					{
					     E T9S, T40, TgB, Taj, T46, T9U;
					     T9S = T3T - T3Z;
					     T40 = T3T + T3Z;
					     TgB = Tag + Tai;
					     Taj = Tag - Tai;
					     T46 = FNMS(T44, T45, T43);
					     T9U = FMA(T41, T45, T9T);
					     {
						  E Tak, T4j, T9X, TgC;
						  Tak = T46 - T4i;
						  T4j = T46 + T4i;
						  T9X = T9U - T9W;
						  TgC = T9U + T9W;
						  T4k = T40 + T4j;
						  Tgw = T40 - T4j;
						  T9Y = T9S - T9X;
						  Tec = T9S + T9X;
						  TgD = TgB - TgC;
						  TiT = TgB + TgC;
						  Tef = Taj - Tak;
						  Tal = Taj + Tak;
					     }
					}
				   }
			      }
			      {
				   E T2O, T9H, T2M, T9n, T37, T2P, T2R, T2Y, T2V;
				   {
					E T2L, T2J, T36, T32;
					T2L = iio[-WS(ios, 61)];
					T2J = rio[WS(ios, 2)];
					T36 = iio[-WS(ios, 13)];
					T32 = rio[WS(ios, 50)];
					T2O = rio[WS(ios, 34)];
					{
					     E T9G, T2K, T9m, T33;
					     T9G = Tt * T2J;
					     T2K = Tr * T2J;
					     T9m = T35 * T32;
					     T33 = T31 * T32;
					     T9H = FMA(Tr, T2L, T9G);
					     T2M = FNMS(Tt, T2L, T2K);
					     T9n = FMA(T31, T36, T9m);
					     T37 = FNMS(T35, T36, T33);
					     T2P = T2N * T2O;
					}
					T2R = iio[-WS(ios, 29)];
					T2Y = iio[-WS(ios, 45)];
					T2V = rio[WS(ios, 18)];
				   }
				   {
					E T1D, T1E, T1F, T20, T24, T8N, T21, T8X;
					{
					     E T2T, T9j, T9J, T9l, T2Z, Tgq, T9K;
					     T1D = rio[WS(ios, 4)];
					     {
						  E T2S, T9I, T9k, T2W;
						  T2S = FMA(T2Q, T2R, T2P);
						  T9I = T2N * T2R;
						  T9k = T2X * T2V;
						  T2W = T2U * T2V;
						  T2T = T2M + T2S;
						  T9j = T2M - T2S;
						  T9J = FNMS(T2Q, T2O, T9I);
						  T9l = FMA(T2U, T2Y, T9k);
						  T2Z = FNMS(T2X, T2Y, T2W);
						  T1E = T9 * T1D;
					     }
					     Tgq = T9H + T9J;
					     T9K = T9H - T9J;
					     {
						  E Tgr, T9o, T38, T9L;
						  Tgr = T9l + T9n;
						  T9o = T9l - T9n;
						  T38 = T2Z + T37;
						  T9L = T2Z - T37;
						  T9p = T9j - T9o;
						  Te5 = T9j + T9o;
						  Tgs = Tgq - Tgr;
						  TiN = Tgq + Tgr;
						  T9M = T9K + T9L;
						  Te8 = T9K - T9L;
						  T39 = T2T + T38;
						  Tgl = T2T - T38;
						  T1F = iio[-WS(ios, 59)];
					     }
					}
					T20 = rio[WS(ios, 52)];
					T24 = iio[-WS(ios, 11)];
					T1J = rio[WS(ios, 36)];
					T1G = FMA(Td, T1F, T1E);
					T8N = T9 * T1F;
					T21 = T1Z * T20;
					T8X = T1Z * T24;
					T1K = T1I * T1J;
					T8O = FNMS(Td, T1D, T8N);
					T25 = FMA(T23, T24, T21);
					T8Y = FNMS(T23, T20, T8X);
					T1N = iio[-WS(ios, 27)];
					T1W = iio[-WS(ios, 43)];
					T1S = rio[WS(ios, 20)];
				   }
			      }
			 }
			 {
			      E T4t, Ta7, T4o, T4J, Ta3, T4q, T4C, T4y;
			      {
				   E T78, T7a, TbR, T7b, TbN, T7t, T7d, T7i, T7m, T9N, T9O, Tbt, Tbu;
				   {
					E T3f, T9y, T9u, T3I, T3x, T3B, T3t, T9A, T3y, T9r, T9C, T3u;
					{
					     E T3b, T3c, T3e, T3H, T3E, T3s, T3p, T9z;
					     {
						  E T1P, T8U, T8Q, T8W, T1X, Tga, T8R;
						  T3b = rio[WS(ios, 10)];
						  {
						       E T1O, T8P, T8V, T1T;
						       T1O = FMA(T1M, T1N, T1K);
						       T8P = T1I * T1N;
						       T8V = T1V * T1S;
						       T1T = T1R * T1S;
						       T1P = T1G + T1O;
						       T8U = T1G - T1O;
						       T8Q = FNMS(T1M, T1J, T8P);
						       T8W = FMA(T1R, T1W, T8V);
						       T1X = FNMS(T1V, T1W, T1T);
						       T3c = T3a * T3b;
						  }
						  Tga = T8O + T8Q;
						  T8R = T8O - T8Q;
						  {
						       E Tgb, T8Z, T26, T8S;
						       Tgb = T8W + T8Y;
						       T8Z = T8W - T8Y;
						       T26 = T1X + T25;
						       T8S = T1X - T25;
						       T90 = T8U - T8Z;
						       TdX = T8U + T8Z;
						       Tgc = Tga - Tgb;
						       TiJ = Tga + Tgb;
						       T8T = T8R + T8S;
						       TdY = T8R - T8S;
						       T27 = T1P + T26;
						       Tg9 = T1P - T26;
						       T3e = iio[-WS(ios, 53)];
						  }
					     }
					     T3H = iio[-WS(ios, 37)];
					     T3E = rio[WS(ios, 26)];
					     T3s = iio[-WS(ios, 21)];
					     T3f = FMA(T3d, T3e, T3c);
					     {
						  E T9x, T9t, T3F, T3o;
						  T9x = T3a * T3e;
						  T9t = T3G * T3E;
						  T3F = T3D * T3E;
						  T3o = rio[WS(ios, 42)];
						  T9y = FNMS(T3d, T3b, T9x);
						  T9u = FMA(T3D, T3H, T9t);
						  T3I = FNMS(T3G, T3H, T3F);
						  T3p = T3n * T3o;
						  T9z = T3r * T3o;
					     }
					     T3x = rio[WS(ios, 58)];
					     T3B = iio[-WS(ios, 5)];
					     T3t = FNMS(T3r, T3s, T3p);
					     T9A = FMA(T3n, T3s, T9z);
					     T3y = T3w * T3x;
					     T9r = T3w * T3B;
					}
					T9C = T3f - T3t;
					T3u = T3f + T3t;
					{
					     E Tgm, T9B, T3C, T9s, T9D, T9w;
					     Tgm = T9y + T9A;
					     T9B = T9y - T9A;
					     T3C = FMA(T3A, T3B, T3y);
					     T9s = FNMS(T3A, T3x, T9r);
					     {
						  E T3J, T9q, Tgn, T9v;
						  T3J = T3C + T3I;
						  T9q = T3C - T3I;
						  Tgn = T9s + T9u;
						  T9v = T9s - T9u;
						  T9N = T9C + T9B;
						  T9D = T9B - T9C;
						  Tgt = T3J - T3u;
						  T3K = T3u + T3J;
						  T9O = T9q - T9v;
						  T9w = T9q + T9v;
						  TiO = Tgm + Tgn;
						  Tgo = Tgm - Tgn;
					     }
					     Te9 = T9D + T9w;
					     T9E = T9w - T9D;
					}
				   }
				   {
					E T75, T7s, T7p, TbQ, T77, T76, TbM, T7q;
					T75 = rio[WS(ios, 7)];
					Te6 = T9N + T9O;
					T9P = T9N - T9O;
					T77 = iio[-WS(ios, 56)];
					T76 = T1i * T75;
					T7s = iio[-WS(ios, 40)];
					T7p = rio[WS(ios, 23)];
					TbQ = T1i * T77;
					T78 = FMA(T1k, T77, T76);
					T7a = rio[WS(ios, 39)];
					TbM = T7r * T7p;
					T7q = T7o * T7p;
					TbR = FNMS(T1k, T75, TbQ);
					T7b = T79 * T7a;
					TbN = FMA(T7o, T7s, TbM);
					T7t = FNMS(T7r, T7s, T7q);
					T7d = iio[-WS(ios, 24)];
					T7i = rio[WS(ios, 55)];
					T7m = iio[-WS(ios, 8)];
				   }
				   {
					E T5h, TaI, T5G, TaE, T5t, T5x, T5p, TaK, T5u, TaB, TaM, T5q;
					{
					     E T5e, T5f, T5g, T5B, T5F, T5o, T5l, TaJ;