hb_64.c

最新的FFT程序

					T49 = T41 + T48;
					T4i = T4a - T4h;
					T4v = T4a + T4h;
					T6q = T49 - T4i;
					T4j = T49 + T4i;
					{
					     E Thp, TcY, Tff, ThA, Tdm, Tdn, Tfc, Tdl, Tfd, Tdd;
					     {
						  E Tdk, T6n, T4w, TcX, Td5, Tdc;
						  Thp = TcT - TcW;
						  TcX = TcT + TcW;
						  T6r = FMA(KP707106781, T6q, T6p);
						  T7d = FNMS(KP707106781, T6q, T6p);
						  T5s = FNMS(KP707106781, T4j, T40);
						  T4k = FMA(KP707106781, T4j, T40);
						  T6n = T4v - T4u;
						  T4w = T4u + T4v;
						  TcY = FMA(KP707106781, TcX, TcQ);
						  Tff = FNMS(KP707106781, TcX, TcQ);
						  T6o = FMA(KP707106781, T6n, T6m);
						  T7c = FNMS(KP707106781, T6n, T6m);
						  T5t = FNMS(KP707106781, T4w, T4t);
						  T4x = FMA(KP707106781, T4w, T4t);
						  Tdk = Tdi - Tdj;
						  ThA = Tdi + Tdj;
						  Tdm = FMA(KP414213562, Td1, Td4);
						  Td5 = FNMS(KP414213562, Td4, Td1);
						  Tdc = FMA(KP414213562, Tdb, Td8);
						  Tdn = FNMS(KP414213562, Td8, Tdb);
						  Tfc = FNMS(KP707106781, Tdk, Tdh);
						  Tdl = FMA(KP707106781, Tdk, Tdh);
						  Tfd = Tdc - Td5;
						  Tdd = Td5 + Tdc;
					     }
					     {
						  E Tfg, Tdo, Tht, Thw;
						  Tfe = FMA(KP923879532, Tfd, Tfc);
						  Tg3 = FNMS(KP923879532, Tfd, Tfc);
						  Tej = FNMS(KP923879532, Tdd, TcY);
						  Tde = FMA(KP923879532, Tdd, TcY);
						  Tfg = Tdm - Tdn;
						  Tdo = Tdm + Tdn;
						  Thq = FMA(KP707106781, Thp, Tho);
						  Tjr = FNMS(KP707106781, Thp, Tho);
						  Tfh = FMA(KP923879532, Tfg, Tff);
						  Tg2 = FNMS(KP923879532, Tfg, Tff);
						  Tei = FNMS(KP923879532, Tdo, Tdl);
						  Tdp = FMA(KP923879532, Tdo, Tdl);
						  ThC = FNMS(KP414213562, Thr, Ths);
						  Tht = FMA(KP414213562, Ths, Thr);
						  Thw = FNMS(KP414213562, Thv, Thu);
						  ThD = FMA(KP414213562, Thu, Thv);
						  Tjo = FNMS(KP707106781, ThA, Thz);
						  ThB = FMA(KP707106781, ThA, Thz);
						  Tjp = Tht - Thw;
						  Thx = Tht + Thw;
					     }
					}
				   }
			      }
			 }
		    }
	       }
	       {
		    E TiX, TiU, TiT, TeL, TeI, TeH;
		    {
			 E Tjq, Tkf, Tiv, Thy, Tjt, Tke, Tiu, ThF, T9F, T9C, T9x, T9G, T9B;
			 {
			      E Tav, Tas, Tan, Taw, Tar;
			      {
				   E Taa, TaY, Tb3, Tb0, TaX, Tb4;
				   {
					E T22, TaL, Tb2, T11, TaM, TaS, TaI, TaO, Tb1;
					{
					     E Tjs, ThE, TaH, TaE;
					     T22 = T1w + T21;
					     TaL = T21 - T1w;
					     Tjq = FMA(KP923879532, Tjp, Tjo);
					     Tkf = FNMS(KP923879532, Tjp, Tjo);
					     Tiv = FNMS(KP923879532, Thx, Thq);
					     Thy = FMA(KP923879532, Thx, Thq);
					     Tjs = ThD - ThC;
					     ThE = ThC + ThD;
					     Tb2 = TaF + TaG;
					     TaH = TaF - TaG;
					     TaE = Tv - T10;
					     T11 = Tv + T10;
					     Tjt = FMA(KP923879532, Tjs, Tjr);
					     Tke = FNMS(KP923879532, Tjs, Tjr);
					     Tiu = FNMS(KP923879532, ThE, ThB);
					     ThF = FMA(KP923879532, ThE, ThB);
					     Taa = Ta8 - Ta9;
					     TaM = Ta9 + Ta8;
					     TaY = T11 - T22;
					     TaS = TaE - TaH;
					     TaI = TaE + TaH;
					}
					rio[0] = T11 + T22;
					TaO = TaM - TaN;
					Tb1 = TaN + TaM;
					{
					     E TaK, TaD, TaP, TaV, TaQ, TaJ;
					     TaK = W[95];
					     TaD = W[94];
					     iio[-WS(ios, 63)] = Tb2 + Tb1;
					     TaP = TaL + TaO;
					     TaV = TaO - TaL;
					     TaQ = TaK * TaI;
					     TaJ = TaD * TaI;
					     {
						  E TaU, TaR, TaW, TaT;
						  TaU = W[31];
						  iio[-WS(ios, 15)] = FMA(TaD, TaP, TaQ);
						  rio[WS(ios, 48)] = FNMS(TaK, TaP, TaJ);
						  TaR = W[30];
						  TaW = TaU * TaS;
						  Tb3 = Tb1 - Tb2;
						  Tb0 = W[63];
						  TaT = TaR * TaS;
						  iio[-WS(ios, 47)] = FMA(TaR, TaV, TaW);
						  TaX = W[62];
						  Tb4 = Tb0 * TaY;
						  rio[WS(ios, 16)] = FNMS(TaU, TaV, TaT);
					     }
					}
				   }
				   {
					E Tao, T9S, T9X, Tab, Tat, Ta2, Tap, Tae, TaZ, Tac, Tad;
					TaZ = TaX * TaY;
					iio[-WS(ios, 31)] = FMA(TaX, Tb3, Tb4);
					Tao = T9O - T9R;
					T9S = T9O + T9R;
					rio[WS(ios, 32)] = FNMS(Tb0, Tb3, TaZ);
					Tac = T9W - T9T;
					T9X = T9T + T9W;
					Tab = Ta7 + Taa;
					Tat = Taa - Ta7;
					Ta2 = T9Y - Ta1;
					Tad = T9Y + Ta1;
					Tap = Tad - Tac;
					Tae = Tac + Tad;
					{
					     E Tal, Taq, Tai, TaB, TaA, Taz;
					     {
						  E Ta6, Tau, Taf, Tay, Ta4, T9N, Ta3;
						  Ta6 = W[111];
						  Tau = T9X - Ta2;
						  Ta3 = T9X + Ta2;
						  Taf = FMA(KP707106781, Tae, Tab);
						  Tal = FNMS(KP707106781, Tae, Tab);
						  Tay = FNMS(KP707106781, Tap, Tao);
						  Taq = FMA(KP707106781, Tap, Tao);
						  Ta4 = FMA(KP707106781, Ta3, T9S);
						  Tai = FNMS(KP707106781, Ta3, T9S);
						  T9N = W[110];
						  {
						       E Tax, TaC, Tag, Ta5;
						       TaB = FNMS(KP707106781, Tau, Tat);
						       Tav = FMA(KP707106781, Tau, Tat);
						       TaA = W[79];
						       Tag = Ta6 * Ta4;
						       Ta5 = T9N * Ta4;
						       Tax = W[78];
						       TaC = TaA * Tay;
						       iio[-WS(ios, 7)] = FMA(T9N, Taf, Tag);
						       rio[WS(ios, 56)] = FNMS(Ta6, Taf, Ta5);
						       Taz = Tax * Tay;
						       iio[-WS(ios, 23)] = FMA(Tax, TaB, TaC);
						  }
					     }
					     {
						  E Tak, Tah, Tam, Taj;
						  Tak = W[47];
						  rio[WS(ios, 40)] = FNMS(TaA, TaB, Taz);
						  Tah = W[46];
						  Tam = Tak * Tai;
						  Tas = W[15];
						  Taj = Tah * Tai;
						  Tan = W[14];
						  iio[-WS(ios, 39)] = FMA(Tah, Tal, Tam);
						  Taw = Tas * Taq;
						  rio[WS(ios, 24)] = FNMS(Tak, Tal, Taj);
						  Tar = Tan * Taq;
					     }
					}
				   }
			      }
			      {
				   E T9j, T97, T9k, T8X, T8U, T8P, T8Y, T8T;
				   {
					E T8E, T8f, T8D, T8V, T84, T8Q, T8q, T8F, T8z, T8C, T83;
					T9j = T8y - T8v;
					T8z = T8v + T8y;
					iio[-WS(ios, 55)] = FMA(Tan, Tav, Taw);
					T97 = T8B - T8A;
					T8C = T8A + T8B;
					rio[WS(ios, 8)] = FNMS(Tas, Tav, Tar);
					T9k = T7X - T82;
					T83 = T7X + T82;
					T8E = FNMS(KP414213562, T89, T8e);
					T8f = FMA(KP414213562, T8e, T89);
					T8D = FMA(KP707106781, T8C, T8z);
					T8V = FNMS(KP707106781, T8C, T8z);
					T84 = FMA(KP707106781, T83, T7S);
					T8Q = FNMS(KP707106781, T83, T7S);
					T8q = FNMS(KP414213562, T8p, T8k);
					T8F = FMA(KP414213562, T8k, T8p);
					{
					     E T8N, T8S, T8K, T93, T92, T91;
					     {
						  E T8u, T8W, T8H, T90, T8s, T7N, T8R, T8G, T8r;
						  T8u = W[119];
						  T8R = T8F - T8E;
						  T8G = T8E + T8F;
						  T8r = T8f + T8q;
						  T8W = T8f - T8q;
						  T8N = FNMS(KP923879532, T8G, T8D);
						  T8H = FMA(KP923879532, T8G, T8D);
						  T8S = FMA(KP923879532, T8R, T8Q);
						  T90 = FNMS(KP923879532, T8R, T8Q);
						  T8K = FNMS(KP923879532, T8r, T84);
						  T8s = FMA(KP923879532, T8r, T84);
						  T7N = W[118];
						  {
						       E T8Z, T94, T8I, T8t;
						       T8X = FMA(KP923879532, T8W, T8V);
						       T93 = FNMS(KP923879532, T8W, T8V);
						       T92 = W[87];
						       T8I = T8u * T8s;
						       T8t = T7N * T8s;
						       T8Z = W[86];
						       T94 = T92 * T90;
						       iio[-WS(ios, 3)] = FMA(T7N, T8H, T8I);
						       rio[WS(ios, 60)] = FNMS(T8u, T8H, T8t);
						       T91 = T8Z * T90;
						       iio[-WS(ios, 19)] = FMA(T8Z, T93, T94);
						  }
					     }
					     {
						  E T8M, T8J, T8O, T8L;
						  T8M = W[55];
						  rio[WS(ios, 44)] = FNMS(T92, T93, T91);
						  T8J = W[54];
						  T8O = T8M * T8K;
						  T8U = W[23];
						  T8L = T8J * T8K;
						  T8P = W[22];
						  iio[-WS(ios, 35)] = FMA(T8J, T8N, T8O);
						  T8Y = T8U * T8S;
						  rio[WS(ios, 28)] = FNMS(T8M, T8N, T8L);
						  T8T = T8P * T8S;
					     }
					}
				   }
				   {
					E T9b, T9m, T9D, T9l, T9y, T98, T9n, T9e;
					iio[-WS(ios, 51)] = FMA(T8P, T8X, T8Y);
					rio[WS(ios, 12)] = FNMS(T8U, T8X, T8T);
					T9b = FNMS(KP414213562, T9a, T99);
					T9m = FMA(KP414213562, T99, T9a);
					T9D = FNMS(KP707106781, T9k, T9j);
					T9l = FMA(KP707106781, T9k, T9j);
					T9y = FNMS(KP707106781, T97, T96);
					T98 = FMA(KP707106781, T97, T96);
					T9n = FNMS(KP414213562, T9c, T9d);
					T9e = FMA(KP414213562, T9d, T9c);
					{
					     E T9v, T9A, T9s, T9L, T9K, T9J;
					     {
						  E T9i, T9E, T9p, T9I, T9g, T95, T9o, T9z, T9f;
						  T9i = W[7];
						  T9o = T9m + T9n;
						  T9z = T9m - T9n;
						  T9E = T9e - T9b;
						  T9f = T9b + T9e;
						  T9p = FMA(KP923879532, T9o, T9l);
						  T9v = FNMS(KP923879532, T9o, T9l);
						  T9I = FNMS(KP923879532, T9z, T9y);
						  T9A = FMA(KP923879532, T9z, T9y);
						  T9g = FMA(KP923879532, T9f, T98);
						  T9s = FNMS(KP923879532, T9f, T98);
						  T95 = W[6];
						  {
						       E T9H, T9M, T9q, T9h;
						       T9L = FNMS(KP923879532, T9E, T9D);
						       T9F = FMA(KP923879532, T9E, T9D);
						       T9K = W[39];
						       T9q = T9i * T9g;
						       T9h = T95 * T9g;
						       T9H = W[38];
						       T9M = T9K * T9I;
						       iio[-WS(ios, 59)] = FMA(T95, T9p, T9q);
						       rio[WS(ios, 4)] = FNMS(T9i, T9p, T9h);
						       T9J = T9H * T9I;
						       iio[-WS(ios, 43)] = FMA(T9H, T9L, T9M);
						  }
					     }
					     {
						  E T9u, T9r, T9w, T9t;
						  T9u = W[71];
						  rio[WS(ios, 20)] = FNMS(T9K, T9L, T9J);
						  T9r = W[70];
						  T9w = T9u * T9s;
						  T9C = W[103];
						  T9t = T9r * T9s;
						  T9x = W[102];
						  iio[-WS(ios, 27)] = FMA(T9r, T9v, T9w);
						  T9G = T9C * T9A;
						  rio[WS(ios, 36)] = FNMS(T9u, T9v, T9t);
						  T9B = T9x * T9A;
					     }
					}
				   }
			      }
			 }
			 {
			      E TjB, Tjf, Tj8, TjE, TiB, Tip, Tio, TiC, Tif, Tic, Tib, TkH, TkE, TkD, TjZ;
			      E TjW, TjV;
			      {
				   E TkK, Tkx, Tki, TkC, TkN, Tku, Tkr;
				   {
					E Tkd, Tkg, Tka, TkF, TkG, Tkq, Tkl, Tkm;
					{
					     E Tko, Tkp, Tk8, Tk9;
					     TjB = FMA(KP923879532, TjA, Tjz);
					     Tk8 = FNMS(KP923879532, TjA, Tjz);
					     iio[-WS(ios, 11)] = FMA(T9x, T9F, T9G);
					     Tjf = Tjb + Tje;
					     Tk9 = Tje - Tjb;
					     rio[WS(ios, 52)] = FNMS(T9C, T9F, T9B);
					     Tkd = FNMS(KP534511135, Tkc, Tkb);
					     Tko = FMA(KP534511135, Tkb, Tkc);
					     Tkp = FMA(KP534511135, Tke, Tkf);
					     Tkg = FNMS(KP534511135, Tkf, Tke);
					     Tka = FMA(KP831469612, Tk9, Tk8);
					     TkF = FNMS(KP831469612, Tk9, Tk8);
					     TkG = Tko - Tkp;
					     Tkq = Tko + Tkp;
					     Tj8 = FNMS(KP923879532, Tj7, Tj6);
					     Tkl = FMA(KP923879532, Tj7, Tj6);
					     Tkm = TjC + TjD;
					     TjE = TjC - TjD;
					}
					{
					     E TkB, Tkh, TkA, Tkn;
					     TkH = FMA(KP881921264, TkG, TkF);
					     TkK = FNMS(KP881921264, TkG, TkF);
					     TkB = Tkd + Tkg;
					     Tkh = Tkd - Tkg;
					     TkA = FNMS(KP831469612, Tkm, Tkl);
					     Tkn = FMA(KP831469612, Tkm, Tkl);
					     Tkx = FNMS(KP881921264, Tkh, Tka);
					     Tki = FMA(KP881921264, Tkh, Tka);
					     TkC = FNMS(KP881921264, TkB, TkA);
					     TkN = FMA(KP881921264, TkB, TkA);
					     Tku = FNMS(KP881921264, Tkq, Tkn);
					     Tkr = FMA(KP881921264, Tkq, Tkn);
					}
				   }
				   {
					E TkM, TkL, Tk7, Tkk, Tkt, Tkw;
					Tk7 = W[116];
					Tkk = W[117];
					{
					     E TkJ, Tks, Tkj, TkO;
					     TkJ = W[84];
					     TkM = W[85];
					     Tks = Tk7 * Tkr;
					     Tkj = Tk7 * Tki;
					     TkO = TkJ * TkN;
					     TkL = TkJ * TkK;
					     rio[WS(ios, 59)] = FNMS(Tkk, Tki, Tks);
					     iio[-WS(ios, 4)] = FMA(Tkk, Tkr, Tkj);
					     rio[WS(ios, 43)] = FNMS(TkM, TkK, TkO);
					}
					iio[-WS(ios, 20)] = FMA(TkM, TkN, TkL);
					Tkt = W[52];
					Tkw = W[53];
					{
					     E Tkz, Tky, Tkv, TkI;
					     Tkz = W[20];
					     TkE = W[21];
					     Tky = Tkt * Tkx;
					     Tkv = Tkt * Tku;
					     TkI = Tkz * TkH;
					     TkD = Tkz * TkC;
					     iio[-WS(ios, 36)] = FMA(Tkw, Tku, Tky);
					     rio[WS(ios, 27)] = FNMS(Tkw, Tkx, Tkv);
					     iio[-WS(ios, 52)] = FMA(TkE, TkC, TkI);
					}
				   }
			      }
			      rio[WS(ios, 11)] = FNMS(TkE, TkH, TkD);
			      {
				   E Tii, Ti5, ThI, Tia, Til, Ti2, ThZ;
				   {
					E Thn, ThG, Tid, Th4, Tie, ThY, ThR, ThU;
					{
					     E TgO, Th3, ThW, ThX;
					     TiB = FNMS(KP923879532, TgN, TgG);
					     TgO = FMA(KP923879532, TgN, TgG);
					     Th3 = TgV - Th2;
					     Tip = TgV + Th2;
					     Thn = FNMS(KP098491403, Thm, Thf);
					     ThW = FMA(KP098491403, Thf, Thm);
					     ThX = FMA(KP098491403, Thy, ThF);
					     ThG = FNMS(KP098491403, ThF, Thy);
					     Tid = FNMS(KP980785280, Th3, TgO);
					     Th4 = FMA(KP980785280, Th3, TgO);
					     Tie = ThW - ThX;
					     ThY = ThW + ThX;
					     Tio = FNMS(KP923879532, ThQ, ThN);
					     ThR = FMA(KP923879532, ThQ, ThN);
					     ThU = ThS + ThT;
					     TiC = ThS - ThT;
					}
					{
					     E Ti9, ThH, Ti8, ThV;
					     Tif = FMA(KP995184726, Tie, Tid);
					     Tii = FNMS(KP995184726, Tie, Tid);
					     Ti9 = Thn + ThG;
					     ThH = Thn - ThG;
					     Ti8 = FNMS(KP980785280, ThU, ThR);
					     ThV = FMA(KP980785280, ThU, ThR);
					     Ti5 = FNMS(KP995184726, ThH, Th4);
					     ThI = FMA(KP995184726, ThH, Th4);
					     Tia = FNMS(KP995184726, Ti9, Ti8);
					     Til = FMA(KP995184726, Ti9, Ti8);
					     Ti2 = FNMS(KP995184726, ThY, ThV);
					     ThZ = FMA(KP995184726, ThY, ThV);
					}
				   }
				   {
					E Tik, Tij, TgD, ThK, Ti1, Ti4;
					TgD = W[124];
					ThK = W[125];
					{
					     E Tih, Ti0, ThJ, Tim;
					     Tih = W[92];
					     Tik = W[93];
					     Ti0 = TgD * ThZ;
					     ThJ = TgD * ThI;
					     Tim = Tih * Til;