hf_64.c

快速fft变换

					     T9I = FNMS(T4z, T4x, T9H);
					     T4B = FMA(T4z, T4A, T4y);
					     T9z = T4D * T4H;
					     T4F = T4D * T4E;
					}
				   }
				   {
					E T9J, Tfx, T4C, T9y, T9A, T4I;
					T9J = T9G - T9I;
					Tfx = T9G + T9I;
					T4C = T4v + T4B;
					T9y = T4v - T4B;
					T9A = FNMS(T4G, T4E, T9z);
					T4I = FMA(T4G, T4H, T4F);
					{
					     E Tfy, T9D, T9K, T4P;
					     Tfy = T9A + T9C;
					     T9D = T9A - T9C;
					     T9K = T4I - T4O;
					     T4P = T4I + T4O;
					     T9E = T9y - T9D;
					     Td6 = T9y + T9D;
					     Tfz = Tfx - Tfy;
					     ThJ = Tfx + Tfy;
					     Tfw = T4C - T4P;
					     T4Q = T4C + T4P;
					     Td7 = T9J - T9K;
					     T9L = T9J + T9K;
					}
				   }
			      }
			      {
				   E T9V, T4W, T9R, T5f, T55, T58, T57, T9X, T52, T9O, T56;
				   {
					E T5b, T5e, T5d, T9Q, T5c;
					{
					     E T4S, T4V, T4R, T4U, T9U, T4T, T5a;
					     T4S = cr[WS(rs, 61)];
					     TfN = Tfw + Tfz;
					     TfA = Tfw - Tfz;
					     Taf = FMA(KP414213562, T9E, T9L);
					     T9M = FNMS(KP414213562, T9L, T9E);
					     Td8 = FMA(KP414213562, Td7, Td6);
					     Tdh = FNMS(KP414213562, Td6, Td7);
					     T4V = ci[WS(rs, 61)];
					     T4R = W[120];
					     T4U = W[121];
					     T5b = cr[WS(rs, 45)];
					     T5e = ci[WS(rs, 45)];
					     T9U = T4R * T4V;
					     T4T = T4R * T4S;
					     T5a = W[88];
					     T5d = W[89];
					     T9V = FNMS(T4U, T4S, T9U);
					     T4W = FMA(T4U, T4V, T4T);
					     T9Q = T5a * T5e;
					     T5c = T5a * T5b;
					}
					{
					     E T4Y, T51, T4X, T50, T9W, T4Z, T54;
					     T4Y = cr[WS(rs, 29)];
					     T51 = ci[WS(rs, 29)];
					     T9R = FNMS(T5d, T5b, T9Q);
					     T5f = FMA(T5d, T5e, T5c);
					     T4X = W[56];
					     T50 = W[57];
					     T55 = cr[WS(rs, 13)];
					     T58 = ci[WS(rs, 13)];
					     T9W = T4X * T51;
					     T4Z = T4X * T4Y;
					     T54 = W[24];
					     T57 = W[25];
					     T9X = FNMS(T50, T4Y, T9W);
					     T52 = FMA(T50, T51, T4Z);
					     T9O = T54 * T58;
					     T56 = T54 * T55;
					}
				   }
				   {
					E T9Y, TfC, T53, T9N, T9P, T59;
					T9Y = T9V - T9X;
					TfC = T9V + T9X;
					T53 = T4W + T52;
					T9N = T4W - T52;
					T9P = FNMS(T57, T55, T9O);
					T59 = FMA(T57, T58, T56);
					{
					     E TfD, T9S, T9Z, T5g;
					     TfD = T9P + T9R;
					     T9S = T9P - T9R;
					     T9Z = T59 - T5f;
					     T5g = T59 + T5f;
					     T9T = T9N - T9S;
					     Td9 = T9N + T9S;
					     TfE = TfC - TfD;
					     ThI = TfC + TfD;
					     TfB = T53 - T5g;
					     T5h = T53 + T5g;
					     Tda = T9Y - T9Z;
					     Ta0 = T9Y + T9Z;
					}
				   }
			      }
			 }
			 {
			      E TaN, Tdp, Tg0, TfX, Tdq, TaU;
			      {
				   E TaQ, T6i, TaL, T6B, T6r, T6u, T6t, TaS, T6o, TaI, T6s;
				   {
					E T6x, T6A, T6z, TaK, T6y;
					{
					     E T6e, T6h, T6d, T6g, TaP, T6f, T6w;
					     T6e = cr[WS(rs, 3)];
					     TfO = TfE - TfB;
					     TfF = TfB + TfE;
					     Tag = FNMS(KP414213562, T9T, Ta0);
					     Ta1 = FMA(KP414213562, Ta0, T9T);
					     Tdb = FNMS(KP414213562, Tda, Td9);
					     Tdi = FMA(KP414213562, Td9, Tda);
					     T6h = ci[WS(rs, 3)];
					     T6d = W[4];
					     T6g = W[5];
					     T6x = cr[WS(rs, 51)];
					     T6A = ci[WS(rs, 51)];
					     TaP = T6d * T6h;
					     T6f = T6d * T6e;
					     T6w = W[100];
					     T6z = W[101];
					     TaQ = FNMS(T6g, T6e, TaP);
					     T6i = FMA(T6g, T6h, T6f);
					     TaK = T6w * T6A;
					     T6y = T6w * T6x;
					}
					{
					     E T6k, T6n, T6j, T6m, TaR, T6l, T6q;
					     T6k = cr[WS(rs, 35)];
					     T6n = ci[WS(rs, 35)];
					     TaL = FNMS(T6z, T6x, TaK);
					     T6B = FMA(T6z, T6A, T6y);
					     T6j = W[68];
					     T6m = W[69];
					     T6r = cr[WS(rs, 19)];
					     T6u = ci[WS(rs, 19)];
					     TaR = T6j * T6n;
					     T6l = T6j * T6k;
					     T6q = W[36];
					     T6t = W[37];
					     TaS = FNMS(T6m, T6k, TaR);
					     T6o = FMA(T6m, T6n, T6l);
					     TaI = T6q * T6u;
					     T6s = T6q * T6r;
					}
				   }
				   {
					E TaT, TfY, T6p, TaH, TaJ, T6v;
					TaT = TaQ - TaS;
					TfY = TaQ + TaS;
					T6p = T6i + T6o;
					TaH = T6i - T6o;
					TaJ = FNMS(T6t, T6r, TaI);
					T6v = FMA(T6t, T6u, T6s);
					{
					     E TfZ, TaM, T6C, TaO;
					     TfZ = TaJ + TaL;
					     TaM = TaJ - TaL;
					     T6C = T6v + T6B;
					     TaO = T6B - T6v;
					     TaN = TaH - TaM;
					     Tdp = TaH + TaM;
					     Tg0 = TfY - TfZ;
					     ThU = TfY + TfZ;
					     TfX = T6p - T6C;
					     T6D = T6p + T6C;
					     Tdq = TaT + TaO;
					     TaU = TaO - TaT;
					}
				   }
			      }
			      {
				   E Tb5, T6J, Tb0, T72, T6S, T6V, T6U, Tb7, T6P, TaX, T6T;
				   {
					E T6Y, T71, T70, TaZ, T6Z;
					{
					     E T6F, T6I, T6E, T6H, Tb4, T6G, T6X;
					     T6F = cr[WS(rs, 59)];
					     Tgf = TfX + Tg0;
					     Tg1 = TfX - Tg0;
					     Tbo = FNMS(KP414213562, TaN, TaU);
					     TaV = FMA(KP414213562, TaU, TaN);
					     Tdr = FMA(KP414213562, Tdq, Tdp);
					     TdA = FNMS(KP414213562, Tdp, Tdq);
					     T6I = ci[WS(rs, 59)];
					     T6E = W[116];
					     T6H = W[117];
					     T6Y = cr[WS(rs, 43)];
					     T71 = ci[WS(rs, 43)];
					     Tb4 = T6E * T6I;
					     T6G = T6E * T6F;
					     T6X = W[84];
					     T70 = W[85];
					     Tb5 = FNMS(T6H, T6F, Tb4);
					     T6J = FMA(T6H, T6I, T6G);
					     TaZ = T6X * T71;
					     T6Z = T6X * T6Y;
					}
					{
					     E T6L, T6O, T6K, T6N, Tb6, T6M, T6R;
					     T6L = cr[WS(rs, 27)];
					     T6O = ci[WS(rs, 27)];
					     Tb0 = FNMS(T70, T6Y, TaZ);
					     T72 = FMA(T70, T71, T6Z);
					     T6K = W[52];
					     T6N = W[53];
					     T6S = cr[WS(rs, 11)];
					     T6V = ci[WS(rs, 11)];
					     Tb6 = T6K * T6O;
					     T6M = T6K * T6L;
					     T6R = W[20];
					     T6U = W[21];
					     Tb7 = FNMS(T6N, T6L, Tb6);
					     T6P = FMA(T6N, T6O, T6M);
					     TaX = T6R * T6V;
					     T6T = T6R * T6S;
					}
				   }
				   {
					E Tb8, Tg3, T6Q, TaW, TaY, T6W;
					Tb8 = Tb5 - Tb7;
					Tg3 = Tb5 + Tb7;
					T6Q = T6J + T6P;
					TaW = T6J - T6P;
					TaY = FNMS(T6U, T6S, TaX);
					T6W = FMA(T6U, T6V, T6T);
					{
					     E Tg4, Tb1, T73, Tb3;
					     Tg4 = TaY + Tb0;
					     Tb1 = TaY - Tb0;
					     T73 = T6W + T72;
					     Tb3 = T72 - T6W;
					     Tb2 = TaW - Tb1;
					     Tds = TaW + Tb1;
					     Tg5 = Tg3 - Tg4;
					     ThT = Tg3 + Tg4;
					     Tg2 = T6Q - T73;
					     T74 = T6Q + T73;
					     Tdt = Tb8 + Tb3;
					     Tb9 = Tb3 - Tb8;
					}
				   }
			      }
			 }
		    }
	       }
	       {
		    E Thq, Tge, Tg6, Tdu, TdB, Tj7, Thv, ThA, Tht, Tj8, ThD, Thy, ThS, Ti0, ThZ;
		    E ThV, ThH, ThP, ThO, ThK, Tkm, TcD, Tk0, Tk4, TjZ, Tk3, Tik, Tin;
		    {
			 E Tbp, Tba, TiI, TiL;
			 {
			      E Tio, T1I, Tj1, T3v, Tj2, TiX, TiN, Tir, T76, TiJ, TiC, TiG, T5j, Tit, Tiw;
			      E TiK;
			      {
				   E TiO, TiW, Tip, Tiq;
				   {
					E TO, T1H, T2B, T3u;
					Thq = Tm - TN;
					TO = Tm + TN;
					Tge = Tg2 - Tg5;
					Tg6 = Tg2 + Tg5;
					Tbp = FMA(KP414213562, Tb2, Tb9);
					Tba = FNMS(KP414213562, Tb9, Tb2);
					Tdu = FNMS(KP414213562, Tdt, Tds);
					TdB = FMA(KP414213562, Tds, Tdt);
					T1H = T1f + T1G;
					Tj7 = T1f - T1G;
					Thv = T29 - T2A;
					T2B = T29 + T2A;
					T3u = T32 + T3t;
					ThA = T32 - T3t;
					Tht = Thr - Ths;
					TiO = Ths + Thr;
					Tio = TO - T1H;
					T1I = TO + T1H;
					Tj1 = T2B - T3u;
					T3v = T2B + T3u;
					TiW = TiP + TiV;
					Tj8 = TiV - TiP;
				   }
				   ThD = ThB - ThC;
				   Tip = ThB + ThC;
				   Tiq = Thw + Thx;
				   Thy = Thw - Thx;
				   {
					E T6c, T75, Tiz, TiA;
					ThS = T5K - T6b;
					T6c = T5K + T6b;
					Tj2 = TiW - TiO;
					TiX = TiO + TiW;
					TiN = Tiq + Tip;
					Tir = Tip - Tiq;
					T75 = T6D + T74;
					Ti0 = T74 - T6D;
					ThZ = ThX - ThY;
					Tiz = ThX + ThY;
					TiA = ThU + ThT;
					ThV = ThT - ThU;
					{
					     E T4p, Tiy, TiB, T5i, Tiu, Tiv;
					     ThH = T3X - T4o;
					     T4p = T3X + T4o;
					     T76 = T6c + T75;
					     Tiy = T6c - T75;
					     TiJ = Tiz + TiA;
					     TiB = Tiz - TiA;
					     T5i = T4Q + T5h;
					     ThP = T4Q - T5h;
					     ThO = ThM - ThN;
					     Tiu = ThM + ThN;
					     Tiv = ThJ + ThI;
					     ThK = ThI - ThJ;
					     TiC = Tiy - TiB;
					     TiG = Tiy + TiB;
					     T5j = T4p + T5i;
					     Tit = T4p - T5i;
					     Tiw = Tiu - Tiv;
					     TiK = Tiu + Tiv;
					}
				   }
			      }
			      {
				   E TiZ, TiD, TiH, TiE, Tis, TiM, TiY, Tj0;
				   {
					E T3w, TiF, Tix, T77, Tj5, Tj3, Tj6, Tj4;
					TiI = T1I - T3v;
					T3w = T1I + T3v;
					TiF = Tit - Tiw;
					Tix = Tit + Tiw;
					T77 = T5j + T76;
					TiZ = T76 - T5j;
					Tj5 = Tj2 - Tj1;
					Tj3 = Tj1 + Tj2;
					TiD = Tix + TiC;
					Tj4 = TiC - Tix;
					cr[0] = T3w + T77;
					ci[WS(rs, 31)] = T3w - T77;
					Tj6 = TiG - TiF;
					TiH = TiF + TiG;
					ci[WS(rs, 39)] = FMA(KP707106781, Tj4, Tj3);
					cr[WS(rs, 56)] = FMS(KP707106781, Tj4, Tj3);
					TiE = Tio + Tir;
					Tis = Tio - Tir;
					ci[WS(rs, 55)] = FMA(KP707106781, Tj6, Tj5);
					cr[WS(rs, 40)] = FMS(KP707106781, Tj6, Tj5);
				   }
				   TiL = TiJ - TiK;
				   TiM = TiK + TiJ;
				   cr[WS(rs, 8)] = FMA(KP707106781, TiD, Tis);
				   ci[WS(rs, 23)] = FNMS(KP707106781, TiD, Tis);
				   ci[WS(rs, 7)] = FMA(KP707106781, TiH, TiE);
				   cr[WS(rs, 24)] = FNMS(KP707106781, TiH, TiE);
				   TiY = TiN + TiX;
				   Tj0 = TiX - TiN;
				   ci[WS(rs, 63)] = TiM + TiY;
				   cr[WS(rs, 32)] = TiM - TiY;
				   ci[WS(rs, 47)] = TiZ + Tj0;
				   cr[WS(rs, 48)] = TiZ - Tj0;
			      }
			 }
			 {
			      E TjW, TbB, Tk2, T99, TbF, TbL, Tbv, Taj, Tcu, Tcy, Tci, Tce, Tcr, Tcx, Tch;
			      E Tc7, Tcn, Tkg, Tka, TbZ, TbP, T7J, TbO, T7u, Tk7, TjT, TbI, TbM, Tbw, Tbs;
			      E T7Y, TbQ;
			      {
				   E TbX, TbW, TbU, TbT, Tc1, Tc5, Tc4, Tc2, TaG, Tbq, Tbn, Tcb, Tcs, Tca, Tcc;
				   E Tbb, Tcm, TbV;
				   {
					E T8W, Tbz, T8z, T97, T8n, T8y;
					TbX = FNMS(KP707106781, T8m, T87);
					T8n = FMA(KP707106781, T8m, T87);
					T8y = FMA(KP707106781, T8x, T8u);
					TbW = FNMS(KP707106781, T8x, T8u);
					TbU = FNMS(KP707106781, T8V, T8G);
					T8W = FMA(KP707106781, T8V, T8G);
					ci[WS(rs, 15)] = TiI + TiL;
					cr[WS(rs, 16)] = TiI - TiL;
					Tbz = FMA(KP198912367, T8n, T8y);
					T8z = FNMS(KP198912367, T8y, T8n);
					T97 = FMA(KP707106781, T96, T93);
					TbT = FNMS(KP707106781, T96, T93);
					{
					     E Tae, TbD, Ta3, Tah;
					     {
						  E T9x, Ta2, TbA, T98;
						  Tc1 = FNMS(KP707106781, T9w, T9h);
						  T9x = FMA(KP707106781, T9w, T9h);
						  Ta2 = T9M + Ta1;
						  Tc5 = Ta1 - T9M;
						  Tc4 = FNMS(KP707106781, Tad, Taa);
						  Tae = FMA(KP707106781, Tad, Taa);
						  TbA = FNMS(KP198912367, T8W, T97);
						  T98 = FMA(KP198912367, T97, T8W);
						  TbD = FNMS(KP923879532, Ta2, T9x);
						  Ta3 = FMA(KP923879532, Ta2, T9x);
						  TjW = Tbz + TbA;
						  TbB = Tbz - TbA;
						  Tk2 = T98 - T8z;
						  T99 = T8z + T98;
						  Tah = Taf + Tag;
						  Tc2 = Taf - Tag;
					     }
					     {
						  E Tc8, Tc9, TbE, Tai;
						  TaG = FMA(KP707106781, TaF, Taq);
						  Tc8 = FNMS(KP707106781, TaF, Taq);
						  Tc9 = Tbp - Tbo;
						  Tbq = Tbo + Tbp;
						  Tbn = FMA(KP707106781, Tbm, Tbj);
						  Tcb = FNMS(KP707106781, Tbm, Tbj);
						  TbE = FNMS(KP923879532, Tah, Tae);
						  Tai = FMA(KP923879532, Tah, Tae);
						  Tcs = FMA(KP923879532, Tc9, Tc8);
						  Tca = FNMS(KP923879532, Tc9, Tc8);
						  TbF = FMA(KP820678790, TbE, TbD);
						  TbL = FNMS(KP820678790, TbD, TbE);
						  Tbv = FMA(KP098491403, Ta3, Tai);
						  Taj = FNMS(KP098491403, Tai, Ta3);
						  Tcc = Tba - TaV;
						  Tbb = TaV + Tba;
					     }
					}
				   }
				   {
					E Tcp, Tc3, Tct, Tcd, Tcq, Tc6;
					Tct = FNMS(KP923879532, Tcc, Tcb);
					Tcd = FMA(KP923879532, Tcc, Tcb);
					Tcp = FMA(KP923879532, Tc2, Tc1);
					Tc3 = FNMS(KP923879532, Tc2, Tc1);
					Tcu = FMA(KP303346683, Tct, Tcs);
					Tcy = FNMS(KP303346683, Tcs, Tct);
					Tci = FMA(KP534511135, Tca, Tcd);
					Tce = FNMS(KP534511135, Tcd, Tca);
					Tcq = FMA(KP923879532, Tc5, Tc4);
					Tc6 = FNMS(KP923879532, Tc5, Tc4);
					Tcm = FNMS(KP668178637, TbT, TbU);
					TbV = FMA(KP668178637, TbU, TbT);
					Tcr = FMA(KP303346683, Tcq, Tcp);
					Tcx = FNMS(KP303346683, Tcp, Tcq);
					Tch = FMA(KP534511135, Tc3, Tc6);
					Tc7 = FNMS(KP534511135, Tc6, Tc3);
				   }