hf2_64.c

最新的FFT程序

			      E TiI, Tkp, TiQ, TiS, TiL, Tkq, TiP, TiV, Tjf, Tja, Tjd, Tji, Tj4, Tj2, Tj1;
			      E Tj7, Tk0, Tkh, Tk3, Tki;
			      {
				   E TjG, T2I, Tkj, T4N, Tkk, Tkf, Tk5, TjJ, T8o, Tk2, TjU, TjY, T6D, TjL, TjO;
				   E Tk1;
				   {
					E Tk6, Tke, TjH, TjI;
					{
					     E T1C, T2H, T3L, T4M;
					     TiI = TY - T1B;
					     T1C = TY + T1B;
					     Th6 = TgT + TgW;
					     TgX = TgT - TgW;
					     Tbx = FMA(KP414213562, Tbb, Tbi);
					     Tbj = FNMS(KP414213562, Tbi, Tbb);
					     Tet = FNMS(KP414213562, Tes, Ter);
					     TeA = FMA(KP414213562, Ter, Tes);
					     T2H = T27 + T2G;
					     Tkp = T2G - T27;
					     TiQ = T39 - T3K;
					     T3L = T39 + T3K;
					     T4M = T4k + T4L;
					     TiS = T4k - T4L;
					     TiL = TiJ - TiK;
					     Tk6 = TiJ + TiK;
					     TjG = T1C - T2H;
					     T2I = T1C + T2H;
					     Tkj = T4M - T3L;
					     T4N = T3L + T4M;
					     Tke = Tk7 + Tkd;
					     Tkq = Tkd - Tk7;
					}
					TiP = TiN - TiO;
					TjH = TiN + TiO;
					TjI = TiT + TiU;
					TiV = TiT - TiU;
					{
					     E T7w, T8n, TjR, TjS;
					     Tjf = T74 - T7v;
					     T7w = T74 + T7v;
					     Tkk = Tke - Tk6;
					     Tkf = Tk6 + Tke;
					     Tk5 = TjH + TjI;
					     TjJ = TjH - TjI;
					     T8n = T7V + T8m;
					     Tja = T8m - T7V;
					     Tjd = Tjb - Tjc;
					     TjR = Tjb + Tjc;
					     TjS = Tjg + Tjh;
					     Tji = Tjg - Tjh;
					     {
						  E T5J, TjQ, TjT, T6C, TjM, TjN;
						  Tj4 = T5d - T5I;
						  T5J = T5d + T5I;
						  T8o = T7w + T8n;
						  TjQ = T7w - T8n;
						  Tk2 = TjR + TjS;
						  TjT = TjR - TjS;
						  T6C = T68 + T6B;
						  Tj2 = T6B - T68;
						  Tj1 = TiZ - Tj0;
						  TjM = TiZ + Tj0;
						  TjN = Tj5 + Tj6;
						  Tj7 = Tj5 - Tj6;
						  TjU = TjQ - TjT;
						  TjY = TjQ + TjT;
						  T6D = T5J + T6C;
						  TjL = T5J - T6C;
						  TjO = TjM - TjN;
						  Tk1 = TjM + TjN;
					     }
					}
				   }
				   {
					E TjW, TjK, TjV, Tko, Tkn, Tkl, Tkm, TjZ;
					{
					     E T4O, TjX, TjP, T8p, Tk4, Tkg;
					     Tk0 = T2I - T4N;
					     T4O = T2I + T4N;
					     TjX = TjO - TjL;
					     TjP = TjL + TjO;
					     T8p = T6D + T8o;
					     Tkh = T8o - T6D;
					     Tk3 = Tk1 - Tk2;
					     Tk4 = Tk1 + Tk2;
					     Tkg = Tk5 + Tkf;
					     Tki = Tkf - Tk5;
					     TjW = TjG - TjJ;
					     TjK = TjG + TjJ;
					     rio[0] = T4O + T8p;
					     iio[-WS(ios, 32)] = T4O - T8p;
					     iio[0] = Tk4 + Tkg;
					     rio[WS(ios, 32)] = Tk4 - Tkg;
					     TjV = TjP + TjU;
					     Tko = TjU - TjP;
					     Tkn = Tkk - Tkj;
					     Tkl = Tkj + Tkk;
					     Tkm = TjX + TjY;
					     TjZ = TjX - TjY;
					}
					rio[WS(ios, 8)] = FMA(KP707106781, TjV, TjK);
					iio[-WS(ios, 40)] = FNMS(KP707106781, TjV, TjK);
					iio[-WS(ios, 8)] = FMA(KP707106781, Tkm, Tkl);
					rio[WS(ios, 40)] = FMS(KP707106781, Tkm, Tkl);
					rio[WS(ios, 24)] = FMA(KP707106781, TjZ, TjW);
					iio[-WS(ios, 56)] = FNMS(KP707106781, TjZ, TjW);
					iio[-WS(ios, 24)] = FMA(KP707106781, Tko, Tkn);
					rio[WS(ios, 56)] = FMS(KP707106781, Tko, Tkn);
				   }
			      }
			      {
				   E Tjq, TiM, Tkx, Tkr, Tjt, Tky, Tks, TiX, Tjz, Tje, Tjx, TjD, Tjn, Tj9;
				   rio[WS(ios, 16)] = Tk0 + Tk3;
				   iio[-WS(ios, 48)] = Tk0 - Tk3;
				   iio[-WS(ios, 16)] = Tkh + Tki;
				   rio[WS(ios, 48)] = Tkh - Tki;
				   Tjq = TiI + TiL;
				   TiM = TiI - TiL;
				   Tkx = Tkq - Tkp;
				   Tkr = Tkp + Tkq;
				   {
					E Tjw, Tj3, Tjv, Tj8;
					{
					     E Tjr, TiR, Tjs, TiW;
					     Tjr = TiQ + TiP;
					     TiR = TiP - TiQ;
					     Tjs = TiS - TiV;
					     TiW = TiS + TiV;
					     Tjt = Tjr + Tjs;
					     Tky = Tjs - Tjr;
					     Tks = TiR + TiW;
					     TiX = TiR - TiW;
					     Tjw = Tj1 + Tj2;
					     Tj3 = Tj1 - Tj2;
					}
					Tjv = Tj4 + Tj7;
					Tj8 = Tj4 - Tj7;
					Tjz = Tjd + Tja;
					Tje = Tja - 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 = FNMS(KP414213562, Tje, Tjj);
						  Tjk = FMA(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 = Tjk - Tj9;
						       Tjl = Tj9 + Tjk;
						       rio[WS(ios, 4)] = FMA(KP923879532, TjB, Tju);
						       iio[-WS(ios, 36)] = FNMS(KP923879532, TjB, Tju);
						       iio[-WS(ios, 60)] = FMA(KP923879532, Tjp, Tjm);
						       rio[WS(ios, 28)] = FNMS(KP923879532, Tjp, Tjm);
						       iio[-WS(ios, 12)] = FMA(KP923879532, TkA, Tkz);
						       rio[WS(ios, 44)] = FMS(KP923879532, TkA, Tkz);
						       iio[-WS(ios, 28)] = FMA(KP923879532, TkC, TkB);
						       rio[WS(ios, 60)] = FMS(KP923879532, TkC, TkB);
						       rio[WS(ios, 12)] = FMA(KP923879532, Tjl, TiY);
						       iio[-WS(ios, 44)] = FNMS(KP923879532, Tjl, TiY);
						       Tku = TjD + TjE;
						       TjF = TjD - TjE;
						  }
					     }
					     iio[-WS(ios, 4)] = FMA(KP923879532, Tku, Tkt);
					     rio[WS(ios, 36)] = FMS(KP923879532, Tku, Tkt);
					     rio[WS(ios, 20)] = FMA(KP923879532, TjF, TjC);
					     iio[-WS(ios, 52)] = FNMS(KP923879532, TjF, TjC);
					}
				   }
			      }
			 }
			 {
			      E TkV, Tl1, ThG, Tgk, TkH, TkN, Tis, Ti0, Thy, Thv, TkI, ThJ, TkO, TgH, TiC;
			      E TiG, Tiq, Tim, ThN, ThT, ThD, Th9, TkW, Tiv, Tl2, Ti7, ThO, Thq, Tiz, TiF;
			      E Tip, Tif;
			      {
				   E Ti2, Ti1, Ti5, Ti4, The, Tij, TiB, Tii, Tik, Thp;
				   {
					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;
					iio[-WS(ios, 20)] = FMA(KP923879532, Tkw, Tkv);
					rio[WS(ios, 52)] = FMS(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;
						  Ti2 = Tgl - Tgo;
						  Tgp = Tgl + Tgo;
						  TkV = FMA(KP707106781, TkU, TkT);
						  Tl1 = FNMS(KP707106781, TkU, TkT);
						  ThG = FNMS(KP707106781, Tgj, Tg8);
						  Tgk = FMA(KP707106781, Tgj, Tg8);
						  TkG = ThX + ThY;
						  ThZ = ThX - ThY;
						  Tgu = Tgs + Tgt;
						  Ti1 = Tgs - Tgt;
						  Ti5 = Tgw - Tgz;
						  TgA = Tgw + Tgz;
						  TkH = FMA(KP707106781, TkG, TkF);
						  TkN = FNMS(KP707106781, TkG, TkF);
						  Tis = FNMS(KP707106781, ThZ, ThW);
						  Ti0 = FMA(KP707106781, ThZ, ThW);
						  ThH = FNMS(KP414213562, Tgp, Tgu);
						  Tgv = FMA(KP414213562, Tgu, Tgp);
						  TgF = TgD + TgE;
						  Ti4 = TgD - TgE;
					     }
					     {
						  E Tig, Tih, ThI, TgG;
						  The = Tha + Thd;
						  Tig = Tha - Thd;
						  Tih = Thx - Thw;
						  Thy = Thw + Thx;
						  Thv = Tht + Thu;
						  Tij = Thu - Tht;
						  ThI = FMA(KP414213562, TgA, TgF);
						  TgG = FNMS(KP414213562, TgF, TgA);
						  TiB = FMA(KP707106781, Tih, Tig);
						  Tii = FNMS(KP707106781, Tih, Tig);
						  TkI = ThH + ThI;
						  ThJ = ThH - ThI;
						  TkO = TgG - Tgv;
						  TgH = Tgv + TgG;
						  Tik = Tho - Thj;
						  Thp = Thj + Tho;
					     }
					}
				   }
				   {
					E Ti9, Tid, Tic, Tia, Tit, Ti3;
					{
					     E Th4, ThL, TgZ, Th7, ThM, Th8;
					     {
						  E TgN, TgY, TiA, Til;
						  Ti9 = TgJ - TgM;
						  TgN = TgJ + TgM;
						  TgY = TgS + TgX;
						  Tid = TgS - TgX;
						  Tic = Th2 - Th3;
						  Th4 = Th2 + Th3;
						  TiA = FMA(KP707106781, Tik, Tij);
						  Til = FNMS(KP707106781, Tik, Tij);
						  ThL = FNMS(KP707106781, TgY, TgN);
						  TgZ = FMA(KP707106781, TgY, TgN);
						  TiC = FMA(KP198912367, TiB, TiA);
						  TiG = FNMS(KP198912367, TiA, TiB);
						  Tiq = FNMS(KP668178637, Tii, Til);
						  Tim = FMA(KP668178637, Til, Tii);
						  Th7 = Th5 + Th6;
						  Tia = Th6 - Th5;
					     }
					     ThM = FNMS(KP707106781, Th7, Th4);
					     Th8 = FMA(KP707106781, Th7, Th4);
					     Tit = FNMS(KP414213562, Ti1, Ti2);
					     Ti3 = FMA(KP414213562, Ti2, Ti1);
					     ThN = FNMS(KP668178637, ThM, ThL);
					     ThT = FMA(KP668178637, ThL, ThM);
					     ThD = FNMS(KP198912367, TgZ, Th8);
					     Th9 = FMA(KP198912367, 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);
					     ThO = FNMS(KP707106781, Thp, The);
					     Thq = FMA(KP707106781, Thp, The);
					     Tiz = FMA(KP198912367, Tiy, Tix);
					     TiF = FNMS(KP198912367, Tix, Tiy);
					     Tip = FNMS(KP668178637, Tib, Tie);
					     Tif = FMA(KP668178637, Tie, Tib);
					}
				   }
			      }
			      {
				   E TkQ, TkP, Tl0, TkZ;
				   {
					E ThC, TgI, TkJ, TkL, ThP, Thz;
					ThC = FNMS(KP923879532, TgH, Tgk);
					TgI = FMA(KP923879532, TgH, Tgk);
					TkJ = FMA(KP923879532, TkI, TkH);
					TkL = FNMS(KP923879532, TkI, TkH);
					ThP = FNMS(KP707106781, Thy, Thv);
					Thz = FMA(KP707106781, Thy, Thv);
					{
					     E ThS, TkR, TkS, ThV;
					     {
						  E ThK, ThU, ThE, ThA, ThR, ThQ;
						  ThS = FMA(KP923879532, ThJ, ThG);
						  ThK = FNMS(KP923879532, ThJ, ThG);
						  ThQ = FMA(KP668178637, ThP, ThO);
						  ThU = FNMS(KP668178637, ThO, ThP);
						  ThE = FMA(KP198912367, Thq, Thz);
						  ThA = FNMS(KP198912367, Thz, Thq);
						  ThR = ThN + ThQ;
						  TkQ = ThQ - ThN;
						  TkP = FMA(KP923879532, TkO, TkN);
						  TkR = FNMS(KP923879532, TkO, TkN);
						  {
						       E ThF, TkK, TkM, ThB;
						       ThF = ThD - ThE;
						       TkK = ThD + ThE;
						       TkM = ThA - Th9;
						       ThB = Th9 + ThA;
						       iio[-WS(ios, 58)] = FMA(KP831469612, ThR, ThK);
						       rio[WS(ios, 26)] = FNMS(KP831469612, ThR, ThK);
						       rio[WS(ios, 18)] = FMA(KP980785280, ThF, ThC);
						       iio[-WS(ios, 50)] = FNMS(KP980785280, ThF, ThC);
						       iio[-WS(ios, 2)] = FMA(KP980785280, TkK, TkJ);
						       rio[WS(ios, 34)] = FMS(KP980785280, TkK, TkJ);
						       iio[-WS(ios, 18)] = FMA(KP980785280, TkM, TkL);
						       rio[WS(ios, 50)] = FMS(KP980785280, TkM, TkL);
						       rio[WS(ios, 2)] = FMA(KP980785280, ThB, TgI);
						       iio[-WS(ios, 34)] = FNMS(KP980785280, ThB, TgI);
						       TkS = ThT + ThU;
						       ThV = ThT - ThU;
						  }
					     }
					     iio[-WS(ios, 26)] = FNMS(KP831469612, TkS, TkR);
					     rio[WS(ios, 58)] = -(FMA(KP831469612, TkS, TkR));
					     rio[WS(ios, 10)] = FMA(KP831469612, ThV, ThS);
					     iio[-WS(ios, 42)] = FNMS(KP831469612, 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);
					iio[-WS(ios, 10)] = FMA(KP831469612, TkQ, TkP);
					rio[WS(ios, 42)] = FMS(KP831469612, TkQ, TkP);
					rio[WS(ios, 6)] = FMA(KP831469612, Tin, Ti8);
					iio[-WS(ios, 38)] = FNMS(KP831469612, Tin, Ti8);
					TkY = Tip - Tiq;
					Tir = Tip + Tiq;
					iio[-WS(ios, 6)] = FMA(KP831469612, TkY, TkX);
					rio[WS(ios, 38)] = FMS(KP831469612, TkY, TkX);
					rio[WS(ios, 22)] = FMA(KP831469612, Tir, Tio);
					iio[-WS(ios, 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 = TiC - Tiz;
					Tl5 = FMA(KP923879532, Tl2, Tl1);
					Tl3 = FNMS(KP923879532, Tl2, Tl1);
					iio[-WS(ios, 22)] = FMA(KP831469612, Tl0, TkZ);
					rio[WS(ios, 54)] = FMS(KP831469612, Tl0, TkZ);
					rio[WS(ios, 14)] = FMA(KP980785280, TiD, Tiw);
					iio[-WS(ios, 46)] = FNMS(KP980785280, TiD, Tiw);
					Tl4 = TiG - TiF;
					TiH = TiF + TiG;
					iio[-WS(ios, 14)] = FMA(KP980785280, Tl4, Tl3);
					rio[WS(ios, 46)] = FMS(KP980785280, Tl4, Tl3);
					iio[-WS(ios, 62)] = FMA(KP980785280, TiH, TiE);
					rio[WS(ios, 30)] = FNMS(KP980785280, TiH, TiE);
				   }
			      }
			 }
			 {
			      E Tla, TdV, TdO, Tm6, Tm5, TdR;
			      {
				   E TlI, TcT, TlO, 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 Tdc, Tdb, Tdf, Tde, Tdj, Tdn, Tdm, Tdk, TbY, TcI, TcF, Tdt, TdL, Tds, Tdu;
					E Tct, TdD, Tdd;
					{
					     E Tae, TcR, T9R, Tap, T9F, T9Q;
					     Tdc = FMA(KP707106781, T9E, T9p);
					     T9F = FNMS(KP707106781, T9E, T9p);
					     T9Q = FNMS(KP707106781, T9P, T9M);
					     Tdb = FMA(KP707106781, T9P, T9M);
					     Tdf = FMA(KP707106781, Tad, T9Y);
					     Tae = FNMS(KP707106781, Tad, T9Y);
					     iio[-WS(ios, 30)] = FMA(KP980785280, Tl6, Tl5);