t2_64.c

fftw-3.0.1

			 Tif = Thu + Thv;
			 ii[WS(ios, 4)] = Tif + Tik;
			 ii[WS(ios, 36)] = Tik - Tif;
			 Thw = Thu - Thv;
			 ri[WS(ios, 52)] = Tht - Thw;
			 ri[WS(ios, 20)] = Tht + Thw;
			 Til = Thr - Tho;
			 ii[WS(ios, 20)] = Til + Tim;
			 ii[WS(ios, 52)] = Tim - Til;
		    }
	       }
	       {
		    E Tez, Tft, TiU, TiW, Tf0, Tfu, Tfr, Tfv;
		    {
			 E Teb, Tey, TiS, TiT;
			 Teb = TdZ - Tea;
			 Tey = Tem - Tex;
			 Tez = Teb + Tey;
			 Tft = Teb - Tey;
			 TiS = Tfz - Tfy;
			 TiT = TiM - TiL;
			 TiU = TiS + TiT;
			 TiW = TiT - TiS;
		    }
		    {
			 E TeQ, TeZ, Tfh, Tfq;
			 TeQ = TeE - TeP;
			 TeZ = TeV - TeY;
			 Tf0 = FMA(KP980785280, TeQ, KP195090322 * TeZ);
			 Tfu = FNMS(KP980785280, TeZ, KP195090322 * TeQ);
			 Tfh = Tf5 - Tfg;
			 Tfq = Tfm - Tfp;
			 Tfr = FNMS(KP980785280, Tfq, KP195090322 * Tfh);
			 Tfv = FMA(KP195090322, Tfq, KP980785280 * Tfh);
		    }
		    {
			 E Tfs, TiR, Tfw, TiV;
			 Tfs = Tf0 + Tfr;
			 ri[WS(ios, 46)] = Tez - Tfs;
			 ri[WS(ios, 14)] = Tez + Tfs;
			 TiR = Tfu + Tfv;
			 ii[WS(ios, 14)] = TiR + TiU;
			 ii[WS(ios, 46)] = TiU - TiR;
			 Tfw = Tfu - Tfv;
			 ri[WS(ios, 62)] = Tft - Tfw;
			 ri[WS(ios, 30)] = Tft + Tfw;
			 TiV = Tfr - Tf0;
			 ii[WS(ios, 30)] = TiV + TiW;
			 ii[WS(ios, 62)] = TiW - TiV;
		    }
	       }
	       {
		    E TfZ, Tgf, TiG, TiI, Tg6, Tgg, Tgd, Tgh;
		    {
			 E TfR, TfY, TiE, TiF;
			 TfR = TfN - TfQ;
			 TfY = TfU - TfX;
			 TfZ = TfR + TfY;
			 Tgf = TfR - TfY;
			 TiE = Tgl - Tgk;
			 TiF = Tiy - Tiv;
			 TiG = TiE + TiF;
			 TiI = TiF - TiE;
		    }
		    {
			 E Tg2, Tg5, Tg9, Tgc;
			 Tg2 = Tg0 - Tg1;
			 Tg5 = Tg3 - Tg4;
			 Tg6 = FMA(KP555570233, Tg2, KP831469612 * Tg5);
			 Tgg = FNMS(KP831469612, Tg2, KP555570233 * Tg5);
			 Tg9 = Tg7 - Tg8;
			 Tgc = Tga - Tgb;
			 Tgd = FNMS(KP831469612, Tgc, KP555570233 * Tg9);
			 Tgh = FMA(KP831469612, Tg9, KP555570233 * Tgc);
		    }
		    {
			 E Tge, TiD, Tgi, TiH;
			 Tge = Tg6 + Tgd;
			 ri[WS(ios, 42)] = TfZ - Tge;
			 ri[WS(ios, 10)] = TfZ + Tge;
			 TiD = Tgg + Tgh;
			 ii[WS(ios, 10)] = TiD + TiG;
			 ii[WS(ios, 42)] = TiG - TiD;
			 Tgi = Tgg - Tgh;
			 ri[WS(ios, 58)] = Tgf - Tgi;
			 ri[WS(ios, 26)] = Tgf + Tgi;
			 TiH = Tgd - Tg6;
			 ii[WS(ios, 26)] = TiH + TiI;
			 ii[WS(ios, 58)] = TiI - TiH;
		    }
	       }
	       {
		    E Tgn, Tgv, TiA, TiC, Tgq, Tgw, Tgt, Tgx;
		    {
			 E Tgj, Tgm, Tiu, Tiz;
			 Tgj = TfN + TfQ;
			 Tgm = Tgk + Tgl;
			 Tgn = Tgj + Tgm;
			 Tgv = Tgj - Tgm;
			 Tiu = TfU + TfX;
			 Tiz = Tiv + Tiy;
			 TiA = Tiu + Tiz;
			 TiC = Tiz - Tiu;
		    }
		    {
			 E Tgo, Tgp, Tgr, Tgs;
			 Tgo = Tg0 + Tg1;
			 Tgp = Tg3 + Tg4;
			 Tgq = FMA(KP980785280, Tgo, KP195090322 * Tgp);
			 Tgw = FNMS(KP195090322, Tgo, KP980785280 * Tgp);
			 Tgr = Tg7 + Tg8;
			 Tgs = Tga + Tgb;
			 Tgt = FNMS(KP195090322, Tgs, KP980785280 * Tgr);
			 Tgx = FMA(KP195090322, Tgr, KP980785280 * Tgs);
		    }
		    {
			 E Tgu, Tit, Tgy, TiB;
			 Tgu = Tgq + Tgt;
			 ri[WS(ios, 34)] = Tgn - Tgu;
			 ri[WS(ios, 2)] = Tgn + Tgu;
			 Tit = Tgw + Tgx;
			 ii[WS(ios, 2)] = Tit + TiA;
			 ii[WS(ios, 34)] = TiA - Tit;
			 Tgy = Tgw - Tgx;
			 ri[WS(ios, 50)] = Tgv - Tgy;
			 ri[WS(ios, 18)] = Tgv + Tgy;
			 TiB = Tgt - Tgq;
			 ii[WS(ios, 18)] = TiB + TiC;
			 ii[WS(ios, 50)] = TiC - TiB;
		    }
	       }
	       {
		    E T7V, TjN, TjT, TaH, T8O, TjK, TaK, TjS, TaO, TaU, T9I, TaE, TaR, TaV, TaB;
		    E TaF, T8N;
		    T7V = T7x - T7U;
		    TjN = TjL + TjM;
		    TjT = TjM - TjL;
		    TaH = T7x + T7U;
		    T8N = FMA(KP195090322, T8D, KP980785280 * T8M);
		    T8O = T8m - T8N;
		    TjK = T8m + T8N;
		    {
			 E TaJ, TaM, TaN, T9u, T9H;
			 TaJ = FNMS(KP980785280, T8D, KP195090322 * T8M);
			 TaK = TaI + TaJ;
			 TjS = TaJ - TaI;
			 TaM = T96 + T9t;
			 TaN = T9D + T9G;
			 TaO = FMA(KP634393284, TaM, KP773010453 * TaN);
			 TaU = FNMS(KP634393284, TaN, KP773010453 * TaM);
			 T9u = T96 - T9t;
			 T9H = T9D - T9G;
			 T9I = FMA(KP995184726, T9u, KP098017140 * T9H);
			 TaE = FNMS(KP995184726, T9H, KP098017140 * T9u);
			 {
			      E TaP, TaQ, Tan, TaA;
			      TaP = T9Z + Tam;
			      TaQ = Taw + Taz;
			      TaR = FNMS(KP634393284, TaQ, KP773010453 * TaP);
			      TaV = FMA(KP773010453, TaQ, KP634393284 * TaP);
			      Tan = T9Z - Tam;
			      TaA = Taw - Taz;
			      TaB = FNMS(KP995184726, TaA, KP098017140 * Tan);
			      TaF = FMA(KP098017140, TaA, KP995184726 * Tan);
			 }
		    }
		    {
			 E T8P, TaC, TjR, TjU;
			 T8P = T7V + T8O;
			 TaC = T9I + TaB;
			 ri[WS(ios, 47)] = T8P - TaC;
			 ri[WS(ios, 15)] = T8P + TaC;
			 TjR = TaE + TaF;
			 TjU = TjS + TjT;
			 ii[WS(ios, 15)] = TjR + TjU;
			 ii[WS(ios, 47)] = TjU - TjR;
		    }
		    {
			 E TaD, TaG, TjV, TjW;
			 TaD = T7V - T8O;
			 TaG = TaE - TaF;
			 ri[WS(ios, 63)] = TaD - TaG;
			 ri[WS(ios, 31)] = TaD + TaG;
			 TjV = TaB - T9I;
			 TjW = TjT - TjS;
			 ii[WS(ios, 31)] = TjV + TjW;
			 ii[WS(ios, 63)] = TjW - TjV;
		    }
		    {
			 E TaL, TaS, TjJ, TjO;
			 TaL = TaH + TaK;
			 TaS = TaO + TaR;
			 ri[WS(ios, 39)] = TaL - TaS;
			 ri[WS(ios, 7)] = TaL + TaS;
			 TjJ = TaU + TaV;
			 TjO = TjK + TjN;
			 ii[WS(ios, 7)] = TjJ + TjO;
			 ii[WS(ios, 39)] = TjO - TjJ;
		    }
		    {
			 E TaT, TaW, TjP, TjQ;
			 TaT = TaH - TaK;
			 TaW = TaU - TaV;
			 ri[WS(ios, 55)] = TaT - TaW;
			 ri[WS(ios, 23)] = TaT + TaW;
			 TjP = TaR - TaO;
			 TjQ = TjN - TjK;
			 ii[WS(ios, 23)] = TjP + TjQ;
			 ii[WS(ios, 55)] = TjQ - TjP;
		    }
	       }
	       {
		    E TbV, Tjj, Tjp, TcT, Tca, Tjg, TcW, Tjo, Td0, Td6, Tcu, TcQ, Td3, Td7, TcN;
		    E TcR, Tc9;
		    TbV = TbN - TbU;
		    Tjj = Tjh + Tji;
		    Tjp = Tji - Tjh;
		    TcT = TbN + TbU;
		    Tc9 = FMA(KP831469612, Tc5, KP555570233 * Tc8);
		    Tca = Tc2 - Tc9;
		    Tjg = Tc2 + Tc9;
		    {
			 E TcV, TcY, TcZ, Tcm, Tct;
			 TcV = FNMS(KP831469612, Tc8, KP555570233 * Tc5);
			 TcW = TcU + TcV;
			 Tjo = TcV - TcU;
			 TcY = Tce + Tcl;
			 TcZ = Tcp + Tcs;
			 Td0 = FMA(KP471396736, TcY, KP881921264 * TcZ);
			 Td6 = FNMS(KP471396736, TcZ, KP881921264 * TcY);
			 Tcm = Tce - Tcl;
			 Tct = Tcp - Tcs;
			 Tcu = FMA(KP956940335, Tcm, KP290284677 * Tct);
			 TcQ = FNMS(KP956940335, Tct, KP290284677 * Tcm);
			 {
			      E Td1, Td2, TcF, TcM;
			      Td1 = Tcx + TcE;
			      Td2 = TcI + TcL;
			      Td3 = FNMS(KP471396736, Td2, KP881921264 * Td1);
			      Td7 = FMA(KP881921264, Td2, KP471396736 * Td1);
			      TcF = Tcx - TcE;
			      TcM = TcI - TcL;
			      TcN = FNMS(KP956940335, TcM, KP290284677 * TcF);
			      TcR = FMA(KP290284677, TcM, KP956940335 * TcF);
			 }
		    }
		    {
			 E Tcb, TcO, Tjn, Tjq;
			 Tcb = TbV + Tca;
			 TcO = Tcu + TcN;
			 ri[WS(ios, 45)] = Tcb - TcO;
			 ri[WS(ios, 13)] = Tcb + TcO;
			 Tjn = TcQ + TcR;
			 Tjq = Tjo + Tjp;
			 ii[WS(ios, 13)] = Tjn + Tjq;
			 ii[WS(ios, 45)] = Tjq - Tjn;
		    }
		    {
			 E TcP, TcS, Tjr, Tjs;
			 TcP = TbV - Tca;
			 TcS = TcQ - TcR;
			 ri[WS(ios, 61)] = TcP - TcS;
			 ri[WS(ios, 29)] = TcP + TcS;
			 Tjr = TcN - Tcu;
			 Tjs = Tjp - Tjo;
			 ii[WS(ios, 29)] = Tjr + Tjs;
			 ii[WS(ios, 61)] = Tjs - Tjr;
		    }
		    {
			 E TcX, Td4, Tjf, Tjk;
			 TcX = TcT + TcW;
			 Td4 = Td0 + Td3;
			 ri[WS(ios, 37)] = TcX - Td4;
			 ri[WS(ios, 5)] = TcX + Td4;
			 Tjf = Td6 + Td7;
			 Tjk = Tjg + Tjj;
			 ii[WS(ios, 5)] = Tjf + Tjk;
			 ii[WS(ios, 37)] = Tjk - Tjf;
		    }
		    {
			 E Td5, Td8, Tjl, Tjm;
			 Td5 = TcT - TcW;
			 Td8 = Td6 - Td7;
			 ri[WS(ios, 53)] = Td5 - Td8;
			 ri[WS(ios, 21)] = Td5 + Td8;
			 Tjl = Td3 - Td0;
			 Tjm = Tjj - Tjg;
			 ii[WS(ios, 21)] = Tjl + Tjm;
			 ii[WS(ios, 53)] = Tjm - Tjl;
		    }
	       }
	       {
		    E Tb1, Tjz, TjF, Tbt, Tb8, Tju, Tbw, TjE, TbA, TbG, Tbg, Tbq, TbD, TbH, Tbn;
		    E Tbr, Tb7;
		    Tb1 = TaX - Tb0;
		    Tjz = Tjv + Tjy;
		    TjF = Tjy - Tjv;
		    Tbt = TaX + Tb0;
		    Tb7 = FMA(KP831469612, Tb5, KP555570233 * Tb6);
		    Tb8 = Tb4 - Tb7;
		    Tju = Tb4 + Tb7;
		    {
			 E Tbv, Tby, Tbz, Tbc, Tbf;
			 Tbv = FNMS(KP555570233, Tb5, KP831469612 * Tb6);
			 Tbw = Tbu + Tbv;
			 TjE = Tbv - Tbu;
			 Tby = Tba + Tbb;
			 Tbz = Tbd + Tbe;
			 TbA = FMA(KP956940335, Tby, KP290284677 * Tbz);
			 TbG = FNMS(KP290284677, Tby, KP956940335 * Tbz);
			 Tbc = Tba - Tbb;
			 Tbf = Tbd - Tbe;
			 Tbg = FMA(KP471396736, Tbc, KP881921264 * Tbf);
			 Tbq = FNMS(KP881921264, Tbc, KP471396736 * Tbf);
			 {
			      E TbB, TbC, Tbj, Tbm;
			      TbB = Tbh + Tbi;
			      TbC = Tbk + Tbl;
			      TbD = FNMS(KP290284677, TbC, KP956940335 * TbB);
			      TbH = FMA(KP290284677, TbB, KP956940335 * TbC);
			      Tbj = Tbh - Tbi;
			      Tbm = Tbk - Tbl;
			      Tbn = FNMS(KP881921264, Tbm, KP471396736 * Tbj);
			      Tbr = FMA(KP881921264, Tbj, KP471396736 * Tbm);
			 }
		    }
		    {
			 E Tb9, Tbo, TjD, TjG;
			 Tb9 = Tb1 + Tb8;
			 Tbo = Tbg + Tbn;
			 ri[WS(ios, 43)] = Tb9 - Tbo;
			 ri[WS(ios, 11)] = Tb9 + Tbo;
			 TjD = Tbq + Tbr;
			 TjG = TjE + TjF;
			 ii[WS(ios, 11)] = TjD + TjG;
			 ii[WS(ios, 43)] = TjG - TjD;
		    }
		    {
			 E Tbp, Tbs, TjH, TjI;
			 Tbp = Tb1 - Tb8;
			 Tbs = Tbq - Tbr;
			 ri[WS(ios, 59)] = Tbp - Tbs;
			 ri[WS(ios, 27)] = Tbp + Tbs;
			 TjH = Tbn - Tbg;
			 TjI = TjF - TjE;
			 ii[WS(ios, 27)] = TjH + TjI;
			 ii[WS(ios, 59)] = TjI - TjH;
		    }
		    {
			 E Tbx, TbE, Tjt, TjA;
			 Tbx = Tbt + Tbw;
			 TbE = TbA + TbD;
			 ri[WS(ios, 35)] = Tbx - TbE;
			 ri[WS(ios, 3)] = Tbx + TbE;
			 Tjt = TbG + TbH;
			 TjA = Tju + Tjz;
			 ii[WS(ios, 3)] = Tjt + TjA;
			 ii[WS(ios, 35)] = TjA - Tjt;
		    }
		    {
			 E TbF, TbI, TjB, TjC;
			 TbF = Tbt - Tbw;
			 TbI = TbG - TbH;
			 ri[WS(ios, 51)] = TbF - TbI;
			 ri[WS(ios, 19)] = TbF + TbI;
			 TjB = TbD - TbA;
			 TjC = Tjz - Tju;
			 ii[WS(ios, 19)] = TjB + TjC;
			 ii[WS(ios, 51)] = TjC - TjB;
		    }
	       }
	       {
		    E Tdd, Tj5, Tjb, TdF, Tdk, TiY, TdI, Tja, TdM, TdS, Tds, TdC, TdP, TdT, Tdz;
		    E TdD, Tdj;
		    Tdd = Td9 - Tdc;
		    Tj5 = TiZ + Tj4;
		    Tjb = Tj4 - TiZ;
		    TdF = Td9 + Tdc;
		    Tdj = FMA(KP195090322, Tdh, KP980785280 * Tdi);
		    Tdk = Tdg - Tdj;
		    TiY = Tdg + Tdj;
		    {
			 E TdH, TdK, TdL, Tdo, Tdr;
			 TdH = FNMS(KP195090322, Tdi, KP980785280 * Tdh);
			 TdI = TdG + TdH;
			 Tja = TdH - TdG;
			 TdK = Tdm + Tdn;
			 TdL = Tdp + Tdq;
			 TdM = FMA(KP995184726, TdK, KP098017140 * TdL);
			 TdS = FNMS(KP098017140, TdK, KP995184726 * TdL);
			 Tdo = Tdm - Tdn;
			 Tdr = Tdp - Tdq;
			 Tds = FMA(KP634393284, Tdo, KP773010453 * Tdr);
			 TdC = FNMS(KP773010453, Tdo, KP634393284 * Tdr);
			 {
			      E TdN, TdO, Tdv, Tdy;
			      TdN = Tdt + Tdu;
			      TdO = Tdw + Tdx;
			      TdP = FNMS(KP098017140, TdO, KP995184726 * TdN);
			      TdT = FMA(KP098017140, TdN, KP995184726 * TdO);
			      Tdv = Tdt - Tdu;
			      Tdy = Tdw - Tdx;
			      Tdz = FNMS(KP773010453, Tdy, KP634393284 * Tdv);
			      TdD = FMA(KP773010453, Tdv, KP634393284 * Tdy);
			 }
		    }
		    {
			 E Tdl, TdA, Tj9, Tjc;
			 Tdl = Tdd + Tdk;
			 TdA = Tds + Tdz;
			 ri[WS(ios, 41)] = Tdl - TdA;
			 ri[WS(ios, 9)] = Tdl + TdA;
			 Tj9 = TdC + TdD;
			 Tjc = Tja + Tjb;
			 ii[WS(ios, 9)] = Tj9 + Tjc;
			 ii[WS(ios, 41)] = Tjc - Tj9;
		    }
		    {
			 E TdB, TdE, Tjd, Tje;
			 TdB = Tdd - Tdk;
			 TdE = TdC - TdD;
			 ri[WS(ios, 57)] = TdB - TdE;
			 ri[WS(ios, 25)] = TdB + TdE;
			 Tjd = Tdz - Tds;
			 Tje = Tjb - Tja;
			 ii[WS(ios, 25)] = Tjd + Tje;
			 ii[WS(ios, 57)] = Tje - Tjd;
		    }
		    {
			 E TdJ, TdQ, TiX, Tj6;
			 TdJ = TdF + TdI;
			 TdQ = TdM + TdP;
			 ri[WS(ios, 33)] = TdJ - TdQ;
			 ri[WS(ios, 1)] = TdJ + TdQ;
			 TiX = TdS + TdT;
			 Tj6 = TiY + Tj5;
			 ii[WS(ios, 1)] = TiX + Tj6;
			 ii[WS(ios, 33)] = Tj6 - TiX;
		    }
		    {
			 E TdR, TdU, Tj7, Tj8;
			 TdR = TdF - TdI;
			 TdU = TdS - TdT;
			 ri[WS(ios, 49)] = TdR - TdU;
			 ri[WS(ios, 17)] = TdR + TdU;
			 Tj7 = TdP - TdM;
			 Tj8 = Tj5 - TiY;
			 ii[WS(ios, 17)] = Tj7 + Tj8;
			 ii[WS(ios, 49)] = Tj8 - Tj7;
		    }
	       }
	  }
     }
     return W;
}

static const tw_instr twinstr[] = {
     {TW_COS, 0, 1},
     {TW_SIN, 0, 1},
     {TW_COS, 0, 3},
     {TW_SIN, 0, 3},
     {TW_COS, 0, 9},
     {TW_SIN, 0, 9},
     {TW_COS, 0, 27},
     {TW_SIN, 0, 27},
     {TW_COS, 0, 63},
     {TW_SIN, 0, 63},
     {TW_NEXT, 1, 0}
};

static const ct_desc desc = { 64, "t2_64", twinstr, {880, 386, 274, 0}, &GENUS, 0, 0, 0 };

void X(codelet_t2_64) (planner *p) {
     X(kdft_dit_register) (p, t2_64, &desc);
}