n1_64.c

快速fft变换

				   TaA = TN - TQ;
				   TR = TN + TQ;
				   TY = TU + TX;
				   Tax = TX - TU;
				   T44 = T3T - T3S;
				   T3U = T3S + T3T;
				   Tdf = TR - TY;
				   TZ = TR + TY;
				   Tay = T2Q - T2T;
				   T2U = T2Q + T2T;
				   TaB = T2X - T30;
				   T31 = T2X + T30;
				   Tch = Tay - Tax;
				   Taz = Tax + Tay;
				   Tci = TaA - TaB;
				   TaC = TaA + TaB;
			      }
			      {
				   E T7T, T7X, T7U, T3R, T3Y, T7W;
				   T7T = T3Q - T3P;
				   T3R = T3P + T3Q;
				   T3Y = T3U - T3X;
				   T7X = T3U + T3X;
				   Tdg = T2U - T31;
				   T32 = T2U + T31;
				   T6J = FMA(KP707106781, T3Y, T3R);
				   T3Z = FNMS(KP707106781, T3Y, T3R);
				   T7U = T44 + T43;
				   T45 = T43 - T44;
				   T42 = T40 - T41;
				   T7W = T40 + T41;
				   T9n = FNMS(KP707106781, T7U, T7T);
				   T7V = FMA(KP707106781, T7U, T7T);
				   T9o = FNMS(KP707106781, T7X, T7W);
				   T7Y = FMA(KP707106781, T7X, T7W);
			      }
			 }
		    }
		    {
			 E T4P, T14, TaH, T4d, T4a, T17, TaI, T4S, T1e, TaZ, T4f, T1b, T4U, T4o, TaY;
			 E T4i;
			 {
			      E T4b, T4c, T15, T16, T12, T13, T4Q, T4R;
			      T12 = ri[WS(is, 1)];
			      T13 = ri[WS(is, 33)];
			      T4b = ii[WS(is, 1)];
			      T6K = FMA(KP707106781, T45, T42);
			      T46 = FNMS(KP707106781, T45, T42);
			      T4P = T12 - T13;
			      T14 = T12 + T13;
			      T4c = ii[WS(is, 33)];
			      T15 = ri[WS(is, 17)];
			      T16 = ri[WS(is, 49)];
			      T4Q = ii[WS(is, 17)];
			      TaH = T4b + T4c;
			      T4d = T4b - T4c;
			      T4a = T15 - T16;
			      T17 = T15 + T16;
			      T4R = ii[WS(is, 49)];
			      {
				   E T4l, T4k, T4m, T1c, T1d;
				   T1c = ri[WS(is, 57)];
				   T1d = ri[WS(is, 25)];
				   T4l = ii[WS(is, 57)];
				   TaI = T4Q + T4R;
				   T4S = T4Q - T4R;
				   T4k = T1c - T1d;
				   T1e = T1c + T1d;
				   T4m = ii[WS(is, 25)];
				   {
					E T4g, T4h, T19, T1a, T4n;
					T19 = ri[WS(is, 9)];
					T1a = ri[WS(is, 41)];
					T4n = T4l - T4m;
					TaZ = T4l + T4m;
					T4g = ii[WS(is, 9)];
					T4f = T19 - T1a;
					T1b = T19 + T1a;
					T4h = ii[WS(is, 41)];
					T4U = T4k + T4n;
					T4o = T4k - T4n;
					TaY = T4g + T4h;
					T4i = T4g - T4h;
				   }
			      }
			 }
			 {
			      E TaG, T4V, TaJ, T4e, T4p, T8e, T8d, T4T, T4W;
			      {
				   E TaX, T4j, Tb0, Tdl, T18, T1f, Tdk;
				   TaX = T14 - T17;
				   T18 = T14 + T17;
				   T1f = T1b + T1e;
				   TaG = T1e - T1b;
				   T4j = T4f + T4i;
				   T4V = T4i - T4f;
				   Tdp = T18 - T1f;
				   T1g = T18 + T1f;
				   Tb0 = TaY - TaZ;
				   Tdl = TaY + TaZ;
				   TaJ = TaH - TaI;
				   Tdk = TaH + TaI;
				   T82 = T4d - T4a;
				   T4e = T4a + T4d;
				   Tcm = TaX - Tb0;
				   Tb1 = TaX + Tb0;
				   Tej = Tdk + Tdl;
				   Tdm = Tdk - Tdl;
				   T4p = T4j - T4o;
				   T8e = T4j + T4o;
			      }
			      T8d = T4P + T4S;
			      T4T = T4P - T4S;
			      T4W = T4U - T4V;
			      T83 = T4V + T4U;
			      T6R = FMA(KP707106781, T4p, T4e);
			      T4q = FNMS(KP707106781, T4p, T4e);
			      T6O = FMA(KP707106781, T4W, T4T);
			      T4X = FNMS(KP707106781, T4W, T4T);
			      T9s = FNMS(KP707106781, T8e, T8d);
			      T8f = FMA(KP707106781, T8e, T8d);
			      Tcp = TaJ - TaG;
			      TaK = TaG + TaJ;
			 }
		    }
		    {
			 E T85, T4L, TaO, T1n, Tdq, TaN, T86, T4G, T4r, T1q, T4s, TaR, T4z, T4w, T1t;
			 E T4t;
			 {
			      E T4C, T1j, T4D, TaL, T4K, T4H, T1m, T4E;
			      {
				   E T4I, T4J, T1h, T1i, T1k, T1l;
				   T1h = ri[WS(is, 5)];
				   T1i = ri[WS(is, 37)];
				   T4I = ii[WS(is, 5)];
				   T9v = FNMS(KP707106781, T83, T82);
				   T84 = FMA(KP707106781, T83, T82);
				   T4C = T1h - T1i;
				   T1j = T1h + T1i;
				   T4J = ii[WS(is, 37)];
				   T1k = ri[WS(is, 21)];
				   T1l = ri[WS(is, 53)];
				   T4D = ii[WS(is, 21)];
				   TaL = T4I + T4J;
				   T4K = T4I - T4J;
				   T4H = T1k - T1l;
				   T1m = T1k + T1l;
				   T4E = ii[WS(is, 53)];
			      }
			      {
				   E T4x, T4y, T1r, T1s;
				   {
					E T1o, T4F, TaM, T1p;
					T1o = ri[WS(is, 61)];
					T85 = T4K - T4H;
					T4L = T4H + T4K;
					TaO = T1j - T1m;
					T1n = T1j + T1m;
					T4F = T4D - T4E;
					TaM = T4D + T4E;
					T1p = ri[WS(is, 29)];
					T4x = ii[WS(is, 61)];
					Tdq = TaL + TaM;
					TaN = TaL - TaM;
					T86 = T4C + T4F;
					T4G = T4C - T4F;
					T4r = T1o - T1p;
					T1q = T1o + T1p;
					T4y = ii[WS(is, 29)];
				   }
				   T1r = ri[WS(is, 13)];
				   T1s = ri[WS(is, 45)];
				   T4s = ii[WS(is, 13)];
				   TaR = T4x + T4y;
				   T4z = T4x - T4y;
				   T4w = T1r - T1s;
				   T1t = T1r + T1s;
				   T4t = ii[WS(is, 45)];
			      }
			 }
			 {
			      E T88, TaP, T89, TaU, T4Z, T4B, T4M, T4Y, T8g, T8h;
			      {
				   E T4A, Tb2, Tdr, T4v, Tb3;
				   {
					E TaQ, T1u, T4u, TaS, TaT;
					T88 = T4z - T4w;
					T4A = T4w + T4z;
					TaQ = T1q - T1t;
					T1u = T1q + T1t;
					T4u = T4s - T4t;
					TaS = T4s + T4t;
					Tb2 = TaO + TaN;
					TaP = TaN - TaO;
					Tdr = TaR + TaS;
					TaT = TaR - TaS;
					T89 = T4r + T4u;
					T4v = T4r - T4u;
					Tdn = T1u - T1n;
					T1v = T1n + T1u;
					Tb3 = TaQ - TaT;
					TaU = TaQ + TaT;
				   }
				   T4Z = FNMS(KP414213562, T4v, T4A);
				   T4B = FMA(KP414213562, T4A, T4v);
				   Tcq = Tb2 - Tb3;
				   Tb4 = Tb2 + Tb3;
				   Tek = Tdq + Tdr;
				   Tds = Tdq - Tdr;
				   T4M = FNMS(KP414213562, T4L, T4G);
				   T4Y = FMA(KP414213562, T4G, T4L);
			      }
			      T87 = FNMS(KP414213562, T86, T85);
			      T8g = FMA(KP414213562, T85, T86);
			      T6P = T4M + T4B;
			      T4N = T4B - T4M;
			      T6S = T4Y + T4Z;
			      T50 = T4Y - T4Z;
			      T8h = FNMS(KP414213562, T88, T89);
			      T8a = FMA(KP414213562, T89, T88);
			      T9w = T8g - T8h;
			      T8i = T8g + T8h;
			      Tcn = TaU - TaP;
			      TaV = TaP + TaU;
			 }
		    }
		    {
			 E T8o, T5E, Tbf, T1S, TdB, Tbe, T8p, T5z, T5k, T1V, T5l, Tbi, T5s, T5p, T1Y;
			 E T5m;
			 {
			      E T5v, T1O, T5w, Tbc, T5D, T5A, T1R, T5x;
			      {
				   E T5B, T5C, T1M, T1N, T1P, T1Q;
				   T1M = ri[WS(is, 3)];
				   T1N = ri[WS(is, 35)];
				   T5B = ii[WS(is, 3)];
				   T9t = T8a - T87;
				   T8b = T87 + T8a;
				   T5v = T1M - T1N;
				   T1O = T1M + T1N;
				   T5C = ii[WS(is, 35)];
				   T1P = ri[WS(is, 19)];
				   T1Q = ri[WS(is, 51)];
				   T5w = ii[WS(is, 19)];
				   Tbc = T5B + T5C;
				   T5D = T5B - T5C;
				   T5A = T1P - T1Q;
				   T1R = T1P + T1Q;
				   T5x = ii[WS(is, 51)];
			      }
			      {
				   E T5q, T5r, T1W, T1X;
				   {
					E T1T, T5y, Tbd, T1U;
					T1T = ri[WS(is, 59)];
					T8o = T5D - T5A;
					T5E = T5A + T5D;
					Tbf = T1O - T1R;
					T1S = T1O + T1R;
					T5y = T5w - T5x;
					Tbd = T5w + T5x;
					T1U = ri[WS(is, 27)];
					T5q = ii[WS(is, 59)];
					TdB = Tbc + Tbd;
					Tbe = Tbc - Tbd;
					T8p = T5v + T5y;
					T5z = T5v - T5y;
					T5k = T1T - T1U;
					T1V = T1T + T1U;
					T5r = ii[WS(is, 27)];
				   }
				   T1W = ri[WS(is, 11)];
				   T1X = ri[WS(is, 43)];
				   T5l = ii[WS(is, 11)];
				   Tbi = T5q + T5r;
				   T5s = T5q - T5r;
				   T5p = T1W - T1X;
				   T1Y = T1W + T1X;
				   T5m = ii[WS(is, 43)];
			      }
			 }
			 {
			      E T8r, Tbg, T8s, Tbl, T5S, T5u, T5F, T5R, T8z, T8A;
			      {
				   E T5t, Tbt, TdC, T5o, Tbu;
				   {
					E Tbh, T1Z, T5n, Tbj, Tbk;
					T8r = T5s - T5p;
					T5t = T5p + T5s;
					Tbh = T1V - T1Y;
					T1Z = T1V + T1Y;
					T5n = T5l - T5m;
					Tbj = T5l + T5m;
					Tbt = Tbf + Tbe;
					Tbg = Tbe - Tbf;
					TdC = Tbi + Tbj;
					Tbk = Tbi - Tbj;
					T8s = T5k + T5n;
					T5o = T5k - T5n;
					Tdy = T1Z - T1S;
					T20 = T1S + T1Z;
					Tbu = Tbh - Tbk;
					Tbl = Tbh + Tbk;
				   }
				   T5S = FNMS(KP414213562, T5o, T5t);
				   T5u = FMA(KP414213562, T5t, T5o);
				   Tcx = Tbt - Tbu;
				   Tbv = Tbt + Tbu;
				   Tep = TdB + TdC;
				   TdD = TdB - TdC;
				   T5F = FNMS(KP414213562, T5E, T5z);
				   T5R = FMA(KP414213562, T5z, T5E);
			      }
			      T8q = FNMS(KP414213562, T8p, T8o);
			      T8z = FMA(KP414213562, T8o, T8p);
			      T6W = T5F + T5u;
			      T5G = T5u - T5F;
			      T6Z = T5R + T5S;
			      T5T = T5R - T5S;
			      T8A = FNMS(KP414213562, T8r, T8s);
			      T8t = FMA(KP414213562, T8s, T8r);
			      T9D = T8z - T8A;
			      T8B = T8z + T8A;
			      Tcu = Tbl - Tbg;
			      Tbm = Tbg + Tbl;
			 }
		    }
	       }
	       {
		    E T9A, T8u, TbE, Tao, Td7, Td8;
		    {
			 E Teq, Ten, Tex, Teh, TeB, Tev, Tey, Tem, Te9, Tec;
			 {
			      E Tef, Teu, Tel, T11, Tei, Tet, T2y, TeI, T23, T22, T33, Teg, TeD, TeG, T34;
			      E TeH;
			      {
				   E TeE, TeF, Tv, T10, T1w, T21;
				   Tef = Tf - Tu;
				   Tv = Tf + Tu;
				   T10 = TK + TZ;
				   Teu = TZ - TK;
				   Tel = Tej - Tek;
				   TeE = Tej + Tek;
				   T9A = T8t - T8q;
				   T8u = T8q + T8t;
				   TeD = Tv - T10;
				   T11 = Tv + T10;
				   TeF = Teo + Tep;
				   Teq = Teo - Tep;
				   Tei = T1g - T1v;
				   T1w = T1g + T1v;
				   T21 = T1L + T20;
				   Ten = T1L - T20;
				   Tet = T2i - T2x;
				   T2y = T2i + T2x;
				   TeI = TeE + TeF;
				   TeG = TeE - TeF;
				   T23 = T21 - T1w;
				   T22 = T1w + T21;
				   T33 = T2N + T32;
				   Teg = T2N - T32;
			      }
			      ro[WS(os, 16)] = TeD + TeG;
			      ro[WS(os, 48)] = TeD - TeG;
			      ro[0] = T11 + T22;
			      ro[WS(os, 32)] = T11 - T22;
			      T34 = T2y - T33;
			      TeH = T2y + T33;
			      io[0] = TeH + TeI;
			      io[WS(os, 32)] = TeH - TeI;
			      io[WS(os, 48)] = T34 - T23;
			      io[WS(os, 16)] = T23 + T34;
			      Tex = Tef - Teg;
			      Teh = Tef + Teg;
			      TeB = Teu + Tet;
			      Tev = Tet - Teu;
			      Tey = Tel - Tei;
			      Tem = Tei + Tel;
			 }
			 {
			      E TdV, Tdb, TdJ, Te5, TdE, Tdz, Te6, Tdi, Teb, Te3, TdZ, TdY, TdW, TdM, TdR;
			      E Tdu;
			      {
				   E TdL, Tde, Tdh, TdK, Tez, Ter;
				   TdV = Td9 + Tda;
				   Tdb = Td9 - Tda;
				   TdJ = TdH - TdI;
				   Te5 = TdI + TdH;
				   Tez = Ten + Teq;
				   Ter = Ten - Teq;
				   TdL = Tdd + Tdc;
				   Tde = Tdc - Tdd;
				   {
					E TeA, TeC, Tew, Tes;
					TeA = Tey - Tez;
					TeC = Tey + Tez;
					Tew = Ter - Tem;
					Tes = Tem + Ter;
					ro[WS(os, 24)] = FMA(KP707106781, TeA, Tex);
					ro[WS(os, 56)] = FNMS(KP707106781, TeA, Tex);
					io[WS(os, 8)] = FMA(KP707106781, TeC, TeB);
					io[WS(os, 40)] = FNMS(KP707106781, TeC, TeB);
					io[WS(os, 24)] = FMA(KP707106781, Tew, Tev);
					io[WS(os, 56)] = FNMS(KP707106781, Tew, Tev);
					ro[WS(os, 8)] = FMA(KP707106781, Tes, Teh);
					ro[WS(os, 40)] = FNMS(KP707106781, Tes, Teh);
					Tdh = Tdf + Tdg;
					TdK = Tdf - Tdg;
				   }
				   {
					E Te1, Te2, Tdo, Tdt;
					TdE = TdA - TdD;
					Te1 = TdA + TdD;
					Te2 = Tdy + Tdx;
					Tdz = Tdx - Tdy;
					Te6 = Tde + Tdh;
					Tdi = Tde - Tdh;
					Teb = FMA(KP414213562, Te1, Te2);
					Te3 = FNMS(KP414213562, Te2, Te1);
					TdZ = Tdn + Tdm;
					Tdo = Tdm - Tdn;
					Tdt = Tdp - Tds;
					TdY = Tdp + Tds;
					TdW = TdL + TdK;
					TdM = TdK - TdL;
					TdR = FNMS(KP414213562, Tdo, Tdt);
					Tdu = FMA(KP414213562, Tdt, Tdo);
				   }
			      }
			      {
				   E TdT, Tea, Te0, TdU;
				   {
					E Tdj, TdQ, TdF, TdP, TdN, TdS, TdO, TdG;
					TdT = FNMS(KP707106781, Tdi, Tdb);
					Tdj = FMA(KP707106781, Tdi, Tdb);
					Tea = FNMS(KP414213562, TdY, TdZ);
					Te0 = FMA(KP414213562, TdZ, TdY);
					TdQ = FMA(KP414213562, Tdz, TdE);
					TdF = FNMS(KP414213562, TdE, Tdz);
					TdP = FMA(KP707106781, TdM, TdJ);
					TdN = FNMS(KP707106781, TdM, TdJ);
					TdS = TdQ - TdR;
					TdU = TdR + TdQ;
					TdO = Tdu + TdF;
					TdG = Tdu - TdF;
					io[WS(os, 12)] = FMA(KP923879532, TdS, TdP);
					io[WS(os, 44)] = FNMS(KP923879532, TdS, TdP);
					ro[WS(os, 12)] = FMA(KP923879532, TdG, Tdj);
					ro[WS(os, 44)] = FNMS(KP923879532, TdG, Tdj);
					io[WS(os, 60)] = FMA(KP923879532, TdO, TdN);
					io[WS(os, 28)] = FNMS(KP923879532, TdO, TdN);
				   }
				   {
					E Te8, Te7, Ted, Tee, TdX, Te4;
					Te9 = FNMS(KP707106781, TdW, TdV);
					TdX = FMA(KP707106781, TdW, TdV);
					Te4 = Te0 + Te3;
					Te8 = Te3 - Te0;
					Te7 = FNMS(KP707106781, Te6, Te5);
					Ted = FMA(KP707106781, Te6, Te5);
					ro[WS(os, 60)] = FMA(KP923879532, TdU, TdT);
					ro[WS(os, 28)] = FNMS(KP923879532, TdU, TdT);
					ro[WS(os, 4)] = FMA(KP923879532, Te4, TdX);
					ro[WS(os, 36)] = FNMS(KP923879532, Te4, TdX);
					Tee = Tea + Teb;
					Tec = Tea - Teb;