hb_64.c

快速fft变换

					     Tjh = ThR - ThS;
					     Tjl = ThV - ThU;
					     ThW = ThU + ThV;
					     Tkb = Tji + Tjh;
					     Tjj = Tjh - Tji;
					     Tkc = Tjk + Tjl;
					     Tjm = Tjk - Tjl;
					}
					{
					     E TfZ, Tdz, Tg2, TdQ, TdS, TdR, Tg3, TdK, TdE, TdJ;
					     TfZ = Tdy + Tdx;
					     Tdz = Tdx - Tdy;
					     Tg2 = TdM + TdP;
					     TdQ = TdM - TdP;
					     TdS = TdD - TdA;
					     TdE = TdA + TdD;
					     TiL = ThT + ThW;
					     ThX = ThT - ThW;
					     TdJ = TdF - TdI;
					     TdR = TdF + TdI;
					     Tg3 = TdJ + TdE;
					     TdK = TdE - TdJ;
					     {
						  E T9A, T3T, T7f, T9p, T4r, T4s, T7c, T4q, T7d, T4g;
						  {
						       E T4p, Tg0, TdT, T3S, T44, T4f;
						       T9A = T3M + T3R;
						       T3S = T3M - T3R;
						       Tg4 = FNMS(KP707106781, Tg3, Tg2);
						       TgW = FMA(KP707106781, Tg3, Tg2);
						       Tf9 = FMA(KP707106781, TdK, Tdz);
						       TdL = FNMS(KP707106781, TdK, Tdz);
						       Tg0 = TdR + TdS;
						       TdT = TdR - TdS;
						       T3T = FMA(KP707106781, T3S, T3H);
						       T7f = FNMS(KP707106781, T3S, T3H);
						       Tg1 = FNMS(KP707106781, Tg0, TfZ);
						       TgX = FMA(KP707106781, Tg0, TfZ);
						       Tfa = FMA(KP707106781, TdT, TdQ);
						       TdU = FNMS(KP707106781, TdT, TdQ);
						       T4p = T4n + T4o;
						       T9p = T4n - T4o;
						       T4r = FNMS(KP414213562, T3Y, T43);
						       T44 = FMA(KP414213562, T43, T3Y);
						       T4f = FNMS(KP414213562, T4e, T49);
						       T4s = FMA(KP414213562, T49, T4e);
						       T7c = FNMS(KP707106781, T4p, T4m);
						       T4q = FMA(KP707106781, T4p, T4m);
						       T7d = T44 - T4f;
						       T4g = T44 + T4f;
						  }
						  {
						       E T7g, T4t, T9t, T9w;
						       T7e = FNMS(KP923879532, T7d, T7c);
						       T80 = FMA(KP923879532, T7d, T7c);
						       T6f = FMA(KP923879532, T4g, T3T);
						       T4h = FNMS(KP923879532, T4g, T3T);
						       T7g = T4s - T4r;
						       T4t = T4r + T4s;
						       T9q = FNMS(KP707106781, T9p, T9o);
						       Tbr = FMA(KP707106781, T9p, T9o);
						       T7h = FNMS(KP923879532, T7g, T7f);
						       T7Z = FMA(KP923879532, T7g, T7f);
						       T6g = FMA(KP923879532, T4t, T4q);
						       T4u = FNMS(KP923879532, T4t, T4q);
						       T9D = FNMS(KP414213562, T9r, T9s);
						       T9t = FMA(KP414213562, T9s, T9r);
						       T9w = FNMS(KP414213562, T9v, T9u);
						       T9C = FMA(KP414213562, T9u, T9v);
						       Tbo = FMA(KP707106781, T9A, T9z);
						       T9B = FNMS(KP707106781, T9A, T9z);
						       Tbp = T9w + T9t;
						       T9x = T9t - T9w;
						  }
					     }
					}
				   }
			      }
			 }
		    }
	       }
	       {
		    E Tbq, Tcf, Tav, T9y, Tbt, Tce, Tau, T9F, T6p, T6d, T6c, T6q, Thf, The, Thd;
		    {
			 E Tk9, Tkm, TjP, TjO, TjN;
			 {
			      E Tj0, TiS, TiU, Tj3, Tj1, Tj4, TiY, Tj2;
			      {
				   E TiQ, TiW, TiV, TiR, TiD, TiG, TiN, TiF, TiO;
				   {
					E T11, T22, TiJ, TiE, TiM, Tbs, T9E;
					TiQ = Tv - T10;
					T11 = Tv + T10;
					Tbq = FNMS(KP923879532, Tbp, Tbo);
					Tcf = FMA(KP923879532, Tbp, Tbo);
					Tav = FMA(KP923879532, T9x, T9q);
					T9y = FNMS(KP923879532, T9x, T9q);
					Tbs = T9C + T9D;
					T9E = T9C - T9D;
					T22 = T1w + T21;
					TiW = T1w - T21;
					TiV = TiH - TiI;
					TiJ = TiH + TiI;
					Tbt = FNMS(KP923879532, Tbs, Tbr);
					Tce = FMA(KP923879532, Tbs, Tbr);
					Tau = FMA(KP923879532, T9E, T9B);
					T9F = FNMS(KP923879532, T9E, T9B);
					TiE = T11 - T22;
					TiR = TiL - TiK;
					TiM = TiK + TiL;
					cr[0] = T11 + T22;
					TiD = W[62];
					TiG = W[63];
					ci[0] = TiJ + TiM;
					TiN = TiJ - TiM;
					TiF = TiD * TiE;
					TiO = TiG * TiE;
				   }
				   cr[WS(rs, 32)] = FNMS(TiG, TiN, TiF);
				   ci[WS(rs, 32)] = FMA(TiD, TiN, TiO);
				   Tj0 = TiQ + TiR;
				   TiS = TiQ - TiR;
				   {
					E TiP, TiX, TiT, TiZ;
					TiP = W[94];
					TiU = W[95];
					TiZ = W[30];
					Tj3 = TiW + TiV;
					TiX = TiV - TiW;
					TiT = TiP * TiS;
					Tj1 = TiZ * Tj0;
					Tj4 = TiZ * Tj3;
					TiY = TiP * TiX;
					cr[WS(rs, 48)] = FNMS(TiU, TiX, TiT);
					Tj2 = W[31];
				   }
			      }
			      {
				   E Tii, Til, Tik, Tih, Tim;
				   {
					E Tib, Tit, Tio, ThG, ThP, ThY, Tie, Tip, Tic, Tid;
					Tib = Ti3 + Tia;
					Tit = Tia - Ti3;
					ci[WS(rs, 48)] = FMA(TiU, TiS, TiY);
					Tio = Thy - ThF;
					ThG = Thy + ThF;
					ci[WS(rs, 16)] = FMA(Tj2, Tj0, Tj4);
					cr[WS(rs, 16)] = FNMS(Tj2, Tj3, Tj1);
					ThP = ThH - ThO;
					Tic = ThH + ThO;
					Tid = ThX - ThQ;
					ThY = ThQ + ThX;
					Tie = Tic + Tid;
					Tip = Tid - Tic;
					{
					     E Tiy, TiB, Ti0, Tiz, TiC, TiA;
					     {
						  E Tin, Tis, Tiq, ThZ, Tiu, Tir, Tiw, Tix, Tiv;
						  Tin = W[110];
						  Tis = W[111];
						  Tiy = FMA(KP707106781, Tip, Tio);
						  Tiq = FNMS(KP707106781, Tip, Tio);
						  ThZ = ThP + ThY;
						  Tiu = ThP - ThY;
						  Tir = Tin * Tiq;
						  Tix = W[46];
						  TiB = FMA(KP707106781, Tiu, Tit);
						  Tiv = FNMS(KP707106781, Tiu, Tit);
						  Ti0 = FNMS(KP707106781, ThZ, ThG);
						  Tii = FMA(KP707106781, ThZ, ThG);
						  cr[WS(rs, 56)] = FNMS(Tis, Tiv, Tir);
						  Tiw = Tin * Tiv;
						  Tiz = Tix * Tiy;
						  TiC = Tix * TiB;
						  TiA = W[47];
						  ci[WS(rs, 56)] = FMA(Tis, Tiq, Tiw);
					     }
					     {
						  E Tif, Ti2, Thx, Tig, Ti1;
						  Til = FMA(KP707106781, Tie, Tib);
						  Tif = FNMS(KP707106781, Tie, Tib);
						  Ti2 = W[79];
						  ci[WS(rs, 24)] = FMA(TiA, Tiy, TiC);
						  cr[WS(rs, 24)] = FNMS(TiA, TiB, Tiz);
						  Thx = W[78];
						  Tig = Ti2 * Ti0;
						  Tik = W[15];
						  Ti1 = Thx * Ti0;
						  ci[WS(rs, 40)] = FMA(Thx, Tif, Tig);
						  Tih = W[14];
						  Tim = Tik * Tii;
						  cr[WS(rs, 40)] = FNMS(Ti2, Tif, Ti1);
					     }
					}
				   }
				   {
					E TjF, TjI, TjU, Tk2, TjZ, Tk5, Tjw, TjM;
					{
					     E TjX, TjG, Tju, Tjg, TjS, Tjn, TjH, Tjf, TjE, Tij, TjT, Tjv, TjY;
					     TjE = TjC - TjD;
					     Tk9 = TjC + TjD;
					     Tij = Tih * Tii;
					     ci[WS(rs, 8)] = FMA(Tih, Til, Tim);
					     Tkm = Tje + Tjb;
					     Tjf = Tjb - Tje;
					     TjX = FNMS(KP707106781, TjE, TjB);
					     TjF = FMA(KP707106781, TjE, TjB);
					     cr[WS(rs, 8)] = FNMS(Tik, Til, Tij);
					     TjG = FMA(KP414213562, Tjq, Tjt);
					     Tju = FNMS(KP414213562, Tjt, Tjq);
					     Tjg = FMA(KP707106781, Tjf, Tj8);
					     TjS = FNMS(KP707106781, Tjf, Tj8);
					     Tjn = FMA(KP414213562, Tjm, Tjj);
					     TjH = FNMS(KP414213562, Tjj, Tjm);
					     TjI = TjG - TjH;
					     TjT = TjG + TjH;
					     Tjv = Tjn - Tju;
					     TjY = Tju + Tjn;
					     TjU = FNMS(KP923879532, TjT, TjS);
					     Tk2 = FMA(KP923879532, TjT, TjS);
					     TjZ = FNMS(KP923879532, TjY, TjX);
					     Tk5 = FMA(KP923879532, TjY, TjX);
					     Tjw = FNMS(KP923879532, Tjv, Tjg);
					     TjM = FMA(KP923879532, Tjv, Tjg);
					}
					{
					     E Tk4, Tk3, TjR, TjW, TjJ, Tjy, Tj5;
					     TjR = W[54];
					     TjW = W[55];
					     {
						  E Tk1, Tk0, TjV, Tk6;
						  Tk1 = W[118];
						  Tk4 = W[119];
						  Tk0 = TjR * TjZ;
						  TjV = TjR * TjU;
						  Tk6 = Tk1 * Tk5;
						  Tk3 = Tk1 * Tk2;
						  ci[WS(rs, 28)] = FMA(TjW, TjU, Tk0);
						  cr[WS(rs, 28)] = FNMS(TjW, TjZ, TjV);
						  ci[WS(rs, 60)] = FMA(Tk4, Tk2, Tk6);
					     }
					     cr[WS(rs, 60)] = FNMS(Tk4, Tk5, Tk3);
					     TjP = FMA(KP923879532, TjI, TjF);
					     TjJ = FNMS(KP923879532, TjI, TjF);
					     Tjy = W[87];
					     Tj5 = W[86];
					     {
						  E TjL, TjQ, TjK, Tjx;
						  TjO = W[23];
						  TjK = Tjy * Tjw;
						  Tjx = Tj5 * Tjw;
						  TjL = W[22];
						  TjQ = TjO * TjM;
						  ci[WS(rs, 44)] = FMA(Tj5, TjJ, TjK);
						  cr[WS(rs, 44)] = FNMS(Tjy, TjJ, Tjx);
						  TjN = TjL * TjM;
						  ci[WS(rs, 12)] = FMA(TjL, TjP, TjQ);
					     }
					}
				   }
			      }
			 }
			 {
			      E T5T, T5S, T5R, Tkx, Tkw, Tkv;
			      {
				   E Tkn, Tkq, TkC, TkK, TkH, TkN, Tki, Tku;
				   {
					E Tkg, Tko, TkF, Tka, TkA, Tkd, Tkp, TkB, Tkh, TkG;
					cr[WS(rs, 12)] = FNMS(TjO, TjP, TjN);
					Tkg = FMA(KP414213562, Tkf, Tke);
					Tko = FNMS(KP414213562, Tke, Tkf);
					TkF = FMA(KP707106781, Tkm, Tkl);
					Tkn = FNMS(KP707106781, Tkm, Tkl);
					Tka = FNMS(KP707106781, Tk9, Tk8);
					TkA = FMA(KP707106781, Tk9, Tk8);
					Tkd = FNMS(KP414213562, Tkc, Tkb);
					Tkp = FMA(KP414213562, Tkb, Tkc);
					Tkq = Tko - Tkp;
					TkB = Tko + Tkp;
					Tkh = Tkd - Tkg;
					TkG = Tkg + Tkd;
					TkC = FNMS(KP923879532, TkB, TkA);
					TkK = FMA(KP923879532, TkB, TkA);
					TkH = FNMS(KP923879532, TkG, TkF);
					TkN = FMA(KP923879532, TkG, TkF);
					Tki = FNMS(KP923879532, Tkh, Tka);
					Tku = FMA(KP923879532, Tkh, Tka);
				   }
				   {
					E TkM, TkL, Tkz, TkE, Tkr, Tkk, Tk7;
					Tkz = W[70];
					TkE = W[71];
					{
					     E TkJ, TkI, TkD, TkO;
					     TkJ = W[6];
					     TkM = W[7];
					     TkI = Tkz * TkH;
					     TkD = Tkz * TkC;
					     TkO = TkJ * TkN;
					     TkL = TkJ * TkK;
					     ci[WS(rs, 36)] = FMA(TkE, TkC, TkI);
					     cr[WS(rs, 36)] = FNMS(TkE, TkH, TkD);
					     ci[WS(rs, 4)] = FMA(TkM, TkK, TkO);
					}
					cr[WS(rs, 4)] = FNMS(TkM, TkN, TkL);
					Tkx = FMA(KP923879532, Tkq, Tkn);
					Tkr = FNMS(KP923879532, Tkq, Tkn);
					Tkk = W[103];
					Tk7 = W[102];
					{
					     E Tkt, Tky, Tks, Tkj;
					     Tkw = W[39];
					     Tks = Tkk * Tki;
					     Tkj = Tk7 * Tki;
					     Tkt = W[38];
					     Tky = Tkw * Tku;
					     ci[WS(rs, 52)] = FMA(Tk7, Tkr, Tks);
					     cr[WS(rs, 52)] = FNMS(Tkk, Tkr, Tkj);
					     Tkv = Tkt * Tku;
					     ci[WS(rs, 20)] = FMA(Tkt, Tkx, Tky);
					}
				   }
			      }
			      {
				   E T5J, T5M, T66, T5Y, T69, T63, T5Q, T5q;
				   {
					E T5o, T4v, T61, T5X, T3C, T5W, T62, T5p;
					{
					     E T5K, T5L, T5F, T5I, T2I, T3B;
					     T5F = FNMS(KP923879532, T5E, T5B);
					     T6p = FMA(KP923879532, T5E, T5B);
					     T6d = T5G + T5H;
					     T5I = T5G - T5H;
					     cr[WS(rs, 20)] = FNMS(Tkw, Tkx, Tkv);
					     T5o = FNMS(KP820678790, T5n, T5a);
					     T5K = FMA(KP820678790, T5a, T5n);
					     T5L = FNMS(KP820678790, T4h, T4u);
					     T4v = FMA(KP820678790, T4u, T4h);
					     T5J = FMA(KP980785280, T5I, T5F);
					     T61 = FNMS(KP980785280, T5I, T5F);
					     T2I = FNMS(KP923879532, T2H, T2k);
					     T6c = FMA(KP923879532, T2H, T2k);
					     T6q = T3A + T39;
					     T3B = T39 - T3A;
					     T5X = T5K + T5L;
					     T5M = T5K - T5L;
					     T3C = FMA(KP980785280, T3B, T2I);
					     T5W = FNMS(KP980785280, T3B, T2I);
					}
					T62 = T5o + T4v;
					T5p = T4v - T5o;
					T66 = FMA(KP773010453, T5X, T5W);
					T5Y = FNMS(KP773010453, T5X, T5W);
					T69 = FMA(KP773010453, T62, T61);
					T63 = FNMS(KP773010453, T62, T61);
					T5Q = FMA(KP773010453, T5p, T3C);
					T5q = FNMS(KP773010453, T5p, T3C);
				   }
				   {
					E T68, T67, T5V, T60, T5N, T5s, T23;
					T5V = W[48];
					T60 = W[49];
					{
					     E T65, T64, T5Z, T6a;
					     T65 = W[112];
					     T68 = W[113];
					     T64 = T5V * T63;
					     T5Z = T5V * T5Y;
					     T6a = T65 * T69;
					     T67 = T65 * T66;
					     ci[WS(rs, 25)] = FMA(T60, T5Y, T64);
					     cr[WS(rs, 25)] = FNMS(T60, T63, T5Z);
					     ci[WS(rs, 57)] = FMA(T68, T66, T6a);
					}
					cr[WS(rs, 57)] = FNMS(T68, T69, T67);
					T5T = FMA(KP773010453, T5M, T5J);
					T5N = FNMS(KP773010453, T5M, T5J);
					T5s = W[81];
					T23 = W[80];
					{
					     E T5P, T5U, T5O, T5r;
					     T5S = W[17];
					     T5O = T5s * T5q;
					     T5r = T23 * T5q;
					     T5P = W[16];
					     T5U = T5S * T5Q;
					     ci[WS(rs, 41)] = FMA(T23, T5N, T5O);
					     cr[WS(rs, 41)] = FNMS(T5s, T5N, T5r);
					     T5R = T5P * T5Q;
					     ci[WS(rs, 9)] = FMA(T5P, T5T, T5U);
					}
				   }
			      }
			      {
				   E Th3, TgR, TgQ, Th4, TgN, TgM, TgL;
				   {
					E TgG, TgF, Tge, Tgu, TgK, TgC, Tgx, Tgr;
					{
					     E Tgp, Tgo, Tgd, Tgn, TfY, TgA, TgB, Tgq;
					     {
						  E Tgj, Tgm, Tg5, Tgc, TfQ, TfX;
						  Tg5 = FMA(KP668178637, Tg4, Tg1);
						  Tgp = FNMS(KP668178637, Tg1, Tg4);
						  Tgo = FMA(KP668178637, Tg8, Tgb);
						  Tgc = FNMS(KP668178637, Tgb, Tg8);
						  cr[WS(rs, 9)] = FNMS(T5S, T5T, T5R);
						  Th3 = FMA(KP707106781, Tgi, Tgh);
						  Tgj = FNMS(KP707106781, Tgi, Tgh);
						  Tgm = Tgk - Tgl;
						  TgR = Tgk + Tgl;
						  TgG = Tgc + Tg5;
						  Tgd = Tg5 - Tgc;
						  TfQ = FNMS(KP707106781, TfP, TfO);
						  TgQ = FMA(KP707106781, TfP, TfO);
						  Th4 = TfW + TfT;
						  TfX = TfT - TfW;
						  Tgn = FMA(KP923879532, Tgm, Tgj);
						  TgF = FNMS(KP923879532, Tgm, Tgj);
						  TfY = FMA(KP923879532, TfX, TfQ);
						  TgA = FNMS(KP923879532, TfX, TfQ);
					     }
					     TgB = Tgo + Tgp;
					     Tgq = Tgo - Tgp;
					     Tge = FNMS(KP831469612, Tgd, TfY);
					     Tgu = FMA(KP831469612, Tgd, TfY);
					     TgK = FMA(KP831469612, TgB, TgA);
					     TgC = FNMS(KP831469612, TgB, TgA);
					     Tgx = FMA(KP831469612, Tgq, Tgn);
					     Tgr = FNMS(KP831469612, Tgq, Tgn);
					}
					{
					     E Tgw, Tgv, TfN, Tgg, TgH, TgE, Tgz;
					     TfN = W[82];
					     Tgg = W[83];
					     {
						  E Tgt, Tgs, Tgf, Tgy;
						  Tgt = W[18];
						  Tgw = W[19];
						  Tgs = TfN * Tgr;
						  Tgf = TfN * Tge;
						  Tgy = Tgt * Tgx;