t1_64.c

快速fft变换

					ri[WS(rs, 20)] = FMA(KP923879532, Tin, Tik);
					ri[WS(rs, 52)] = FNMS(KP923879532, Tin, Tik);
				   }
			      }
			 }
		    }
		    {
			 E TjD, TjJ, Tgo, Tf2, Tjp, Tjv, Tha, TgI, Tgd, Tgr, Tjw, Tjq, Tfp, Tgg, Thk;
			 E Tho, Th8, Th4, Tgv, TgB, Tgl, TfR, TjE, Thd, TjK, TgP, Tgx, Tg8, Thh, Thn;
			 E Th7, TgX;
			 {
			      E TgJ, TgK, TgM, TgN, Tg7, TfW, Th1, Thj, Th0, Th2;
			      {
				   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;
				   ii[WS(rs, 52)] = FNMS(KP923879532, Tje, Tjd);
				   ii[WS(rs, 20)] = FMA(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;
					     TgJ = Tf5 - Tf6;
					     Tf7 = Tf5 + Tf6;
					     TjD = FMA(KP707106781, TjC, TjB);
					     TjJ = FNMS(KP707106781, TjC, TjB);
					     Tgo = FMA(KP707106781, Tf1, TeQ);
					     Tf2 = FNMS(KP707106781, Tf1, TeQ);
					     Tjo = TgF + TgG;
					     TgH = TgF - TgG;
					     Tfc = Tf8 + Tfb;
					     TgK = Tf8 - Tfb;
					     TgM = Tfg - Tfh;
					     Tfi = Tfg + Tfh;
					     Tjp = FMA(KP707106781, Tjo, Tjn);
					     Tjv = FNMS(KP707106781, Tjo, Tjn);
					     Tha = FNMS(KP707106781, TgH, TgE);
					     TgI = FMA(KP707106781, TgH, TgE);
					     Tgp = FMA(KP414213562, Tf7, Tfc);
					     Tfd = FNMS(KP414213562, Tfc, Tf7);
					     Tfn = Tfj + Tfm;
					     TgN = Tfj - Tfm;
					}
					{
					     E TgY, TgZ, Tgq, Tfo;
					     Tgd = Tg9 + Tgc;
					     TgY = Tg9 - Tgc;
					     TgZ = Tg6 - Tg1;
					     Tg7 = Tg1 + Tg6;
					     TfW = TfU + TfV;
					     Th1 = TfU - TfV;
					     Tgq = FNMS(KP414213562, Tfi, Tfn);
					     Tfo = FMA(KP414213562, Tfn, Tfi);
					     Thj = FMA(KP707106781, TgZ, TgY);
					     Th0 = FNMS(KP707106781, TgZ, TgY);
					     Tgr = Tgp + Tgq;
					     Tjw = Tgq - Tgp;
					     Tjq = Tfd + Tfo;
					     Tfp = Tfd - Tfo;
					     Th2 = Tge - Tgf;
					     Tgg = Tge + Tgf;
					}
				   }
			      }
			      {
				   E TgU, TgS, TgR, TgV, Thb, TgL;
				   {
					E TfM, Tgu, TfH, TfP, Tgt, TfQ;
					{
					     E Tfv, TfG, Thi, Th3;
					     TgU = Tft - Tfu;
					     Tfv = Tft + Tfu;
					     TfG = TfA + TfF;
					     TgS = TfF - TfA;
					     TgR = TfI - TfL;
					     TfM = TfI + TfL;
					     Thi = FMA(KP707106781, Th2, Th1);
					     Th3 = FNMS(KP707106781, Th2, Th1);
					     Tgu = FMA(KP707106781, TfG, Tfv);
					     TfH = FNMS(KP707106781, TfG, Tfv);
					     Thk = FNMS(KP198912367, Thj, Thi);
					     Tho = FMA(KP198912367, Thi, Thj);
					     Th8 = FMA(KP668178637, Th0, Th3);
					     Th4 = FNMS(KP668178637, Th3, Th0);
					     TfP = TfN + TfO;
					     TgV = TfN - TfO;
					}
					Tgt = FMA(KP707106781, TfP, TfM);
					TfQ = FNMS(KP707106781, TfP, TfM);
					Thb = FNMS(KP414213562, TgJ, TgK);
					TgL = FMA(KP414213562, TgK, TgJ);
					Tgv = FMA(KP198912367, Tgu, Tgt);
					TgB = FNMS(KP198912367, Tgt, Tgu);
					Tgl = FNMS(KP668178637, TfH, TfQ);
					TfR = FMA(KP668178637, 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);
					Tgx = FMA(KP707106781, Tg7, TfW);
					Tg8 = FNMS(KP707106781, Tg7, TfW);
					Thh = FMA(KP198912367, Thg, Thf);
					Thn = FNMS(KP198912367, Thf, Thg);
					Th7 = FNMS(KP668178637, TgT, TgW);
					TgX = FMA(KP668178637, TgW, TgT);
				   }
			      }
			 }
			 {
			      E Tju, Tjt, TjI, TjH;
			      {
				   E Tgk, Tfq, Tjx, Tjz, Tgw, Tgh;
				   Tgk = FNMS(KP923879532, Tfp, Tf2);
				   Tfq = FMA(KP923879532, Tfp, Tf2);
				   Tjx = FMA(KP923879532, Tjw, Tjv);
				   Tjz = FNMS(KP923879532, Tjw, Tjv);
				   Tgw = FMA(KP707106781, Tgg, Tgd);
				   Tgh = FNMS(KP707106781, Tgg, Tgd);
				   {
					E TgA, Tjr, Tjs, TgD;
					{
					     E Tgs, TgC, Tgm, Tgi, Tgz, Tgy;
					     TgA = FNMS(KP923879532, Tgr, Tgo);
					     Tgs = FMA(KP923879532, Tgr, Tgo);
					     Tgy = FNMS(KP198912367, Tgx, Tgw);
					     TgC = FMA(KP198912367, Tgw, Tgx);
					     Tgm = FMA(KP668178637, Tg8, Tgh);
					     Tgi = FNMS(KP668178637, Tgh, Tg8);
					     Tgz = Tgv + Tgy;
					     Tju = Tgy - Tgv;
					     Tjt = FNMS(KP923879532, Tjq, Tjp);
					     Tjr = FMA(KP923879532, Tjq, Tjp);
					     {
						  E Tgn, Tjy, TjA, Tgj;
						  Tgn = Tgl + Tgm;
						  Tjy = Tgm - Tgl;
						  TjA = TfR + Tgi;
						  Tgj = TfR - Tgi;
						  ri[WS(rs, 2)] = FMA(KP980785280, Tgz, Tgs);
						  ri[WS(rs, 34)] = FNMS(KP980785280, Tgz, Tgs);
						  ri[WS(rs, 58)] = FMA(KP831469612, Tgn, Tgk);
						  ri[WS(rs, 26)] = FNMS(KP831469612, Tgn, Tgk);
						  ii[WS(rs, 42)] = FNMS(KP831469612, Tjy, Tjx);
						  ii[WS(rs, 10)] = FMA(KP831469612, Tjy, Tjx);
						  ii[WS(rs, 58)] = FMA(KP831469612, TjA, Tjz);
						  ii[WS(rs, 26)] = FNMS(KP831469612, TjA, Tjz);
						  ri[WS(rs, 10)] = FMA(KP831469612, Tgj, Tfq);
						  ri[WS(rs, 42)] = FNMS(KP831469612, Tgj, Tfq);
						  Tjs = TgB + TgC;
						  TgD = TgB - TgC;
					     }
					}
					ii[WS(rs, 34)] = FNMS(KP980785280, Tjs, Tjr);
					ii[WS(rs, 2)] = FMA(KP980785280, Tjs, Tjr);
					ri[WS(rs, 18)] = FMA(KP980785280, TgD, TgA);
					ri[WS(rs, 50)] = FNMS(KP980785280, 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);
				   ii[WS(rs, 50)] = FNMS(KP980785280, Tju, Tjt);
				   ii[WS(rs, 18)] = FMA(KP980785280, Tju, Tjt);
				   ri[WS(rs, 6)] = FMA(KP831469612, Th5, TgQ);
				   ri[WS(rs, 38)] = FNMS(KP831469612, Th5, TgQ);
				   TjG = Th7 + Th8;
				   Th9 = Th7 - Th8;
				   ii[WS(rs, 38)] = FNMS(KP831469612, TjG, TjF);
				   ii[WS(rs, 6)] = FMA(KP831469612, TjG, TjF);
				   ri[WS(rs, 22)] = FMA(KP831469612, Th9, Th6);
				   ri[WS(rs, 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 = Thh + Thk;
				   TjN = FMA(KP923879532, TjK, TjJ);
				   TjL = FNMS(KP923879532, TjK, TjJ);
				   ii[WS(rs, 54)] = FNMS(KP831469612, TjI, TjH);
				   ii[WS(rs, 22)] = FMA(KP831469612, TjI, TjH);
				   ri[WS(rs, 14)] = FMA(KP980785280, Thl, The);
				   ri[WS(rs, 46)] = FNMS(KP980785280, Thl, The);
				   TjM = Tho - Thn;
				   Thp = Thn + Tho;
				   ii[WS(rs, 46)] = FNMS(KP980785280, TjM, TjL);
				   ii[WS(rs, 14)] = FMA(KP980785280, TjM, TjL);
				   ri[WS(rs, 62)] = FMA(KP980785280, Thp, Thm);
				   ri[WS(rs, 30)] = FNMS(KP980785280, Thp, Thm);
			      }
			 }
		    }
		    {
			 E TjS, TcD, Tcw, TkO, TkN, Tcz;
			 {
			      E TbB, Tkw, Tkq, 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 TbT, TbU, TbW, TbX, Tc4, Tc2, Tc1, Tc5, Tbn, Tbb, TaG, Tcb, Tct, Tca, Tcc;
				   E Tbq, Tcl, TbV;
				   {
					E T8W, Tbz, T8z, T97, T8n, T8y;
					TbT = FMA(KP707106781, T8m, T87);
					T8n = FNMS(KP707106781, T8m, T87);
					T8y = FNMS(KP707106781, T8x, T8u);
					TbU = FMA(KP707106781, T8x, T8u);
					TbW = FMA(KP707106781, T8V, T8G);
					T8W = FNMS(KP707106781, T8V, T8G);
					ii[WS(rs, 62)] = FMA(KP980785280, TjO, TjN);
					ii[WS(rs, 30)] = FNMS(KP980785280, TjO, TjN);
					Tbz = FMA(KP668178637, T8n, T8y);
					T8z = FNMS(KP668178637, T8y, T8n);
					T97 = FNMS(KP707106781, T96, T93);
					TbX = FMA(KP707106781, T96, T93);
					{
					     E Tae, TbE, Ta3, Tah;
					     {
						  E T9x, Ta2, TbA, T98;
						  Tc4 = FMA(KP707106781, T9w, T9h);
						  T9x = FNMS(KP707106781, T9w, T9h);
						  Ta2 = T9M - Ta1;
						  Tc2 = Ta1 + T9M;
						  Tc1 = FMA(KP707106781, Tad, Taa);
						  Tae = FNMS(KP707106781, Tad, Taa);
						  TbA = FNMS(KP668178637, T8W, T97);
						  T98 = FMA(KP668178637, T97, T8W);
						  TbE = FMA(KP923879532, Ta2, T9x);
						  Ta3 = FNMS(KP923879532, Ta2, T9x);
						  TbB = Tbz + TbA;
						  Tkw = TbA - Tbz;
						  Tkq = T8z + T98;
						  T99 = T8z - T98;
						  Tah = Taf - Tag;
						  Tc5 = Taf + Tag;
					     }
					     {
						  E Tc8, Tc9, TbD, Tai;
						  Tbn = FNMS(KP707106781, Tbm, Tbj);
						  Tc8 = FMA(KP707106781, Tbm, Tbj);
						  Tc9 = Tba + TaV;
						  Tbb = TaV - Tba;
						  TaG = FNMS(KP707106781, TaF, Taq);
						  Tcb = FMA(KP707106781, TaF, Taq);
						  TbD = FMA(KP923879532, Tah, Tae);
						  Tai = FNMS(KP923879532, Tah, Tae);
						  Tct = FMA(KP923879532, Tc9, Tc8);
						  Tca = FNMS(KP923879532, Tc9, Tc8);
						  TbF = FMA(KP303346683, TbE, TbD);
						  TbL = FNMS(KP303346683, TbD, TbE);
						  Tbv = FNMS(KP534511135, Ta3, Tai);
						  Taj = FMA(KP534511135, Tai, Ta3);
						  Tcc = Tbo + Tbp;
						  Tbq = Tbo - Tbp;
					     }
					}
				   }
				   {
					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 = FNMS(KP098491403, Tct, Tcs);
					Tcy = FMA(KP098491403, Tcs, Tct);
					Tci = FMA(KP820678790, Tca, Tcd);
					Tce = FNMS(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 TbH, Tbc, Tcm, TbY;
					Tcm = FMA(KP198912367, TbW, TbX);
					TbY = FNMS(KP198912367, TbX, TbW);
					TbH = FMA(KP923879532, Tbb, TaG);
					Tbc = FNMS(KP923879532, Tbb, TaG);
					TkE = Tcm - Tcl;
					Tcn = Tcl + Tcm;
					TkK = TbV + TbY;
					TbZ = TbV - TbY;
					{
					     E T7t, Tkm, TbG, Tbr;
					     TjS = T7l + T7s;
					     T7t = T7l - T7s;
					     Tkm = TcC - TcB;
					     TcD = TcB + TcC;
					     TbP = FNMS(KP414213562, T7B, T7I);
					     T7J = FMA(KP414213562, T7I, T7B);
					     TbG = 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(KP303346683, TbH, TbG);
					     TbM = FMA(KP303346683, TbG, TbH);
					     Tbw = FMA(KP534511135, Tbc, Tbr);
					     Tbs = FNMS(KP534511135, Tbr, Tbc);
					     T7Y = FNMS(KP414213562, T7X, T7Q);
					     TbQ = FMA(KP414213562, T7Q, T7X);
					}
				   }
			      }
			      {
				   E TkJ, TkD, Tck, TbS, TbK, Tku, Tkt, TbN;
				   {
					E TkA, Tby, Tkp, Tbu, Tkz, Tbx;
					{
					     E Tbt, T9a, Tkx, Tky, Tkv;
					     TkA = Taj + Tbs;
					     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 = FMA(KP923879532, T7Z, T7u);
						  T80 = FNMS(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);
					     }
					     Tkz = FNMS(KP831469612, Tkw, Tkv);
					     Tkx = FMA(KP831469612, Tkw, Tkv);
					     Tky = Tbw - Tbv;
					     Tbx = Tbv + Tbw;
					     ri[WS(rs, 11)] = FMA(KP881921264, Tbt, T9a);
					     ri[WS(rs, 43)] = FNMS(KP881921264, Tbt, T9a);
					     ii[WS(rs, 43)] = FNMS(KP881921264, Tky, Tkx);
					     ii[WS(rs, 11)] = FMA(KP881921264, Tky, Tkx);
					}
					{
					     E TbC, TbJ, Tkr, Tks;
					     TbK = FNMS(KP831469612, TbB, Tby);
					     TbC = FMA(KP831469612, TbB, Tby);
					     ri[WS(rs, 59)] = FMA(KP881921264, Tbx, Tbu);
					     ri[WS(rs, 27)] = FNMS(KP881921264, Tbx, Tbu);
					     ii[WS(rs, 59)] = FMA(KP881921264, TkA, Tkz);
					     ii[WS(rs, 27)] = FNMS(KP881921264, TkA, Tkz);
					     TbJ = TbF + TbI;
					     Tku = TbI - TbF;
					     Tkt = FNMS(KP831469612, Tkq, Tkp);
					     Tkr = FMA(KP831469612, Tkq, Tkp);
					     Tks = TbL + TbM;
					     TbN = TbL - TbM;
					     ri[WS(rs, 3)] = FMA(KP956940335, TbJ, TbC);
					     ri[WS(rs, 35)] = FNMS(KP956940335, TbJ, TbC);
					     ii[WS(rs, 35)] = FNMS(KP956940335, Tks, Tkr);
					     ii[WS(rs, 3)] = FMA(KP956940335, Tks, Tkr);
					}
				   }
				   {
					E Tcg, TkI, TkH, Tcj;
					{
					     E Tc0, Tcf, TkF, TkG;
					     Tcg = FNMS(KP980785280, TbZ, TbS);
					     Tc0 = FMA(KP980785280, TbZ, TbS);
					     ri[WS(rs, 19)] = FMA(KP956940335, TbN, TbK);
					     ri[WS(rs, 51)] = FNMS(KP956940335, TbN, TbK);
					     ii[WS(rs, 51)] = FNMS(KP956940335, Tku, Tkt);
					     ii[WS(rs, 19)] = FMA(KP956940335, Tku, Tkt);
					     Tcf = Tc7 + Tce;
					     TkI = Tce - Tc7;
					     TkH = FNMS(KP980785280, TkE, TkD);
					     TkF = FMA(KP980785280, TkE, TkD);
					     TkG = Tch + Tci;
					     Tcj = Tch - Tci;
					     ri[WS(rs, 7)] = FMA(KP773010453, Tcf, Tc0);
					     ri[WS(rs, 39)] = FNMS(KP773010453, Tcf, Tc0);
					     ii[WS(rs, 39)] = FNMS(KP773010453, TkG, TkF);
					     ii[WS(rs, 7)] = FMA(KP773010453, TkG, TkF);
					}
					{
					     E Tco, Tcv, TkL, TkM;
					     Tcw = FMA(KP980785280, Tcn, Tck);
					     Tco = FNMS(KP980785280, Tcn, Tck);
					     ri[WS(rs, 23)] = FMA(KP773010453, Tcj, Tcg);
					     ri[WS(rs, 55)] = FNMS(KP773010453, Tcj, Tcg);
					     ii[WS(rs, 55)] = FNMS(KP773010453, TkI, TkH);
					     ii[WS(rs, 23)] = FMA(KP773010453, TkI, TkH);
					     Tcv = Tcr - Tcu;
					     TkO = Tcr + Tcu;
					     TkN = FMA(KP980785280, TkK, TkJ);
					     TkL = FNMS(KP980785280, TkK, TkJ);
					     TkM = Tcy - Tcx;
					     Tcz = Tcx + Tcy;
					     ri[WS(rs, 15)] = FMA(KP995184726, Tcv, Tco);
					     ri[WS(rs, 47)] = FNMS(KP995184726, Tcv, Tco);
					     ii[WS(rs, 47)] = FNMS(KP995184726, TkM, TkL);
					     ii[WS(rs, 15)] = FMA(KP995184726, TkM, TkL);
					}
				   }
			      }
			 }
			 {
			      E TdN, Tk2, TjW, 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 Te5, Te6, Te8, Te9, Teg, Tee, Ted, Teh, Tdz, Tdv, Tdo, Ten, TeF, Tem, Teo;
				   E TdC, Tex, Te7;
				   {
					E TcP, TcS, TcW, TcZ;
					Te5 = FNMS(KP707106781, TcO, TcN);
					TcP = FMA(KP707106781, TcO, TcN);
					ri[WS(rs, 63)] = FMA(KP995184726, Tcz, Tcw);
					ri[WS(rs, 31)] = FNMS(KP995184726, Tcz, Tcw);
					ii[WS(rs, 63)] = FMA(KP995184726, TkO, TkN);
					ii[WS(rs, 31)] = FNMS(KP995184726, TkO, TkN);
					TcS = FMA(KP707106781, TcR, TcQ);
					Te6 = FNMS(KP707106781, TcR, TcQ);
					Te8 = FNMS(KP707106781, TcV, TcU);
					TcW = FMA(KP707106781, TcV, TcU);