r2cf_128.c

快速fft变换

					     {
						  E T8i, T8h, T5u, T5p;
						  T8i = FNMS(KP707106781, T5t, T5s);
						  T5u = FMA(KP707106781, T5t, T5s);
						  T5p = FMA(KP707106781, T5o, T5h);
						  T8h = FNMS(KP707106781, T5o, T5h);
						  TeM = TdM - TdL;
						  TdN = TdL + TdM;
						  T5I = FMA(KP198912367, T5p, T5u);
						  T5v = FNMS(KP198912367, T5u, T5p);
						  T8q = FNMS(KP668178637, T8h, T8i);
						  T8j = FMA(KP668178637, T8i, T8h);
						  TeP = TdF - TdC;
						  TdG = TdC + TdF;
					     }
					}
				   }
			      }
			 }
		    }
	       }
	       {
		    E T8p, T8g, TcH, TdW, TdT, TcO, Tfp, Tfk, Tfj, Tfq;
		    {
			 E Tbj, Tba, Tcy, Tco, TcB, Tcl, Tcx, Tcv, Tcz, Tcr;
			 {
			      E Tch, Tct, Tcp, Tcq, Tci, T1w, TbV, T11, Tcf, Tc9, T21, Tcj, Tcm, TbY, Tc1;
			      E Tcn, Tcu, Tck;
			      {
				   E Tv, T10, Tc5, Tc8;
				   Tch = Tf - Tu;
				   Tv = Tf + Tu;
				   T8p = FMA(KP668178637, T8e, T8f);
				   T8g = FNMS(KP668178637, T8f, T8e);
				   Tbj = Tb9 - Tb2;
				   Tba = Tb2 + Tb9;
				   T10 = TK + TZ;
				   Tct = TZ - TK;
				   Tcp = Tc3 - Tc4;
				   Tc5 = Tc3 + Tc4;
				   Tc8 = Tc6 + Tc7;
				   Tcq = Tc7 - Tc6;
				   Tci = T1g - T1v;
				   T1w = T1g + T1v;
				   TbV = Tv - T10;
				   T11 = Tv + T10;
				   Tcf = Tc5 + Tc8;
				   Tc9 = Tc5 - Tc8;
				   T21 = T1L + T20;
				   Tcj = T1L - T20;
				   Tcm = TbW - TbX;
				   TbY = TbW + TbX;
				   Tc1 = TbZ + Tc0;
				   Tcn = Tc0 - TbZ;
			      }
			      {
				   E Tcb, T22, Tce, Tc2;
				   Tcb = T21 - T1w;
				   T22 = T1w + T21;
				   Tce = TbY + Tc1;
				   Tc2 = TbY - Tc1;
				   {
					E Tcd, Tcg, Tca, Tcc;
					Tcd = T11 + T22;
					Cr[WS(csr, 32)] = T11 - T22;
					Tcg = Tce + Tcf;
					Ci[WS(csi, 32)] = Tcf - Tce;
					Tca = Tc2 + Tc9;
					Tcc = Tc9 - Tc2;
					Cr[0] = Tcd + Tcg;
					Cr[WS(csr, 64)] = Tcd - Tcg;
					Ci[WS(csi, 48)] = FMS(KP707106781, Tcc, Tcb);
					Ci[WS(csi, 16)] = FMA(KP707106781, Tcc, Tcb);
					Cr[WS(csr, 16)] = FMA(KP707106781, Tca, TbV);
					Cr[WS(csr, 48)] = FNMS(KP707106781, Tca, TbV);
					Tcu = Tcj - Tci;
					Tck = Tci + Tcj;
					Tcy = FNMS(KP414213562, Tcm, Tcn);
					Tco = FMA(KP414213562, Tcn, Tcm);
				   }
			      }
			      TcB = FNMS(KP707106781, Tck, Tch);
			      Tcl = FMA(KP707106781, Tck, Tch);
			      Tcx = FMA(KP707106781, Tcu, Tct);
			      Tcv = FNMS(KP707106781, Tcu, Tct);
			      Tcz = FMA(KP414213562, Tcp, Tcq);
			      Tcr = FNMS(KP414213562, Tcq, Tcp);
			 }
			 {
			      E TbT, TbO, TbN, TbU;
			      {
				   E Ta9, TbB, Tbb, TbL, Tbp, TbM, Tag, Tbk, TbR, TbJ, Tbw, TaO, TbC, Tbs, TbQ;
				   E TbG;
				   {
					E Tbq, Tbr, TbH, TbI;
					{
					     E Tbo, Ta8, Tac, Taf;
					     Tbo = Ta7 - Ta6;
					     Ta8 = Ta6 + Ta7;
					     {
						  E TcC, TcA, Tcw, Tcs;
						  TcC = Tcz - Tcy;
						  TcA = Tcy + Tcz;
						  Tcw = Tcr - Tco;
						  Tcs = Tco + Tcr;
						  Cr[WS(csr, 24)] = FMA(KP923879532, TcC, TcB);
						  Cr[WS(csr, 40)] = FNMS(KP923879532, TcC, TcB);
						  Ci[WS(csi, 56)] = FMS(KP923879532, TcA, Tcx);
						  Ci[WS(csi, 8)] = FMA(KP923879532, TcA, Tcx);
						  Ci[WS(csi, 40)] = FMA(KP923879532, Tcw, Tcv);
						  Ci[WS(csi, 24)] = FMS(KP923879532, Tcw, Tcv);
						  Cr[WS(csr, 8)] = FMA(KP923879532, Tcs, Tcl);
						  Cr[WS(csr, 56)] = FNMS(KP923879532, Tcs, Tcl);
						  Ta9 = FMA(KP707106781, Ta8, Ta5);
						  TbB = FNMS(KP707106781, Ta8, Ta5);
					     }
					     Tbq = FNMS(KP414213562, Taa, Tab);
					     Tac = FMA(KP414213562, Tab, Taa);
					     Taf = FNMS(KP414213562, Tae, Tad);
					     Tbr = FMA(KP414213562, Tad, Tae);
					     Tbb = FMA(KP707106781, Tba, TaV);
					     TbH = FNMS(KP707106781, Tba, TaV);
					     TbL = FNMS(KP707106781, Tbo, Tbn);
					     Tbp = FMA(KP707106781, Tbo, Tbn);
					     TbM = Taf - Tac;
					     Tag = Tac + Taf;
					     TbI = FNMS(KP707106781, Tbj, Tbi);
					     Tbk = FMA(KP707106781, Tbj, Tbi);
					}
					{
					     E TbE, TbF, TaE, TaN;
					     TbE = FNMS(KP707106781, TaD, Tao);
					     TaE = FMA(KP707106781, TaD, Tao);
					     TaN = FMA(KP707106781, TaM, TaL);
					     TbF = FNMS(KP707106781, TaM, TaL);
					     TbR = FNMS(KP668178637, TbH, TbI);
					     TbJ = FMA(KP668178637, TbI, TbH);
					     Tbw = FNMS(KP198912367, TaE, TaN);
					     TaO = FMA(KP198912367, TaN, TaE);
					     TbC = Tbr - Tbq;
					     Tbs = Tbq + Tbr;
					     TbQ = FMA(KP668178637, TbE, TbF);
					     TbG = FNMS(KP668178637, TbF, TbE);
					}
				   }
				   {
					E Tbz, Tah, Tbv, Tbt, Tbx, Tbl;
					Tbz = FNMS(KP923879532, Tag, Ta9);
					Tah = FMA(KP923879532, Tag, Ta9);
					Tbv = FMA(KP923879532, Tbs, Tbp);
					Tbt = FNMS(KP923879532, Tbs, Tbp);
					Tbx = FMA(KP198912367, Tbb, Tbk);
					Tbl = FNMS(KP198912367, Tbk, Tbb);
					{
					     E TbD, TbK, TbP, TbS;
					     TbT = FNMS(KP923879532, TbC, TbB);
					     TbD = FMA(KP923879532, TbC, TbB);
					     {
						  E TbA, Tby, Tbu, Tbm;
						  TbA = Tbx - Tbw;
						  Tby = Tbw + Tbx;
						  Tbu = Tbl - TaO;
						  Tbm = TaO + Tbl;
						  Cr[WS(csr, 28)] = FMA(KP980785280, TbA, Tbz);
						  Cr[WS(csr, 36)] = FNMS(KP980785280, TbA, Tbz);
						  Ci[WS(csi, 60)] = FMS(KP980785280, Tby, Tbv);
						  Ci[WS(csi, 4)] = FMA(KP980785280, Tby, Tbv);
						  Ci[WS(csi, 36)] = FMA(KP980785280, Tbu, Tbt);
						  Ci[WS(csi, 28)] = FMS(KP980785280, Tbu, Tbt);
						  Cr[WS(csr, 4)] = FMA(KP980785280, Tbm, Tah);
						  Cr[WS(csr, 60)] = FNMS(KP980785280, Tbm, Tah);
						  TbK = TbG + TbJ;
						  TbO = TbJ - TbG;
					     }
					     TbN = FMA(KP923879532, TbM, TbL);
					     TbP = FNMS(KP923879532, TbM, TbL);
					     TbS = TbQ + TbR;
					     TbU = TbQ - TbR;
					     Cr[WS(csr, 12)] = FMA(KP831469612, TbK, TbD);
					     Cr[WS(csr, 52)] = FNMS(KP831469612, TbK, TbD);
					     Ci[WS(csi, 52)] = FNMS(KP831469612, TbS, TbP);
					     Ci[WS(csi, 12)] = -(FMA(KP831469612, TbS, TbP));
					}
				   }
			      }
			      {
				   E TeN, Tf7, Tev, Tfm, Tfc, TeQ, TeX, TeW, Tfn, Tff, Tfi, TeC, Tf2, TeK, Tfh;
				   E TeV, Tf8;
				   {
					E TeG, TeJ, Tfd, Tfe, Tey, TeB, TeT, TeU;
					{
					     E Tet, Teu, Tfa, Tfb;
					     TcH = FMA(KP707106781, TcG, TcD);
					     Tet = FNMS(KP707106781, TcG, TcD);
					     Ci[WS(csi, 44)] = FMS(KP831469612, TbO, TbN);
					     Ci[WS(csi, 20)] = FMA(KP831469612, TbO, TbN);
					     Cr[WS(csr, 20)] = FMA(KP831469612, TbU, TbT);
					     Cr[WS(csr, 44)] = FNMS(KP831469612, TbU, TbT);
					     Teu = TdV - TdU;
					     TdW = TdU + TdV;
					     TeG = FNMS(KP923879532, TeF, TeE);
					     Tfa = FMA(KP923879532, TeF, TeE);
					     Tfb = FMA(KP923879532, TeI, TeH);
					     TeJ = FNMS(KP923879532, TeI, TeH);
					     TeN = FNMS(KP923879532, TeM, TeL);
					     Tfd = FMA(KP923879532, TeM, TeL);
					     Tf7 = FMA(KP923879532, Teu, Tet);
					     Tev = FNMS(KP923879532, Teu, Tet);
					     Tfm = FMA(KP303346683, Tfa, Tfb);
					     Tfc = FNMS(KP303346683, Tfb, Tfa);
					     Tfe = FNMS(KP923879532, TeP, TeO);
					     TeQ = FMA(KP923879532, TeP, TeO);
					     TeX = FNMS(KP668178637, Tew, Tex);
					     Tey = FMA(KP668178637, Tex, Tew);
					     TeB = FNMS(KP668178637, TeA, Tez);
					     TeW = FMA(KP668178637, Tez, TeA);
					}
					Tfn = FNMS(KP303346683, Tfd, Tfe);
					Tff = FMA(KP303346683, Tfe, Tfd);
					Tfi = Tey + TeB;
					TeC = Tey - TeB;
					TdT = FMA(KP707106781, TdS, TdR);
					TeT = FNMS(KP707106781, TdS, TdR);
					TeU = TcN - TcK;
					TcO = TcK + TcN;
					Tf2 = FNMS(KP534511135, TeG, TeJ);
					TeK = FMA(KP534511135, TeJ, TeG);
					Tfh = FNMS(KP923879532, TeU, TeT);
					TeV = FMA(KP923879532, TeU, TeT);
				   }
				   {
					E Tf5, TeD, TeY, Tf3, TeR;
					Tf5 = FNMS(KP831469612, TeC, Tev);
					TeD = FMA(KP831469612, TeC, Tev);
					Tf8 = TeX + TeW;
					TeY = TeW - TeX;
					Tf3 = FMA(KP534511135, TeN, TeQ);
					TeR = FNMS(KP534511135, TeQ, TeN);
					{
					     E Tf1, TeZ, Tf6, Tf4, Tf0, TeS;
					     Tf1 = FMA(KP831469612, TeY, TeV);
					     TeZ = FNMS(KP831469612, TeY, TeV);
					     Tf6 = Tf3 - Tf2;
					     Tf4 = Tf2 + Tf3;
					     Tf0 = TeR - TeK;
					     TeS = TeK + TeR;
					     Ci[WS(csi, 54)] = FMS(KP881921264, Tf4, Tf1);
					     Ci[WS(csi, 10)] = FMA(KP881921264, Tf4, Tf1);
					     Ci[WS(csi, 42)] = FMA(KP881921264, Tf0, TeZ);
					     Ci[WS(csi, 22)] = FMS(KP881921264, Tf0, TeZ);
					     Cr[WS(csr, 10)] = FMA(KP881921264, TeS, TeD);
					     Cr[WS(csr, 54)] = FNMS(KP881921264, TeS, TeD);
					     Cr[WS(csr, 42)] = FNMS(KP881921264, Tf6, Tf5);
					     Cr[WS(csr, 22)] = FMA(KP881921264, Tf6, Tf5);
					}
				   }
				   {
					E Tf9, Tfg, Tfl, Tfo;
					Tfp = FNMS(KP831469612, Tf8, Tf7);
					Tf9 = FMA(KP831469612, Tf8, Tf7);
					Tfg = Tfc + Tff;
					Tfk = Tff - Tfc;
					Tfj = FNMS(KP831469612, Tfi, Tfh);
					Tfl = FMA(KP831469612, Tfi, Tfh);
					Tfo = Tfm + Tfn;
					Tfq = Tfm - Tfn;
					Cr[WS(csr, 6)] = FMA(KP956940335, Tfg, Tf9);
					Cr[WS(csr, 58)] = FNMS(KP956940335, Tfg, Tf9);
					Ci[WS(csi, 58)] = FNMS(KP956940335, Tfo, Tfl);
					Ci[WS(csi, 6)] = -(FMA(KP956940335, Tfo, Tfl));
				   }
			      }
			 }
		    }
		    {
			 E T2f, T5W, T5T, T2y, T5J, T5w, T4u, T4h, T7p, T7q;
			 {
			      E Ter, Tem, Tel, Tes;
			      {
				   E TdH, Te9, TcP, Teo, Tee, TdO, TdY, TdZ, Tep, Teh, Tek, Td8, Te4, Tdu, Tej;
				   E TdX, Tea;
				   {
					E Tdm, Tdt, Tef, Teg, TcY, Td7, Tec, Ted;
					Ci[WS(csi, 38)] = FMS(KP956940335, Tfk, Tfj);
					Ci[WS(csi, 26)] = FMA(KP956940335, Tfk, Tfj);
					Cr[WS(csr, 26)] = FMA(KP956940335, Tfq, Tfp);
					Cr[WS(csr, 38)] = FNMS(KP956940335, Tfq, Tfp);
					Tdm = FMA(KP923879532, Tdl, Tde);
					Tec = FNMS(KP923879532, Tdl, Tde);
					Ted = FNMS(KP923879532, Tds, Tdp);
					Tdt = FMA(KP923879532, Tds, Tdp);
					TdH = FMA(KP923879532, TdG, Tdz);
					Tef = FNMS(KP923879532, TdG, Tdz);
					Te9 = FNMS(KP923879532, TcO, TcH);
					TcP = FMA(KP923879532, TcO, TcH);
					Teo = FMA(KP820678790, Tec, Ted);
					Tee = FNMS(KP820678790, Ted, Tec);
					Teg = FNMS(KP923879532, TdN, TdK);
					TdO = FMA(KP923879532, TdN, TdK);
					TdY = FNMS(KP198912367, TcU, TcX);
					TcY = FMA(KP198912367, TcX, TcU);
					Td7 = FNMS(KP198912367, Td6, Td3);
					TdZ = FMA(KP198912367, Td3, Td6);
					Tep = FNMS(KP820678790, Tef, Teg);
					Teh = FMA(KP820678790, Teg, Tef);
					Tek = Td7 - TcY;
					Td8 = TcY + Td7;
					Te4 = FNMS(KP098491403, Tdm, Tdt);
					Tdu = FMA(KP098491403, Tdt, Tdm);
					Tej = FNMS(KP923879532, TdW, TdT);
					TdX = FMA(KP923879532, TdW, TdT);
				   }
				   {
					E Te7, Td9, Te0, Te5, TdP;
					Te7 = FNMS(KP980785280, Td8, TcP);
					Td9 = FMA(KP980785280, Td8, TcP);
					Tea = TdZ - TdY;
					Te0 = TdY + TdZ;
					Te5 = FMA(KP098491403, TdH, TdO);
					TdP = FNMS(KP098491403, TdO, TdH);
					{
					     E Te3, Te1, Te8, Te6, Te2, TdQ;
					     Te3 = FMA(KP980785280, Te0, TdX);
					     Te1 = FNMS(KP980785280, Te0, TdX);
					     Te8 = Te5 - Te4;
					     Te6 = Te4 + Te5;
					     Te2 = TdP - Tdu;
					     TdQ = Tdu + TdP;
					     Ci[WS(csi, 62)] = FMS(KP995184726, Te6, Te3);
					     Ci[WS(csi, 2)] = FMA(KP995184726, Te6, Te3);
					     Ci[WS(csi, 34)] = FMA(KP995184726, Te2, Te1);
					     Ci[WS(csi, 30)] = FMS(KP995184726, Te2, Te1);
					     Cr[WS(csr, 2)] = FMA(KP995184726, TdQ, Td9);
					     Cr[WS(csr, 62)] = FNMS(KP995184726, TdQ, Td9);
					     Cr[WS(csr, 34)] = FNMS(KP995184726, Te8, Te7);
					     Cr[WS(csr, 30)] = FMA(KP995184726, Te8, Te7);
					}
				   }
				   {
					E Teb, Tei, Ten, Teq;
					Ter = FNMS(KP980785280, Tea, Te9);
					Teb = FMA(KP980785280, Tea, Te9);
					Tei = Tee + Teh;
					Tem = Teh - Tee;
					Tel = FMA(KP980785280, Tek, Tej);
					Ten = FNMS(KP980785280, Tek, Tej);
					Teq = Teo + Tep;
					Tes = Teo - Tep;
					Cr[WS(csr, 14)] = FMA(KP773010453, Tei, Teb);
					Cr[WS(csr, 50)] = FNMS(KP773010453, Tei, Teb);
					Ci[WS(csi, 50)] = FNMS(KP773010453, Teq, Ten);
					Ci[WS(csi, 14)] = -(FMA(KP773010453, Teq, Ten));
				   }
			      }
			      {
				   E T77, T6v, T7i, T6C, T78, T6Y, T7h, T6V, T6N, T7d, T6P, T6F, T6I;
				   {
					E T6W, T6X, T6T, T6U, T6M;
					{
					     E T6t, T6u, T6y, T6B;
					     T2f = FMA(KP923879532, T2e, T27);
					     T6t = FNMS(KP923879532, T2e, T27);
					     Ci[WS(csi, 46)] = FMS(KP773010453, Tem, Tel);
					     Ci[WS(csi, 18)] = FMA(KP773010453, Tem, Tel);
					     Cr[WS(csr, 18)] = FMA(KP773010453, Tes, Ter);
					     Cr[WS(csr, 46)] = FNMS(KP773010453, Tes, Ter);
					     T6u = T5U - T5V;
					     T5W = T5U + T5V;
					     T6W = FNMS(KP820678790, T6w, T6x);
					     T6y = FMA(KP820678790, T6x, T6w);
					     T6B = FNMS(KP820678790, T6A, T6z);
					     T6X = FMA(KP820678790, T6z, T6A);
					     T77 = FMA(KP980785280, T6u, T6t);
					     T6v = FNMS(KP980785280, T6u, T6t);
					     T7i = T6B + T6y;
					     T6C = T6y - T6B;
					}
					T5T = FMA(KP923879532, T5S, T5P);
					T6T = FNMS(KP923879532, T5S, T5P);
					T6U = T2x - T2o;
					T2y = T2o + T2x;
					T5J = T5H + T5I;
					T6M = T5I - T5H;
					T78 = T6X + T6W;
					T6Y = T6W - T6X;
					T7h = FMA(KP980785280, T6U, T6T);
					T6V = FNMS(KP980785280, T6U, T6T);
					T6N = FNMS(KP980785280, T6M, T6L);
					T7d = FMA(KP980785280, T6M, T6L);
					T6P = T5v - T5e;
					T5w = T5e + T5v;
					T4u = T4s + T4t;
					T6F = T4s - T4t;
					T6I = T4g - T3Z;
					T4h = T3Z + T4g;
				   }
				   {
					E T75, T7f, T7n, T7c, T7m, T76;
					{
					     E T6D, T72, T6R, T73, T6K, T71, T6Z, T7e, T6Q, T74, T70, T6S;
					     T75 = FNMS(KP773010453, T6C, T6v);
					     T6D = FMA(KP773010453, T6C, T6v);
					     T7e = FNMS(KP980785280, T6P, T6O);
					     T6Q = FMA(KP980785280, T6P, T6O);
					     {
						  E T7a, T6G, T7b, T6J;
						  T7a = FMA(KP980785280, T6F, T6E);
						  T6G = FNMS(KP980785280, T6F, T6E);
						  T7b = FMA(KP980785280, T6I, T6H);
						  T6J = FNMS(KP980785280, T6I, T6H);
						  T7f = FMA(KP357805721, T7e, T7d);
						  T7n = FNMS(KP357805721, T7d, T7e);
						  T72 = FMA(KP472964775, T6N, T6Q);
						  T6R = FNMS(KP472964775, T6Q, T6N);
						  T7c = FMA(KP357805721, T7b, T7a);
						  T7m = FNMS(KP357805721, T7a, T7b);
						  T73 = FMA(KP472964775, T6G, T6J);