t2_64.c

快速fft变换

							    TeL = TbZ + Tc4;
							    Thn = Thl - Thm;
							    Tjh = Thl + Thm;
							    TeT = FNMS(KP414213562, TeK, TeL);
							    TeM = FMA(KP414213562, TeL, TeK);
							    Tcd = FMA(KP414213562, Tcc, Tc5);
							    TcH = FNMS(KP414213562, Tc5, Tcc);
							    Tho = Thk + Thn;
							    Thx = Thk - Thn;
							    TjG = T1C - T2H;
							    T2I = T1C + T2H;
						       }
						  }
					     }
					}
					TiQ = T39 - T3K;
					T3L = T39 + T3K;
					T4M = T4k + T4L;
					TiS = T4k - T4L;
					TiL = TiJ - TiK;
					Tk6 = TiJ + TiK;
					Tke = Tk7 + Tkd;
					Tkq = Tkd - Tk7;
					TiP = TiN - TiO;
					TjH = TiN + TiO;
					Tkj = T4M - T3L;
					T4N = T3L + T4M;
					Tkk = Tke - Tk6;
					Tkf = Tk6 + Tke;
					TjI = TiT + TiU;
					TiV = TiT - TiU;
					{
					     E TjR, TjQ, TjS, T7w, T8n;
					     Tjf = T74 - T7v;
					     T7w = T74 + T7v;
					     T8n = T7V + T8m;
					     Tjd = T8m - T7V;
					     Tjc = Tja - Tjb;
					     TjR = Tja + Tjb;
					     Tk5 = TjH + TjI;
					     TjJ = TjH - TjI;
					     TjQ = T7w - T8n;
					     T8o = T7w + T8n;
					     Tji = Tjg - Tjh;
					     TjS = Tjg + Tjh;
					     {
						  E TjM, TjN, T5J, T6C, TjT;
						  Tj4 = T5d - T5I;
						  T5J = T5d + T5I;
						  T6C = T68 + T6B;
						  Tj2 = T6B - T68;
						  TjT = TjR - TjS;
						  Tk2 = TjR + TjS;
						  Tj1 = TiZ - Tj0;
						  TjM = TiZ + Tj0;
						  TjL = T5J - T6C;
						  T6D = T5J + T6C;
						  Tj7 = Tj5 - Tj6;
						  TjN = Tj5 + Tj6;
						  TjY = TjQ + TjT;
						  TjU = TjQ - TjT;
						  Tk1 = TjM + TjN;
						  TjO = TjM - TjN;
					     }
					}
				   }
				   {
					E Tk0, Tk3, TjW, Tko, Tkn, Tkl, Tkm, TjZ;
					{
					     E TjP, TjX, Tk4, Tkg, T4O, T8p, TjK, TjV;
					     Tk0 = T2I - T4N;
					     T4O = T2I + T4N;
					     T8p = T6D + T8o;
					     Tkh = T8o - T6D;
					     TjP = TjL + TjO;
					     TjX = TjO - TjL;
					     Tk3 = Tk1 - Tk2;
					     Tk4 = Tk1 + Tk2;
					     ri[0] = T4O + T8p;
					     ri[WS(rs, 32)] = T4O - T8p;
					     Tkg = Tk5 + Tkf;
					     Tki = Tkf - Tk5;
					     TjW = TjG - TjJ;
					     TjK = TjG + TjJ;
					     TjV = TjP + TjU;
					     Tko = TjU - TjP;
					     Tkn = Tkk - Tkj;
					     Tkl = Tkj + Tkk;
					     ii[WS(rs, 32)] = Tkg - Tk4;
					     ii[0] = Tk4 + Tkg;
					     ri[WS(rs, 8)] = FMA(KP707106781, TjV, TjK);
					     ri[WS(rs, 40)] = FNMS(KP707106781, TjV, TjK);
					     Tkm = TjX + TjY;
					     TjZ = TjX - TjY;
					}
					ii[WS(rs, 40)] = FNMS(KP707106781, Tkm, Tkl);
					ii[WS(rs, 8)] = FMA(KP707106781, Tkm, Tkl);
					ri[WS(rs, 24)] = FMA(KP707106781, TjZ, TjW);
					ri[WS(rs, 56)] = FNMS(KP707106781, TjZ, TjW);
					ii[WS(rs, 56)] = FNMS(KP707106781, Tko, Tkn);
					ii[WS(rs, 24)] = FMA(KP707106781, Tko, Tkn);
					ri[WS(rs, 16)] = Tk0 + Tk3;
					ri[WS(rs, 48)] = Tk0 - Tk3;
				   }
			      }
			      {
				   E Tjq, TiM, Tkx, Tkr, Tjt, Tky, Tks, TiX, Tjz, Tje, Tjx, TjD, Tjn, Tj9, Tjr;
				   E TiR;
				   ii[WS(rs, 48)] = Tki - Tkh;
				   ii[WS(rs, 16)] = Tkh + Tki;
				   Tjq = TiI + TiL;
				   TiM = TiI - TiL;
				   Tkx = Tkq - Tkp;
				   Tkr = Tkp + Tkq;
				   Tjr = TiQ + TiP;
				   TiR = TiP - TiQ;
				   {
					E Tjw, Tj3, Tjs, TiW, Tjv, Tj8;
					Tjs = TiS - TiV;
					TiW = TiS + TiV;
					Tjw = Tj1 + Tj2;
					Tj3 = Tj1 - Tj2;
					Tjt = Tjr + Tjs;
					Tky = Tjs - Tjr;
					Tks = TiR + TiW;
					TiX = TiR - TiW;
					Tjv = Tj4 + Tj7;
					Tj8 = Tj4 - Tj7;
					Tjz = Tjc + Tjd;
					Tje = Tjc - Tjd;
					Tjx = FMA(KP414213562, Tjw, Tjv);
					TjD = FNMS(KP414213562, Tjv, Tjw);
					Tjn = FNMS(KP414213562, Tj3, Tj8);
					Tj9 = FMA(KP414213562, Tj8, Tj3);
				   }
				   {
					E Tjm, TiY, Tkz, TkB, Tjy, Tjj;
					Tjm = FNMS(KP707106781, TiX, TiM);
					TiY = FMA(KP707106781, TiX, TiM);
					Tkz = FMA(KP707106781, Tky, Tkx);
					TkB = FNMS(KP707106781, Tky, Tkx);
					Tjy = Tjf + Tji;
					Tjj = Tjf - Tji;
					{
					     E TjC, Tkt, Tku, TjF;
					     {
						  E Tju, TjE, Tjo, Tjk, TjB, TjA;
						  TjC = FNMS(KP707106781, Tjt, Tjq);
						  Tju = FMA(KP707106781, Tjt, Tjq);
						  TjA = FNMS(KP414213562, Tjz, Tjy);
						  TjE = FMA(KP414213562, Tjy, Tjz);
						  Tjo = FMA(KP414213562, Tje, Tjj);
						  Tjk = FNMS(KP414213562, Tjj, Tje);
						  TjB = Tjx + TjA;
						  Tkw = TjA - Tjx;
						  Tkv = FNMS(KP707106781, Tks, Tkr);
						  Tkt = FMA(KP707106781, Tks, Tkr);
						  {
						       E Tjp, TkA, TkC, Tjl;
						       Tjp = Tjn + Tjo;
						       TkA = Tjo - Tjn;
						       TkC = Tj9 + Tjk;
						       Tjl = Tj9 - Tjk;
						       ri[WS(rs, 4)] = FMA(KP923879532, TjB, Tju);
						       ri[WS(rs, 36)] = FNMS(KP923879532, TjB, Tju);
						       ri[WS(rs, 60)] = FMA(KP923879532, Tjp, Tjm);
						       ri[WS(rs, 28)] = FNMS(KP923879532, Tjp, Tjm);
						       ii[WS(rs, 44)] = FNMS(KP923879532, TkA, Tkz);
						       ii[WS(rs, 12)] = FMA(KP923879532, TkA, Tkz);
						       ii[WS(rs, 60)] = FMA(KP923879532, TkC, TkB);
						       ii[WS(rs, 28)] = FNMS(KP923879532, TkC, TkB);
						       ri[WS(rs, 12)] = FMA(KP923879532, Tjl, TiY);
						       ri[WS(rs, 44)] = FNMS(KP923879532, Tjl, TiY);
						       Tku = TjD + TjE;
						       TjF = TjD - TjE;
						  }
					     }
					     ii[WS(rs, 36)] = FNMS(KP923879532, Tku, Tkt);
					     ii[WS(rs, 4)] = FMA(KP923879532, Tku, Tkt);
					     ri[WS(rs, 20)] = FMA(KP923879532, TjF, TjC);
					     ri[WS(rs, 52)] = FNMS(KP923879532, TjF, TjC);
					}
				   }
			      }
			 }
			 {
			      E TkV, Tl1, ThG, Tgk, TkH, TkN, Tis, Ti0, Thv, ThJ, TkO, TkI, TgH, Thy, TiC;
			      E TiG, Tiq, Tim, ThN, ThT, ThD, Th9, TkW, Tiv, Tl2, Ti7, ThP, Thq, Tiz, TiF;
			      E Tip, Tif;
			      {
				   E Ti1, Ti2, Ti4, Ti5, Thp, The, Tij, TiB, Tii, Tik;
				   {
					E ThW, Tg8, TkT, TkF, ThX, ThY, TkU, Tgj, Tgd, Tgi;
					ThW = Tg4 - Tg7;
					Tg8 = Tg4 + Tg7;
					TkT = TkE - TkD;
					TkF = TkD + TkE;
					ThX = Tgc - Tg9;
					Tgd = Tg9 + Tgc;
					ii[WS(rs, 52)] = FNMS(KP923879532, Tkw, Tkv);
					ii[WS(rs, 20)] = FMA(KP923879532, Tkw, Tkv);
					Tgi = Tge - Tgh;
					ThY = Tge + Tgh;
					TkU = Tgi - Tgd;
					Tgj = Tgd + Tgi;
					{
					     E TgA, ThH, Tgv, TgF;
					     {
						  E Tgp, TkG, ThZ, Tgu;
						  Ti1 = Tgn - Tgo;
						  Tgp = Tgn + Tgo;
						  TkV = FMA(KP707106781, TkU, TkT);
						  Tl1 = FNMS(KP707106781, TkU, TkT);
						  ThG = FMA(KP707106781, Tgj, Tg8);
						  Tgk = FNMS(KP707106781, Tgj, Tg8);
						  TkG = ThX + ThY;
						  ThZ = ThX - ThY;
						  Tgu = Tgq + Tgt;
						  Ti2 = Tgq - Tgt;
						  Ti4 = Tgy - Tgz;
						  TgA = Tgy + Tgz;
						  TkH = FMA(KP707106781, TkG, TkF);
						  TkN = FNMS(KP707106781, TkG, TkF);
						  Tis = FNMS(KP707106781, ThZ, ThW);
						  Ti0 = FMA(KP707106781, ThZ, ThW);
						  ThH = FMA(KP414213562, Tgp, Tgu);
						  Tgv = FNMS(KP414213562, Tgu, Tgp);
						  TgF = TgB + TgE;
						  Ti5 = TgB - TgE;
					     }
					     {
						  E Tig, Tih, ThI, TgG;
						  Thv = Thr + Thu;
						  Tig = Thr - Thu;
						  Tih = Tho - Thj;
						  Thp = Thj + Tho;
						  The = Thc + Thd;
						  Tij = Thc - Thd;
						  ThI = FNMS(KP414213562, TgA, TgF);
						  TgG = FMA(KP414213562, TgF, TgA);
						  TiB = FMA(KP707106781, Tih, Tig);
						  Tii = FNMS(KP707106781, Tih, Tig);
						  ThJ = ThH + ThI;
						  TkO = ThI - ThH;
						  TkI = Tgv + TgG;
						  TgH = Tgv - TgG;
						  Tik = Thw - Thx;
						  Thy = Thw + Thx;
					     }
					}
				   }
				   {
					E Tic, Tia, Ti9, Tid, Tit, Ti3;
					{
					     E Th4, ThM, TgZ, Th7, ThL, Th8;
					     {
						  E TgN, TgY, TiA, Til;
						  Tic = TgL - TgM;
						  TgN = TgL + TgM;
						  TgY = TgS + TgX;
						  Tia = TgX - TgS;
						  Ti9 = Th0 - Th3;
						  Th4 = Th0 + Th3;
						  TiA = FMA(KP707106781, Tik, Tij);
						  Til = FNMS(KP707106781, Tik, Tij);
						  ThM = FMA(KP707106781, TgY, TgN);
						  TgZ = FNMS(KP707106781, TgY, TgN);
						  TiC = FNMS(KP198912367, TiB, TiA);
						  TiG = FMA(KP198912367, TiA, TiB);
						  Tiq = FMA(KP668178637, Tii, Til);
						  Tim = FNMS(KP668178637, Til, Tii);
						  Th7 = Th5 + Th6;
						  Tid = Th5 - Th6;
					     }
					     ThL = FMA(KP707106781, Th7, Th4);
					     Th8 = FNMS(KP707106781, Th7, Th4);
					     Tit = FNMS(KP414213562, Ti1, Ti2);
					     Ti3 = FMA(KP414213562, Ti2, Ti1);
					     ThN = FMA(KP198912367, ThM, ThL);
					     ThT = FNMS(KP198912367, ThL, ThM);
					     ThD = FNMS(KP668178637, TgZ, Th8);
					     Th9 = FMA(KP668178637, Th8, TgZ);
					}
					{
					     E Tiy, Tib, Tiu, Ti6, Tix, Tie;
					     Tiu = FMA(KP414213562, Ti4, Ti5);
					     Ti6 = FNMS(KP414213562, Ti5, Ti4);
					     Tiy = FMA(KP707106781, Tia, Ti9);
					     Tib = FNMS(KP707106781, Tia, Ti9);
					     TkW = Tiu - Tit;
					     Tiv = Tit + Tiu;
					     Tl2 = Ti3 + Ti6;
					     Ti7 = Ti3 - Ti6;
					     Tix = FMA(KP707106781, Tid, Tic);
					     Tie = FNMS(KP707106781, Tid, Tic);
					     ThP = FMA(KP707106781, Thp, The);
					     Thq = FNMS(KP707106781, Thp, The);
					     Tiz = FMA(KP198912367, Tiy, Tix);
					     TiF = FNMS(KP198912367, Tix, Tiy);
					     Tip = FNMS(KP668178637, Tib, Tie);
					     Tif = FMA(KP668178637, Tie, Tib);
					}
				   }
			      }
			      {
				   E TkM, TkL, Tl0, TkZ;
				   {
					E ThC, TgI, TkP, TkR, ThO, Thz;
					ThC = FNMS(KP923879532, TgH, Tgk);
					TgI = FMA(KP923879532, TgH, Tgk);
					TkP = FMA(KP923879532, TkO, TkN);
					TkR = FNMS(KP923879532, TkO, TkN);
					ThO = FMA(KP707106781, Thy, Thv);
					Thz = FNMS(KP707106781, Thy, Thv);
					{
					     E ThS, TkJ, TkK, ThV;
					     {
						  E ThK, ThU, ThE, ThA, ThR, ThQ;
						  ThS = FNMS(KP923879532, ThJ, ThG);
						  ThK = FMA(KP923879532, ThJ, ThG);
						  ThQ = FNMS(KP198912367, ThP, ThO);
						  ThU = FMA(KP198912367, ThO, ThP);
						  ThE = FMA(KP668178637, Thq, Thz);
						  ThA = FNMS(KP668178637, Thz, Thq);
						  ThR = ThN + ThQ;
						  TkM = ThQ - ThN;
						  TkL = FNMS(KP923879532, TkI, TkH);
						  TkJ = FMA(KP923879532, TkI, TkH);
						  {
						       E ThF, TkQ, TkS, ThB;
						       ThF = ThD + ThE;
						       TkQ = ThE - ThD;
						       TkS = Th9 + ThA;
						       ThB = Th9 - ThA;
						       ri[WS(rs, 2)] = FMA(KP980785280, ThR, ThK);
						       ri[WS(rs, 34)] = FNMS(KP980785280, ThR, ThK);
						       ri[WS(rs, 58)] = FMA(KP831469612, ThF, ThC);
						       ri[WS(rs, 26)] = FNMS(KP831469612, ThF, ThC);
						       ii[WS(rs, 42)] = FNMS(KP831469612, TkQ, TkP);
						       ii[WS(rs, 10)] = FMA(KP831469612, TkQ, TkP);
						       ii[WS(rs, 58)] = FMA(KP831469612, TkS, TkR);
						       ii[WS(rs, 26)] = FNMS(KP831469612, TkS, TkR);
						       ri[WS(rs, 10)] = FMA(KP831469612, ThB, TgI);
						       ri[WS(rs, 42)] = FNMS(KP831469612, ThB, TgI);
						       TkK = ThT + ThU;
						       ThV = ThT - ThU;
						  }
					     }
					     ii[WS(rs, 34)] = FNMS(KP980785280, TkK, TkJ);
					     ii[WS(rs, 2)] = FMA(KP980785280, TkK, TkJ);
					     ri[WS(rs, 18)] = FMA(KP980785280, ThV, ThS);
					     ri[WS(rs, 50)] = FNMS(KP980785280, ThV, ThS);
					}
				   }
				   {
					E Tio, TkX, TkY, Tir, Ti8, Tin;
					Tio = FNMS(KP923879532, Ti7, Ti0);
					Ti8 = FMA(KP923879532, Ti7, Ti0);
					Tin = Tif + Tim;
					Tl0 = Tim - Tif;
					TkZ = FNMS(KP923879532, TkW, TkV);
					TkX = FMA(KP923879532, TkW, TkV);
					ii[WS(rs, 50)] = FNMS(KP980785280, TkM, TkL);
					ii[WS(rs, 18)] = FMA(KP980785280, TkM, TkL);
					ri[WS(rs, 6)] = FMA(KP831469612, Tin, Ti8);
					ri[WS(rs, 38)] = FNMS(KP831469612, Tin, Ti8);
					TkY = Tip + Tiq;
					Tir = Tip - Tiq;
					ii[WS(rs, 38)] = FNMS(KP831469612, TkY, TkX);
					ii[WS(rs, 6)] = FMA(KP831469612, TkY, TkX);
					ri[WS(rs, 22)] = FMA(KP831469612, Tir, Tio);
					ri[WS(rs, 54)] = FNMS(KP831469612, Tir, Tio);
				   }
				   {
					E TiE, Tl3, Tl4, TiH, Tiw, TiD;
					TiE = FMA(KP923879532, Tiv, Tis);
					Tiw = FNMS(KP923879532, Tiv, Tis);
					TiD = Tiz - TiC;
					Tl6 = Tiz + TiC;
					Tl5 = FMA(KP923879532, Tl2, Tl1);
					Tl3 = FNMS(KP923879532, Tl2, Tl1);
					ii[WS(rs, 54)] = FNMS(KP831469612, Tl0, TkZ);
					ii[WS(rs, 22)] = FMA(KP831469612, Tl0, TkZ);
					ri[WS(rs, 14)] = FMA(KP980785280, TiD, Tiw);
					ri[WS(rs, 46)] = FNMS(KP980785280, TiD, Tiw);
					Tl4 = TiG - TiF;
					TiH = TiF + TiG;
					ii[WS(rs, 46)] = FNMS(KP980785280, Tl4, Tl3);
					ii[WS(rs, 14)] = FMA(KP980785280, Tl4, Tl3);
					ri[WS(rs, 62)] = FMA(KP980785280, TiH, TiE);
					ri[WS(rs, 30)] = FNMS(KP980785280, TiH, TiE);
				   }
			      }
			 }
			 {
			      E Tla, TdV, TdO, Tm6, Tm5, TdR;
			      {
				   E TcT, TlO, TlI, Tar, TcX, Td3, TcN, TbB, TdM, TdQ, TdA, Tdw, TdJ, TdP, Tdz;
				   E Tdp, TlW, TdF, Tm2, Tdh, Td7, T91, Td6, T8M, TlT, TlF, Td0, Td4, TcO, TcK;
				   E T9g, Td8;
				   {
					E Tdb, Tdc, Tde, Tdf, Tdm, Tdk, Tdj, Tdn, TcF, Tct, TbY, Tdt, TdL, Tds, Tdu;
					E TcI, TdD, Tdd;
					{
					     E Tae, TcR, T9R, Tap, T9F, T9Q;
					     Tdb = FMA(KP707106781, T9E, T9p);
					     T9F = FNMS(KP707106781, T9E, T9p);
					     T9Q = FNMS(KP707106781, T9P, T9M);
					     Tdc = FMA(KP707106781, T9P, T9M);
					     Tde = FMA(KP707106781, Tad, T9Y);
					     Tae = FNMS(KP707106781, Tad, T9Y);
					     ii[WS(rs, 62)] = FMA(KP980785280, Tl6, Tl5);
					     ii[WS(rs, 30)] = FNMS(KP980785280, Tl6, Tl5);
					     TcR = FMA(KP668178637, T9F, T9Q);
					     T9R = FNMS(KP668178637, T9Q, T9F);
					     Tap = FNMS(KP707106781, Tao, Tal);
					     Tdf = FMA(KP707106781, Tao, Tal);
					     {
						  E Tbw, TcW, Tbl, Tbz;
						  {
						       E TaP, Tbk, TcS, Taq;
						       Tdm = FMA(KP707106781, TaO, Taz);
						       TaP = FNMS(KP707106781, TaO, Taz);
						       Tbk = Tb4 - Tbj;