r2cf_128.c

快速fft变换

				   TaR = TaP + TaQ;
				   Tdv = TaP - TaQ;
				   TdI = TaT - TaS;
				   TaU = TaS + TaT;
				   T4U = T4S - T4T;
				   Tbd = T4T + T4S;
				   T8m = FNMS(KP707106781, T5B, T5A);
				   T5C = FMA(KP707106781, T5B, T5A);
				   {
					E Tbe, Tdx, T5E, T4V;
					Tbe = Tbc + Tbd;
					Tdx = Tbc - Tbd;
					T5E = FMA(KP414213562, T4R, T4U);
					T4V = FNMS(KP414213562, T4U, T4R);
					{
					     E Tdy, TdJ, T5F, T8c, T4W;
					     Tdy = Tdw + Tdx;
					     TdJ = Tdx - Tdw;
					     T5F = T5D + T5E;
					     T8c = T5E - T5D;
					     T8n = T4V - T4O;
					     T4W = T4O + T4V;
					     TeO = FNMS(KP707106781, TdJ, TdI);
					     TdK = FMA(KP707106781, TdJ, TdI);
					     TeL = FNMS(KP707106781, Tdy, Tdv);
					     Tdz = FMA(KP707106781, Tdy, Tdv);
					     T9p = FNMS(KP923879532, T8c, T8b);
					     T8d = FMA(KP923879532, T8c, T8b);
					     T6O = FNMS(KP923879532, T5F, T5C);
					     T5G = FMA(KP923879532, T5F, T5C);
					     T6L = FNMS(KP923879532, T4W, T4H);
					     T4X = FMA(KP923879532, T4W, T4H);
					}
					Tc3 = TaR + TaU;
					TaV = TaR - TaU;
					Tc4 = Tbh + Tbe;
					Tbi = Tbe - Tbh;
				   }
			      }
			 }
			 {
			      E Tai, T3k, Taj, T4l, Tal, Tam, T4m, T3r, T3D, TaF, T3C, Tdb, TaK, T3z, T4o;
			      E T3E;
			      {
				   E T4j, T4k, T3i, T3j;
				   T3i = R1[0];
				   T3j = R1[WS(rs, 32)];
				   T4j = R1[WS(rs, 16)];
				   T9s = FMA(KP923879532, T8n, T8m);
				   T8o = FNMS(KP923879532, T8n, T8m);
				   Tai = T3i + T3j;
				   T3k = T3i - T3j;
				   T4k = R1[WS(rs, 48)];
				   {
					E T3o, T3n, T3p, T3l, T3m;
					T3l = R1[WS(rs, 8)];
					T3m = R1[WS(rs, 40)];
					T3o = R1[WS(rs, 56)];
					Taj = T4j + T4k;
					T4l = T4j - T4k;
					T3n = T3l - T3m;
					Tal = T3l + T3m;
					T3p = R1[WS(rs, 24)];
					{
					     E T3w, TaI, T3v, T3x, T3t, T3u, T3q;
					     T3t = R1[WS(rs, 4)];
					     T3u = R1[WS(rs, 36)];
					     T3q = T3o - T3p;
					     Tam = T3o + T3p;
					     T3w = R1[WS(rs, 20)];
					     TaI = T3t + T3u;
					     T3v = T3t - T3u;
					     T4m = T3n - T3q;
					     T3r = T3n + T3q;
					     T3x = R1[WS(rs, 52)];
					     {
						  E T3A, T3B, TaJ, T3y;
						  T3A = R1[WS(rs, 60)];
						  T3B = R1[WS(rs, 28)];
						  T3D = R1[WS(rs, 12)];
						  TaJ = T3w + T3x;
						  T3y = T3w - T3x;
						  TaF = T3A + T3B;
						  T3C = T3A - T3B;
						  Tdb = TaI - TaJ;
						  TaK = TaI + TaJ;
						  T3z = FNMS(KP414213562, T3y, T3v);
						  T4o = FMA(KP414213562, T3v, T3y);
						  T3E = R1[WS(rs, 44)];
					     }
					}
				   }
			      }
			      {
				   E T3s, T7S, Tak, Tda, Tdn, Tan, T3F, TaG, T4n;
				   T3s = FMA(KP707106781, T3r, T3k);
				   T7S = FNMS(KP707106781, T3r, T3k);
				   Tak = Tai + Taj;
				   Tda = Tai - Taj;
				   Tdn = Tam - Tal;
				   Tan = Tal + Tam;
				   T3F = T3D - T3E;
				   TaG = T3D + T3E;
				   T83 = FNMS(KP707106781, T4m, T4l);
				   T4n = FMA(KP707106781, T4m, T4l);
				   {
					E TaH, Tdc, T4p, T3G;
					TaH = TaF + TaG;
					Tdc = TaF - TaG;
					T4p = FNMS(KP414213562, T3C, T3F);
					T3G = FMA(KP414213562, T3F, T3C);
					{
					     E Tdd, Tdo, T4q, T7T, T3H;
					     Tdd = Tdb + Tdc;
					     Tdo = Tdc - Tdb;
					     T4q = T4o + T4p;
					     T7T = T4o - T4p;
					     T84 = T3G - T3z;
					     T3H = T3z + T3G;
					     TeH = FNMS(KP707106781, Tdo, Tdn);
					     Tdp = FMA(KP707106781, Tdo, Tdn);
					     TeE = FNMS(KP707106781, Tdd, Tda);
					     Tde = FMA(KP707106781, Tdd, Tda);
					     T9i = FNMS(KP923879532, T7T, T7S);
					     T7U = FMA(KP923879532, T7T, T7S);
					     T6H = FNMS(KP923879532, T4q, T4n);
					     T4r = FMA(KP923879532, T4q, T4n);
					     T6E = FNMS(KP923879532, T3H, T3s);
					     T3I = FMA(KP923879532, T3H, T3s);
					}
					TbW = Tak + Tan;
					Tao = Tak - Tan;
					TbX = TaK + TaH;
					TaL = TaH - TaK;
				   }
			      }
			 }
			 {
			      E T7N, T3d, T3a, T7M;
			      {
				   E T2V, T1z, T38, T1C, T1G, T1J, T2Y, T39, T1W, T33, T1V, T3b, T32, Td0, T1S;
				   E T1X;
				   {
					E T1A, T1B, T1x, T1y;
					T1x = R0[WS(rs, 63)];
					T1y = R0[WS(rs, 31)];
					T1A = R0[WS(rs, 15)];
					T9l = FNMS(KP923879532, T84, T83);
					T85 = FMA(KP923879532, T84, T83);
					T2V = T1x - T1y;
					T1z = T1x + T1y;
					T1B = R0[WS(rs, 47)];
					{
					     E T1H, T2W, T1I, T1E, T1F;
					     T1E = R0[WS(rs, 7)];
					     T1F = R0[WS(rs, 39)];
					     T1H = R0[WS(rs, 55)];
					     T38 = T1A - T1B;
					     T1C = T1A + T1B;
					     T2W = T1E - T1F;
					     T1G = T1E + T1F;
					     T1I = R0[WS(rs, 23)];
					     {
						  E T1P, T30, T1O, T1Q, T1M, T1N, T2X;
						  T1M = R0[WS(rs, 3)];
						  T1N = R0[WS(rs, 35)];
						  T2X = T1H - T1I;
						  T1J = T1H + T1I;
						  T1P = R0[WS(rs, 19)];
						  T30 = T1M - T1N;
						  T1O = T1M + T1N;
						  T2Y = T2W + T2X;
						  T39 = T2W - T2X;
						  T1Q = R0[WS(rs, 51)];
						  {
						       E T1T, T1U, T31, T1R;
						       T1T = R0[WS(rs, 59)];
						       T1U = R0[WS(rs, 27)];
						       T1W = R0[WS(rs, 11)];
						       T31 = T1P - T1Q;
						       T1R = T1P + T1Q;
						       T33 = T1T - T1U;
						       T1V = T1T + T1U;
						       T3b = FMA(KP414213562, T30, T31);
						       T32 = FNMS(KP414213562, T31, T30);
						       Td0 = T1O - T1R;
						       T1S = T1O + T1R;
						       T1X = R0[WS(rs, 43)];
						  }
					     }
					}
				   }
				   {
					E TcZ, Td4, T35, T3c, T1Z, Td2, Td5, Td1;
					{
					     E T1D, T1K, T1Y, T34;
					     T1D = T1z + T1C;
					     TcZ = T1z - T1C;
					     Td4 = T1J - T1G;
					     T1K = T1G + T1J;
					     T1Y = T1W + T1X;
					     T34 = T1W - T1X;
					     T1L = T1D + T1K;
					     Tad = T1D - T1K;
					     T35 = FMA(KP414213562, T34, T33);
					     T3c = FNMS(KP414213562, T33, T34);
					     T1Z = T1V + T1Y;
					     Td1 = T1V - T1Y;
					}
					Td2 = Td0 + Td1;
					Td5 = Td1 - Td0;
					{
					     E T7J, T2Z, T36, T7K;
					     T7J = FNMS(KP707106781, T2Y, T2V);
					     T2Z = FMA(KP707106781, T2Y, T2V);
					     Tae = T1Z - T1S;
					     T20 = T1S + T1Z;
					     Tez = FNMS(KP707106781, Td5, Td4);
					     Td6 = FMA(KP707106781, Td5, Td4);
					     TeA = FNMS(KP707106781, Td2, TcZ);
					     Td3 = FMA(KP707106781, Td2, TcZ);
					     T36 = T32 + T35;
					     T7N = T35 - T32;
					     T7K = T3b - T3c;
					     T3d = T3b + T3c;
					     T3a = FMA(KP707106781, T39, T38);
					     T7M = FNMS(KP707106781, T39, T38);
					     T6x = FNMS(KP923879532, T36, T2Z);
					     T37 = FMA(KP923879532, T36, T2Z);
					     T9e = FNMS(KP923879532, T7K, T7J);
					     T7L = FMA(KP923879532, T7K, T7J);
					}
				   }
			      }
			      {
				   E Tav, T7V, T7W, TaC;
				   {
					E T3L, T3W, Tdf, Tar, T42, T4d, Tay, Tdi, T46, Tau, Tdg, T3X, T3S, Taz, T45;
					E T47, Taw, Tax;
					{
					     E T3J, T3K, T3U, T3V;
					     T3J = R1[WS(rs, 2)];
					     T9d = FNMS(KP923879532, T7N, T7M);
					     T7O = FMA(KP923879532, T7N, T7M);
					     T6w = FNMS(KP923879532, T3d, T3a);
					     T3e = FMA(KP923879532, T3d, T3a);
					     T3K = R1[WS(rs, 34)];
					     T3U = R1[WS(rs, 18)];
					     T3V = R1[WS(rs, 50)];
					     {
						  E T40, Tap, Taq, T41, T4b, T4c;
						  T40 = R1[WS(rs, 62)];
						  T3L = T3J - T3K;
						  Tap = T3J + T3K;
						  T3W = T3U - T3V;
						  Taq = T3U + T3V;
						  T41 = R1[WS(rs, 30)];
						  T4b = R1[WS(rs, 14)];
						  T4c = R1[WS(rs, 46)];
						  Tdf = Tap - Taq;
						  Tar = Tap + Taq;
						  T42 = T40 - T41;
						  Taw = T40 + T41;
						  Tax = T4b + T4c;
						  T4d = T4b - T4c;
					     }
					}
					{
					     E T3M, T3N, T3P, T3Q;
					     T3M = R1[WS(rs, 10)];
					     Tay = Taw + Tax;
					     Tdi = Taw - Tax;
					     T3N = R1[WS(rs, 42)];
					     T3P = R1[WS(rs, 58)];
					     T3Q = R1[WS(rs, 26)];
					     {
						  E T43, Tas, T3O, Tat, T3R, T44;
						  T43 = R1[WS(rs, 6)];
						  Tas = T3M + T3N;
						  T3O = T3M - T3N;
						  Tat = T3P + T3Q;
						  T3R = T3P - T3Q;
						  T44 = R1[WS(rs, 38)];
						  T46 = R1[WS(rs, 54)];
						  Tau = Tas + Tat;
						  Tdg = Tat - Tas;
						  T3X = T3O - T3R;
						  T3S = T3O + T3R;
						  Taz = T43 + T44;
						  T45 = T43 - T44;
						  T47 = R1[WS(rs, 22)];
					     }
					}
					{
					     E Tdq, Tdh, T49, T4e, Tdr, Tdk;
					     Tav = Tar - Tau;
					     TbZ = Tar + Tau;
					     {
						  E T3T, T3Y, TaA, T48, Tdj, TaB;
						  T3T = FMA(KP707106781, T3S, T3L);
						  T7V = FNMS(KP707106781, T3S, T3L);
						  T7W = FNMS(KP707106781, T3X, T3W);
						  T3Y = FMA(KP707106781, T3X, T3W);
						  TaA = T46 + T47;
						  T48 = T46 - T47;
						  Tdq = FNMS(KP414213562, Tdf, Tdg);
						  Tdh = FMA(KP414213562, Tdg, Tdf);
						  T3Z = FNMS(KP198912367, T3Y, T3T);
						  T4s = FMA(KP198912367, T3T, T3Y);
						  Tdj = TaA - Taz;
						  TaB = Taz + TaA;
						  T49 = T45 + T48;
						  T4e = T45 - T48;
						  TaC = Tay - TaB;
						  Tc0 = Tay + TaB;
						  Tdr = FMA(KP414213562, Tdi, Tdj);
						  Tdk = FNMS(KP414213562, Tdj, Tdi);
					     }
					     {
						  E T7Z, T7Y, T4f, T4a;
						  T7Z = FNMS(KP707106781, T4e, T4d);
						  T4f = FMA(KP707106781, T4e, T4d);
						  T4a = FMA(KP707106781, T49, T42);
						  T7Y = FNMS(KP707106781, T49, T42);
						  TeF = Tdr - Tdq;
						  Tds = Tdq + Tdr;
						  T4t = FNMS(KP198912367, T4a, T4f);
						  T4g = FMA(KP198912367, T4f, T4a);
						  T87 = FMA(KP668178637, T7Y, T7Z);
						  T80 = FNMS(KP668178637, T7Z, T7Y);
						  TeI = Tdh - Tdk;
						  Tdl = Tdh + Tdk;
					     }
					}
				   }
				   {
					E T50, T5b, TdA, TaY, T5h, T5s, Tb5, TdD, T5l, Tb1, TdB, T5c, T57, Tb6, T5k;
					E T5m, Tb3, Tb4;
					{
					     E T4Y, T4Z, T59, T5a;
					     T4Y = R1[WS(rs, 1)];
					     T86 = FNMS(KP668178637, T7V, T7W);
					     T7X = FMA(KP668178637, T7W, T7V);
					     TaM = TaC - Tav;
					     TaD = Tav + TaC;
					     T4Z = R1[WS(rs, 33)];
					     T59 = R1[WS(rs, 49)];
					     T5a = R1[WS(rs, 17)];
					     {
						  E T5f, TaW, TaX, T5g, T5q, T5r;
						  T5f = R1[WS(rs, 61)];
						  T50 = T4Y - T4Z;
						  TaW = T4Y + T4Z;
						  T5b = T59 - T5a;
						  TaX = T5a + T59;
						  T5g = R1[WS(rs, 29)];
						  T5q = R1[WS(rs, 45)];
						  T5r = R1[WS(rs, 13)];
						  TdA = TaW - TaX;
						  TaY = TaW + TaX;
						  T5h = T5f - T5g;
						  Tb3 = T5f + T5g;
						  Tb4 = T5r + T5q;
						  T5s = T5q - T5r;
					     }
					}
					{
					     E T51, T52, T54, T55;
					     T51 = R1[WS(rs, 9)];
					     Tb5 = Tb3 + Tb4;
					     TdD = Tb3 - Tb4;
					     T52 = R1[WS(rs, 41)];
					     T54 = R1[WS(rs, 57)];
					     T55 = R1[WS(rs, 25)];
					     {
						  E T5i, TaZ, T53, Tb0, T56, T5j;
						  T5i = R1[WS(rs, 5)];
						  TaZ = T51 + T52;
						  T53 = T51 - T52;
						  Tb0 = T54 + T55;
						  T56 = T54 - T55;
						  T5j = R1[WS(rs, 37)];
						  T5l = R1[WS(rs, 53)];
						  Tb1 = TaZ + Tb0;
						  TdB = Tb0 - TaZ;
						  T5c = T56 - T53;
						  T57 = T53 + T56;
						  Tb6 = T5i + T5j;
						  T5k = T5i - T5j;
						  T5m = R1[WS(rs, 21)];
					     }
					}
					{
					     E TdL, TdC, T5o, T5t, TdM, TdF;
					     Tb2 = TaY - Tb1;
					     Tc6 = TaY + Tb1;
					     {
						  E T58, T5d, Tb7, T5n, TdE, Tb8;
						  T58 = FMA(KP707106781, T57, T50);
						  T8e = FNMS(KP707106781, T57, T50);
						  T8f = FNMS(KP707106781, T5c, T5b);
						  T5d = FMA(KP707106781, T5c, T5b);
						  Tb7 = T5l + T5m;
						  T5n = T5l - T5m;
						  TdL = FNMS(KP414213562, TdA, TdB);
						  TdC = FMA(KP414213562, TdB, TdA);
						  T5e = FMA(KP198912367, T5d, T58);
						  T5H = FNMS(KP198912367, T58, T5d);
						  TdE = Tb7 - Tb6;
						  Tb8 = Tb6 + Tb7;
						  T5o = T5k + T5n;
						  T5t = T5n - T5k;
						  Tb9 = Tb5 - Tb8;
						  Tc7 = Tb5 + Tb8;
						  TdM = FMA(KP414213562, TdD, TdE);
						  TdF = FNMS(KP414213562, TdE, TdD);
					     }