hf_64.c

快速fft变换

				   {
					E TbG, Tbc, Tcl, TbY;
					Tcl = FMA(KP668178637, TbW, TbX);
					TbY = FNMS(KP668178637, TbX, TbW);
					TbG = FNMS(KP923879532, Tbb, TaG);
					Tbc = FMA(KP923879532, Tbb, TaG);
					Tcn = Tcl + Tcm;
					Tkg = Tcl - Tcm;
					Tka = TbY + TbV;
					TbZ = TbV - TbY;
					{
					     E T7t, TjS, TbH, Tbr;
					     Tkm = T7s - T7l;
					     T7t = T7l + T7s;
					     TjS = TcB - TcC;
					     TcD = TcB + TcC;
					     TbP = FMA(KP414213562, T7B, T7I);
					     T7J = FNMS(KP414213562, T7I, T7B);
					     TbH = FNMS(KP923879532, Tbq, Tbn);
					     Tbr = FMA(KP923879532, Tbq, Tbn);
					     TbO = FNMS(KP707106781, T7t, T7e);
					     T7u = FMA(KP707106781, T7t, T7e);
					     Tk7 = FNMS(KP707106781, TjS, TjR);
					     TjT = FMA(KP707106781, TjS, TjR);
					     TbI = FMA(KP820678790, TbH, TbG);
					     TbM = FNMS(KP820678790, TbG, TbH);
					     Tbw = FMA(KP098491403, Tbc, Tbr);
					     Tbs = FNMS(KP098491403, Tbr, Tbc);
					     T7Y = FMA(KP414213562, T7X, T7Q);
					     TbQ = FNMS(KP414213562, T7Q, T7X);
					}
				   }
			      }
			      {
				   E Tk1, TjV, Tck, TbS, Tkd, Tcz, Tkh, Tcf, TjY, Tk6, Tke, Tcv, Tki, Tcj;
				   {
					E Tbu, TbC, Tkb, Tkc, Tkj, Tkk, Tbx, TbJ;
					{
					     E Tbt, Tkf, Tk9, T9a, TbK, TbN, Tby;
					     Tk0 = Tbs - Taj;
					     Tbt = Taj + Tbs;
					     {
						  E Tk8, T7Z, TjU, TbR, T80;
						  Tk8 = T7Y - T7J;
						  T7Z = T7J + T7Y;
						  TjU = TbP + TbQ;
						  TbR = TbP - TbQ;
						  Tkf = FNMS(KP923879532, Tk8, Tk7);
						  Tk9 = FMA(KP923879532, Tk8, Tk7);
						  Tby = FNMS(KP923879532, T7Z, T7u);
						  T80 = FMA(KP923879532, T7Z, T7u);
						  Tk1 = FNMS(KP923879532, TjU, TjT);
						  TjV = FMA(KP923879532, TjU, TjT);
						  Tck = FMA(KP923879532, TbR, TbO);
						  TbS = FNMS(KP923879532, TbR, TbO);
						  T9a = FMA(KP980785280, T99, T80);
						  Tbu = FNMS(KP980785280, T99, T80);
					     }
					     TbC = FMA(KP980785280, TbB, Tby);
					     TbK = FNMS(KP980785280, TbB, Tby);
					     TbN = TbL + TbM;
					     Tk4 = TbL - TbM;
					     Tkd = FNMS(KP831469612, Tka, Tk9);
					     Tkb = FMA(KP831469612, Tka, Tk9);
					     ci[0] = FMA(KP995184726, Tbt, T9a);
					     cr[WS(rs, 31)] = FNMS(KP995184726, Tbt, T9a);
					     ci[WS(rs, 8)] = FNMS(KP773010453, TbN, TbK);
					     cr[WS(rs, 23)] = FMA(KP773010453, TbN, TbK);
					     Tkc = Tcx - Tcy;
					     Tcz = Tcx + Tcy;
					     Tkh = FMA(KP831469612, Tkg, Tkf);
					     Tkj = FNMS(KP831469612, Tkg, Tkf);
					     Tkk = Tce - Tc7;
					     Tcf = Tc7 + Tce;
					}
					ci[WS(rs, 60)] = FMA(KP956940335, Tkc, Tkb);
					cr[WS(rs, 35)] = FMS(KP956940335, Tkc, Tkb);
					ci[WS(rs, 52)] = FMA(KP881921264, Tkk, Tkj);
					cr[WS(rs, 43)] = FMS(KP881921264, Tkk, Tkj);
					Tbx = Tbv + Tbw;
					TjY = Tbw - Tbv;
					TbJ = TbF + TbI;
					Tk6 = TbI - TbF;
					cr[WS(rs, 15)] = FMA(KP995184726, Tbx, Tbu);
					ci[WS(rs, 16)] = FNMS(KP995184726, Tbx, Tbu);
					cr[WS(rs, 7)] = FMA(KP773010453, TbJ, TbC);
					ci[WS(rs, 24)] = FNMS(KP773010453, TbJ, TbC);
					Tke = Tcu - Tcr;
					Tcv = Tcr + Tcu;
					Tki = Tci - Tch;
					Tcj = Tch + Tci;
				   }
				   {
					E Tcg, Tco, TjX, Tk5, Tc0, Tcw;
					Tcg = FNMS(KP831469612, TbZ, TbS);
					Tc0 = FMA(KP831469612, TbZ, TbS);
					ci[WS(rs, 44)] = FMA(KP956940335, Tke, Tkd);
					cr[WS(rs, 51)] = FMS(KP956940335, Tke, Tkd);
					ci[WS(rs, 36)] = FMA(KP881921264, Tki, Tkh);
					cr[WS(rs, 59)] = FMS(KP881921264, Tki, Tkh);
					Tco = FMA(KP831469612, Tcn, Tck);
					Tcw = FNMS(KP831469612, Tcn, Tck);
					TjZ = FNMS(KP980785280, TjW, TjV);
					TjX = FMA(KP980785280, TjW, TjV);
					ci[WS(rs, 4)] = FMA(KP881921264, Tcf, Tc0);
					cr[WS(rs, 27)] = FNMS(KP881921264, Tcf, Tc0);
					ci[WS(rs, 12)] = FNMS(KP956940335, Tcz, Tcw);
					cr[WS(rs, 19)] = FMA(KP956940335, Tcz, Tcw);
					Tk3 = FMA(KP980785280, Tk2, Tk1);
					Tk5 = FNMS(KP980785280, Tk2, Tk1);
					ci[WS(rs, 32)] = FMA(KP995184726, TjY, TjX);
					cr[WS(rs, 63)] = FMS(KP995184726, TjY, TjX);
					ci[WS(rs, 40)] = FMA(KP773010453, Tk6, Tk5);
					cr[WS(rs, 55)] = FMS(KP773010453, Tk6, Tk5);
					cr[WS(rs, 11)] = FMA(KP881921264, Tcj, Tcg);
					ci[WS(rs, 20)] = FNMS(KP881921264, Tcj, Tcg);
					cr[WS(rs, 3)] = FMA(KP956940335, Tcv, Tco);
					ci[WS(rs, 28)] = FNMS(KP956940335, Tcv, Tco);
				   }
			      }
			 }
		    }
		    {
			 E Ti8, Thu, Tjf, Tj9, Tib, Tjg, Tja, ThF, Tig, ThW, Tif, Til, Ti6, ThR;
			 ci[WS(rs, 48)] = FMA(KP995184726, Tk0, TjZ);
			 cr[WS(rs, 47)] = FMS(KP995184726, Tk0, TjZ);
			 ci[WS(rs, 56)] = FMA(KP773010453, Tk4, Tk3);
			 cr[WS(rs, 39)] = FMS(KP773010453, Tk4, Tk3);
			 Ti8 = Thq + Tht;
			 Thu = Thq - Tht;
			 Tjf = Tj8 - Tj7;
			 Tj9 = Tj7 + Tj8;
			 {
			      E Tid, ThL, Tie, ThQ;
			      {
				   E Ti9, Thz, Tia, ThE;
				   Ti9 = Thv - Thy;
				   Thz = Thv + Thy;
				   Tia = ThA + ThD;
				   ThE = ThA - ThD;
				   Tib = Ti9 + Tia;
				   Tjg = Tia - Ti9;
				   Tja = Thz - ThE;
				   ThF = Thz + ThE;
				   Tid = ThH + ThK;
				   ThL = ThH - ThK;
			      }
			      Tie = ThO + ThP;
			      ThQ = ThO - ThP;
			      Tig = ThS + ThV;
			      ThW = ThS - ThV;
			      Tif = FNMS(KP414213562, Tie, Tid);
			      Til = FMA(KP414213562, Tid, Tie);
			      Ti6 = FNMS(KP414213562, ThL, ThQ);
			      ThR = FMA(KP414213562, ThQ, ThL);
			 }
			 {
			      E Ti4, ThG, Tjh, Tjj, Tih, Ti1;
			      Ti4 = FNMS(KP707106781, ThF, Thu);
			      ThG = FMA(KP707106781, ThF, Thu);
			      Tjh = FMA(KP707106781, Tjg, Tjf);
			      Tjj = FNMS(KP707106781, Tjg, Tjf);
			      Tih = Ti0 - ThZ;
			      Ti1 = ThZ + Ti0;
			      {
				   E Tje, Tjd, Tjb, Tjc;
				   {
					E Tic, Tim, Ti5, Ti2, Tij, Tii;
					Tik = FNMS(KP707106781, Tib, Ti8);
					Tic = FMA(KP707106781, Tib, Ti8);
					Tii = FNMS(KP414213562, Tih, Tig);
					Tim = FMA(KP414213562, Tig, Tih);
					Ti5 = FMA(KP414213562, ThW, Ti1);
					Ti2 = FNMS(KP414213562, Ti1, ThW);
					Tij = Tif + Tii;
					Tje = Tii - Tif;
					Tjd = FNMS(KP707106781, Tja, Tj9);
					Tjb = FMA(KP707106781, Tja, Tj9);
					{
					     E Ti7, Tji, Tjk, Ti3;
					     Ti7 = Ti5 - Ti6;
					     Tji = Ti6 + Ti5;
					     Tjk = Ti2 - ThR;
					     Ti3 = ThR + Ti2;
					     ci[WS(rs, 3)] = FMA(KP923879532, Tij, Tic);
					     cr[WS(rs, 28)] = FNMS(KP923879532, Tij, Tic);
					     ci[WS(rs, 11)] = FMA(KP923879532, Ti7, Ti4);
					     cr[WS(rs, 20)] = FNMS(KP923879532, Ti7, Ti4);
					     ci[WS(rs, 59)] = FMA(KP923879532, Tji, Tjh);
					     cr[WS(rs, 36)] = FMS(KP923879532, Tji, Tjh);
					     ci[WS(rs, 43)] = FMA(KP923879532, Tjk, Tjj);
					     cr[WS(rs, 52)] = FMS(KP923879532, Tjk, Tjj);
					     cr[WS(rs, 4)] = FMA(KP923879532, Ti3, ThG);
					     ci[WS(rs, 27)] = FNMS(KP923879532, Ti3, ThG);
					     Tjc = Tim - Til;
					     Tin = Til + Tim;
					}
				   }
				   ci[WS(rs, 35)] = FMA(KP923879532, Tjc, Tjb);
				   cr[WS(rs, 60)] = FMS(KP923879532, Tjc, Tjb);
				   ci[WS(rs, 51)] = FMA(KP923879532, Tje, Tjd);
				   cr[WS(rs, 44)] = FMS(KP923879532, Tje, Tjd);
			      }
			 }
		    }
		    {
			 E Tjy, Tju, Tjt, Tjx;
			 {
			      E TjD, TjJ, Tgo, Tf2, Tjp, Tjv, Tha, TgI, Tgg, Tgd, Tgr, Tjw, Tjq, Tfp, Thk;
			      E Tho, Th7, Th4, Tgv, TgB, Tgl, TfR, TjE, Thd, TjK, TgP, Tgw, Tg8, Thh, Thn;
			      E Th8, TgX;
			      {
				   E TgK, TgJ, TgN, TgM, TfW, Th1, Thi, 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 = TeR + TeU;
					TeV = TeR - TeU;
					cr[WS(rs, 12)] = FMA(KP923879532, Tin, Tik);
					ci[WS(rs, 19)] = FNMS(KP923879532, Tin, Tik);
					Tf0 = TeW + TeZ;
					TgG = TeW - TeZ;
					TjC = Tf0 - TeV;
					Tf1 = TeV + Tf0;
					{
					     E Tfi, Tgp, Tfd, Tfn;
					     {
						  E Tf7, Tjo, TgH, Tfc;
						  TgK = 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;
						  TgJ = Tf8 - Tfb;
						  TgN = 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 = FNMS(KP414213562, Tf7, Tfc);
						  Tfd = FMA(KP414213562, Tfc, Tf7);
						  Tfn = Tfj + Tfm;
						  TgM = Tfj - Tfm;
					     }
					     {
						  E TgY, TgZ, Tgq, Tfo;
						  TfW = TfS + TfV;
						  TgY = TfS - TfV;
						  TgZ = Tgf + Tge;
						  Tgg = Tge - Tgf;
						  Tgd = Tg9 - Tgc;
						  Th1 = Tgc + Tg9;
						  Tgq = FMA(KP414213562, Tfi, Tfn);
						  Tfo = FNMS(KP414213562, Tfn, Tfi);
						  Thi = FNMS(KP707106781, TgZ, TgY);
						  Th0 = FMA(KP707106781, TgZ, TgY);
						  Tgr = Tgp + Tgq;
						  Tjw = Tgq - Tgp;
						  Tjq = Tfd + Tfo;
						  Tfp = Tfd - Tfo;
						  Th2 = Tg6 - Tg1;
						  Tg7 = Tg1 + Tg6;
					     }
					}
				   }
				   {
					E TgR, TgV, TgU, TgS, Thc, TgL;
					{
					     E TfM, Tgt, TfH, TfP, Tgu, TfQ;
					     {
						  E Tfv, TfG, Thj, Th3;
						  TgR = Tfr - Tfu;
						  Tfv = Tfr + Tfu;
						  TfG = TfA + TfF;
						  TgV = TfF - TfA;
						  TgU = TfK - TfL;
						  TfM = TfK + TfL;
						  Thj = FNMS(KP707106781, Th2, Th1);
						  Th3 = FMA(KP707106781, Th2, Th1);
						  Tgt = FMA(KP707106781, TfG, Tfv);
						  TfH = FNMS(KP707106781, TfG, Tfv);
						  Thk = FMA(KP668178637, Thj, Thi);
						  Tho = FNMS(KP668178637, Thi, Thj);
						  Th7 = FMA(KP198912367, Th0, Th3);
						  Th4 = FNMS(KP198912367, Th3, Th0);
						  TfP = TfN + TfO;
						  TgS = TfN - TfO;
					     }
					     Tgu = FMA(KP707106781, TfP, TfM);
					     TfQ = FNMS(KP707106781, TfP, TfM);
					     Thc = FNMS(KP414213562, TgJ, TgK);
					     TgL = FMA(KP414213562, TgK, TgJ);
					     Tgv = FNMS(KP198912367, Tgu, Tgt);
					     TgB = FMA(KP198912367, Tgt, Tgu);
					     Tgl = FNMS(KP668178637, TfH, TfQ);
					     TfR = FMA(KP668178637, TfQ, TfH);
					}
					{
					     E Thf, TgT, Thb, TgO, Thg, TgW;
					     Thb = FMA(KP414213562, TgM, TgN);
					     TgO = FNMS(KP414213562, TgN, TgM);
					     Thf = FNMS(KP707106781, TgS, TgR);
					     TgT = FMA(KP707106781, TgS, TgR);
					     TjE = Thc + Thb;
					     Thd = Thb - Thc;
					     TjK = TgL - TgO;
					     TgP = TgL + TgO;
					     Thg = FNMS(KP707106781, TgV, TgU);
					     TgW = FMA(KP707106781, TgV, TgU);
					     Tgw = FMA(KP707106781, Tg7, TfW);
					     Tg8 = FNMS(KP707106781, Tg7, TfW);
					     Thh = FNMS(KP668178637, Thg, Thf);
					     Thn = FMA(KP668178637, Thf, Thg);
					     Th8 = FNMS(KP198912367, TgT, TgW);
					     TgX = FMA(KP198912367, TgW, TgT);
					}
				   }
			      }
			      {
				   E TjH, Th9, TjL, Tjs, TjA, Thl, TjI, Th5, TjM, Thp;
				   {
					E Tgk, Tfq, TgA, Tgs, TjN, Tgy, Tgm, TgD, Tgj, TjO, Tgn, Tgz;
					Tgk = FNMS(KP923879532, Tfp, Tf2);
					Tfq = FMA(KP923879532, Tfp, Tf2);
					TgA = FNMS(KP923879532, Tgr, Tgo);
					Tgs = FMA(KP923879532, Tgr, Tgo);
					{
					     E TjF, Tgx, Tgh, TjG, TgC, Tgi;
					     TjH = FNMS(KP923879532, TjE, TjD);
					     TjF = FMA(KP923879532, TjE, TjD);
					     Tgx = FMA(KP707106781, Tgg, Tgd);
					     Tgh = FNMS(KP707106781, Tgg, Tgd);
					     TjG = Th8 + Th7;
					     Th9 = Th7 - Th8;
					     TjL = FMA(KP923879532, TjK, TjJ);
					     TjN = FNMS(KP923879532, TjK, TjJ);
					     Tgy = FNMS(KP198912367, Tgx, Tgw);
					     TgC = FMA(KP198912367, Tgw, Tgx);
					     Tgm = FNMS(KP668178637, Tg8, Tgh);
					     Tgi = FMA(KP668178637, Tgh, Tg8);
					     ci[WS(rs, 61)] = FMA(KP980785280, TjG, TjF);
					     cr[WS(rs, 34)] = FMS(KP980785280, TjG, TjF);
					     TgD = TgB + TgC;
					     Tjs = TgC - TgB;
					     TjA = Tgi - TfR;
					     Tgj = TfR + Tgi;
					     TjO = Thk - Thh;
					     Thl = Thh + Thk;
					}
					cr[WS(rs, 14)] = FMA(KP980785280, TgD, TgA);
					ci[WS(rs, 17)] = FNMS(KP980785280, TgD, TgA);
					cr[WS(rs, 6)] = FMA(KP831469612, Tgj, Tfq);
					ci[WS(rs, 25)] = FNMS(KP831469612, Tgj, Tfq);
					ci[WS(rs, 53)] = FMA(KP831469612, TjO, TjN);
					cr[WS(rs, 42)] = FMS(KP831469612, TjO, TjN);
					Tgn = Tgl + Tgm;
					Tjy = Tgl - Tgm;
					Tgz = Tgv + Tgy;
					Tju = Tgy - Tgv;
					ci[WS(rs, 9)] = FNMS(KP831469612, Tgn, Tgk);
					cr[WS(rs, 22)] = FMA(KP831469612, Tgn, Tgk);
					ci[WS(rs, 1)] = FMA(KP980785280, Tgz, Tgs);
					cr[WS(rs, 30)] = FNMS(KP980785280, Tgz, Tgs);
					TjI = Th4 - TgX;
					Th5 = TgX + Th4;
					TjM = Thn + Tho;
					Thp = Thn - Tho;
				   }
				   {
					E Th6, The, Tjr, Tjz, TgQ, Thm;
					Th6 = FNMS(KP923879532, TgP, TgI);
					TgQ = FMA(KP923879532, TgP, TgI);
					ci[WS(rs, 45)] = FMA(KP980785280, TjI, TjH);
					cr[WS(rs, 50)] = FMS(KP980785280, TjI, TjH);
					ci[WS(rs, 37)] = FNMS(KP831469612, TjM, TjL);
					cr[WS(rs, 58)] = -(FMA(KP831469612, TjM, TjL));
					The = FMA(KP923879532, Thd, Tha);
					Thm = FNMS(KP923879532, Thd, Tha);
					Tjt = FNMS(KP923879532, Tjq, Tjp);
					Tjr = FMA(KP923879532, Tjq, Tjp);
					cr[WS(rs, 2)] = FMA(KP980785280, Th5, TgQ);
					ci[WS(rs, 29)] = FNMS(KP980785280, Th5, TgQ);
					cr[WS(rs, 10)] = FMA(KP831469612, Thp, Thm);
					ci[WS(rs, 21)] = FNMS(KP831469612, Thp, Thm);
					Tjx = FMA(KP923879532, Tjw, Tjv);
					Tjz = FNMS(KP923879532, Tjw, Tjv);
					ci[WS(rs, 33)] = FMA(KP980785280, Tjs, Tjr);
					cr[WS(rs, 62)] = FMS(KP980785280, Tjs, Tjr);
					ci[WS(rs, 41)] = FMA(KP831469612, TjA, Tjz);
					cr[WS(rs, 54)] = FMS(KP831469612, TjA, Tjz);
					ci[WS(rs, 13)] = FMA(KP980785280, Th9, Th6);
					cr[WS(rs, 18)] = FNMS(KP980785280, Th9, Th6);
					ci[WS(rs, 5)] = FMA(KP831469612, Thl, The);
					cr[WS(rs, 26)] = FNMS(KP831469612, Thl, The);
				   }
			      }
			 }
			 {
			      E Tkq, TdN, Tkw, Td1, TdR, TdX, TdI, Tdl, TeG, TeK, Tet, Teq, TeD, TeJ, Teu;
			      E Tej, Tez, TkK, TkE, Teb, Te2, TcH, Te0, TcE, TkB, Tkn, TdU, TdY, TdH, TdE;
			      E TcK, Te1;
			      {
				   E Te6, Te5, Te9, Te8, Ted, Teh, Teg, Tee, Tdo, TdC, Tdz, Ten, TeE, Tem, Teo;
				   E Tdv, Tex, Te7;
				   {
					E TcP, TcS, TcW, TcZ;
					Te6 = FNMS(KP707106781, TcO, TcN);
					TcP = FMA(KP707106781, TcO, TcN);
					ci[WS(rs, 49)] = FMA(KP980785280, Tju, Tjt);
					cr[WS(rs, 46)] = FMS(KP980785280, Tju, Tjt);
					ci[WS(rs, 57)] = FMA(KP831469612, Tjy, Tjx);
					cr[WS(rs, 38)] = FMS(KP831469612, Tjy, Tjx);
					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, TdM, TcT, TdL, Td0, Tdc;
						  Ted = FNMS(KP707106781, Td4, Td3);
						  Td5 = FMA(KP707106781, Td4, Td3);
						  TdM = FNMS(KP198912367, TcP, TcS);
						  TcT = FMA(KP198912367, TcS, TcP);
						  TdL = FMA(KP198912367, TcW, TcZ);
						  Td0 = FNMS(KP198912367, TcZ, TcW);
						  Tdc = Td8 + Tdb;
						  Teh = Td8 - Tdb;
						  Teg = FNMS(KP707106781, Tdf, Tde);
						  Tdg = FMA(KP707106781, Tdf, Tde);
						  Tkq = TdM + TdL;
						  TdN = TdL - TdM;
						  Tkw = TcT - Td0;
						  Td1 = TcT + Td0;
						  TdP = FNMS(KP923879532, Tdc, Td5);
						  Tdd = FMA(KP923879532, Tdc, Td5);
						  Tdj = Tdh + Tdi;
						  Tee = Tdi - Tdh;
					     }
					     {