hf_64.c

最新的FFT程序

		    }
		    {
			 E TjD, TjJ, Tgo, Tf2, Tjp, Tjv, Tha, TgI, Tgg, Tgd, Tjq, Tgr, Tjw, Tfp, Thk;
			 E Tho, Th8, Th4, Tgv, TgB, Tgl, TfR, TjE, Thd, TjK, TgP, Tgw, Tg8, Thh, Thn;
			 E Th7, TgX;
			 {
			      E TgK, TgJ, TgN, TgM, TfW, Th1, Thj, Th0, Th2, Tg7;
			      {
				   E TgE, TeQ, TjB, Tjn, TgF, TgG, TjC, Tf1, TeV, Tf0;
				   TgE = TeM - TeP;
				   TeQ = TeM + TeP;
				   TjB = Tjm - Tjl;
				   Tjn = Tjl + Tjm;
				   TgF = TeU - TeR;
				   TeV = TeR + TeU;
				   iio[-WS(ios, 20)] = FMA(KP923879532, Tje, Tjd);
				   rio[WS(ios, 52)] = FMS(KP923879532, Tje, Tjd);
				   Tf0 = TeW - TeZ;
				   TgG = TeW + TeZ;
				   TjC = Tf0 - TeV;
				   Tf1 = TeV + Tf0;
				   {
					E Tfi, Tgp, Tfd, Tfn;
					{
					     E Tf7, Tjo, TgH, Tfc;
					     TgK = Tf3 - Tf6;
					     Tf7 = Tf3 + Tf6;
					     TjD = FMA(KP707106781, TjC, TjB);
					     TjJ = FNMS(KP707106781, TjC, TjB);
					     Tgo = FNMS(KP707106781, Tf1, TeQ);
					     Tf2 = FMA(KP707106781, Tf1, TeQ);
					     Tjo = TgF + TgG;
					     TgH = TgF - TgG;
					     Tfc = Tfa + Tfb;
					     TgJ = Tfa - Tfb;
					     TgN = Tfe - Tfh;
					     Tfi = Tfe + Tfh;
					     Tjp = FMA(KP707106781, Tjo, Tjn);
					     Tjv = FNMS(KP707106781, Tjo, Tjn);
					     Tha = FNMS(KP707106781, TgH, TgE);
					     TgI = FMA(KP707106781, TgH, TgE);
					     Tgp = FNMS(KP414213562, Tf7, Tfc);
					     Tfd = FMA(KP414213562, Tfc, Tf7);
					     Tfn = Tfl + Tfm;
					     TgM = Tfl - Tfm;
					}
					{
					     E TgY, TgZ, Tgq, Tfo;
					     TfW = TfS + TfV;
					     TgY = TfS - TfV;
					     TgZ = Tgf - Tge;
					     Tgg = Tge + Tgf;
					     Tgd = Tgb + Tgc;
					     Th1 = Tgc - Tgb;
					     Tgq = FMA(KP414213562, Tfi, Tfn);
					     Tfo = FNMS(KP414213562, Tfn, Tfi);
					     Thj = FMA(KP707106781, TgZ, TgY);
					     Th0 = FNMS(KP707106781, TgZ, TgY);
					     Tjq = Tgp + Tgq;
					     Tgr = Tgp - Tgq;
					     Tjw = Tfo - Tfd;
					     Tfp = Tfd + Tfo;
					     Th2 = Tg6 - Tg1;
					     Tg7 = Tg1 + Tg6;
					}
				   }
			      }
			      {
				   E TgR, TgV, TgU, TgS, Thb, TgL;
				   {
					E TfM, Tgt, TfH, TfP, Tgu, TfQ;
					{
					     E Tfv, TfG, Thi, Th3;
					     TgR = Tfr - Tfu;
					     Tfv = Tfr + Tfu;
					     TfG = TfA + TfF;
					     TgV = TfA - TfF;
					     TgU = TfK - TfL;
					     TfM = TfK + TfL;
					     Thi = FMA(KP707106781, Th2, Th1);
					     Th3 = FNMS(KP707106781, Th2, Th1);
					     Tgt = FNMS(KP707106781, TfG, Tfv);
					     TfH = FMA(KP707106781, TfG, Tfv);
					     Thk = FMA(KP198912367, Thj, Thi);
					     Tho = FNMS(KP198912367, Thi, Thj);
					     Th8 = FNMS(KP668178637, Th0, Th3);
					     Th4 = FMA(KP668178637, Th3, Th0);
					     TfP = TfN + TfO;
					     TgS = TfO - TfN;
					}
					Tgu = FNMS(KP707106781, TfP, TfM);
					TfQ = FMA(KP707106781, TfP, TfM);
					Thb = FNMS(KP414213562, TgJ, TgK);
					TgL = FMA(KP414213562, TgK, TgJ);
					Tgv = FNMS(KP668178637, Tgu, Tgt);
					TgB = FMA(KP668178637, Tgt, Tgu);
					Tgl = FNMS(KP198912367, TfH, TfQ);
					TfR = FMA(KP198912367, TfQ, TfH);
				   }
				   {
					E Thg, TgT, Thc, TgO, Thf, TgW;
					Thc = FMA(KP414213562, TgM, TgN);
					TgO = FNMS(KP414213562, TgN, TgM);
					Thg = FMA(KP707106781, TgS, TgR);
					TgT = FNMS(KP707106781, TgS, TgR);
					TjE = Thc - Thb;
					Thd = Thb + Thc;
					TjK = TgL + TgO;
					TgP = TgL - TgO;
					Thf = FMA(KP707106781, TgV, TgU);
					TgW = FNMS(KP707106781, TgV, TgU);
					Tgw = FNMS(KP707106781, Tg7, TfW);
					Tg8 = FMA(KP707106781, Tg7, TfW);
					Thh = FMA(KP198912367, Thg, Thf);
					Thn = FNMS(KP198912367, Thf, Thg);
					Th7 = FNMS(KP668178637, TgT, TgW);
					TgX = FMA(KP668178637, TgW, TgT);
				   }
			      }
			 }
			 {
			      E Tjy, Tjx, TjI, TjH;
			      {
				   E Tgk, Tfq, Tjr, Tjt, Tgx, Tgh;
				   Tgk = FNMS(KP923879532, Tfp, Tf2);
				   Tfq = FMA(KP923879532, Tfp, Tf2);
				   Tjr = FMA(KP923879532, Tjq, Tjp);
				   Tjt = FNMS(KP923879532, Tjq, Tjp);
				   Tgx = FNMS(KP707106781, Tgg, Tgd);
				   Tgh = FMA(KP707106781, Tgg, Tgd);
				   {
					E TgA, Tjz, TjA, TgD;
					{
					     E Tgs, TgC, Tgm, Tgi, Tgz, Tgy;
					     TgA = FMA(KP923879532, Tgr, Tgo);
					     Tgs = FNMS(KP923879532, Tgr, Tgo);
					     Tgy = FMA(KP668178637, Tgx, Tgw);
					     TgC = FNMS(KP668178637, Tgw, Tgx);
					     Tgm = FMA(KP198912367, Tg8, Tgh);
					     Tgi = FNMS(KP198912367, Tgh, Tg8);
					     Tgz = Tgv + Tgy;
					     Tjy = Tgy - Tgv;
					     Tjx = FMA(KP923879532, Tjw, Tjv);
					     Tjz = FNMS(KP923879532, Tjw, Tjv);
					     {
						  E Tgn, Tjs, Tju, Tgj;
						  Tgn = Tgl - Tgm;
						  Tjs = Tgl + Tgm;
						  Tju = Tgi - TfR;
						  Tgj = TfR + Tgi;
						  iio[-WS(ios, 58)] = FMA(KP831469612, Tgz, Tgs);
						  rio[WS(ios, 26)] = FNMS(KP831469612, Tgz, Tgs);
						  rio[WS(ios, 18)] = FMA(KP980785280, Tgn, Tgk);
						  iio[-WS(ios, 50)] = FNMS(KP980785280, Tgn, Tgk);
						  iio[-WS(ios, 2)] = FMA(KP980785280, Tjs, Tjr);
						  rio[WS(ios, 34)] = FMS(KP980785280, Tjs, Tjr);
						  iio[-WS(ios, 18)] = FMA(KP980785280, Tju, Tjt);
						  rio[WS(ios, 50)] = FMS(KP980785280, Tju, Tjt);
						  rio[WS(ios, 2)] = FMA(KP980785280, Tgj, Tfq);
						  iio[-WS(ios, 34)] = FNMS(KP980785280, Tgj, Tfq);
						  TjA = TgB + TgC;
						  TgD = TgB - TgC;
					     }
					}
					iio[-WS(ios, 26)] = FNMS(KP831469612, TjA, Tjz);
					rio[WS(ios, 58)] = -(FMA(KP831469612, TjA, Tjz));
					rio[WS(ios, 10)] = FMA(KP831469612, TgD, TgA);
					iio[-WS(ios, 42)] = FNMS(KP831469612, TgD, TgA);
				   }
			      }
			      {
				   E Th6, TjF, TjG, Th9, TgQ, Th5;
				   Th6 = FNMS(KP923879532, TgP, TgI);
				   TgQ = FMA(KP923879532, TgP, TgI);
				   Th5 = TgX + Th4;
				   TjI = Th4 - TgX;
				   TjH = FNMS(KP923879532, TjE, TjD);
				   TjF = FMA(KP923879532, TjE, TjD);
				   iio[-WS(ios, 10)] = FMA(KP831469612, Tjy, Tjx);
				   rio[WS(ios, 42)] = FMS(KP831469612, Tjy, Tjx);
				   rio[WS(ios, 6)] = FMA(KP831469612, Th5, TgQ);
				   iio[-WS(ios, 38)] = FNMS(KP831469612, Th5, TgQ);
				   TjG = Th7 - Th8;
				   Th9 = Th7 + Th8;
				   iio[-WS(ios, 6)] = FMA(KP831469612, TjG, TjF);
				   rio[WS(ios, 38)] = FMS(KP831469612, TjG, TjF);
				   rio[WS(ios, 22)] = FMA(KP831469612, Th9, Th6);
				   iio[-WS(ios, 54)] = FNMS(KP831469612, Th9, Th6);
			      }
			      {
				   E Thm, TjL, TjM, Thp, The, Thl;
				   Thm = FMA(KP923879532, Thd, Tha);
				   The = FNMS(KP923879532, Thd, Tha);
				   Thl = Thh + Thk;
				   TjO = Thk - Thh;
				   TjN = FMA(KP923879532, TjK, TjJ);
				   TjL = FNMS(KP923879532, TjK, TjJ);
				   iio[-WS(ios, 22)] = FMA(KP831469612, TjI, TjH);
				   rio[WS(ios, 54)] = FMS(KP831469612, TjI, TjH);
				   rio[WS(ios, 14)] = FMA(KP980785280, Thl, The);
				   iio[-WS(ios, 46)] = FNMS(KP980785280, Thl, The);
				   TjM = Tho - Thn;
				   Thp = Thn + Tho;
				   iio[-WS(ios, 14)] = FMA(KP980785280, TjM, TjL);
				   rio[WS(ios, 46)] = FMS(KP980785280, TjM, TjL);
				   iio[-WS(ios, 62)] = FMA(KP980785280, Thp, Thm);
				   rio[WS(ios, 30)] = FNMS(KP980785280, Thp, Thm);
			      }
			 }
		    }
		    {
			 E TjS, TcD, Tcw, TkO, TkN, Tcz;
			 {
			      E Tkq, TbB, Tkw, T99, TbF, TbL, Tbv, Taj, Tcu, Tcy, Tci, Tce, Tcr, Tcx, Tch;
			      E Tc7, TkE, Tcn, TkK, TbZ, TbP, T7J, TbO, T7u, TkB, Tkn, TbI, TbM, Tbw, Tbs;
			      E T7Y, TbQ;
			      {
				   E TbU, TbT, TbX, TbW, Tc1, Tc5, Tc4, Tc2, TaG, Tbq, Tbn, Tcb, Tct, Tca, Tcc;
				   E Tbb, Tcl, TbV;
				   {
					E T8W, Tbz, T8z, T97, T8n, T8y;
					TbU = FMA(KP707106781, T8m, T87);
					T8n = FNMS(KP707106781, T8m, T87);
					T8y = FNMS(KP707106781, T8x, T8u);
					TbT = FMA(KP707106781, T8x, T8u);
					TbX = FMA(KP707106781, T8V, T8G);
					T8W = FNMS(KP707106781, T8V, T8G);
					iio[-WS(ios, 30)] = FMA(KP980785280, TjO, TjN);
					rio[WS(ios, 62)] = FMS(KP980785280, TjO, TjN);
					Tbz = FNMS(KP668178637, T8n, T8y);
					T8z = FMA(KP668178637, T8y, T8n);
					T97 = FNMS(KP707106781, T96, T93);
					TbW = FMA(KP707106781, T96, T93);
					{
					     E Tae, TbD, Ta3, Tah;
					     {
						  E T9x, Ta2, TbA, T98;
						  Tc1 = FMA(KP707106781, T9w, T9h);
						  T9x = FNMS(KP707106781, T9w, T9h);
						  Ta2 = T9M - Ta1;
						  Tc5 = T9M + Ta1;
						  Tc4 = FMA(KP707106781, Tad, Taa);
						  Tae = FNMS(KP707106781, Tad, Taa);
						  TbA = FMA(KP668178637, T8W, T97);
						  T98 = FNMS(KP668178637, T97, T8W);
						  TbD = FNMS(KP923879532, Ta2, T9x);
						  Ta3 = FMA(KP923879532, Ta2, T9x);
						  Tkq = Tbz + TbA;
						  TbB = Tbz - TbA;
						  Tkw = T98 - T8z;
						  T99 = T8z + T98;
						  Tah = Taf - Tag;
						  Tc2 = Tag + Taf;
					     }
					     {
						  E Tc8, Tc9, TbE, Tai;
						  TaG = FNMS(KP707106781, TaF, Taq);
						  Tc8 = FMA(KP707106781, TaF, Taq);
						  Tc9 = Tbp + Tbo;
						  Tbq = Tbo - Tbp;
						  Tbn = FNMS(KP707106781, Tbm, Tbj);
						  Tcb = FMA(KP707106781, Tbm, Tbj);
						  TbE = FNMS(KP923879532, Tah, Tae);
						  Tai = FMA(KP923879532, Tah, Tae);
						  Tct = FMA(KP923879532, Tc9, Tc8);
						  Tca = FNMS(KP923879532, Tc9, Tc8);
						  TbF = FNMS(KP534511135, TbE, TbD);
						  TbL = FMA(KP534511135, TbD, TbE);
						  Tbv = FNMS(KP303346683, Ta3, Tai);
						  Taj = FMA(KP303346683, Tai, Ta3);
						  Tcc = Tba + TaV;
						  Tbb = TaV - Tba;
					     }
					}
				   }
				   {
					E Tcq, Tc3, Tcs, Tcd, Tcp, Tc6;
					Tcs = FMA(KP923879532, Tcc, Tcb);
					Tcd = FNMS(KP923879532, Tcc, Tcb);
					Tcq = FMA(KP923879532, Tc2, Tc1);
					Tc3 = FNMS(KP923879532, Tc2, Tc1);
					Tcu = FMA(KP098491403, Tct, Tcs);
					Tcy = FNMS(KP098491403, Tcs, Tct);
					Tci = FNMS(KP820678790, Tca, Tcd);
					Tce = FMA(KP820678790, Tcd, Tca);
					Tcp = FMA(KP923879532, Tc5, Tc4);
					Tc6 = FNMS(KP923879532, Tc5, Tc4);
					Tcl = FNMS(KP198912367, TbT, TbU);
					TbV = FMA(KP198912367, TbU, TbT);
					Tcr = FMA(KP098491403, Tcq, Tcp);
					Tcx = FNMS(KP098491403, Tcp, Tcq);
					Tch = FNMS(KP820678790, Tc3, Tc6);
					Tc7 = FMA(KP820678790, Tc6, Tc3);
				   }
				   {
					E TbG, Tbc, Tcm, TbY;
					Tcm = FMA(KP198912367, TbW, TbX);
					TbY = FNMS(KP198912367, TbX, TbW);
					TbG = FNMS(KP923879532, Tbb, TaG);
					Tbc = FMA(KP923879532, Tbb, TaG);
					TkE = Tcm - Tcl;
					Tcn = Tcl + Tcm;
					TkK = TbV + TbY;
					TbZ = TbV - TbY;
					{
					     E T7t, Tkm, TbH, Tbr;
					     TjS = T7l + T7s;
					     T7t = T7l - T7s;
					     Tkm = TcC - TcB;
					     TcD = TcB + TcC;
					     TbP = FNMS(KP414213562, T7B, T7I);
					     T7J = FMA(KP414213562, T7I, T7B);
					     TbH = FMA(KP923879532, Tbq, Tbn);
					     Tbr = FNMS(KP923879532, Tbq, Tbn);
					     TbO = FNMS(KP707106781, T7t, T7e);
					     T7u = FMA(KP707106781, T7t, T7e);
					     TkB = FNMS(KP707106781, Tkm, Tkl);
					     Tkn = FMA(KP707106781, Tkm, Tkl);
					     TbI = FNMS(KP534511135, TbH, TbG);
					     TbM = FMA(KP534511135, TbG, TbH);
					     Tbw = FNMS(KP303346683, Tbc, Tbr);
					     Tbs = FMA(KP303346683, Tbr, Tbc);
					     T7Y = FNMS(KP414213562, T7X, T7Q);
					     TbQ = FMA(KP414213562, T7Q, T7X);
					}
				   }
			      }
			      {
				   E TkJ, TkD, Tck, TbS, TbK, Tky, Tkx, TbN;
				   {
					E Tku, Tby, Tkv, Tbu, Tkt, Tbx;
					{
					     E Tbt, T9a, Tkr, Tks, Tkp;
					     Tku = Tbs - Taj;
					     Tbt = Taj + Tbs;
					     {
						  E TkC, T7Z, Tko, TbR, T80;
						  TkC = T7J + T7Y;
						  T7Z = T7J - T7Y;
						  Tko = TbQ - TbP;
						  TbR = TbP + TbQ;
						  TkJ = FMA(KP923879532, TkC, TkB);
						  TkD = FNMS(KP923879532, TkC, TkB);
						  Tby = FNMS(KP923879532, T7Z, T7u);
						  T80 = FMA(KP923879532, T7Z, T7u);
						  Tkv = FNMS(KP923879532, Tko, Tkn);
						  Tkp = FMA(KP923879532, Tko, Tkn);
						  Tck = FMA(KP923879532, TbR, TbO);
						  TbS = FNMS(KP923879532, TbR, TbO);
						  T9a = FMA(KP831469612, T99, T80);
						  Tbu = FNMS(KP831469612, T99, T80);
					     }
					     Tkt = FNMS(KP831469612, Tkq, Tkp);
					     Tkr = FMA(KP831469612, Tkq, Tkp);
					     Tks = Tbv - Tbw;
					     Tbx = Tbv + Tbw;
					     rio[WS(ios, 3)] = FMA(KP956940335, Tbt, T9a);
					     iio[-WS(ios, 35)] = FNMS(KP956940335, Tbt, T9a);
					     iio[-WS(ios, 3)] = FMA(KP956940335, Tks, Tkr);
					     rio[WS(ios, 35)] = FMS(KP956940335, Tks, Tkr);
					}
					{
					     E TbC, TbJ, Tkz, TkA;
					     TbK = FMA(KP831469612, TbB, Tby);
					     TbC = FNMS(KP831469612, TbB, Tby);
					     rio[WS(ios, 19)] = FMA(KP956940335, Tbx, Tbu);
					     iio[-WS(ios, 51)] = FNMS(KP956940335, Tbx, Tbu);
					     iio[-WS(ios, 19)] = FMA(KP956940335, Tku, Tkt);
					     rio[WS(ios, 51)] = FMS(KP956940335, Tku, Tkt);
					     TbJ = TbF + TbI;
					     Tky = TbI - TbF;
					     Tkx = FMA(KP831469612, Tkw, Tkv);
					     Tkz = FNMS(KP831469612, Tkw, Tkv);
					     TkA = TbM - TbL;
					     TbN = TbL + TbM;
					     iio[-WS(ios, 59)] = FMA(KP881921264, TbJ, TbC);
					     rio[WS(ios, 27)] = FNMS(KP881921264, TbJ, TbC);
					     iio[-WS(ios, 27)] = FMA(KP881921264, TkA, Tkz);
					     rio[WS(ios, 59)] = FMS(KP881921264, TkA, Tkz);
					}
				   }
				   {
					E Tcg, TkI, TkH, Tcj;
					{
					     E Tc0, Tcf, TkF, TkG;
					     Tcg = FNMS(KP980785280, TbZ, TbS);
					     Tc0 = FMA(KP980785280, TbZ, TbS);
					     rio[WS(ios, 11)] = FMA(KP881921264, TbN, TbK);
					     iio[-WS(ios, 43)] = FNMS(KP881921264, TbN, TbK);
					     iio[-WS(ios, 11)] = FMA(KP881921264, Tky, Tkx);
					     rio[WS(ios, 43)] = FMS(KP881921264, Tky, Tkx);
					     Tcf = Tc7 + Tce;
					     TkI = Tce - Tc7;
					     TkH = FNMS(KP980785280, TkE, TkD);
					     TkF = FMA(KP980785280, TkE, TkD);
					     TkG = Tch - Tci;
					     Tcj = Tch + Tci;
					     rio[WS(ios, 7)] = FMA(KP773010453, Tcf, Tc0);
					     iio[-WS(ios, 39)] = FNMS(KP773010453, Tcf, Tc0);
					     iio[-WS(ios, 7)] = FMA(KP773010453, TkG, TkF);
					     rio[WS(ios, 39)] = FMS(KP773010453, TkG, TkF);
					}
					{
					     E Tco, Tcv, TkL, TkM;
					     Tcw = FMA(KP980785280, Tcn, Tck);
					     Tco = FNMS(KP980785280, Tcn, Tck);
					     rio[WS(ios, 23)] = FMA(KP773010453, Tcj, Tcg);
					     iio[-WS(ios, 55)] = FNMS(KP773010453, Tcj, Tcg);
					     iio[-WS(ios, 23)] = FMA(KP773010453, TkI, TkH);
					     rio[WS(ios, 55)] = FMS(KP773010453, TkI, TkH);
					     Tcv = Tcr + Tcu;
					     TkO = Tcu - Tcr;
					     TkN = FMA(KP980785280, TkK, TkJ);
					     TkL = FNMS(KP980785280, TkK, TkJ);
					     TkM = Tcy - Tcx;
					     Tcz = Tcx + Tcy;
					     rio[WS(ios, 15)] = FMA(KP995184726, Tcv, Tco);
					     iio[-WS(ios, 47)] = FNMS(KP995184726, Tcv, Tco);
					     iio[-WS(ios, 15)] = FMA(KP995184726, TkM, TkL);
					     rio[WS(ios, 47)] = FMS(KP995184726, TkM, TkL);
					}
				   }
			      }
			 }
			 {
			      E TjW, TdN, Tk2, Td1, TdR, TdX, TdH, Tdl, TeG, TeK, Teu, Teq, TeD, TeJ, Tet;
			      E Tej, Tka, Tez, Tkg, Teb, Te1, TcH, Te0, TcE, Tk7, TjT, TdU, TdY, TdI, TdE;
			      E TcK, Te2;
			      {
				   E Te6, Te5, Te9, Te8, Ted, Teh, Teg, Tee, Tdo, TdC, Tdz, Ten, TeF, Tem, Teo;
				   E Tdv, Tex, Te7;
				   {
					E TcP, TcS, TcW, TcZ;
					Te6 = FNMS(KP707106781, TcO, TcN);
					TcP = FMA(KP707106781, TcO, TcN);
					iio[-WS(ios, 63)] = FMA(KP995184726, Tcz, Tcw);
					rio[WS(ios, 31)] = FNMS(KP995184726, Tcz, Tcw);
					iio[-WS(ios, 31)] = FMA(KP995184726, TkO, TkN);
					rio[WS(ios, 63)] = FMS(KP995184726, TkO, TkN);
					TcS = FMA(KP707106781, TcR, TcQ);
					Te5 = FNMS(KP707106781, TcR, TcQ);
					Te9 = FNMS(KP707106781, TcV, TcU);
					TcW = FMA(KP707106781, TcV, TcU);
					TcZ = FMA(KP707106781, TcY, TcX);
					Te8 = FNMS(KP707106781, TcY, TcX);
					{
					     E Tdg, TdP, Tdd, Tdj;
					     {
						  E Td5, TdL, TcT, TdM, Td0, Tdc;
						  Ted = FNMS(KP707106781, Td4, Td3);
						  Td5 = FMA(KP707106781, Td4, Td3);
						  TdL = FNMS(KP198912367, TcP, TcS);
						  TcT = FMA(KP198912367, TcS, TcP);