r2cb_128.c

快速fft变换

						  T3n = T3l + T3m;
						  T7C = T3p + T3s;
						  T3t = T3p - T3s;
						  TaC = TW - TZ;
						  T10 = TW + TZ;
						  TdB = TaE + TaD;
						  TaF = TaD - TaE;
						  T7D = T3k + T3n;
						  T3o = T3k - T3n;
					     }
					     TaB = Tax - TaA;
					     TaO = Tax + TaA;
					     TaP = TaF - TaC;
					     TaG = TaC + TaF;
					     TdF = TT - T10;
					     T11 = TT + T10;
					     Tct = TaO - TaP;
					     TaQ = TaO + TaP;
					     T3G = FMA(KP414213562, T3d, T3i);
					     T3j = FNMS(KP414213562, T3i, T3d);
					     T3u = FMA(KP414213562, T3t, T3o);
					     T3H = FNMS(KP414213562, T3o, T3t);
					}
					TeQ = TdC + TdB;
					TdD = TdB - TdC;
					T6m = T3u - T3j;
					T3v = T3j + T3u;
					T7B = FNMS(KP414213562, T7A, T7z);
					T7K = FMA(KP414213562, T7z, T7A);
					T7L = FMA(KP414213562, T7C, T7D);
					T7E = FNMS(KP414213562, T7D, T7C);
					T6j = T3G - T3H;
					T3I = T3G + T3H;
					T99 = T7K + T7L;
					T7M = T7K - T7L;
				   }
			      }
			 }
		    }
	       }
	       {
		    E Tcw, T9c, T7F, Tev, Teu, TeD, Tep, TeG, Tez, TeE, Tes;
		    {
			 E TbX, TbY, Tc7, TbP, Tar, Tc5, Tc1, Tc0, Tc4, Tba, TbS, TbL, TbQ, TaS, Tbx;
			 E Tc8;
			 {
			      E TeO, TaH, TeR, TeL, TeU, TeZ, Tf2, TeX, Tfh, Tfn, Tfo, Tfm;
			      {
				   E T12, Tfg, Tfj, Tx, Tff, T24, Tfi, Tfk, Th, T1y, T23;
				   TeO = TM - T11;
				   T12 = TM + T11;
				   Tcw = TaG - TaB;
				   TaH = TaB + TaG;
				   T9c = T7B - T7E;
				   T7F = T7B + T7E;
				   Tfg = TeQ + TeP;
				   TeR = TeP - TeQ;
				   TeL = FNMS(KP2_000000000, Tg, T9);
				   Th = FMA(KP2_000000000, Tg, T9);
				   T1y = T1i + T1x;
				   TeU = T1i - T1x;
				   TeZ = T1N - T22;
				   T23 = T1N + T22;
				   Tfj = Tf1 + Tf0;
				   Tf2 = Tf0 - Tf1;
				   Tx = FMA(KP2_000000000, Tw, Th);
				   Tff = FNMS(KP2_000000000, Tw, Th);
				   T24 = T1y + T23;
				   Tfi = T1y - T23;
				   TeX = TeV - TeW;
				   Tfk = TeW + TeV;
				   {
					E T13, Tfp, Tfl, Tfq;
					T13 = FMA(KP2_000000000, T12, Tx);
					Tfp = FNMS(KP2_000000000, T12, Tx);
					Tfh = FNMS(KP2_000000000, Tfg, Tff);
					Tfn = FMA(KP2_000000000, Tfg, Tff);
					Tfl = Tfj - Tfk;
					Tfq = Tfk + Tfj;
					R0[0] = FMA(KP2_000000000, T24, T13);
					R0[WS(rs, 32)] = FNMS(KP2_000000000, T24, T13);
					R0[WS(rs, 48)] = FMA(KP2_000000000, Tfq, Tfp);
					R0[WS(rs, 16)] = FNMS(KP2_000000000, Tfq, Tfp);
					Tfo = Tfi + Tfl;
					Tfm = Tfi - Tfl;
				   }
			      }
			      {
				   E Tf7, TeN, Tfa, Tf3, Tf8, TeS;
				   R0[WS(rs, 8)] = FMA(KP1_414213562, Tfm, Tfh);
				   R0[WS(rs, 40)] = FNMS(KP1_414213562, Tfm, Tfh);
				   R0[WS(rs, 56)] = FMA(KP1_414213562, Tfo, Tfn);
				   R0[WS(rs, 24)] = FNMS(KP1_414213562, Tfo, Tfn);
				   Tf7 = FMA(KP2_000000000, TeM, TeL);
				   TeN = FNMS(KP2_000000000, TeM, TeL);
				   Tfa = Tf2 - TeZ;
				   Tf3 = TeZ + Tf2;
				   Tf8 = TeO + TeR;
				   TeS = TeO - TeR;
				   {
					E TbH, TbK, TaI, TaR;
					{
					     E Taf, Tf9, Tfd, Tf5, TeT, Tfb, TeY, Taq;
					     TbX = FNMS(KP1_414213562, Tae, Ta9);
					     Taf = FMA(KP1_414213562, Tae, Ta9);
					     Tf9 = FNMS(KP1_414213562, Tf8, Tf7);
					     Tfd = FMA(KP1_414213562, Tf8, Tf7);
					     Tf5 = FNMS(KP1_414213562, TeS, TeN);
					     TeT = FMA(KP1_414213562, TeS, TeN);
					     Tfb = TeU - TeX;
					     TeY = TeU + TeX;
					     Taq = FNMS(KP414213562, Tap, Tak);
					     TbY = FMA(KP414213562, Tak, Tap);
					     Tc7 = FNMS(KP707106781, TbG, TbD);
					     TbH = FMA(KP707106781, TbG, TbD);
					     {
						  E Tfc, Tfe, Tf6, Tf4;
						  Tfc = FNMS(KP414213562, Tfb, Tfa);
						  Tfe = FMA(KP414213562, Tfa, Tfb);
						  Tf6 = FMA(KP414213562, TeY, Tf3);
						  Tf4 = FNMS(KP414213562, Tf3, TeY);
						  TbP = FNMS(KP1_847759065, Taq, Taf);
						  Tar = FMA(KP1_847759065, Taq, Taf);
						  R0[WS(rs, 44)] = FMA(KP1_847759065, Tfc, Tf9);
						  R0[WS(rs, 12)] = FNMS(KP1_847759065, Tfc, Tf9);
						  R0[WS(rs, 60)] = FMA(KP1_847759065, Tfe, Tfd);
						  R0[WS(rs, 28)] = FNMS(KP1_847759065, Tfe, Tfd);
						  R0[WS(rs, 52)] = FMA(KP1_847759065, Tf6, Tf5);
						  R0[WS(rs, 20)] = FNMS(KP1_847759065, Tf6, Tf5);
						  R0[WS(rs, 4)] = FMA(KP1_847759065, Tf4, TeT);
						  R0[WS(rs, 36)] = FNMS(KP1_847759065, Tf4, TeT);
						  TbK = TbI + TbJ;
						  Tc5 = TbI - TbJ;
					     }
					}
					Tc1 = FNMS(KP707106781, TaH, Taw);
					TaI = FMA(KP707106781, TaH, Taw);
					TaR = FMA(KP707106781, TaQ, TaN);
					Tc0 = FNMS(KP707106781, TaQ, TaN);
					Tc4 = FNMS(KP707106781, Tb9, TaY);
					Tba = FMA(KP707106781, Tb9, TaY);
					TbS = FNMS(KP923879532, TbK, TbH);
					TbL = FMA(KP923879532, TbK, TbH);
					TbQ = FMA(KP198912367, TaI, TaR);
					TaS = FNMS(KP198912367, TaR, TaI);
					Tbx = Tbl + Tbw;
					Tc8 = Tbw - Tbl;
				   }
			      }
			 }
			 {
			      E Ten, Teo, Tex, Tef, Tdz, Ter, Teq, TdQ, Tei, Teb, Teg, TdK, Te1, Tey;
			      {
				   E Te7, Tea, TdE, TdJ;
				   {
					E Tdt, TbR, TbV, TbN, TaT, TbT, Tby, Tdy;
					Ten = FMA(KP2_000000000, Tds, Tdr);
					Tdt = FNMS(KP2_000000000, Tds, Tdr);
					TbR = FNMS(KP1_961570560, TbQ, TbP);
					TbV = FMA(KP1_961570560, TbQ, TbP);
					TbN = FNMS(KP1_961570560, TaS, Tar);
					TaT = FMA(KP1_961570560, TaS, Tar);
					TbT = FNMS(KP923879532, Tbx, Tba);
					Tby = FMA(KP923879532, Tbx, Tba);
					Tdy = Tdu - Tdx;
					Teo = Tdu + Tdx;
					Tex = Te6 - Te3;
					Te7 = Te3 + Te6;
					{
					     E TbU, TbW, TbO, TbM;
					     TbU = FNMS(KP820678790, TbT, TbS);
					     TbW = FMA(KP820678790, TbS, TbT);
					     TbO = FMA(KP098491403, Tby, TbL);
					     TbM = FNMS(KP098491403, TbL, Tby);
					     Tef = FNMS(KP1_414213562, Tdy, Tdt);
					     Tdz = FMA(KP1_414213562, Tdy, Tdt);
					     R0[WS(rs, 41)] = FMA(KP1_546020906, TbU, TbR);
					     R0[WS(rs, 9)] = FNMS(KP1_546020906, TbU, TbR);
					     R0[WS(rs, 57)] = FMA(KP1_546020906, TbW, TbV);
					     R0[WS(rs, 25)] = FNMS(KP1_546020906, TbW, TbV);
					     R0[WS(rs, 49)] = FMA(KP1_990369453, TbO, TbN);
					     R0[WS(rs, 17)] = FNMS(KP1_990369453, TbO, TbN);
					     R0[WS(rs, 1)] = FMA(KP1_990369453, TbM, TaT);
					     R0[WS(rs, 33)] = FNMS(KP1_990369453, TbM, TaT);
					     Tea = Te8 + Te9;
					     Tev = Te8 - Te9;
					}
				   }
				   Ter = TdA - TdD;
				   TdE = TdA + TdD;
				   TdJ = TdF + TdI;
				   Teq = TdI - TdF;
				   Teu = TdM - TdP;
				   TdQ = TdM + TdP;
				   Tei = FNMS(KP707106781, Tea, Te7);
				   Teb = FMA(KP707106781, Tea, Te7);
				   Teg = FMA(KP414213562, TdE, TdJ);
				   TdK = FNMS(KP414213562, TdJ, TdE);
				   Te1 = TdV + Te0;
				   Tey = Te0 - TdV;
			      }
			      {
				   E Tcd, TbZ, Tcg, Tc9, Tce, Tc2;
				   {
					E Teh, Tel, Ted, TdL, Tej, Te2;
					Teh = FNMS(KP1_847759065, Teg, Tef);
					Tel = FMA(KP1_847759065, Teg, Tef);
					Ted = FNMS(KP1_847759065, TdK, Tdz);
					TdL = FMA(KP1_847759065, TdK, Tdz);
					Tej = FNMS(KP707106781, Te1, TdQ);
					Te2 = FMA(KP707106781, Te1, TdQ);
					{
					     E Tek, Tem, Tee, Tec;
					     Tek = FNMS(KP668178637, Tej, Tei);
					     Tem = FMA(KP668178637, Tei, Tej);
					     Tee = FMA(KP198912367, Te2, Teb);
					     Tec = FNMS(KP198912367, Teb, Te2);
					     Tcd = FMA(KP1_847759065, TbY, TbX);
					     TbZ = FNMS(KP1_847759065, TbY, TbX);
					     R0[WS(rs, 42)] = FMA(KP1_662939224, Tek, Teh);
					     R0[WS(rs, 10)] = FNMS(KP1_662939224, Tek, Teh);
					     R0[WS(rs, 58)] = FMA(KP1_662939224, Tem, Tel);
					     R0[WS(rs, 26)] = FNMS(KP1_662939224, Tem, Tel);
					     R0[WS(rs, 50)] = FMA(KP1_961570560, Tee, Ted);
					     R0[WS(rs, 18)] = FNMS(KP1_961570560, Tee, Ted);
					     R0[WS(rs, 2)] = FMA(KP1_961570560, Tec, TdL);
					     R0[WS(rs, 34)] = FNMS(KP1_961570560, Tec, TdL);
					}
				   }
				   Tcg = FMA(KP923879532, Tc8, Tc7);
				   Tc9 = FNMS(KP923879532, Tc8, Tc7);
				   Tce = FMA(KP668178637, Tc0, Tc1);
				   Tc2 = FNMS(KP668178637, Tc1, Tc0);
				   {
					E Tcf, Tcj, Tcb, Tc3, Tch, Tc6;
					Tcf = FNMS(KP1_662939224, Tce, Tcd);
					Tcj = FMA(KP1_662939224, Tce, Tcd);
					Tcb = FMA(KP1_662939224, Tc2, TbZ);
					Tc3 = FNMS(KP1_662939224, Tc2, TbZ);
					Tch = FMA(KP923879532, Tc5, Tc4);
					Tc6 = FNMS(KP923879532, Tc5, Tc4);
					{
					     E Tci, Tck, Tcc, Tca;
					     Tci = FNMS(KP303346683, Tch, Tcg);
					     Tck = FMA(KP303346683, Tcg, Tch);
					     Tcc = FMA(KP534511135, Tc6, Tc9);
					     Tca = FNMS(KP534511135, Tc9, Tc6);
					     TeD = FMA(KP1_414213562, Teo, Ten);
					     Tep = FNMS(KP1_414213562, Teo, Ten);
					     R0[WS(rs, 45)] = FMA(KP1_913880671, Tci, Tcf);
					     R0[WS(rs, 13)] = FNMS(KP1_913880671, Tci, Tcf);
					     R0[WS(rs, 61)] = FMA(KP1_913880671, Tck, Tcj);
					     R0[WS(rs, 29)] = FNMS(KP1_913880671, Tck, Tcj);
					     R0[WS(rs, 53)] = FMA(KP1_763842528, Tcc, Tcb);
					     R0[WS(rs, 21)] = FNMS(KP1_763842528, Tcc, Tcb);
					     R0[WS(rs, 5)] = FMA(KP1_763842528, Tca, Tc3);
					     R0[WS(rs, 37)] = FNMS(KP1_763842528, Tca, Tc3);
					}
				   }
				   TeG = FMA(KP707106781, Tey, Tex);
				   Tez = FNMS(KP707106781, Tey, Tex);
				   TeE = FMA(KP414213562, Teq, Ter);
				   Tes = FNMS(KP414213562, Ter, Teq);
			      }
			 }
		    }
		    {
			 E T5L, T5M, T61, T62;
			 {
			      E Td3, Td4, Tdd, TcV, Tcr, Tdb, Td7, Td6, Tda, TcC, TcY, TcR, TcW, Tcy, TcJ;
			      E Tde;
			      {
				   E TcN, TcQ, Tcu, Tcx;
				   {
					E Tcn, TeF, TeJ, TeB, Tet, TeH, Tew, Tcq;
					Td3 = FMA(KP1_414213562, Tcm, Tcl);
					Tcn = FNMS(KP1_414213562, Tcm, Tcl);
					TeF = FNMS(KP1_847759065, TeE, TeD);
					TeJ = FMA(KP1_847759065, TeE, TeD);
					TeB = FMA(KP1_847759065, Tes, Tep);
					Tet = FNMS(KP1_847759065, Tes, Tep);
					TeH = FMA(KP707106781, Tev, Teu);
					Tew = FNMS(KP707106781, Tev, Teu);
					Tcq = FNMS(KP414213562, Tcp, Tco);
					Td4 = FMA(KP414213562, Tco, Tcp);
					Tdd = FMA(KP707106781, TcM, TcL);
					TcN = FNMS(KP707106781, TcM, TcL);
					{
					     E TeI, TeK, TeC, TeA;
					     TeI = FNMS(KP198912367, TeH, TeG);
					     TeK = FMA(KP198912367, TeG, TeH);
					     TeC = FMA(KP668178637, Tew, Tez);
					     TeA = FNMS(KP668178637, Tez, Tew);
					     TcV = FMA(KP1_847759065, Tcq, Tcn);
					     Tcr = FNMS(KP1_847759065, Tcq, Tcn);
					     R0[WS(rs, 46)] = FMA(KP1_961570560, TeI, TeF);
					     R0[WS(rs, 14)] = FNMS(KP1_961570560, TeI, TeF);
					     R0[WS(rs, 62)] = FMA(KP1_961570560, TeK, TeJ);
					     R0[WS(rs, 30)] = FNMS(KP1_961570560, TeK, TeJ);
					     R0[WS(rs, 54)] = FMA(KP1_662939224, TeC, TeB);
					     R0[WS(rs, 22)] = FNMS(KP1_662939224, TeC, TeB);
					     R0[WS(rs, 6)] = FMA(KP1_662939224, TeA, Tet);
					     R0[WS(rs, 38)] = FNMS(KP1_662939224, TeA, Tet);
					     TcQ = TcO - TcP;
					     Tdb = TcO + TcP;
					}
				   }
				   Td7 = FMA(KP707106781, Tct, Tcs);
				   Tcu = FNMS(KP707106781, Tct, Tcs);
				   Tcx = FNMS(KP707106781, Tcw, Tcv);
				   Td6 = FMA(KP707106781, Tcw, Tcv);
				   Tda = FMA(KP707106781, TcB, TcA);
				   TcC = FNMS(KP707106781, TcB, TcA);
				   TcY = FNMS(KP923879532, TcQ, TcN);
				   TcR = FMA(KP923879532, TcQ, TcN);
				   TcW = FMA(KP668178637, Tcu, Tcx);
				   Tcy = FNMS(KP668178637, Tcx, Tcu);
				   TcJ = TcF - TcI;
				   Tde = TcI + TcF;
			      }
			      {
				   E Tdj, Td5, Tdm, Tdf, Tdk, Td8;
				   {
					E TcX, Td1, TcT, Tcz, TcZ, TcK;
					TcX = FNMS(KP1_662939224, TcW, TcV);
					Td1 = FMA(KP1_662939224, TcW, TcV);
					TcT = FNMS(KP1_662939224, Tcy, Tcr);
					Tcz = FMA(KP1_662939224, Tcy, Tcr);
					TcZ = FNMS(KP923879532, TcJ, TcC);
					TcK = FMA(KP923879532, TcJ, TcC);
					{
					     E Td0, Td2, TcU, TcS;
					     Td0 = FNMS(KP534511135, TcZ, TcY);
					     Td2 = FMA(KP534511135, TcY, TcZ);
					     TcU = FMA(KP303346683, TcK, TcR);
					     TcS = FNMS(KP303346683, TcR, TcK);
					     Tdj = FMA(KP1_847759065, Td4, Td3);
					     Td5 = FNMS(KP1_847759065, Td4, Td3);
					     R0[WS(rs, 43)] = FMA(KP1_763842528, Td0, TcX);
					     R0[WS(rs, 11)] = FNMS(KP1_763842528, Td0, TcX);
					     R0[WS(rs, 59)] = FMA(KP1_763842528, Td2, Td1);
					     R0[WS(rs, 27)] = FNMS(KP1_763842528, Td2, Td1);
					     R0[WS(rs, 51)] = FMA(KP1_913880671, TcU, TcT);
					     R0[WS(rs, 19)] = FNMS(KP1_913880671, TcU, TcT);
					     R0[WS(rs, 3)] = FMA(KP1_913880671, TcS, Tcz);
					     R0[WS(rs, 35)] = FNMS(KP1_913880671, TcS, Tcz);
					}
				   }
				   Tdm = FMA(KP923879532, Tde, Tdd);
				   Tdf = FNMS(KP923879532, Tde, Tdd);
				   Tdk = FMA(KP198912367, Td6, Td7);
				   Td8 = FNMS(KP198912367, Td7, Td6);
				   {
					E T5F, T2R, T5G, T3K, T64, T5S, T5X, T5x, T5U, T4q, T4R, T63, T5P, T5i, T5V;
					E T5A;
					{
					     E T5N, T5O, T5R, T3w, T3J, T5Q, T5y, T5z;
					     {
						  E T2p, Tdl, Tdp, Tdh, Td9, Tdn, Tdc, T2Q;
						  T5N = FNMS(KP1_847759065, T2o, T2d);
						  T2p = FMA(KP1_847759065, T2o, T2d);
						  Tdl = FNMS(KP1_961570560, Tdk, Tdj);
						  Tdp = FMA(KP1_961570560, Tdk, Tdj);
						  Tdh = FMA(KP1_961570560, Td8, Td5);
						  Td9 = FNMS(KP1_961570560, Td8, Td5);
						  Tdn = FMA(KP923879532, Tdb, Tda);
						  Tdc = FNMS(KP923879532, Tdb, Tda);
						  T2Q = FNMS(KP198912367, T2P, T2G);
						  T5O = FMA(KP198912367, T2G, T2P);
						  T5R = FNMS(KP923879532, T3v, T38);
						  T3w = FMA(KP923879532, T3v, T38);
						  {
						       E Tdo, Tdq, Tdi, Tdg;
						       Tdo = FNMS(KP098491403, Tdn, Tdm);
						       Tdq = FMA(KP098491403, Tdm, Tdn);
						       Tdi = FMA(KP820678790, Tdc, Tdf);
						       Tdg = FNMS(KP820678790, Tdf, Tdc);
						       T5F = FNMS(KP1_961570560, T2Q, T2p);
						       T2R = FMA(KP1_961570560, T2Q, T2p);
						       R0[WS(rs, 47)] = FMA(KP1_990369453, Tdo, Tdl);
						       R0[WS(rs, 15)] = FNMS(KP1_990369453, Tdo, Tdl);
						       R0[WS(rs, 63)] = FMA(KP1_990369453, Tdq, Tdp);
						       R0[WS(rs, 31)] = FNMS(KP1_990369453, Tdq, Tdp);
						       R0[WS(rs, 55)] = FMA(KP1_546020906, Tdi, Tdh);
						       R0[WS(rs, 23)] = FNMS(KP1_546020906, Tdi, Tdh);
						       R0[WS(rs, 7)] = FMA(KP1_546020906, Tdg, Td9);
						       R0[WS(rs, 39)] = FNMS(KP1_546020906, Tdg, Td9);
						       T3J = FMA(KP923879532, T3I, T3F);
						       T5Q = FNMS(KP923879532, T3I, T3F);
						  }
					     }
					     T5G = FMA(KP098491403, T3w, T3J);
					     T3K = FNMS(KP098491403, T3J, T3w);
					     T64 = FMA(KP820678790, T5Q, T5R);
					     T5S = FNMS(KP820678790, T5R, T5Q);
					     T5X = FNMS(KP923879532, T5w, T5t);
					     T5x = FMA(KP923879532, T5w, T5t);
					     T5U = FNMS(KP923879532, T4p, T42);
					     T4q = FMA(KP923879532, T4p, T42);
					     T4R = FNMS(KP198912367, T4Q, T4H);
					     T5y = FMA(KP198912367, T4H, T4Q);
					     T63 = FMA(KP1_961570560, T5O, T5N);
					     T5P = FNMS(KP1_961570560, T5O, T5N);
					     T5z = FNMS(KP198912367, T58, T5h);
					     T5i = FMA(KP198912367, T5h, T58);