hb_64.c

快速fft变换

					TjD = Tja + Tj9;
					Tjb = Tj9 - Tja;
					TeB = FNMS(KP414213562, Tdb, Tdg);
					Tdh = FMA(KP414213562, Tdg, Tdb);
					T90 = T2S + T2X;
					T2Y = T2S - T2X;
					Tgl = FMA(KP414213562, TfR, TfS);
					TfT = FNMS(KP414213562, TfS, TfR);
					Tgk = FNMS(KP414213562, TfU, TfV);
					TfW = FMA(KP414213562, TfV, TfU);
					T2Z = FMA(KP707106781, T2Y, T2N);
					T6X = FNMS(KP707106781, T2Y, T2N);
					T37 = T35 + T36;
					T8X = T35 - T36;
					{
					     E T8Q, T8T, T3p, T6Y, T38;
					     T3y = T3w + T3x;
					     T8Q = T3x - T3w;
					     T8T = T3j + T3o;
					     T3p = T3j - T3o;
					     T6Y = FNMS(KP707106781, T37, T34);
					     T38 = FMA(KP707106781, T37, T34);
					     {
						  E Tb9, T8R, Tba, T8U;
						  Tb9 = FMA(KP707106781, T8Q, T8P);
						  T8R = FNMS(KP707106781, T8Q, T8P);
						  Tba = FMA(KP707106781, T8T, T8S);
						  T8U = FNMS(KP707106781, T8T, T8S);
						  T6Z = FMA(KP668178637, T6Y, T6X);
						  T7r = FNMS(KP668178637, T6X, T6Y);
						  T5H = FMA(KP198912367, T2Z, T38);
						  T39 = FNMS(KP198912367, T38, T2Z);
						  Tbb = FNMS(KP198912367, Tba, Tb9);
						  TbC = FMA(KP198912367, Tb9, Tba);
						  T9S = FNMS(KP668178637, T8R, T8U);
						  T8V = FMA(KP668178637, T8U, T8R);
						  T3q = FMA(KP707106781, T3p, T3e);
						  T70 = FNMS(KP707106781, T3p, T3e);
					     }
					}
				   }
			      }
			 }
		    }
		    {
			 E T97, Tbk, T9j, T9k, Tbh, T9i, Tbi, T9e;
			 {
			      E T9g, T5f, T18, Ted, TdY, ThI, T4A, T95, T9b, T57, T1u, Te1, Te4, ThM, T52;
			      E T9c, T5h, T4K, TdZ, T1f, ThJ, Teg, T5g, T4F, T1j, Te8, T98, T4W, T4N, T1m;
			      E Te7, T4Q, T1n, Te6;
			      {
				   E T1q, Te3, T4Y, T1t, Te2, T51;
				   {
					E T15, T5b, T14, TdX, T5e, T16, T4x, T4y;
					{
					     E T12, T13, T5c, T5d, T71, T3z;
					     T12 = cr[WS(rs, 1)];
					     T71 = FNMS(KP707106781, T3y, T3v);
					     T3z = FMA(KP707106781, T3y, T3v);
					     {
						  E Tbc, T8Y, Tbd, T91;
						  Tbc = FMA(KP707106781, T8X, T8W);
						  T8Y = FNMS(KP707106781, T8X, T8W);
						  Tbd = FMA(KP707106781, T90, T8Z);
						  T91 = FNMS(KP707106781, T90, T8Z);
						  T72 = FNMS(KP668178637, T71, T70);
						  T7q = FMA(KP668178637, T70, T71);
						  T5G = FNMS(KP198912367, T3q, T3z);
						  T3A = FMA(KP198912367, T3z, T3q);
						  Tbe = FNMS(KP198912367, Tbd, Tbc);
						  TbD = FMA(KP198912367, Tbc, Tbd);
						  T9T = FNMS(KP668178637, T8Y, T91);
						  T92 = FMA(KP668178637, T91, T8Y);
						  T13 = ci[WS(rs, 30)];
					     }
					     T5c = ci[WS(rs, 46)];
					     T5d = cr[WS(rs, 49)];
					     T15 = cr[WS(rs, 17)];
					     T5b = T12 - T13;
					     T14 = T12 + T13;
					     TdX = T5c - T5d;
					     T5e = T5c + T5d;
					     T16 = ci[WS(rs, 14)];
					     T4x = ci[WS(rs, 62)];
					     T4y = cr[WS(rs, 33)];
					}
					{
					     E T4w, T17, TdW, T4z;
					     T9g = T5b + T5e;
					     T5f = T5b - T5e;
					     T4w = T15 - T16;
					     T17 = T15 + T16;
					     TdW = T4x - T4y;
					     T4z = T4x + T4y;
					     T18 = T14 + T17;
					     Ted = T14 - T17;
					     TdY = TdW - TdX;
					     ThI = TdW + TdX;
					     T4A = T4w + T4z;
					     T95 = T4z - T4w;
					}
				   }
				   {
					E T1r, T53, T56, T1s, T4Z, T50;
					{
					     E T1o, T1p, T54, T55;
					     T1o = ci[WS(rs, 2)];
					     T1p = cr[WS(rs, 29)];
					     T54 = ci[WS(rs, 50)];
					     T55 = cr[WS(rs, 45)];
					     T1r = cr[WS(rs, 13)];
					     T53 = T1o - T1p;
					     T1q = T1o + T1p;
					     Te3 = T54 - T55;
					     T56 = T54 + T55;
					     T1s = ci[WS(rs, 18)];
					     T4Z = ci[WS(rs, 34)];
					     T50 = cr[WS(rs, 61)];
					}
					T9b = T53 + T56;
					T57 = T53 - T56;
					T4Y = T1r - T1s;
					T1t = T1r + T1s;
					Te2 = T4Z - T50;
					T51 = T4Z + T50;
				   }
				   T1u = T1q + T1t;
				   Te1 = T1q - T1t;
				   Te4 = Te2 - Te3;
				   ThM = Te2 + Te3;
				   T52 = T4Y - T51;
				   T9c = T4Y + T51;
				   {
					E T1c, T4B, T1b, Tee, T4J, T1d, T4C, T4D;
					{
					     E T19, T1a, T4H, T4I;
					     T19 = cr[WS(rs, 9)];
					     T1a = ci[WS(rs, 22)];
					     T4H = ci[WS(rs, 38)];
					     T4I = cr[WS(rs, 57)];
					     T1c = ci[WS(rs, 6)];
					     T4B = T19 - T1a;
					     T1b = T19 + T1a;
					     Tee = T4H - T4I;
					     T4J = T4H + T4I;
					     T1d = cr[WS(rs, 25)];
					     T4C = ci[WS(rs, 54)];
					     T4D = cr[WS(rs, 41)];
					}
					{
					     E T1k, T4S, T4V, T1l, T4O, T4P;
					     {
						  E T1h, T1i, T4T, T4U;
						  T1h = cr[WS(rs, 5)];
						  {
						       E T4G, T1e, Tef, T4E;
						       T4G = T1c - T1d;
						       T1e = T1c + T1d;
						       Tef = T4C - T4D;
						       T4E = T4C + T4D;
						       T5h = T4G - T4J;
						       T4K = T4G + T4J;
						       TdZ = T1b - T1e;
						       T1f = T1b + T1e;
						       ThJ = Tef + Tee;
						       Teg = Tee - Tef;
						       T5g = T4B - T4E;
						       T4F = T4B + T4E;
						       T1i = ci[WS(rs, 26)];
						  }
						  T4T = ci[WS(rs, 42)];
						  T4U = cr[WS(rs, 53)];
						  T1k = cr[WS(rs, 21)];
						  T4S = T1h - T1i;
						  T1j = T1h + T1i;
						  Te8 = T4T - T4U;
						  T4V = T4T + T4U;
						  T1l = ci[WS(rs, 10)];
						  T4O = ci[WS(rs, 58)];
						  T4P = cr[WS(rs, 37)];
					     }
					     T98 = T4S + T4V;
					     T4W = T4S - T4V;
					     T4N = T1k - T1l;
					     T1m = T1k + T1l;
					     Te7 = T4O - T4P;
					     T4Q = T4O + T4P;
					}
				   }
			      }
			      T1n = T1j + T1m;
			      Te6 = T1j - T1m;
			      {
				   E Te9, ThL, T4R, T99;
				   Te9 = Te7 - Te8;
				   ThL = Te7 + Te8;
				   T4R = T4N + T4Q;
				   T99 = T4Q - T4N;
				   {
					E Tjr, ThK, Tjs, ThN;
					{
					     E T1g, T1v, Tjp, Tjo;
					     Tjr = T18 - T1f;
					     T1g = T18 + T1f;
					     T1v = T1n + T1u;
					     Tjp = T1n - T1u;
					     ThK = ThI + ThJ;
					     Tjo = ThI - ThJ;
					     ThH = T1g - T1v;
					     T1w = T1g + T1v;
					     Tke = Tjp + Tjo;
					     Tjq = Tjo - Tjp;
					     Tjs = ThM - ThL;
					     ThN = ThL + ThM;
					}
					{
					     E Tg6, Te0, Tg9, Teh, Tej, Tei, Tga, Teb, Te5, Tea;
					     Tg6 = TdZ + TdY;
					     Te0 = TdY - TdZ;
					     Tkf = Tjr + Tjs;
					     Tjt = Tjr - Tjs;
					     TiK = ThK + ThN;
					     ThO = ThK - ThN;
					     Tg9 = Ted + Teg;
					     Teh = Ted - Teg;
					     Tej = Te4 - Te1;
					     Te5 = Te1 + Te4;
					     Tea = Te6 - Te9;
					     Tei = Te6 + Te9;
					     Tga = Tea + Te5;
					     Teb = Te5 - Tea;
					     {
						  E T9h, T4M, T78, T96, T5k, T5l, T75, T5j, T76, T59;
						  {
						       E T5i, Tg7, Tek, T4L, T4X, T58;
						       T9h = T4F + T4K;
						       T4L = T4F - T4K;
						       Tgb = FNMS(KP707106781, Tga, Tg9);
						       TgT = FMA(KP707106781, Tga, Tg9);
						       Tfc = FMA(KP707106781, Teb, Te0);
						       Tec = FNMS(KP707106781, Teb, Te0);
						       Tg7 = Tei + Tej;
						       Tek = Tei - Tej;
						       T4M = FMA(KP707106781, T4L, T4A);
						       T78 = FNMS(KP707106781, T4L, T4A);
						       Tg8 = FNMS(KP707106781, Tg7, Tg6);
						       TgU = FMA(KP707106781, Tg7, Tg6);
						       Tfd = FMA(KP707106781, Tek, Teh);
						       Tel = FNMS(KP707106781, Tek, Teh);
						       T5i = T5g + T5h;
						       T96 = T5h - T5g;
						       T5k = FNMS(KP414213562, T4R, T4W);
						       T4X = FMA(KP414213562, T4W, T4R);
						       T58 = FNMS(KP414213562, T57, T52);
						       T5l = FMA(KP414213562, T52, T57);
						       T75 = FNMS(KP707106781, T5i, T5f);
						       T5j = FMA(KP707106781, T5i, T5f);
						       T76 = T4X - T58;
						       T59 = T4X + T58;
						  }
						  {
						       E T79, T5m, T9a, T9d;
						       T77 = FNMS(KP923879532, T76, T75);
						       T83 = FMA(KP923879532, T76, T75);
						       T6i = FMA(KP923879532, T59, T4M);
						       T5a = FNMS(KP923879532, T59, T4M);
						       T79 = T5l - T5k;
						       T5m = T5k + T5l;
						       T97 = FNMS(KP707106781, T96, T95);
						       Tbk = FMA(KP707106781, T96, T95);
						       T7a = FNMS(KP923879532, T79, T78);
						       T82 = FMA(KP923879532, T79, T78);
						       T6j = FMA(KP923879532, T5m, T5j);
						       T5n = FNMS(KP923879532, T5m, T5j);
						       T9j = FNMS(KP414213562, T98, T99);
						       T9a = FMA(KP414213562, T99, T98);
						       T9d = FMA(KP414213562, T9c, T9b);
						       T9k = FNMS(KP414213562, T9b, T9c);
						       Tbh = FMA(KP707106781, T9h, T9g);
						       T9i = FNMS(KP707106781, T9h, T9g);
						       Tbi = T9a + T9d;
						       T9e = T9a - T9d;
						  }
					     }
					}
				   }
			      }
			 }
			 {
			      E T9z, T4m, T1D, TdM, ThR, Tdx, T3H, T9o, T9r, T4e, T1Z, TdA, TdD, ThV, T49;
			      E T9s, T4o, T3R, Tdy, T1K, ThS, TdP, T4n, T3M, T1O, T3V, TdH, T3U, T1R, T3W;
			      E T9u, T43;
			      {
				   E T1V, T46, TdC, T45, T1Y, T47, T48, TdB;
				   {
					E Tdw, T3D, T3G, Tdv, T4a, T4d;
					{
					     E T4i, T1z, T3E, T4l, T1C, T3F;
					     {
						  E T4j, T4k, T1A, T1B;
						  {
						       E T1x, Tbl, T9l, T1y;
						       T1x = ci[0];
						       Tbj = FNMS(KP923879532, Tbi, Tbh);
						       Tcc = FMA(KP923879532, Tbi, Tbh);
						       Tas = FMA(KP923879532, T9e, T97);
						       T9f = FNMS(KP923879532, T9e, T97);
						       Tbl = T9j - T9k;
						       T9l = T9j + T9k;
						       T1y = cr[WS(rs, 31)];
						       T4j = ci[WS(rs, 48)];
						       Tbm = FNMS(KP923879532, Tbl, Tbk);
						       Tcb = FMA(KP923879532, Tbl, Tbk);
						       Tar = FNMS(KP923879532, T9l, T9i);
						       T9m = FMA(KP923879532, T9l, T9i);
						       T4i = T1x - T1y;
						       T1z = T1x + T1y;
						       T4k = cr[WS(rs, 47)];
						  }
						  T1A = cr[WS(rs, 15)];
						  T1B = ci[WS(rs, 16)];
						  T3E = ci[WS(rs, 32)];
						  Tdw = T4j - T4k;
						  T4l = T4j + T4k;
						  T3D = T1A - T1B;
						  T1C = T1A + T1B;
						  T3F = cr[WS(rs, 63)];
					     }
					     T9z = T4i + T4l;
					     T4m = T4i - T4l;
					     T1D = T1z + T1C;
					     TdM = T1z - T1C;
					     T3G = T3E + T3F;
					     Tdv = T3E - T3F;
					}
					{
					     E T4b, T4c, T1T, T1U, T1W, T1X;
					     T1T = ci[WS(rs, 4)];
					     T1U = cr[WS(rs, 27)];
					     ThR = Tdv + Tdw;
					     Tdx = Tdv - Tdw;
					     T3H = T3D - T3G;
					     T9o = T3D + T3G;
					     T4a = T1T - T1U;
					     T1V = T1T + T1U;
					     T4b = ci[WS(rs, 52)];
					     T4c = cr[WS(rs, 43)];
					     T1W = cr[WS(rs, 11)];
					     T1X = ci[WS(rs, 20)];
					     T46 = ci[WS(rs, 36)];
					     TdC = T4b - T4c;
					     T4d = T4b + T4c;
					     T45 = T1W - T1X;
					     T1Y = T1W + T1X;
					     T47 = cr[WS(rs, 59)];
					}
					T9r = T4a + T4d;
					T4e = T4a - T4d;
				   }
				   T1Z = T1V + T1Y;
				   TdA = T1V - T1Y;
				   T48 = T46 + T47;
				   TdB = T46 - T47;
				   {
					E T3I, T1G, T3J, TdN, T3Q, T3N, T1J, T3K, T3Z, T42;
					{
					     E T3O, T3P, T1E, T1F, T1H, T1I;
					     T1E = cr[WS(rs, 7)];
					     T1F = ci[WS(rs, 24)];
					     TdD = TdB - TdC;
					     ThV = TdB + TdC;
					     T49 = T45 - T48;
					     T9s = T45 + T48;
					     T3I = T1E - T1F;
					     T1G = T1E + T1F;
					     T3O = ci[WS(rs, 40)];
					     T3P = cr[WS(rs, 55)];
					     T1H = ci[WS(rs, 8)];
					     T1I = cr[WS(rs, 23)];
					     T3J = ci[WS(rs, 56)];
					     TdN = T3O - T3P;
					     T3Q = T3O + T3P;
					     T3N = T1H - T1I;
					     T1J = T1H + T1I;
					     T3K = cr[WS(rs, 39)];
					}
					{
					     E T40, T41, T1P, T1Q;
					     {
						  E T1M, TdO, T3L, T1N;
						  T1M = cr[WS(rs, 3)];
						  T4o = T3N - T3Q;
						  T3R = T3N + T3Q;
						  Tdy = T1G - T1J;
						  T1K = T1G + T1J;
						  TdO = T3J - T3K;
						  T3L = T3J + T3K;
						  T1N = ci[WS(rs, 28)];
						  T40 = ci[WS(rs, 44)];
						  ThS = TdO + TdN;
						  TdP = TdN - TdO;
						  T4n = T3I - T3L;
						  T3M = T3I + T3L;
						  T3Z = T1M - T1N;
						  T1O = T1M + T1N;
						  T41 = cr[WS(rs, 51)];
					     }
					     T1P = cr[WS(rs, 19)];
					     T1Q = ci[WS(rs, 12)];
					     T3V = ci[WS(rs, 60)];
					     TdH = T40 - T41;
					     T42 = T40 + T41;
					     T3U = T1P - T1Q;
					     T1R = T1P + T1Q;
					     T3W = cr[WS(rs, 35)];
					}
					T9u = T3Z + T42;
					T43 = T3Z - T42;
				   }
			      }
			      {
				   E T1S, TdF, T3X, TdG;
				   T1S = T1O + T1R;
				   TdF = T1O - T1R;
				   T3X = T3V + T3W;
				   TdG = T3V - T3W;
				   {
					E TdI, T3Y, T9v, ThT, ThW;
					{
					     E Tjk, Tji, ThU, Tjh, T1L, T20, Tjl;
					     Tjk = T1D - T1K;
					     T1L = T1D + T1K;
					     T20 = T1S + T1Z;
					     Tji = T1S - T1Z;
					     TdI = TdG - TdH;
					     ThU = TdG + TdH;
					     T3Y = T3U + T3X;
					     T9v = T3U - T3X;
					     ThQ = T1L - T20;
					     T21 = T1L + T20;
					     ThT = ThR + ThS;