hb_64.c

最新的FFT程序

					     Tij = Tih * Tii;
					     rio[WS(ios, 63)] = FNMS(ThK, ThI, Ti0);
					     iio[0] = FMA(ThK, ThZ, ThJ);
					     rio[WS(ios, 47)] = FNMS(Tik, Tii, Tim);
					}
					iio[-WS(ios, 16)] = FMA(Tik, Til, Tij);
					Ti1 = W[60];
					Ti4 = W[61];
					{
					     E Ti7, Ti6, Ti3, Tig;
					     Ti7 = W[28];
					     Tic = W[29];
					     Ti6 = Ti1 * Ti5;
					     Ti3 = Ti1 * Ti2;
					     Tig = Ti7 * Tif;
					     Tib = Ti7 * Tia;
					     iio[-WS(ios, 32)] = FMA(Ti4, Ti2, Ti6);
					     rio[WS(ios, 31)] = FNMS(Ti4, Ti5, Ti3);
					     iio[-WS(ios, 48)] = FMA(Tic, Tia, Tig);
					}
				   }
			      }
			      rio[WS(ios, 15)] = FNMS(Tic, Tif, Tib);
			      {
				   E Tk2, TjP, Tjw, TjU, Tk5, TjM, TjJ;
				   {
					E Tjn, Tju, TjX, Tjg, TjY, TjI, TjG, TjH;
					Tjn = FNMS(KP303346683, Tjm, Tjj);
					TjG = FMA(KP303346683, Tjj, Tjm);
					TjH = FMA(KP303346683, Tjq, Tjt);
					Tju = FNMS(KP303346683, Tjt, Tjq);
					TjX = FNMS(KP831469612, Tjf, Tj8);
					Tjg = FMA(KP831469612, Tjf, Tj8);
					TjY = TjG + TjH;
					TjI = TjG - TjH;
					{
					     E TjT, Tjv, TjS, TjF;
					     TjZ = FMA(KP956940335, TjY, TjX);
					     Tk2 = FNMS(KP956940335, TjY, TjX);
					     TjT = Tju - Tjn;
					     Tjv = Tjn + Tju;
					     TjS = FNMS(KP831469612, TjE, TjB);
					     TjF = FMA(KP831469612, TjE, TjB);
					     TjP = FNMS(KP956940335, Tjv, Tjg);
					     Tjw = FMA(KP956940335, Tjv, Tjg);
					     TjU = FMA(KP956940335, TjT, TjS);
					     Tk5 = FNMS(KP956940335, TjT, TjS);
					     TjM = FNMS(KP956940335, TjI, TjF);
					     TjJ = FMA(KP956940335, TjI, TjF);
					}
				   }
				   {
					E Tk4, Tk3, Tj5, Tjy, TjL, TjO;
					Tj5 = W[4];
					Tjy = W[5];
					{
					     E Tk1, TjK, Tjx, Tk6;
					     Tk1 = W[36];
					     Tk4 = W[37];
					     TjK = Tj5 * TjJ;
					     Tjx = Tj5 * Tjw;
					     Tk6 = Tk1 * Tk5;
					     Tk3 = Tk1 * Tk2;
					     iio[-WS(ios, 60)] = FMA(Tjy, Tjw, TjK);
					     rio[WS(ios, 3)] = FNMS(Tjy, TjJ, Tjx);
					     iio[-WS(ios, 44)] = FMA(Tk4, Tk2, Tk6);
					}
					rio[WS(ios, 19)] = FNMS(Tk4, Tk5, Tk3);
					TjL = W[68];
					TjO = W[69];
					{
					     E TjR, TjQ, TjN, Tk0;
					     TjR = W[100];
					     TjW = W[101];
					     TjQ = TjL * TjP;
					     TjN = TjL * TjM;
					     Tk0 = TjR * TjZ;
					     TjV = TjR * TjU;
					     rio[WS(ios, 35)] = FNMS(TjO, TjM, TjQ);
					     iio[-WS(ios, 28)] = FMA(TjO, TjP, TjN);
					     rio[WS(ios, 51)] = FNMS(TjW, TjU, Tk0);
					}
				   }
			      }
			      iio[-WS(ios, 12)] = FMA(TjW, TjZ, TjV);
			      {
				   E Tj0, TiN, Tiy, TiS, Tj3, TiK, TiH;
				   {
					E Tit, Tiw, TiV, Tiq, TiW, TiG, TiE, TiF;
					Tit = FNMS(KP820678790, Tis, Tir);
					TiE = FMA(KP820678790, Tir, Tis);
					TiF = FMA(KP820678790, Tiu, Tiv);
					Tiw = FNMS(KP820678790, Tiv, Tiu);
					TiV = FMA(KP980785280, Tip, Tio);
					Tiq = FNMS(KP980785280, Tip, Tio);
					TiW = TiE + TiF;
					TiG = TiE - TiF;
					{
					     E TiR, Tix, TiQ, TiD;
					     TiX = FMA(KP773010453, TiW, TiV);
					     Tj0 = FNMS(KP773010453, TiW, TiV);
					     TiR = Tiw - Tit;
					     Tix = Tit + Tiw;
					     TiQ = FNMS(KP980785280, TiC, TiB);
					     TiD = FMA(KP980785280, TiC, TiB);
					     TiN = FNMS(KP773010453, Tix, Tiq);
					     Tiy = FMA(KP773010453, Tix, Tiq);
					     TiS = FMA(KP773010453, TiR, TiQ);
					     Tj3 = FNMS(KP773010453, TiR, TiQ);
					     TiK = FNMS(KP773010453, TiG, TiD);
					     TiH = FMA(KP773010453, TiG, TiD);
					}
				   }
				   {
					E Tj2, Tj1, Tin, TiA, TiJ, TiM;
					Tin = W[12];
					TiA = W[13];
					{
					     E TiZ, TiI, Tiz, Tj4;
					     TiZ = W[44];
					     Tj2 = W[45];
					     TiI = Tin * TiH;
					     Tiz = Tin * Tiy;
					     Tj4 = TiZ * Tj3;
					     Tj1 = TiZ * Tj0;
					     iio[-WS(ios, 56)] = FMA(TiA, Tiy, TiI);
					     rio[WS(ios, 7)] = FNMS(TiA, TiH, Tiz);
					     iio[-WS(ios, 40)] = FMA(Tj2, Tj0, Tj4);
					}
					rio[WS(ios, 23)] = FNMS(Tj2, Tj3, Tj1);
					TiJ = W[76];
					TiM = W[77];
					{
					     E TiP, TiO, TiL, TiY;
					     TiP = W[108];
					     TiU = W[109];
					     TiO = TiJ * TiN;
					     TiL = TiJ * TiK;
					     TiY = TiP * TiX;
					     TiT = TiP * TiS;
					     rio[WS(ios, 39)] = FNMS(TiM, TiK, TiO);
					     iio[-WS(ios, 24)] = FMA(TiM, TiN, TiL);
					     rio[WS(ios, 55)] = FNMS(TiU, TiS, TiY);
					}
				   }
			      }
			 }
		    }
		    iio[-WS(ios, 8)] = FMA(TiU, TiX, TiT);
		    {
			 E T5z, T5n, T5m, T5A, Tep, Ted, Tec, Teq, Te3, Te0, TdZ;
			 {
			      E T5d, T5a, T55, T5e, T59;
			      {
				   E T5b, T4T, T4z, T5c, T57, T4W, T3a, T56, T4P, T4S;
				   T5z = FNMS(KP707106781, T4O, T4L);
				   T4P = FMA(KP707106781, T4O, T4L);
				   T4S = T4Q + T4R;
				   T5n = T4R - T4Q;
				   {
					E T4U, T4V, T2w, T39, T3R, T4y;
					T4U = FNMS(KP198912367, T3D, T3Q);
					T3R = FMA(KP198912367, T3Q, T3D);
					T4y = FNMS(KP198912367, T4x, T4k);
					T4V = FMA(KP198912367, T4k, T4x);
					T5b = FNMS(KP923879532, T4S, T4P);
					T4T = FMA(KP923879532, T4S, T4P);
					T4z = T3R + T4y;
					T5c = T3R - T4y;
					T2w = FMA(KP707106781, T2v, T2c);
					T5m = FNMS(KP707106781, T2v, T2c);
					T5A = T2P - T38;
					T39 = T2P + T38;
					T57 = T4V - T4U;
					T4W = T4U + T4V;
					T3a = FMA(KP923879532, T39, T2w);
					T56 = FNMS(KP923879532, T39, T2w);
				   }
				   {
					E T53, T58, T50, T5j, T5i, T5h;
					{
					     E T4X, T4C, T5g, T4A, T23;
					     T53 = FNMS(KP980785280, T4W, T4T);
					     T4X = FMA(KP980785280, T4W, T4T);
					     T4C = W[123];
					     T5g = FNMS(KP980785280, T57, T56);
					     T58 = FMA(KP980785280, T57, T56);
					     T4A = FMA(KP980785280, T4z, T3a);
					     T50 = FNMS(KP980785280, T4z, T3a);
					     T23 = W[122];
					     {
						  E T5f, T5k, T4Y, T4B;
						  T5j = FNMS(KP980785280, T5c, T5b);
						  T5d = FMA(KP980785280, T5c, T5b);
						  T5i = W[91];
						  T4Y = T4C * T4A;
						  T4B = T23 * T4A;
						  T5f = W[90];
						  T5k = T5i * T5g;
						  iio[-WS(ios, 1)] = FMA(T23, T4X, T4Y);
						  rio[WS(ios, 62)] = FNMS(T4C, T4X, T4B);
						  T5h = T5f * T5g;
						  iio[-WS(ios, 17)] = FMA(T5f, T5j, T5k);
					     }
					}
					{
					     E T52, T4Z, T54, T51;
					     T52 = W[59];
					     rio[WS(ios, 46)] = FNMS(T5i, T5j, T5h);
					     T4Z = W[58];
					     T54 = T52 * T50;
					     T5a = W[27];
					     T51 = T4Z * T50;
					     T55 = W[26];
					     iio[-WS(ios, 33)] = FMA(T4Z, T53, T54);
					     T5e = T5a * T58;
					     rio[WS(ios, 30)] = FNMS(T52, T53, T51);
					     T59 = T55 * T58;
					}
				   }
			      }
			      {
				   E Te6, TdT, Tds, TdY, Te9, TdQ, TdN;
				   {
					E TcN, Tdq, Tca, Te1, Te2, TdM, TdF, TdI;
					{
					     E TdK, TdL, Tbw, Tc9;
					     Tep = FNMS(KP923879532, Tbv, Tbg);
					     Tbw = FMA(KP923879532, Tbv, Tbg);
					     iio[-WS(ios, 49)] = FMA(T55, T5d, T5e);
					     Ted = Tc8 - TbP;
					     Tc9 = TbP + Tc8;
					     rio[WS(ios, 14)] = FNMS(T5a, T5d, T59);
					     TcN = FNMS(KP098491403, TcM, TcB);
					     TdK = FMA(KP098491403, TcB, TcM);
					     TdL = FNMS(KP098491403, Tde, Tdp);
					     Tdq = FMA(KP098491403, Tdp, Tde);
					     Tca = FMA(KP980785280, Tc9, Tbw);
					     Te1 = FNMS(KP980785280, Tc9, Tbw);
					     Te2 = TdK - TdL;
					     TdM = TdK + TdL;
					     Tec = FNMS(KP923879532, TdE, TdB);
					     TdF = FMA(KP923879532, TdE, TdB);
					     TdI = TdG + TdH;
					     Teq = TdG - TdH;
					}
					{
					     E TdX, Tdr, TdW, TdJ;
					     Te3 = FMA(KP995184726, Te2, Te1);
					     Te6 = FNMS(KP995184726, Te2, Te1);
					     TdX = Tdq - TcN;
					     Tdr = TcN + Tdq;
					     TdW = FNMS(KP980785280, TdI, TdF);
					     TdJ = FMA(KP980785280, TdI, TdF);
					     TdT = FNMS(KP995184726, Tdr, Tca);
					     Tds = FMA(KP995184726, Tdr, Tca);
					     TdY = FMA(KP995184726, TdX, TdW);
					     Te9 = FNMS(KP995184726, TdX, TdW);
					     TdQ = FNMS(KP995184726, TdM, TdJ);
					     TdN = FMA(KP995184726, TdM, TdJ);
					}
				   }
				   {
					E Te8, Te7, Tb5, Tdu, TdP, TdS;
					Tb5 = W[0];
					Tdu = W[1];
					{
					     E Te5, TdO, Tdt, Tea;
					     Te5 = W[32];
					     Te8 = W[33];
					     TdO = Tb5 * TdN;
					     Tdt = Tb5 * Tds;
					     Tea = Te5 * Te9;
					     Te7 = Te5 * Te6;
					     iio[-WS(ios, 62)] = FMA(Tdu, Tds, TdO);
					     rio[WS(ios, 1)] = FNMS(Tdu, TdN, Tdt);
					     iio[-WS(ios, 46)] = FMA(Te8, Te6, Tea);
					}
					rio[WS(ios, 17)] = FNMS(Te8, Te9, Te7);
					TdP = W[64];
					TdS = W[65];
					{
					     E TdV, TdU, TdR, Te4;
					     TdV = W[96];
					     Te0 = W[97];
					     TdU = TdP * TdT;
					     TdR = TdP * TdQ;
					     Te4 = TdV * Te3;
					     TdZ = TdV * TdY;
					     rio[WS(ios, 33)] = FNMS(TdS, TdQ, TdU);
					     iio[-WS(ios, 30)] = FMA(TdS, TdT, TdR);
					     rio[WS(ios, 49)] = FNMS(Te0, TdY, Te4);
					}
				   }
			      }
			 }
			 iio[-WS(ios, 14)] = FMA(Te0, Te3, TdZ);
			 {
			      E T5V, T5S, T5N, T5W, T5R;
			      {
				   E T5T, T5B, T5v, T5U, T5P, T5E, T5o, T5O;
				   {
					E T5C, T5D, T5r, T5u;
					T5C = FMA(KP668178637, T5p, T5q);
					T5r = FNMS(KP668178637, T5q, T5p);
					T5u = FMA(KP668178637, T5t, T5s);
					T5D = FNMS(KP668178637, T5s, T5t);
					T5T = FNMS(KP923879532, T5A, T5z);
					T5B = FMA(KP923879532, T5A, T5z);
					T5v = T5r + T5u;
					T5U = T5u - T5r;
					T5P = T5C - T5D;
					T5E = T5C + T5D;
					T5o = FMA(KP923879532, T5n, T5m);
					T5O = FNMS(KP923879532, T5n, T5m);
				   }
				   {
					E T5L, T5Q, T5I, T61, T60, T5Z;
					{
					     E T5F, T5y, T5Y, T5w, T5l;
					     T5L = FNMS(KP831469612, T5E, T5B);
					     T5F = FMA(KP831469612, T5E, T5B);
					     T5y = W[11];
					     T5Y = FNMS(KP831469612, T5P, T5O);
					     T5Q = FMA(KP831469612, T5P, T5O);
					     T5w = FMA(KP831469612, T5v, T5o);
					     T5I = FNMS(KP831469612, T5v, T5o);
					     T5l = W[10];
					     {
						  E T5X, T62, T5G, T5x;
						  T61 = FNMS(KP831469612, T5U, T5T);
						  T5V = FMA(KP831469612, T5U, T5T);
						  T60 = W[43];
						  T5G = T5y * T5w;
						  T5x = T5l * T5w;
						  T5X = W[42];
						  T62 = T60 * T5Y;
						  iio[-WS(ios, 57)] = FMA(T5l, T5F, T5G);
						  rio[WS(ios, 6)] = FNMS(T5y, T5F, T5x);
						  T5Z = T5X * T5Y;
						  iio[-WS(ios, 41)] = FMA(T5X, T61, T62);
					     }
					}
					{
					     E T5K, T5H, T5M, T5J;
					     T5K = W[75];
					     rio[WS(ios, 22)] = FNMS(T60, T61, T5Z);
					     T5H = W[74];
					     T5M = T5K * T5I;
					     T5S = W[107];
					     T5J = T5H * T5I;
					     T5N = W[106];
					     iio[-WS(ios, 25)] = FMA(T5H, T5L, T5M);
					     T5W = T5S * T5Q;
					     rio[WS(ios, 38)] = FNMS(T5K, T5L, T5J);
					     T5R = T5N * T5Q;
					}
				   }
			      }
			      {
				   E TeO, TeB, Tem, TeG, TeR, Tey, Tev;
				   {
					E Teh, Tek, Tee, TeJ, TeK, Teu, Tes, Tet;
					iio[-WS(ios, 9)] = FMA(T5N, T5V, T5W);
					rio[WS(ios, 54)] = FNMS(T5S, T5V, T5R);
					Teh = FNMS(KP820678790, Teg, Tef);
					Tes = FMA(KP820678790, Tef, Teg);
					Tet = FNMS(KP820678790, Tei, Tej);
					Tek = FMA(KP820678790, Tej, Tei);
					Tee = FMA(KP980785280, Ted, Tec);
					TeJ = FNMS(KP980785280, Ted, Tec);
					TeK = Tes - Tet;
					Teu = Tes + Tet;
					{
					     E TeF, Tel, TeE, Ter;
					     TeL = FMA(KP773010453, TeK, TeJ);
					     TeO = FNMS(KP773010453, TeK, TeJ);
					     TeF = Tek - Teh;
					     Tel = Teh + Tek;
					     TeE = FNMS(KP980785280, Teq, Tep);
					     Ter = FMA(KP980785280, Teq, Tep);
					     TeB = FNMS(KP773010453, Tel, Tee);
					     Tem = FMA(KP773010453, Tel, Tee);
					     TeG = FMA(KP773010453, TeF, TeE);
					     TeR = FNMS(KP773010453, TeF, TeE);
					     Tey = FNMS(KP773010453, Teu, Ter);
					     Tev = FMA(KP773010453, Teu, Ter);
					}
				   }
				   {
					E TeQ, TeP, Teb, Teo, Tex, TeA;
					Teb = W[112];
					Teo = W[113];
					{
					     E TeN, Tew, Ten, TeS;
					     TeN = W[80];
					     TeQ = W[81];
					     Tew = Teb * Tev;
					     Ten = Teb * Tem;
					     TeS = TeN * TeR;
					     TeP = TeN * TeO;
					     rio[WS(ios, 57)] = FNMS(Teo, Tem, Tew);
					     iio[-WS(ios, 6)] = FMA(Teo, Tev, Ten);
					     rio[WS(ios, 41)] = FNMS(TeQ, TeO, TeS);
					}
					iio[-WS(ios, 22)] = FMA(TeQ, TeR, TeP);
					Tex = W[48];
					TeA = W[49];
					{
					     E TeD, TeC, Tez, TeM;
					     TeD = W[16];
					     TeI = W[17];
					     TeC = Tex * TeB;
					     Tez = Tex * Tey;
					     TeM = TeD * TeL;
					     TeH = TeD * TeG;
					     iio[-WS(ios, 38)] = FMA(TeA, Tey, TeC);
					     rio[WS(ios, 25)] = FNMS(TeA, TeB, Tez);
					     iio[-WS(ios, 54)] = FMA(TeI, TeG, TeM);
					}
				   }
			      }
			 }
		    }
		    rio[WS(ios, 9)] = FNMS(TeI, TeL, TeH);
		    {
			 E T7j, T77, T76, T7k, Tg9, TfX, TfW, Tga, TfN, TfK, TfJ;
			 {
			      E T6X, T6U, T6P, T6Y, T6T;
			      {
				   E T6V, T6D, T6t, T6W, T6R, T6G, T6e, T6Q, T6z, T6C;
				   T7j = FNMS(KP707106781, T6y, T6x);
				   T6z = FMA(KP707106781, T6y, T6x);
				   T6C = T6A + T6B;
				   T77 = T6A - T6B;
				   {
					E T6E, T6F, T66, T6d, T6l, T6s;
					T6E = FMA(KP198912367, T6h, T6k);
					T6l = FNMS(KP198912367, T6k, T6h);
					T6s = FMA(KP198912367, T6r, T6o);
					T6F = FNMS(KP198912367, T6o, T6r);
					T6V = FNMS(KP923879532, T6C, T6z);
					T6D = FMA(KP923879532, T6C, T6z);