hb_64.c

最新的FFT程序

				   T6c = FMA(KP414213562, T6b, T6a);
				   TbH = FMA(KP707106781, TbG, Tbz);
				   TeY = FNMS(KP707106781, TbG, Tbz);
				   TbN = TbL - TbM;
				   TgT = TbL + TbM;
				   TgX = TbV - TbY;
				   TbZ = TbV + TbY;
				   TeX = FNMS(KP707106781, TbN, TbK);
				   TbO = FMA(KP707106781, TbN, TbK);
				   {
					E Tja, TgR, Tj9, TgU;
					Tja = FNMS(KP707106781, TgQ, TgP);
					TgR = FMA(KP707106781, TgQ, TgP);
					Tj9 = FNMS(KP707106781, TgT, TgS);
					TgU = FMA(KP707106781, TgT, TgS);
					TeZ = FNMS(KP668178637, TeY, TeX);
					Tfq = FMA(KP668178637, TeX, TeY);
					TdG = FMA(KP198912367, TbH, TbO);
					TbP = FNMS(KP198912367, TbO, TbH);
					Tjb = FNMS(KP668178637, Tja, Tj9);
					TjC = FMA(KP668178637, Tj9, Tja);
					ThS = FMA(KP198912367, TgR, TgU);
					TgV = FNMS(KP198912367, TgU, TgR);
					Tc0 = FMA(KP707106781, TbZ, TbS);
					Tf1 = FNMS(KP707106781, TbZ, TbS);
				   }
				   Tc6 = Tc4 - Tc5;
				   Th0 = Tc4 + Tc5;
			      }
			 }
		    }
		    {
			 E Th7, Tjk, Thj, Thk, Tjh, Thi, Tji, The;
			 {
			      E TcE, Th5, T3b, T18, Tcd, Thg, T8b, T3L, Tcy, Thb, T3t, T1u, Tcv, Thc, T87;
			      E T3A, Tcj, TcF, T1f, T3E, TcG, Tcg, T8c, T3i, T1k, Tcm, T1j, Tcq, T3n, T1l;
			      E T3o, T3p;
			      {
				   E T1r, Tct, T1q, Tcx, T3w, T1s, T3x, T3y;
				   {
					E T15, Tcb, T14, TcD, T3H, T16, T3I, T3J;
					{
					     E T12, T13, T3F, T3G, Tf0, Tc7;
					     T12 = rio[WS(ios, 1)];
					     Tf0 = FNMS(KP707106781, Tc6, Tc3);
					     Tc7 = FMA(KP707106781, Tc6, Tc3);
					     {
						  E Tjd, TgY, Tjc, Th1;
						  Tjd = FNMS(KP707106781, TgX, TgW);
						  TgY = FMA(KP707106781, TgX, TgW);
						  Tjc = FNMS(KP707106781, Th0, TgZ);
						  Th1 = FMA(KP707106781, Th0, TgZ);
						  Tf2 = FMA(KP668178637, Tf1, Tf0);
						  Tfr = FNMS(KP668178637, Tf0, Tf1);
						  TdH = FNMS(KP198912367, Tc0, Tc7);
						  Tc8 = FMA(KP198912367, Tc7, Tc0);
						  Tje = FNMS(KP668178637, Tjd, Tjc);
						  TjD = FMA(KP668178637, Tjc, Tjd);
						  ThT = FMA(KP198912367, TgY, Th1);
						  Th2 = FNMS(KP198912367, Th1, TgY);
						  T13 = iio[-WS(ios, 33)];
					     }
					     T3F = iio[-WS(ios, 1)];
					     T3G = rio[WS(ios, 33)];
					     T15 = rio[WS(ios, 17)];
					     Tcb = T12 - T13;
					     T14 = T12 + T13;
					     TcD = T3F + T3G;
					     T3H = T3F - T3G;
					     T16 = iio[-WS(ios, 49)];
					     T3I = iio[-WS(ios, 17)];
					     T3J = rio[WS(ios, 49)];
					}
					{
					     E T1o, T1p, T3u, T3v;
					     T1o = iio[-WS(ios, 61)];
					     {
						  E TcC, T17, Tcc, T3K;
						  TcC = T15 - T16;
						  T17 = T15 + T16;
						  Tcc = T3I + T3J;
						  T3K = T3I - T3J;
						  TcE = TcC + TcD;
						  Th5 = TcD - TcC;
						  T3b = T14 - T17;
						  T18 = T14 + T17;
						  Tcd = Tcb - Tcc;
						  Thg = Tcb + Tcc;
						  T8b = T3H + T3K;
						  T3L = T3H - T3K;
						  T1p = rio[WS(ios, 29)];
					     }
					     T3u = iio[-WS(ios, 29)];
					     T3v = rio[WS(ios, 61)];
					     T1r = rio[WS(ios, 13)];
					     Tct = T1o - T1p;
					     T1q = T1o + T1p;
					     Tcx = T3v + T3u;
					     T3w = T3u - T3v;
					     T1s = iio[-WS(ios, 45)];
					     T3x = iio[-WS(ios, 13)];
					     T3y = rio[WS(ios, 45)];
					}
				   }
				   {
					E T1c, Tci, T1b, Tch, T3e, T1d, T3f, T3g;
					{
					     E T19, T1a, T3c, T3d;
					     T19 = rio[WS(ios, 9)];
					     {
						  E Tcw, T1t, Tcu, T3z;
						  Tcw = T1r - T1s;
						  T1t = T1r + T1s;
						  Tcu = T3x + T3y;
						  T3z = T3x - T3y;
						  Tcy = Tcw - Tcx;
						  Thb = Tcw + Tcx;
						  T3t = T1q - T1t;
						  T1u = T1q + T1t;
						  Tcv = Tct - Tcu;
						  Thc = Tct + Tcu;
						  T87 = T3w + T3z;
						  T3A = T3w - T3z;
						  T1a = iio[-WS(ios, 41)];
					     }
					     T3c = iio[-WS(ios, 9)];
					     T3d = rio[WS(ios, 41)];
					     T1c = iio[-WS(ios, 57)];
					     Tci = T19 - T1a;
					     T1b = T19 + T1a;
					     Tch = T3c + T3d;
					     T3e = T3c - T3d;
					     T1d = rio[WS(ios, 25)];
					     T3f = iio[-WS(ios, 25)];
					     T3g = rio[WS(ios, 57)];
					}
					{
					     E Tce, T1e, Tcf, T3h;
					     Tcj = Tch - Tci;
					     TcF = Tci + Tch;
					     Tce = T1c - T1d;
					     T1e = T1c + T1d;
					     Tcf = T3g + T3f;
					     T3h = T3f - T3g;
					     T1f = T1b + T1e;
					     T3E = T1e - T1b;
					     TcG = Tce + Tcf;
					     Tcg = Tce - Tcf;
					     T8c = T3e + T3h;
					     T3i = T3e - T3h;
					}
				   }
				   {
					E T1h, T1i, T3l, T3m;
					T1h = rio[WS(ios, 5)];
					T1i = iio[-WS(ios, 37)];
					T3l = iio[-WS(ios, 5)];
					T3m = rio[WS(ios, 37)];
					T1k = rio[WS(ios, 21)];
					Tcm = T1h - T1i;
					T1j = T1h + T1i;
					Tcq = T3l + T3m;
					T3n = T3l - T3m;
					T1l = iio[-WS(ios, 53)];
					T3o = iio[-WS(ios, 21)];
					T3p = rio[WS(ios, 53)];
				   }
			      }
			      {
				   E T85, Tcr, Th8, T3k, Tco, Th9, T3r, T9U, T9V, T88;
				   {
					E T1g, T86, T1v, T8a, T8d;
					T85 = T18 - T1f;
					T1g = T18 + T1f;
					{
					     E Tcp, T1m, Tcn, T3q, T1n;
					     Tcp = T1k - T1l;
					     T1m = T1k + T1l;
					     Tcn = T3o + T3p;
					     T3q = T3o - T3p;
					     Tcr = Tcp + Tcq;
					     Th8 = Tcq - Tcp;
					     T3k = T1j - T1m;
					     T1n = T1j + T1m;
					     Tco = Tcm - Tcn;
					     Th9 = Tcm + Tcn;
					     T86 = T3n + T3q;
					     T3r = T3n - T3q;
					     T1v = T1n + T1u;
					     T8a = T1u - T1n;
					}
					T8d = T8b - T8c;
					T9U = T8c + T8b;
					T9T = T1g - T1v;
					T1w = T1g + T1v;
					T9a = T8d - T8a;
					T8e = T8a + T8d;
					T9V = T86 + T87;
					T88 = T86 - T87;
				   }
				   {
					E T6f, T3j, T6i, T3M, T3N, T3O, T6j, T3C, T3s, T3B;
					T6f = T3b - T3i;
					T3j = T3b + T3i;
					TaF = T9V + T9U;
					T9W = T9U - T9V;
					T99 = T85 - T88;
					T89 = T85 + T88;
					T6i = T3L - T3E;
					T3M = T3E + T3L;
					T3N = T3r - T3k;
					T3s = T3k + T3r;
					T3B = T3t - T3A;
					T3O = T3t + T3A;
					T6j = T3s - T3B;
					T3C = T3s + T3B;
					{
					     E Th6, Tcl, Tf8, Thh, TcJ, TcK, Tf5, TcI, Tf6, TcA;
					     {
						  E TcH, T6g, T3P, Tck, Tcs, Tcz;
						  Th6 = Tcj + Tcg;
						  Tck = Tcg - Tcj;
						  T6k = FMA(KP707106781, T6j, T6i);
						  T7a = FNMS(KP707106781, T6j, T6i);
						  T5p = FNMS(KP707106781, T3C, T3j);
						  T3D = FMA(KP707106781, T3C, T3j);
						  T6g = T3O - T3N;
						  T3P = T3N + T3O;
						  Tcl = FMA(KP707106781, Tck, Tcd);
						  Tf8 = FNMS(KP707106781, Tck, Tcd);
						  T6h = FMA(KP707106781, T6g, T6f);
						  T79 = FNMS(KP707106781, T6g, T6f);
						  T5q = FNMS(KP707106781, T3P, T3M);
						  T3Q = FMA(KP707106781, T3P, T3M);
						  TcH = TcF - TcG;
						  Thh = TcF + TcG;
						  TcJ = FMA(KP414213562, Tco, Tcr);
						  Tcs = FNMS(KP414213562, Tcr, Tco);
						  Tcz = FMA(KP414213562, Tcy, Tcv);
						  TcK = FNMS(KP414213562, Tcv, Tcy);
						  Tf5 = FNMS(KP707106781, TcH, TcE);
						  TcI = FMA(KP707106781, TcH, TcE);
						  Tf6 = Tcz - Tcs;
						  TcA = Tcs + Tcz;
					     }
					     {
						  E Tf9, TcL, Tha, Thd;
						  Tf7 = FMA(KP923879532, Tf6, Tf5);
						  Tg0 = FNMS(KP923879532, Tf6, Tf5);
						  Teg = FNMS(KP923879532, TcA, Tcl);
						  TcB = FMA(KP923879532, TcA, Tcl);
						  Tf9 = TcJ - TcK;
						  TcL = TcJ + TcK;
						  Th7 = FMA(KP707106781, Th6, Th5);
						  Tjk = FNMS(KP707106781, Th6, Th5);
						  Tfa = FMA(KP923879532, Tf9, Tf8);
						  TfZ = FNMS(KP923879532, Tf9, Tf8);
						  Tef = FNMS(KP923879532, TcL, TcI);
						  TcM = FMA(KP923879532, TcL, TcI);
						  Thj = FMA(KP414213562, Th8, Th9);
						  Tha = FNMS(KP414213562, Th9, Th8);
						  Thd = FNMS(KP414213562, Thc, Thb);
						  Thk = FMA(KP414213562, Thb, Thc);
						  Tjh = FNMS(KP707106781, Thh, Thg);
						  Thi = FMA(KP707106781, Thh, Thg);
						  Tji = Tha + Thd;
						  The = Tha - Thd;
					     }
					}
				   }
			      }
			 }
			 {
			      E Tdh, Tho, T3S, T1D, TcQ, Thz, T8m, T4s, Tdb, Thu, T4a, T1Z, Td8, Thv, T8i;
			      E T4h, TcT, Tdi, T1K, T4l, Tdj, TcW, T3Z, T8n, TcZ, T1O, T45, Td3, T44, Td2;
			      E T1R, T46;
			      {
				   E Td6, T1V, T4e, Tda, T4d, Td9, T1Y, T4f;
				   {
					E TcO, T1z, T4p, Tdg, T4o, Tdf, T1C, T4q;
					{
					     E T4m, T4n, T1A, T1B;
					     {
						  E T1x, Tjl, Thl, T1y;
						  T1x = iio[-WS(ios, 63)];
						  Tjj = FNMS(KP923879532, Tji, Tjh);
						  Tkc = FMA(KP923879532, Tji, Tjh);
						  Tis = FNMS(KP923879532, The, Th7);
						  Thf = FMA(KP923879532, The, Th7);
						  Tjl = Thj - Thk;
						  Thl = Thj + Thk;
						  T1y = rio[WS(ios, 31)];
						  T4m = iio[-WS(ios, 31)];
						  Tjm = FMA(KP923879532, Tjl, Tjk);
						  Tkb = FNMS(KP923879532, Tjl, Tjk);
						  Tir = FNMS(KP923879532, Thl, Thi);
						  Thm = FMA(KP923879532, Thl, Thi);
						  TcO = T1x - T1y;
						  T1z = T1x + T1y;
						  T4n = rio[WS(ios, 63)];
					     }
					     T1A = rio[WS(ios, 15)];
					     T1B = iio[-WS(ios, 47)];
					     T4p = iio[-WS(ios, 15)];
					     Tdg = T4n + T4m;
					     T4o = T4m - T4n;
					     Tdf = T1A - T1B;
					     T1C = T1A + T1B;
					     T4q = rio[WS(ios, 47)];
					}
					{
					     E T4b, T4c, T1W, T1X;
					     {
						  E T1T, TcP, T4r, T1U;
						  T1T = iio[-WS(ios, 59)];
						  Tdh = Tdf - Tdg;
						  Tho = Tdf + Tdg;
						  T3S = T1z - T1C;
						  T1D = T1z + T1C;
						  TcP = T4p + T4q;
						  T4r = T4p - T4q;
						  T1U = rio[WS(ios, 27)];
						  T4b = iio[-WS(ios, 27)];
						  TcQ = TcO - TcP;
						  Thz = TcO + TcP;
						  T8m = T4o + T4r;
						  T4s = T4o - T4r;
						  Td6 = T1T - T1U;
						  T1V = T1T + T1U;
						  T4c = rio[WS(ios, 59)];
					     }
					     T1W = rio[WS(ios, 11)];
					     T1X = iio[-WS(ios, 43)];
					     T4e = iio[-WS(ios, 11)];
					     Tda = T4c + T4b;
					     T4d = T4b - T4c;
					     Td9 = T1W - T1X;
					     T1Y = T1W + T1X;
					     T4f = rio[WS(ios, 43)];
					}
				   }
				   {
					E T3V, TcU, T3Y, TcV;
					{
					     E TcR, T1G, T3W, TcS, T1J, T3X;
					     {
						  E T3T, T3U, T1H, T1I;
						  {
						       E T1E, Td7, T4g, T1F;
						       T1E = rio[WS(ios, 7)];
						       Tdb = Td9 - Tda;
						       Thu = Td9 + Tda;
						       T4a = T1V - T1Y;
						       T1Z = T1V + T1Y;
						       Td7 = T4e + T4f;
						       T4g = T4e - T4f;
						       T1F = iio[-WS(ios, 39)];
						       T3T = iio[-WS(ios, 7)];
						       Td8 = Td6 - Td7;
						       Thv = Td6 + Td7;
						       T8i = T4d + T4g;
						       T4h = T4d - T4g;
						       TcR = T1E - T1F;
						       T1G = T1E + T1F;
						       T3U = rio[WS(ios, 39)];
						  }
						  T1H = iio[-WS(ios, 55)];
						  T1I = rio[WS(ios, 23)];
						  T3W = iio[-WS(ios, 23)];
						  TcS = T3T + T3U;
						  T3V = T3T - T3U;
						  TcU = T1H - T1I;
						  T1J = T1H + T1I;
						  T3X = rio[WS(ios, 55)];
					     }
					     TcT = TcR - TcS;
					     Tdi = TcR + TcS;
					     T1K = T1G + T1J;
					     T4l = T1J - T1G;
					     T3Y = T3W - T3X;
					     TcV = T3X + T3W;
					}
					{
					     E T42, T43, T1M, T1N, T1P, T1Q;
					     T1M = rio[WS(ios, 3)];
					     T1N = iio[-WS(ios, 35)];
					     Tdj = TcU + TcV;
					     TcW = TcU - TcV;
					     T3Z = T3V - T3Y;
					     T8n = T3V + T3Y;
					     TcZ = T1M - T1N;
					     T1O = T1M + T1N;
					     T42 = iio[-WS(ios, 3)];
					     T43 = rio[WS(ios, 35)];
					     T1P = rio[WS(ios, 19)];
					     T1Q = iio[-WS(ios, 51)];
					     T45 = iio[-WS(ios, 19)];
					     Td3 = T42 + T43;
					     T44 = T42 - T43;
					     Td2 = T1P - T1Q;
					     T1R = T1P + T1Q;
					     T46 = rio[WS(ios, 51)];
					}
				   }
			      }
			      {
				   E T8g, Td4, Thr, T41, T9Z, Ths, Td1, T48, Ta0, T8j, T8h;
				   {
					E T1L, T1S, Td0, T47, T8o, T20, T8l;
					T8g = T1D - T1K;
					T1L = T1D + T1K;
					Td4 = Td2 + Td3;
					Thr = Td2 - Td3;
					T41 = T1O - T1R;
					T1S = T1O + T1R;
					Td0 = T45 + T46;
					T47 = T45 - T46;
					T8o = T8m - T8n;
					T9Z = T8n + T8m;
					T20 = T1S + T1Z;
					T8l = T1Z - T1S;
					Ths = TcZ + Td0;
					Td1 = TcZ - Td0;
					T48 = T44 - T47;
					T8h = T44 + T47;
					T21 = T1L + T20;
					T9Y = T1L - T20;
					T9d = T8o - T8l;
					T8p = T8l + T8o;
				   }
				   Ta0 = T8h + T8i;
				   T8j = T8h - T8i;
				   {
					E T6m, T40, T6p, T4t, T4u, T4v, T6q, T4j, T49, T4i;
					T6m = T3S - T3Z;
					T40 = T3S + T3Z;
					TaG = Ta0 + T9Z;
					Ta1 = T9Z - Ta0;
					T9c = T8g - T8j;
					T8k = T8g + T8j;
					T6p = T4s - T4l;
					T4t = T4l + T4s;
					T4u = T48 - T41;