hf2_64.c

最新的FFT程序

					     rio[WS(ios, 62)] = FMS(KP980785280, Tl6, Tl5);
					     TcR = FNMS(KP668178637, T9F, T9Q);
					     T9R = FMA(KP668178637, T9Q, T9F);
					     Tap = FNMS(KP707106781, Tao, Tal);
					     Tde = FMA(KP707106781, Tao, Tal);
					     {
						  E Tbw, TcV, Tbl, Tbz;
						  {
						       E TaP, Tbk, TcS, Taq;
						       Tdj = FMA(KP707106781, TaO, Taz);
						       TaP = FNMS(KP707106781, TaO, Taz);
						       Tbk = Tb4 - Tbj;
						       Tdn = Tb4 + Tbj;
						       Tdm = FMA(KP707106781, Tbv, Tbs);
						       Tbw = FNMS(KP707106781, Tbv, Tbs);
						       TcS = FMA(KP668178637, Tae, Tap);
						       Taq = FNMS(KP668178637, Tap, Tae);
						       TcV = FNMS(KP923879532, Tbk, TaP);
						       Tbl = FMA(KP923879532, Tbk, TaP);
						       TlI = TcR + TcS;
						       TcT = TcR - TcS;
						       TlO = Taq - T9R;
						       Tar = T9R + Taq;
						       Tbz = Tbx - Tby;
						       Tdk = Tby + Tbx;
						  }
						  {
						       E Tdq, Tdr, TcW, TbA;
						       TbY = FNMS(KP707106781, TbX, TbI);
						       Tdq = FMA(KP707106781, TbX, TbI);
						       Tdr = TcH + TcG;
						       TcI = TcG - TcH;
						       TcF = FNMS(KP707106781, TcE, TcB);
						       Tdt = FMA(KP707106781, TcE, TcB);
						       TcW = FNMS(KP923879532, Tbz, Tbw);
						       TbA = FMA(KP923879532, Tbz, Tbw);
						       TdL = FMA(KP923879532, Tdr, Tdq);
						       Tds = FNMS(KP923879532, Tdr, Tdq);
						       TcX = FNMS(KP534511135, TcW, TcV);
						       Td3 = FMA(KP534511135, TcV, TcW);
						       TcN = FNMS(KP303346683, Tbl, TbA);
						       TbB = FMA(KP303346683, TbA, Tbl);
						       Tdu = Tcs + Tcd;
						       Tct = Tcd - Tcs;
						  }
					     }
					}
					{
					     E TdI, Tdl, TdK, Tdv, TdH, Tdo;
					     TdK = FMA(KP923879532, Tdu, Tdt);
					     Tdv = FNMS(KP923879532, Tdu, Tdt);
					     TdI = FMA(KP923879532, Tdk, Tdj);
					     Tdl = FNMS(KP923879532, Tdk, Tdj);
					     TdM = FMA(KP098491403, TdL, TdK);
					     TdQ = FNMS(KP098491403, TdK, TdL);
					     TdA = FNMS(KP820678790, Tds, Tdv);
					     Tdw = FMA(KP820678790, Tdv, Tds);
					     TdH = FMA(KP923879532, Tdn, Tdm);
					     Tdo = FNMS(KP923879532, Tdn, Tdm);
					     TdD = FNMS(KP198912367, Tdb, Tdc);
					     Tdd = FMA(KP198912367, Tdc, Tdb);
					     TdJ = FMA(KP098491403, TdI, TdH);
					     TdP = FNMS(KP098491403, TdH, TdI);
					     Tdz = FNMS(KP820678790, Tdl, Tdo);
					     Tdp = FMA(KP820678790, Tdo, Tdl);
					}
					{
					     E TcY, Tcu, TdE, Tdg;
					     TdE = FMA(KP198912367, Tde, Tdf);
					     Tdg = FNMS(KP198912367, Tdf, Tde);
					     TcY = FNMS(KP923879532, Tct, TbY);
					     Tcu = FMA(KP923879532, Tct, TbY);
					     TlW = TdE - TdD;
					     TdF = TdD + TdE;
					     Tm2 = Tdd + Tdg;
					     Tdh = Tdd - Tdg;
					     {
						  E T8L, TlE, TcZ, TcJ;
						  Tla = T8D + T8K;
						  T8L = T8D - T8K;
						  TlE = TdU - TdT;
						  TdV = TdT + TdU;
						  Td7 = FNMS(KP414213562, T8T, T90);
						  T91 = FMA(KP414213562, T90, T8T);
						  TcZ = FMA(KP923879532, TcI, TcF);
						  TcJ = FNMS(KP923879532, TcI, TcF);
						  Td6 = FNMS(KP707106781, T8L, T8w);
						  T8M = FMA(KP707106781, T8L, T8w);
						  TlT = FNMS(KP707106781, TlE, TlD);
						  TlF = FMA(KP707106781, TlE, TlD);
						  Td0 = FNMS(KP534511135, TcZ, TcY);
						  Td4 = FMA(KP534511135, TcY, TcZ);
						  TcO = FNMS(KP303346683, Tcu, TcJ);
						  TcK = FMA(KP303346683, TcJ, Tcu);
						  T9g = FNMS(KP414213562, T9f, T98);
						  Td8 = FMA(KP414213562, T98, T9f);
					     }
					}
				   }
				   {
					E Tm1, TlV, TdC, Tda, Td2, TlQ, TlP, Td5;
					{
					     E TlM, TcQ, TlN, TcM, TlL, TcP;
					     {
						  E TcL, Tas, TlJ, TlK, TlH;
						  TlM = TcK - TbB;
						  TcL = TbB + TcK;
						  {
						       E TlU, T9h, TlG, Td9, T9i;
						       TlU = T91 + T9g;
						       T9h = T91 - T9g;
						       TlG = Td8 - Td7;
						       Td9 = Td7 + Td8;
						       Tm1 = FMA(KP923879532, TlU, TlT);
						       TlV = FNMS(KP923879532, TlU, TlT);
						       TcQ = FNMS(KP923879532, T9h, T8M);
						       T9i = FMA(KP923879532, T9h, T8M);
						       TlN = FNMS(KP923879532, TlG, TlF);
						       TlH = FMA(KP923879532, TlG, TlF);
						       TdC = FMA(KP923879532, Td9, Td6);
						       Tda = FNMS(KP923879532, Td9, Td6);
						       Tas = FMA(KP831469612, Tar, T9i);
						       TcM = FNMS(KP831469612, Tar, T9i);
						  }
						  TlL = FNMS(KP831469612, TlI, TlH);
						  TlJ = FMA(KP831469612, TlI, TlH);
						  TlK = TcN - TcO;
						  TcP = TcN + TcO;
						  rio[WS(ios, 3)] = FMA(KP956940335, TcL, Tas);
						  iio[-WS(ios, 35)] = FNMS(KP956940335, TcL, Tas);
						  iio[-WS(ios, 3)] = FMA(KP956940335, TlK, TlJ);
						  rio[WS(ios, 35)] = FMS(KP956940335, TlK, TlJ);
					     }
					     {
						  E TcU, Td1, TlR, TlS;
						  Td2 = FMA(KP831469612, TcT, TcQ);
						  TcU = FNMS(KP831469612, TcT, TcQ);
						  rio[WS(ios, 19)] = FMA(KP956940335, TcP, TcM);
						  iio[-WS(ios, 51)] = FNMS(KP956940335, TcP, TcM);
						  iio[-WS(ios, 19)] = FMA(KP956940335, TlM, TlL);
						  rio[WS(ios, 51)] = FMS(KP956940335, TlM, TlL);
						  Td1 = TcX + Td0;
						  TlQ = Td0 - TcX;
						  TlP = FMA(KP831469612, TlO, TlN);
						  TlR = FNMS(KP831469612, TlO, TlN);
						  TlS = Td4 - Td3;
						  Td5 = Td3 + Td4;
						  iio[-WS(ios, 59)] = FMA(KP881921264, Td1, TcU);
						  rio[WS(ios, 27)] = FNMS(KP881921264, Td1, TcU);
						  iio[-WS(ios, 27)] = FMA(KP881921264, TlS, TlR);
						  rio[WS(ios, 59)] = FMS(KP881921264, TlS, TlR);
					     }
					}
					{
					     E Tdy, Tm0, TlZ, TdB;
					     {
						  E Tdi, Tdx, TlX, TlY;
						  Tdy = FNMS(KP980785280, Tdh, Tda);
						  Tdi = FMA(KP980785280, Tdh, Tda);
						  rio[WS(ios, 11)] = FMA(KP881921264, Td5, Td2);
						  iio[-WS(ios, 43)] = FNMS(KP881921264, Td5, Td2);
						  iio[-WS(ios, 11)] = FMA(KP881921264, TlQ, TlP);
						  rio[WS(ios, 43)] = FMS(KP881921264, TlQ, TlP);
						  Tdx = Tdp + Tdw;
						  Tm0 = Tdw - Tdp;
						  TlZ = FNMS(KP980785280, TlW, TlV);
						  TlX = FMA(KP980785280, TlW, TlV);
						  TlY = Tdz - TdA;
						  TdB = Tdz + TdA;
						  rio[WS(ios, 7)] = FMA(KP773010453, Tdx, Tdi);
						  iio[-WS(ios, 39)] = FNMS(KP773010453, Tdx, Tdi);
						  iio[-WS(ios, 7)] = FMA(KP773010453, TlY, TlX);
						  rio[WS(ios, 39)] = FMS(KP773010453, TlY, TlX);
					     }
					     {
						  E TdG, TdN, Tm3, Tm4;
						  TdO = FMA(KP980785280, TdF, TdC);
						  TdG = FNMS(KP980785280, TdF, TdC);
						  rio[WS(ios, 23)] = FMA(KP773010453, TdB, Tdy);
						  iio[-WS(ios, 55)] = FNMS(KP773010453, TdB, Tdy);
						  iio[-WS(ios, 23)] = FMA(KP773010453, Tm0, TlZ);
						  rio[WS(ios, 55)] = FMS(KP773010453, Tm0, TlZ);
						  TdN = TdJ + TdM;
						  Tm6 = TdM - TdJ;
						  Tm5 = FMA(KP980785280, Tm2, Tm1);
						  Tm3 = FNMS(KP980785280, Tm2, Tm1);
						  Tm4 = TdQ - TdP;
						  TdR = TdP + TdQ;
						  rio[WS(ios, 15)] = FMA(KP995184726, TdN, TdG);
						  iio[-WS(ios, 47)] = FNMS(KP995184726, TdN, TdG);
						  iio[-WS(ios, 15)] = FMA(KP995184726, Tm4, Tm3);
						  rio[WS(ios, 47)] = FMS(KP995184726, Tm4, Tm3);
					     }
					}
				   }
			      }
			      {
				   E Tle, Tf5, Tlk, Tej, Tf9, Tff, TeZ, TeD, TfY, Tg2, TfM, TfI, TfV, Tg1, TfL;
				   E TfB, Tls, TfR, Tly, Tft, Tfj, TdZ, Tfi, TdW, Tlp, Tlb, Tfc, Tfg, Tf0, TeW;
				   E Te2, Tfk;
				   {
					E Tfo, Tfn, Tfr, Tfq, Tfv, Tfz, Tfy, Tfw, TeG, TeU, TeR, TfF, TfX, TfE, TfG;
					E TeN, TfP, Tfp;
					{
					     E Te7, Tea, Tee, Teh;
					     Tfo = FNMS(KP707106781, Te6, Te5);
					     Te7 = FMA(KP707106781, Te6, Te5);
					     iio[-WS(ios, 63)] = FMA(KP995184726, TdR, TdO);
					     rio[WS(ios, 31)] = FNMS(KP995184726, TdR, TdO);
					     iio[-WS(ios, 31)] = FMA(KP995184726, Tm6, Tm5);
					     rio[WS(ios, 63)] = FMS(KP995184726, Tm6, Tm5);
					     Tea = FMA(KP707106781, Te9, Te8);
					     Tfn = FNMS(KP707106781, Te9, Te8);
					     Tfr = FNMS(KP707106781, Ted, Tec);
					     Tee = FMA(KP707106781, Ted, Tec);
					     Teh = FMA(KP707106781, Teg, Tef);
					     Tfq = FNMS(KP707106781, Teg, Tef);
					     {
						  E Tey, Tf7, Tev, TeB;
						  {
						       E Ten, Tf3, Teb, Tf4, Tei, Teu;
						       Tfv = FNMS(KP707106781, Tem, Tel);
						       Ten = FMA(KP707106781, Tem, Tel);
						       Tf3 = FNMS(KP198912367, Te7, Tea);
						       Teb = FMA(KP198912367, Tea, Te7);
						       Tf4 = FMA(KP198912367, Tee, Teh);
						       Tei = FNMS(KP198912367, Teh, Tee);
						       Teu = Teq + Tet;
						       Tfz = Teq - Tet;
						       Tfy = FNMS(KP707106781, Tex, Tew);
						       Tey = FMA(KP707106781, Tex, Tew);
						       Tle = Tf3 + Tf4;
						       Tf5 = Tf3 - Tf4;
						       Tlk = Tei - Teb;
						       Tej = Teb + Tei;
						       Tf7 = FNMS(KP923879532, Teu, Ten);
						       Tev = FMA(KP923879532, Teu, Ten);
						       TeB = Tez + TeA;
						       Tfw = TeA - Tez;
						  }
						  {
						       E TfC, TfD, Tf8, TeC;
						       TeG = FMA(KP707106781, TeF, TeE);
						       TfC = FNMS(KP707106781, TeF, TeE);
						       TfD = TeT - TeS;
						       TeU = TeS + TeT;
						       TeR = FMA(KP707106781, TeQ, TeP);
						       TfF = FNMS(KP707106781, TeQ, TeP);
						       Tf8 = FNMS(KP923879532, TeB, Tey);
						       TeC = FMA(KP923879532, TeB, Tey);
						       TfX = FMA(KP923879532, TfD, TfC);
						       TfE = FNMS(KP923879532, TfD, TfC);
						       Tf9 = FNMS(KP820678790, Tf8, Tf7);
						       Tff = FMA(KP820678790, Tf7, Tf8);
						       TeZ = FNMS(KP098491403, Tev, TeC);
						       TeD = FMA(KP098491403, TeC, Tev);
						       TfG = TeM - TeJ;
						       TeN = TeJ + TeM;
						  }
					     }
					}
					{
					     E TfU, Tfx, TfW, TfH, TfT, TfA;
					     TfW = FNMS(KP923879532, TfG, TfF);
					     TfH = FMA(KP923879532, TfG, TfF);
					     TfU = FMA(KP923879532, Tfw, Tfv);
					     Tfx = FNMS(KP923879532, Tfw, Tfv);
					     TfY = FNMS(KP303346683, TfX, TfW);
					     Tg2 = FMA(KP303346683, TfW, TfX);
					     TfM = FMA(KP534511135, TfE, TfH);
					     TfI = FNMS(KP534511135, TfH, TfE);
					     TfT = FMA(KP923879532, Tfz, Tfy);
					     TfA = FNMS(KP923879532, Tfz, Tfy);
					     TfP = FNMS(KP668178637, Tfn, Tfo);
					     Tfp = FMA(KP668178637, Tfo, Tfn);
					     TfV = FMA(KP303346683, TfU, TfT);
					     Tg1 = FNMS(KP303346683, TfT, TfU);
					     TfL = FNMS(KP534511135, Tfx, TfA);
					     TfB = FMA(KP534511135, TfA, Tfx);
					}
					{
					     E Tfa, TeO, TfQ, Tfs, Tfb, TeV;
					     TfQ = FMA(KP668178637, Tfq, Tfr);
					     Tfs = FNMS(KP668178637, Tfr, Tfq);
					     Tfa = FNMS(KP923879532, TeN, TeG);
					     TeO = FMA(KP923879532, TeN, TeG);
					     Tls = TfQ - TfP;
					     TfR = TfP + TfQ;
					     Tly = Tfp + Tfs;
					     Tft = Tfp - Tfs;
					     Tfj = FNMS(KP414213562, TdX, TdY);
					     TdZ = FMA(KP414213562, TdY, TdX);
					     Tfb = FNMS(KP923879532, TeU, TeR);
					     TeV = FMA(KP923879532, TeU, TeR);
					     Tfi = FNMS(KP707106781, TdV, TdS);
					     TdW = FMA(KP707106781, TdV, TdS);
					     Tlp = FNMS(KP707106781, Tla, Tl9);
					     Tlb = FMA(KP707106781, Tla, Tl9);
					     Tfc = FMA(KP820678790, Tfb, Tfa);
					     Tfg = FNMS(KP820678790, Tfa, Tfb);
					     Tf0 = FMA(KP098491403, TeO, TeV);
					     TeW = FNMS(KP098491403, TeV, TeO);
					     Te2 = FNMS(KP414213562, Te1, Te0);
					     Tfk = FMA(KP414213562, Te0, Te1);
					}
				   }
				   {
					E Tlx, Tlr, TfO, Tfm, Tfe, Tlm, Tll, Tfh;
					{
					     E Tli, Tf2, Tlj, TeY, Tlh, Tf1;
					     {
						  E TeX, Tek, Tlf, Tlg, Tld;
						  Tli = TeW - TeD;
						  TeX = TeD + TeW;
						  {
						       E Tlq, Te3, Tlc, Tfl, Te4;
						       Tlq = Te2 - TdZ;
						       Te3 = TdZ + Te2;
						       Tlc = Tfj + Tfk;
						       Tfl = Tfj - Tfk;
						       Tlx = FNMS(KP923879532, Tlq, Tlp);
						       Tlr = FMA(KP923879532, Tlq, Tlp);
						       Tf2 = FNMS(KP923879532, Te3, TdW);
						       Te4 = FMA(KP923879532, Te3, TdW);
						       Tlj = FNMS(KP923879532, Tlc, Tlb);
						       Tld = FMA(KP923879532, Tlc, Tlb);
						       TfO = FNMS(KP923879532, Tfl, Tfi);
						       Tfm = FMA(KP923879532, Tfl, Tfi);
						       Tek = FMA(KP980785280, Tej, Te4);
						       TeY = FNMS(KP980785280, Tej, Te4);
						  }
						  Tlh = FNMS(KP980785280, Tle, Tld);
						  Tlf = FMA(KP980785280, Tle, Tld);
						  Tlg = TeZ + Tf0;
						  Tf1 = TeZ - Tf0;
						  rio[WS(ios, 1)] = FMA(KP995184726, TeX, Tek);
						  iio[-WS(ios, 33)] = FNMS(KP995184726, TeX, Tek);
						  iio[-WS(ios, 1)] = FMA(KP995184726, Tlg, Tlf);
						  rio[WS(ios, 33)] = FMS(KP995184726, Tlg, Tlf);
					     }
					     {
						  E Tf6, Tfd, Tln, Tlo;
						  Tfe = FMA(KP980785280, Tf5, Tf2);
						  Tf6 = FNMS(KP980785280, Tf5, Tf2);
						  rio[WS(ios, 17)] = FMA(KP995184726, Tf1, TeY);
						  iio[-WS(ios, 49)] = FNMS(KP995184726, Tf1, TeY);
						  iio[-WS(ios, 17)] = FMA(KP995184726, Tli, Tlh);
						  rio[WS(ios, 49)] = FMS(KP995184726, Tli, Tlh);
						  Tfd = Tf9 + Tfc;
						  Tlm = Tfc - Tf9;
						  Tll = FMA(KP980785280, Tlk, Tlj);
						  Tln = FNMS(KP980785280, Tlk, Tlj);
						  Tlo = Tff + Tfg;
						  Tfh = Tff - Tfg;
						  iio[-WS(ios, 57)] = FMA(KP773010453, Tfd, Tf6);
						  rio[WS(ios, 25)] = FNMS(KP773010453, Tfd, Tf6);
						  iio[-WS(ios, 25)] = FNMS(KP773010453, Tlo, Tln);
						  rio[WS(ios, 57)] = -(FMA(KP773010453, Tlo, Tln));
					     }
					}
					{
					     E TfK, Tlw, Tlv, TfN;
					     {
						  E Tfu, TfJ, Tlt, Tlu;
						  TfK = FNMS(KP831469612, Tft, Tfm);
						  Tfu = FMA(KP831469612, Tft, Tfm);
						  rio[WS(ios, 9)] = FMA(KP773010453, Tfh, Tfe);
						  iio[-WS(ios, 41)] = FNMS(KP773010453, Tfh, Tfe);
						  iio[-WS(ios, 9)] = FMA(KP773010453, Tlm, Tll);
						  rio[WS(ios, 41)] = FMS(KP773010453, Tlm, Tll);
						  TfJ = TfB + TfI;
						  Tlw = TfI - TfB;
						  Tlv = FNMS(KP831469612, Tls, Tlr);
						  Tlt = FMA(KP831469612, Tls, Tlr);
						  Tlu = TfL + TfM;
						  TfN = TfL - TfM;
						  rio[WS(ios, 5)] = FMA(KP881921264, TfJ, Tfu);
						  iio[-WS(ios, 37)] = FNMS(KP881921264, TfJ, Tfu);
						  iio[-WS(ios, 5)] = FMA(KP881921264, Tlu, Tlt);
						  rio[WS(ios, 37)] = FMS(KP881921264, Tlu, Tlt);
					     }
					     {
						  E TfS, TfZ, Tlz, TlA;
						  Tg0 = FMA(KP831469612, TfR, TfO);
						  TfS = FNMS(KP831469612, TfR, TfO);
						  rio[WS(ios, 21)] = FMA(KP881921264, TfN, TfK);
						  iio[-WS(ios, 53)] = FNMS(KP881921264, TfN, TfK);
						  iio[-WS(ios, 21)] = FMA(KP881921264, Tlw, Tlv);
						  rio[WS(ios, 53)] = FMS(KP881921264, Tlw, Tlv);
						  TfZ = TfV - TfY;
						  TlC = TfV + TfY;
						  TlB = FMA(KP831469612, Tly, Tlx);
						  Tlz = FNMS(KP831469612, Tly, Tlx);
						  TlA = Tg2 - Tg1;
						  Tg3 = Tg1 + Tg2;
						  rio[WS(ios, 13)] = FMA(KP956940335, TfZ, TfS);
						  iio[-WS(ios, 45)] = FNMS(KP956940335, TfZ, TfS);
						  iio[-WS(ios, 13)] = FMA(KP956940335, TlA, Tlz);
						  rio[WS(ios, 45)] = FMS(KP956940335, TlA, Tlz);
					     }
					}
				   }
			      }
			 }
		    }
	       }
	  }
	  iio[-WS(ios, 61)] = FMA(KP956940335, Tg3, Tg0);
	  rio[WS(ios, 29)] = FNMS(KP956940335, Tg3, Tg0);
	  iio[-WS(ios, 29)] = FNMS(KP956940335, TlC, TlB);
	  rio[WS(ios, 61)] = -(FMA(KP956940335, TlC, TlB));
     }
     return W;
}

static const tw_instr twinstr[] = {