hb_64.c

快速fft变换

						  Tgv = Tgt * Tgu;
						  ci[WS(rs, 42)] = FMA(Tgg, Tge, Tgs);
						  cr[WS(rs, 42)] = FNMS(Tgg, Tgr, Tgf);
						  ci[WS(rs, 10)] = FMA(Tgw, Tgu, Tgy);
					     }
					     cr[WS(rs, 10)] = FNMS(Tgw, Tgx, Tgv);
					     TgN = FMA(KP831469612, TgG, TgF);
					     TgH = FNMS(KP831469612, TgG, TgF);
					     TgE = W[51];
					     Tgz = W[50];
					     {
						  E TgJ, TgO, TgI, TgD;
						  TgM = W[115];
						  TgI = TgE * TgC;
						  TgD = Tgz * TgC;
						  TgJ = W[114];
						  TgO = TgM * TgK;
						  ci[WS(rs, 26)] = FMA(Tgz, TgH, TgI);
						  cr[WS(rs, 26)] = FNMS(TgE, TgH, TgD);
						  TgL = TgJ * TgK;
						  ci[WS(rs, 58)] = FMA(TgJ, TgN, TgO);
					     }
					}
				   }
				   {
					E Th5, Th8, Ths, Thk, Thv, Thp, Thc, Th0;
					{
					     E TgV, TgY, Thn, Thj, TgS, Thi, Th6, Th7, Tho, TgZ;
					     cr[WS(rs, 58)] = FNMS(TgM, TgN, TgL);
					     TgV = FNMS(KP198912367, TgU, TgT);
					     Th6 = FMA(KP198912367, TgT, TgU);
					     Th7 = FNMS(KP198912367, TgW, TgX);
					     TgY = FMA(KP198912367, TgX, TgW);
					     Th5 = FMA(KP923879532, Th4, Th3);
					     Thn = FNMS(KP923879532, Th4, Th3);
					     Thj = Th7 - Th6;
					     Th8 = Th6 + Th7;
					     TgS = FMA(KP923879532, TgR, TgQ);
					     Thi = FNMS(KP923879532, TgR, TgQ);
					     Tho = TgV - TgY;
					     TgZ = TgV + TgY;
					     Ths = FMA(KP980785280, Thj, Thi);
					     Thk = FNMS(KP980785280, Thj, Thi);
					     Thv = FMA(KP980785280, Tho, Thn);
					     Thp = FNMS(KP980785280, Tho, Thn);
					     Thc = FMA(KP980785280, TgZ, TgS);
					     Th0 = FNMS(KP980785280, TgZ, TgS);
					}
					{
					     E Thu, Tht, Thh, Thm, Th9, Th2, TgP;
					     Thh = W[98];
					     Thm = W[99];
					     {
						  E Thr, Thq, Thl, Thw;
						  Thr = W[34];
						  Thu = W[35];
						  Thq = Thh * Thp;
						  Thl = Thh * Thk;
						  Thw = Thr * Thv;
						  Tht = Thr * Ths;
						  ci[WS(rs, 50)] = FMA(Thm, Thk, Thq);
						  cr[WS(rs, 50)] = FNMS(Thm, Thp, Thl);
						  ci[WS(rs, 18)] = FMA(Thu, Ths, Thw);
					     }
					     cr[WS(rs, 18)] = FNMS(Thu, Thv, Tht);
					     Thf = FMA(KP980785280, Th8, Th5);
					     Th9 = FNMS(KP980785280, Th8, Th5);
					     Th2 = W[67];
					     TgP = W[66];
					     {
						  E Thb, Thg, Tha, Th1;
						  The = W[3];
						  Tha = Th2 * Th0;
						  Th1 = TgP * Th0;
						  Thb = W[2];
						  Thg = The * Thc;
						  ci[WS(rs, 34)] = FMA(TgP, Th9, Tha);
						  cr[WS(rs, 34)] = FNMS(Th2, Th9, Th1);
						  Thd = Thb * Thc;
						  ci[WS(rs, 2)] = FMA(Thb, Thf, Thg);
					     }
					}
				   }
			      }
			 }
		    }
		    {
			 E Tcl, Tc9, Tc8, Tcm, T9R, T93, T8O, T9U, Tez, Tdt, Td6, TeC, Tfv, Tfu, Tft;
			 E T8B, T8A, T8z;
			 {
			      E TbP, TbO, TbN, T6B, T6A, T6z, TaN, TaM, TaL;
			      {
				   E T6r, T6u, T6O, T6G, T6R, T6L, T6y, T6m;
				   {
					E T6k, T6h, T6J, T6F, T6e, T6E, T6s, T6t, T6K, T6l;
					cr[WS(rs, 2)] = FNMS(The, Thf, Thd);
					T6k = FMA(KP098491403, T6j, T6i);
					T6s = FNMS(KP098491403, T6i, T6j);
					T6t = FMA(KP098491403, T6f, T6g);
					T6h = FNMS(KP098491403, T6g, T6f);
					T6r = FNMS(KP980785280, T6q, T6p);
					T6J = FMA(KP980785280, T6q, T6p);
					T6F = T6s + T6t;
					T6u = T6s - T6t;
					T6e = FNMS(KP980785280, T6d, T6c);
					T6E = FMA(KP980785280, T6d, T6c);
					T6K = T6k + T6h;
					T6l = T6h - T6k;
					T6O = FMA(KP995184726, T6F, T6E);
					T6G = FNMS(KP995184726, T6F, T6E);
					T6R = FMA(KP995184726, T6K, T6J);
					T6L = FNMS(KP995184726, T6K, T6J);
					T6y = FMA(KP995184726, T6l, T6e);
					T6m = FNMS(KP995184726, T6l, T6e);
				   }
				   {
					E T6Q, T6P, T6D, T6I, T6v, T6o, T6b;
					T6D = W[64];
					T6I = W[65];
					{
					     E T6N, T6M, T6H, T6S;
					     T6N = W[0];
					     T6Q = W[1];
					     T6M = T6D * T6L;
					     T6H = T6D * T6G;
					     T6S = T6N * T6R;
					     T6P = T6N * T6O;
					     ci[WS(rs, 33)] = FMA(T6I, T6G, T6M);
					     cr[WS(rs, 33)] = FNMS(T6I, T6L, T6H);
					     ci[WS(rs, 1)] = FMA(T6Q, T6O, T6S);
					}
					cr[WS(rs, 1)] = FNMS(T6Q, T6R, T6P);
					T6B = FMA(KP995184726, T6u, T6r);
					T6v = FNMS(KP995184726, T6u, T6r);
					T6o = W[97];
					T6b = W[96];
					{
					     E T6x, T6C, T6w, T6n;
					     T6A = W[33];
					     T6w = T6o * T6m;
					     T6n = T6b * T6m;
					     T6x = W[32];
					     T6C = T6A * T6y;
					     ci[WS(rs, 49)] = FMA(T6b, T6v, T6w);
					     cr[WS(rs, 49)] = FNMS(T6o, T6v, T6n);
					     T6z = T6x * T6y;
					     ci[WS(rs, 17)] = FMA(T6x, T6B, T6C);
					}
				   }
			      }
			      {
				   E TbF, TbI, Tc2, TbU, Tc5, TbZ, TbM, Tbw;
				   {
					E Tbn, Tbu, TbX, TbT, Tbg, TbS, TbY, Tbv;
					{
					     E TbG, TbH, TbB, TbE, Tb8, Tbf;
					     TbB = FMA(KP923879532, TbA, Tbz);
					     Tcl = FNMS(KP923879532, TbA, Tbz);
					     Tc9 = TbC + TbD;
					     TbE = TbC - TbD;
					     cr[WS(rs, 17)] = FNMS(T6A, T6B, T6z);
					     Tbn = FNMS(KP820678790, Tbm, Tbj);
					     TbG = FMA(KP820678790, Tbj, Tbm);
					     TbH = FMA(KP820678790, Tbq, Tbt);
					     Tbu = FNMS(KP820678790, Tbt, Tbq);
					     TbF = FMA(KP980785280, TbE, TbB);
					     TbX = FNMS(KP980785280, TbE, TbB);
					     Tb8 = FNMS(KP923879532, Tb7, Tb6);
					     Tc8 = FMA(KP923879532, Tb7, Tb6);
					     Tcm = Tbe - Tbb;
					     Tbf = Tbb + Tbe;
					     TbT = TbG + TbH;
					     TbI = TbG - TbH;
					     Tbg = FNMS(KP980785280, Tbf, Tb8);
					     TbS = FMA(KP980785280, Tbf, Tb8);
					}
					TbY = Tbn - Tbu;
					Tbv = Tbn + Tbu;
					Tc2 = FMA(KP773010453, TbT, TbS);
					TbU = FNMS(KP773010453, TbT, TbS);
					Tc5 = FNMS(KP773010453, TbY, TbX);
					TbZ = FMA(KP773010453, TbY, TbX);
					TbM = FMA(KP773010453, Tbv, Tbg);
					Tbw = FNMS(KP773010453, Tbv, Tbg);
				   }
				   {
					E Tc4, Tc3, TbR, TbW, TbJ, Tby, Tb5;
					TbR = W[44];
					TbW = W[45];
					{
					     E Tc1, Tc0, TbV, Tc6;
					     Tc1 = W[108];
					     Tc4 = W[109];
					     Tc0 = TbR * TbZ;
					     TbV = TbR * TbU;
					     Tc6 = Tc1 * Tc5;
					     Tc3 = Tc1 * Tc2;
					     ci[WS(rs, 23)] = FMA(TbW, TbU, Tc0);
					     cr[WS(rs, 23)] = FNMS(TbW, TbZ, TbV);
					     ci[WS(rs, 55)] = FMA(Tc4, Tc2, Tc6);
					}
					cr[WS(rs, 55)] = FNMS(Tc4, Tc5, Tc3);
					TbP = FMA(KP773010453, TbI, TbF);
					TbJ = FNMS(KP773010453, TbI, TbF);
					Tby = W[77];
					Tb5 = W[76];
					{
					     E TbL, TbQ, TbK, Tbx;
					     TbO = W[13];
					     TbK = Tby * Tbw;
					     Tbx = Tb5 * Tbw;
					     TbL = W[12];
					     TbQ = TbO * TbM;
					     ci[WS(rs, 39)] = FMA(Tb5, TbJ, TbK);
					     cr[WS(rs, 39)] = FNMS(Tby, TbJ, Tbx);
					     TbN = TbL * TbM;
					     ci[WS(rs, 7)] = FMA(TbL, TbP, TbQ);
					}
				   }
			      }
			      {
				   E TaD, TaG, Tb0, TaS, Tb3, TaX, TaK, Tay;
				   {
					E Tat, Taw, TaV, TaR, Taq, TaQ, TaW, Tax;
					{
					     E TaE, TaF, TaB, TaC, Tao, Tap;
					     TaB = FMA(KP923879532, T9Q, T9N);
					     T9R = FNMS(KP923879532, T9Q, T9N);
					     T93 = T8V + T92;
					     TaC = T8V - T92;
					     cr[WS(rs, 7)] = FNMS(TbO, TbP, TbN);
					     Tat = FNMS(KP303346683, Tas, Tar);
					     TaE = FMA(KP303346683, Tar, Tas);
					     TaF = FMA(KP303346683, Tau, Tav);
					     Taw = FNMS(KP303346683, Tav, Tau);
					     TaD = FMA(KP831469612, TaC, TaB);
					     TaV = FNMS(KP831469612, TaC, TaB);
					     Tao = FNMS(KP923879532, T8N, T8G);
					     T8O = FMA(KP923879532, T8N, T8G);
					     T9U = T9S - T9T;
					     Tap = T9S + T9T;
					     TaR = TaE + TaF;
					     TaG = TaE - TaF;
					     Taq = FMA(KP831469612, Tap, Tao);
					     TaQ = FNMS(KP831469612, Tap, Tao);
					}
					TaW = Tat - Taw;
					Tax = Tat + Taw;
					Tb0 = FMA(KP956940335, TaR, TaQ);
					TaS = FNMS(KP956940335, TaR, TaQ);
					Tb3 = FNMS(KP956940335, TaW, TaV);
					TaX = FMA(KP956940335, TaW, TaV);
					TaK = FMA(KP956940335, Tax, Taq);
					Tay = FNMS(KP956940335, Tax, Taq);
				   }
				   {
					E Tb2, Tb1, TaP, TaU, TaH, TaA, Tan;
					TaP = W[36];
					TaU = W[37];
					{
					     E TaZ, TaY, TaT, Tb4;
					     TaZ = W[100];
					     Tb2 = W[101];
					     TaY = TaP * TaX;
					     TaT = TaP * TaS;
					     Tb4 = TaZ * Tb3;
					     Tb1 = TaZ * Tb0;
					     ci[WS(rs, 19)] = FMA(TaU, TaS, TaY);
					     cr[WS(rs, 19)] = FNMS(TaU, TaX, TaT);
					     ci[WS(rs, 51)] = FMA(Tb2, Tb0, Tb4);
					}
					cr[WS(rs, 51)] = FNMS(Tb2, Tb3, Tb1);
					TaN = FMA(KP956940335, TaG, TaD);
					TaH = FNMS(KP956940335, TaG, TaD);
					TaA = W[69];
					Tan = W[68];
					{
					     E TaJ, TaO, TaI, Taz;
					     TaM = W[5];
					     TaI = TaA * Tay;
					     Taz = Tan * Tay;
					     TaJ = W[4];
					     TaO = TaM * TaK;
					     ci[WS(rs, 35)] = FMA(Tan, TaH, TaI);
					     cr[WS(rs, 35)] = FNMS(TaA, TaH, Taz);
					     TaL = TaJ * TaK;
					     ci[WS(rs, 3)] = FMA(TaJ, TaN, TaO);
					}
				   }
			      }
			      {
				   E Tfl, Tfo, TfI, TfA, TfL, TfF, Tfs, Tfg;
				   {
					E Tfe, Tfb, TfD, Tfz, Tf8, Tfy, TfE, Tff;
					{
					     E Tfm, Tfn, Tfj, Tfk, Tf6, Tf7;
					     Tfj = FNMS(KP707106781, Tey, Tev);
					     Tez = FMA(KP707106781, Tey, Tev);
					     Tdt = Tdh - Tds;
					     Tfk = Tds + Tdh;
					     cr[WS(rs, 3)] = FNMS(TaM, TaN, TaL);
					     Tfe = FNMS(KP198912367, Tfd, Tfc);
					     Tfm = FMA(KP198912367, Tfc, Tfd);
					     Tfn = FNMS(KP198912367, Tf9, Tfa);
					     Tfb = FMA(KP198912367, Tfa, Tf9);
					     Tfl = FNMS(KP923879532, Tfk, Tfj);
					     TfD = FMA(KP923879532, Tfk, Tfj);
					     Tf6 = FNMS(KP707106781, Td5, TcU);
					     Td6 = FMA(KP707106781, Td5, TcU);
					     TeC = TeA - TeB;
					     Tf7 = TeA + TeB;
					     Tfz = Tfm + Tfn;
					     Tfo = Tfm - Tfn;
					     Tf8 = FNMS(KP923879532, Tf7, Tf6);
					     Tfy = FMA(KP923879532, Tf7, Tf6);
					}
					TfE = Tfe + Tfb;
					Tff = Tfb - Tfe;
					TfI = FMA(KP980785280, Tfz, Tfy);
					TfA = FNMS(KP980785280, Tfz, Tfy);
					TfL = FMA(KP980785280, TfE, TfD);
					TfF = FNMS(KP980785280, TfE, TfD);
					Tfs = FMA(KP980785280, Tff, Tf8);
					Tfg = FNMS(KP980785280, Tff, Tf8);
				   }
				   {
					E TfK, TfJ, Tfx, TfC, Tfp, Tfi, Tf5;
					Tfx = W[58];
					TfC = W[59];
					{
					     E TfH, TfG, TfB, TfM;
					     TfH = W[122];
					     TfK = W[123];
					     TfG = Tfx * TfF;
					     TfB = Tfx * TfA;
					     TfM = TfH * TfL;
					     TfJ = TfH * TfI;
					     ci[WS(rs, 30)] = FMA(TfC, TfA, TfG);
					     cr[WS(rs, 30)] = FNMS(TfC, TfF, TfB);
					     ci[WS(rs, 62)] = FMA(TfK, TfI, TfM);
					}
					cr[WS(rs, 62)] = FNMS(TfK, TfL, TfJ);
					Tfv = FMA(KP980785280, Tfo, Tfl);
					Tfp = FNMS(KP980785280, Tfo, Tfl);
					Tfi = W[91];
					Tf5 = W[90];
					{
					     E Tfr, Tfw, Tfq, Tfh;
					     Tfu = W[27];
					     Tfq = Tfi * Tfg;
					     Tfh = Tf5 * Tfg;
					     Tfr = W[26];
					     Tfw = Tfu * Tfs;
					     ci[WS(rs, 46)] = FMA(Tf5, Tfp, Tfq);
					     cr[WS(rs, 46)] = FNMS(Tfi, Tfp, Tfh);
					     Tft = Tfr * Tfs;
					     ci[WS(rs, 14)] = FMA(Tfr, Tfv, Tfw);
					}
				   }
			      }
			 }
			 {
			      E T89, T7X, T7W, T8a, T7D, T7C, T7B;
			      {
				   E T7t, T7w, T7Q, T7I, T7T, T7N, T7A, T7k;
				   {
					E T7b, T7i, T7L, T7H, T74, T7G, T7M, T7j;
					{
					     E T7u, T7v, T7p, T7s, T6W, T73;
					     T7p = FMA(KP923879532, T7o, T7n);
					     T89 = FNMS(KP923879532, T7o, T7n);
					     T7X = T7q + T7r;
					     T7s = T7q - T7r;
					     cr[WS(rs, 14)] = FNMS(Tfu, Tfv, Tft);
					     T7b = FNMS(KP534511135, T7a, T77);
					     T7u = FMA(KP534511135, T77, T7a);
					     T7v = FNMS(KP534511135, T7e, T7h);
					     T7i = FMA(KP534511135, T7h, T7e);
					     T7t = FMA(KP831469612, T7s, T7p);
					     T7L = FNMS(KP831469612, T7s, T7p);
					     T6W = FMA(KP923879532, T6V, T6U);
					     T7W = FNMS(KP923879532, T6V, T6U);
					     T8a = T72 + T6Z;
					     T73 = T6Z - T72;
					     T7H = T7v - T7u;
					     T7w = T7u + T7v;
					     T74 = FMA(KP831469612, T73, T6W);
					     T7G = FNMS(KP831469612, T73, T6W);
					}
					T7M = T7b - T7i;
					T7j = T7b + T7i;
					T7Q = FMA(KP881921264, T7H, T7G);
					T7I = FNMS(KP881921264, T7H, T7G);
					T7T = FMA(KP881921264, T7M, T7L);
					T7N = FNMS(KP881921264, T7M, T7L);
					T7A = FMA(KP881921264, T7j, T74);
					T7k = FNMS(KP881921264, T7j, T74);
				   }
				   {
					E T7S, T7R, T7F, T7K, T7x, T7m, T6T;
					T7F = W[104];
					T7K = W[105];
					{
					     E T7P, T7O, T7J, T7U;
					     T7P = W[40];
					     T7S = W[41];
					     T7O = T7F * T7N;
					     T7J = T7F * T7I;
					     T7U = T7P * T7T;
					     T7R = T7P * T7Q;
					     ci[WS(rs, 53)] = FMA(T7K, T7I, T7O);
					     cr[WS(rs, 53)] = FNMS(T7K, T7N, T7J);
					     ci[WS(rs, 21)] = FMA(T7S, T7Q, T7U);
					}
					cr[WS(rs, 21)] = FNMS(T7S, T7T, T7R);
					T7D = FMA(KP881921264, T7w, T7t);
					T7x = FNMS(KP881921264, T7w, T7t);
					T7m = W[73];
					T6T = W[72];
					{
					     E T7z, T7E, T7y, T7l;
					     T7C = W[9];
					     T7y = T7m * T7k;
					     T7l = T6T * T7k;
					     T7z = W[8];
					     T7E = T7C * T7A;
					     ci[WS(rs, 37)] = FMA(T6T, T7x, T7y);
					     cr[WS(rs, 37)] = FNMS(T7m, T7x, T7l);
					     T7B = T7z * T7A;
					     ci[WS(rs, 5)] = FMA(T7z, T7D, T7E);
					}
				   }
			      }
			      {
				   E T8u, T8t, T86, T8i, T8y, T8q, T8l, T8f;
				   {
					E T8d, T8c, T85, T8b, T7Y, T8o, T81, T84, T8p, T8e;
					T81 = FMA(KP303346683, T80, T7Z);
					T8d = FNMS(KP303346683, T7Z, T80);
					T8c = FMA(KP303346683, T82, T83);
					T84 = FNMS(KP303346683, T83, T82);