hf_64.c

快速fft变换

					TeU = TeS - TeT;
					Ths = TeS + TeT;
					T7I = T7G + T7H;
					TcG = T7G - T7H;
					T1f = T11 + T1e;
					TeR = T11 - T1e;
					T2C = W[122];
				   }
				   T2F = W[123];
				   T2W = cr[WS(rs, 46)];
				   T2Z = ci[WS(rs, 46)];
				   T8X = T2C * T2G;
				   T2E = T2C * T2D;
				   T2V = W[90];
				   T2Y = W[91];
				   T8Y = FNMS(T2F, T2D, T8X);
				   T2H = FMA(T2F, T2G, T2E);
				   T8D = T2V * T2Z;
				   T2X = T2V * T2W;
			      }
			      {
				   E T2J, T2M, T2I, T2L, T8Z, T2K, T2P;
				   T2J = cr[WS(rs, 30)];
				   T2M = ci[WS(rs, 30)];
				   T8E = FNMS(T2Y, T2W, T8D);
				   T30 = FMA(T2Y, T2Z, T2X);
				   T2I = W[58];
				   T2L = W[59];
				   T2Q = cr[WS(rs, 14)];
				   T2T = ci[WS(rs, 14)];
				   T8Z = T2I * T2M;
				   T2K = T2I * T2J;
				   T2P = W[26];
				   T2S = W[27];
				   T90 = FNMS(T2L, T2J, T8Z);
				   T2N = FMA(T2L, T2M, T2K);
				   T8B = T2P * T2T;
				   T2R = T2P * T2Q;
			      }
			 }
			 {
			      E T91, Tfe, T2O, T8A, T8C, T2U;
			      T91 = T8Y - T90;
			      Tfe = T8Y + T90;
			      T2O = T2H + T2N;
			      T8A = T2H - T2N;
			      T8C = FNMS(T2S, T2Q, T8B);
			      T2U = FMA(T2S, T2T, T2R);
			      {
				   E Tff, T8F, T92, T31;
				   Tff = T8C + T8E;
				   T8F = T8C - T8E;
				   T92 = T2U - T30;
				   T31 = T2U + T30;
				   T8G = T8A - T8F;
				   TcU = T8A + T8F;
				   Tfg = Tfe - Tff;
				   ThB = Tfe + Tff;
				   T32 = T2O + T31;
				   Tfj = T2O - T31;
				   TcX = T91 - T92;
				   T93 = T91 + T92;
			      }
			 }
		    }
		    {
			 E Ta5, T3C, T9f, T3V, T3L, T3O, T3N, Ta7, T3I, T9c, T3M;
			 {
			      E T3R, T3U, T3T, T9e, T3S;
			      {
				   E T3y, T3B, T3x, T3A, Ta4, T3z, T3Q;
				   T3y = cr[WS(rs, 1)];
				   T3B = ci[WS(rs, 1)];
				   T3x = W[0];
				   T3A = W[1];
				   T3R = cr[WS(rs, 49)];
				   T3U = ci[WS(rs, 49)];
				   Ta4 = T3x * T3B;
				   T3z = T3x * T3y;
				   T3Q = W[96];
				   T3T = W[97];
				   Ta5 = FNMS(T3A, T3y, Ta4);
				   T3C = FMA(T3A, T3B, T3z);
				   T9e = T3Q * T3U;
				   T3S = T3Q * T3R;
			      }
			      {
				   E T3E, T3H, T3D, T3G, Ta6, T3F, T3K;
				   T3E = cr[WS(rs, 33)];
				   T3H = ci[WS(rs, 33)];
				   T9f = FNMS(T3T, T3R, T9e);
				   T3V = FMA(T3T, T3U, T3S);
				   T3D = W[64];
				   T3G = W[65];
				   T3L = cr[WS(rs, 17)];
				   T3O = ci[WS(rs, 17)];
				   Ta6 = T3D * T3H;
				   T3F = T3D * T3E;
				   T3K = W[32];
				   T3N = W[33];
				   Ta7 = FNMS(T3G, T3E, Ta6);
				   T3I = FMA(T3G, T3H, T3F);
				   T9c = T3K * T3O;
				   T3M = T3K * T3L;
			      }
			 }
			 {
			      E Ta8, TfI, T3J, T9b, T9d, T3P;
			      Ta8 = Ta5 - Ta7;
			      TfI = Ta5 + Ta7;
			      T3J = T3C + T3I;
			      T9b = T3C - T3I;
			      T9d = FNMS(T3N, T3L, T9c);
			      T3P = FMA(T3N, T3O, T3M);
			      {
				   E TfJ, T9g, Ta9, T3W;
				   TfJ = T9d + T9f;
				   T9g = T9d - T9f;
				   Ta9 = T3P - T3V;
				   T3W = T3P + T3V;
				   T9h = T9b - T9g;
				   Td3 = T9b + T9g;
				   TfK = TfI - TfJ;
				   ThM = TfI + TfJ;
				   T3X = T3J + T3W;
				   Tfr = T3J - T3W;
				   Tde = Ta8 - Ta9;
				   Taa = Ta8 + Ta9;
			      }
			 }
		    }
	       }
	       {
		    E TaC, T69, Taw, TfU, T5X, Tar, TaA, T63;
		    {
			 E T8S, T3r, T8M, Tfl, T3f, T8H, T8Q, T3l;
			 {
			      E T8k, T8f, T8v, T8e;
			      {
				   E T8a, T2f, T8j, T2y, T2o, T2r, T2q, T8c, T2l, T8g, T2p;
				   {
					E T2u, T2x, T2w, T8i, T2v;
					{
					     E T2b, T2e, T2a, T2d, T89, T2c, T2t;
					     T2b = cr[WS(rs, 10)];
					     T2e = ci[WS(rs, 10)];
					     T2a = W[18];
					     T2d = W[19];
					     T2u = cr[WS(rs, 26)];
					     T2x = ci[WS(rs, 26)];
					     T89 = T2a * T2e;
					     T2c = T2a * T2b;
					     T2t = W[50];
					     T2w = W[51];
					     T8a = FNMS(T2d, T2b, T89);
					     T2f = FMA(T2d, T2e, T2c);
					     T8i = T2t * T2x;
					     T2v = T2t * T2u;
					}
					{
					     E T2h, T2k, T2g, T2j, T8b, T2i, T2n;
					     T2h = cr[WS(rs, 42)];
					     T2k = ci[WS(rs, 42)];
					     T8j = FNMS(T2w, T2u, T8i);
					     T2y = FMA(T2w, T2x, T2v);
					     T2g = W[82];
					     T2j = W[83];
					     T2o = cr[WS(rs, 58)];
					     T2r = ci[WS(rs, 58)];
					     T8b = T2g * T2k;
					     T2i = T2g * T2h;
					     T2n = W[114];
					     T2q = W[115];
					     T8c = FNMS(T2j, T2h, T8b);
					     T2l = FMA(T2j, T2k, T2i);
					     T8g = T2n * T2r;
					     T2p = T2n * T2o;
					}
				   }
				   {
					E T8d, Tfa, T2m, T88, T8h, T2s, Tf9, T2z;
					T8d = T8a - T8c;
					Tfa = T8a + T8c;
					T2m = T2f + T2l;
					T88 = T2f - T2l;
					T8h = FNMS(T2q, T2o, T8g);
					T2s = FMA(T2q, T2r, T2p);
					T8k = T8h - T8j;
					Tf9 = T8h + T8j;
					T2z = T2s + T2y;
					T8f = T2s - T2y;
					T8v = T88 + T8d;
					T8e = T88 - T8d;
					Thx = Tfa + Tf9;
					Tfb = Tf9 - Tfa;
					Tf6 = T2m - T2z;
					T2A = T2m + T2z;
				   }
			      }
			      {
				   E T38, T8J, T3h, T3k, T8L, T3e, T3g, T3j, T8P, T3i;
				   {
					E T3n, T3q, T3m, T3p;
					{
					     E T34, T37, T33, T8w, T8l, T36, T8I, T35;
					     T34 = cr[WS(rs, 6)];
					     T37 = ci[WS(rs, 6)];
					     T33 = W[10];
					     T8w = T8k - T8f;
					     T8l = T8f + T8k;
					     T36 = W[11];
					     T8I = T33 * T37;
					     T35 = T33 * T34;
					     T8x = T8v + T8w;
					     TcO = T8v - T8w;
					     T8m = T8e + T8l;
					     TcR = T8l - T8e;
					     T38 = FMA(T36, T37, T35);
					     T8J = FNMS(T36, T34, T8I);
					}
					T3n = cr[WS(rs, 22)];
					T3q = ci[WS(rs, 22)];
					T3m = W[42];
					T3p = W[43];
					{
					     E T3a, T3d, T3c, T8K, T3b, T8R, T3o, T39;
					     T3a = cr[WS(rs, 38)];
					     T3d = ci[WS(rs, 38)];
					     T8R = T3m * T3q;
					     T3o = T3m * T3n;
					     T39 = W[74];
					     T3c = W[75];
					     T8S = FNMS(T3p, T3n, T8R);
					     T3r = FMA(T3p, T3q, T3o);
					     T8K = T39 * T3d;
					     T3b = T39 * T3a;
					     T3h = cr[WS(rs, 54)];
					     T3k = ci[WS(rs, 54)];
					     T8L = FNMS(T3c, T3a, T8K);
					     T3e = FMA(T3c, T3d, T3b);
					     T3g = W[106];
					     T3j = W[107];
					}
				   }
				   T8M = T8J - T8L;
				   Tfl = T8J + T8L;
				   T3f = T38 + T3e;
				   T8H = T38 - T3e;
				   T8P = T3g * T3k;
				   T3i = T3g * T3h;
				   T8Q = FNMS(T3j, T3h, T8P);
				   T3l = FMA(T3j, T3k, T3i);
			      }
			 }
			 {
			      E T9u, T9p, Tab, T9o;
			      {
				   E T9k, T43, T9t, T4m, T4c, T4f, T4e, T9m, T49, T9q, T4d;
				   {
					E T4i, T4l, T4k, T9s, T4j;
					{
					     E T3Z, T42, T3Y, T41, T9j, T40, T4h;
					     {
						  E T94, T8N, T8T, Tfk, T8O, T3s, T8U, T95;
						  T3Z = cr[WS(rs, 9)];
						  T94 = T8H + T8M;
						  T8N = T8H - T8M;
						  T8T = T8Q - T8S;
						  Tfk = T8Q + T8S;
						  T8O = T3l - T3r;
						  T3s = T3l + T3r;
						  T42 = ci[WS(rs, 9)];
						  Tfm = Tfk - Tfl;
						  ThC = Tfl + Tfk;
						  T8U = T8O + T8T;
						  T95 = T8T - T8O;
						  T3t = T3f + T3s;
						  Tfh = T3f - T3s;
						  T96 = T94 + T95;
						  TcV = T94 - T95;
						  T8V = T8N + T8U;
						  TcY = T8U - T8N;
						  T3Y = W[16];
					     }
					     T41 = W[17];
					     T4i = cr[WS(rs, 25)];
					     T4l = ci[WS(rs, 25)];
					     T9j = T3Y * T42;
					     T40 = T3Y * T3Z;
					     T4h = W[48];
					     T4k = W[49];
					     T9k = FNMS(T41, T3Z, T9j);
					     T43 = FMA(T41, T42, T40);
					     T9s = T4h * T4l;
					     T4j = T4h * T4i;
					}
					{
					     E T45, T48, T44, T47, T9l, T46, T4b;
					     T45 = cr[WS(rs, 41)];
					     T48 = ci[WS(rs, 41)];
					     T9t = FNMS(T4k, T4i, T9s);
					     T4m = FMA(T4k, T4l, T4j);
					     T44 = W[80];
					     T47 = W[81];
					     T4c = cr[WS(rs, 57)];
					     T4f = ci[WS(rs, 57)];
					     T9l = T44 * T48;
					     T46 = T44 * T45;
					     T4b = W[112];
					     T4e = W[113];
					     T9m = FNMS(T47, T45, T9l);
					     T49 = FMA(T47, T48, T46);
					     T9q = T4b * T4f;
					     T4d = T4b * T4c;
					}
				   }
				   {
					E T9n, Tft, T4a, T9i, T9r, T4g, Tfs, T4n;
					T9n = T9k - T9m;
					Tft = T9k + T9m;
					T4a = T43 + T49;
					T9i = T43 - T49;
					T9r = FNMS(T4e, T4c, T9q);
					T4g = FMA(T4e, T4f, T4d);
					T9u = T9r - T9t;
					Tfs = T9r + T9t;
					T4n = T4g + T4m;
					T9p = T4g - T4m;
					Tab = T9i + T9n;
					T9o = T9i - T9n;
					ThN = Tft + Tfs;
					Tfu = Tfs - Tft;
					TfL = T4a - T4n;
					T4o = T4a + T4n;
				   }
			      }
			      {
				   E T5Q, Tat, T5Z, T62, Tav, T5W, T5Y, T61, Taz, T60;
				   {
					E T65, T68, T64, T67;
					{
					     E T5M, T5P, T5L, Tac, T9v, T5O, Tas, T5N;
					     T5M = cr[WS(rs, 7)];
					     T5P = ci[WS(rs, 7)];
					     T5L = W[12];
					     Tac = T9u - T9p;
					     T9v = T9p + T9u;
					     T5O = W[13];
					     Tas = T5L * T5P;
					     T5N = T5L * T5M;
					     Tad = Tab + Tac;
					     Td4 = Tab - Tac;
					     T9w = T9o + T9v;
					     Tdf = T9v - T9o;
					     T5Q = FMA(T5O, T5P, T5N);
					     Tat = FNMS(T5O, T5M, Tas);
					}
					T65 = cr[WS(rs, 23)];
					T68 = ci[WS(rs, 23)];
					T64 = W[44];
					T67 = W[45];
					{
					     E T5S, T5V, T5U, Tau, T5T, TaB, T66, T5R;
					     T5S = cr[WS(rs, 39)];
					     T5V = ci[WS(rs, 39)];
					     TaB = T64 * T68;
					     T66 = T64 * T65;
					     T5R = W[76];
					     T5U = W[77];
					     TaC = FNMS(T67, T65, TaB);
					     T69 = FMA(T67, T68, T66);
					     Tau = T5R * T5V;
					     T5T = T5R * T5S;
					     T5Z = cr[WS(rs, 55)];
					     T62 = ci[WS(rs, 55)];
					     Tav = FNMS(T5U, T5S, Tau);
					     T5W = FMA(T5U, T5V, T5T);
					     T5Y = W[108];
					     T61 = W[109];
					}
				   }
				   Taw = Tat - Tav;
				   TfU = Tat + Tav;
				   T5X = T5Q + T5W;
				   Tar = T5Q - T5W;
				   Taz = T5Y * T62;
				   T60 = T5Y * T5Z;
				   TaA = FNMS(T61, T5Z, Taz);
				   T63 = FMA(T61, T62, T60);
			      }
			 }
		    }
		    {
			 E T9T, Td9, TfE, TfB, Tda, Ta0;
			 {
			      E T9E, Td6, Tfz, Tfw, Td7, T9L;
			      {
				   E T9G, T4v, T9C, T4O, T4E, T4H, T4G, T9I, T4B, T9z, T4F;
				   {
					E T4K, T4N, T4M, T9B, T4L;
					{
					     E T4r, T4u, T4q, T4t, T9F, T4s, T4J;
					     {
						  E Tbl, Tax, TaD, TfT, Tay, T6a, TaE, Tbk;
						  T4r = cr[WS(rs, 5)];
						  Tbl = Tar + Taw;
						  Tax = Tar - Taw;
						  TaD = TaA - TaC;
						  TfT = TaA + TaC;
						  Tay = T63 - T69;
						  T6a = T63 + T69;
						  T4u = ci[WS(rs, 5)];
						  TfV = TfT - TfU;
						  ThY = TfU + TfT;
						  TaE = Tay + TaD;
						  Tbk = Tay - TaD;
						  T6b = T5X + T6a;
						  Tg9 = T6a - T5X;
						  Tbm = Tbk - Tbl;
						  Tdn = Tbl + Tbk;
						  TaF = Tax + TaE;
						  Tdy = TaE - Tax;
						  T4q = W[8];
					     }
					     T4t = W[9];
					     T4K = cr[WS(rs, 53)];
					     T4N = ci[WS(rs, 53)];
					     T9F = T4q * T4u;
					     T4s = T4q * T4r;
					     T4J = W[104];
					     T4M = W[105];
					     T9G = FNMS(T4t, T4r, T9F);
					     T4v = FMA(T4t, T4u, T4s);
					     T9B = T4J * T4N;
					     T4L = T4J * T4K;
					}
					{
					     E T4x, T4A, T4w, T4z, T9H, T4y, T4D;
					     T4x = cr[WS(rs, 37)];
					     T4A = ci[WS(rs, 37)];
					     T9C = FNMS(T4M, T4K, T9B);
					     T4O = FMA(T4M, T4N, T4L);
					     T4w = W[72];
					     T4z = W[73];
					     T4E = cr[WS(rs, 21)];
					     T4H = ci[WS(rs, 21)];
					     T9H = T4w * T4A;
					     T4y = T4w * T4x;
					     T4D = W[40];
					     T4G = W[41];