hf2_64.c

快速傅立叶变换库函数

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