hb_64.c

fftw-3.0.1

		    TeY = TeW + TeX;
		    Tf2 = Tf0 + Tf1;
		    TeV = W[104];
		    TeZ = W[105];
		    iio[-WS(ios, 10)] = FMA(TeV, TeY, TeZ * Tf2);
		    rio[WS(ios, 53)] = FNMS(TeZ, TeY, TeV * Tf2);
	       }
	       {
		    E Tf4, Tf6, Tf3, Tf5;
		    Tf4 = Tf1 - Tf0;
		    Tf6 = TeW - TeX;
		    Tf3 = W[40];
		    Tf5 = W[41];
		    rio[WS(ios, 21)] = FNMS(Tf5, Tf6, Tf3 * Tf4);
		    iio[-WS(ios, 42)] = FMA(Tf3, Tf6, Tf5 * Tf4);
	       }
	       {
		    E Tei, Tek, Teh, Tej;
		    Tei = Tef - Te8;
		    Tek = TdO - Te3;
		    Teh = W[56];
		    Tej = W[57];
		    rio[WS(ios, 29)] = FNMS(Tej, Tek, Teh * Tei);
		    iio[-WS(ios, 34)] = FMA(Teh, Tek, Tej * Tei);
	       }
	       {
		    E Teo, Tes, Tel, Tep;
		    Teo = Tem + Ten;
		    Tes = Teq + Ter;
		    Tel = W[24];
		    Tep = W[25];
		    rio[WS(ios, 13)] = FNMS(Tep, Tes, Tel * Teo);
		    iio[-WS(ios, 50)] = FMA(Tel, Tes, Tep * Teo);
	       }
	       {
		    E TeI, TeQ, Tex, TeJ;
		    TeI = TeE + TeH;
		    TeQ = TeM + TeP;
		    Tex = W[8];
		    TeJ = W[9];
		    rio[WS(ios, 5)] = FNMS(TeJ, TeQ, Tex * TeI);
		    iio[-WS(ios, 58)] = FMA(Tex, TeQ, TeJ * TeI);
	       }
	       {
		    E TeS, TeU, TeR, TeT;
		    TeS = TeM - TeP;
		    TeU = TeH - TeE;
		    TeR = W[72];
		    TeT = W[73];
		    iio[-WS(ios, 26)] = FMA(TeR, TeS, TeT * TeU);
		    rio[WS(ios, 37)] = FNMS(TeT, TeS, TeR * TeU);
	       }
	       {
		    E Teu, Tew, Tet, Tev;
		    Teu = Teq - Ter;
		    Tew = Ten - Tem;
		    Tet = W[88];
		    Tev = W[89];
		    iio[-WS(ios, 18)] = FMA(Tet, Teu, Tev * Tew);
		    rio[WS(ios, 45)] = FNMS(Tev, Teu, Tet * Tew);
	       }
	  }
	  {
	       E Tcr, Tdw, TcX, Td6, TcI, Tdt, TcS, Tdl, Tbm, TcW, TcL, TcT, Tdd, Tdx, Tdi;
	       E Tds;
	       {
		    E Tcq, Td4, TbZ, Td5, TbF, TbY;
		    Tcq = Tce + Tcp;
		    Td4 = TcA - TcD;
		    TbF = FNMS(KP195090322, TbE, KP980785280 * Tbx);
		    TbY = FMA(KP195090322, TbQ, KP980785280 * TbX);
		    TbZ = TbF + TbY;
		    Td5 = TbY - TbF;
		    Tcr = TbZ + Tcq;
		    Tdw = Td4 - Td5;
		    TcX = Tcq - TbZ;
		    Td6 = Td4 + Td5;
	       }
	       {
		    E TcE, Tdk, TcH, Tdj, TcF, TcG;
		    TcE = TcA + TcD;
		    Tdk = Tcp - Tce;
		    TcF = FMA(KP980785280, TbE, KP195090322 * Tbx);
		    TcG = FNMS(KP195090322, TbX, KP980785280 * TbQ);
		    TcH = TcF + TcG;
		    Tdj = TcF - TcG;
		    TcI = TcE + TcH;
		    Tdt = Tdk - Tdj;
		    TcS = TcE - TcH;
		    Tdl = Tdj + Tdk;
	       }
	       {
		    E TaI, TcJ, Tbl, TcK;
		    {
			 E Taw, TaH, Tb9, Tbk;
			 Taw = Tak + Tav;
			 TaH = TaD + TaG;
			 TaI = FNMS(KP098017140, TaH, KP995184726 * Taw);
			 TcJ = FMA(KP995184726, TaH, KP098017140 * Taw);
			 Tb9 = TaT + Tb8;
			 Tbk = Tbc + Tbj;
			 Tbl = FMA(KP098017140, Tb9, KP995184726 * Tbk);
			 TcK = FNMS(KP098017140, Tbk, KP995184726 * Tb9);
		    }
		    Tbm = TaI + Tbl;
		    TcW = TcJ - TcK;
		    TcL = TcJ + TcK;
		    TcT = Tbl - TaI;
	       }
	       {
		    E Td9, Tdg, Tdc, Tdh;
		    {
			 E Td7, Td8, Tda, Tdb;
			 Td7 = TaD - TaG;
			 Td8 = Tav - Tak;
			 Td9 = FNMS(KP634393284, Td8, KP773010453 * Td7);
			 Tdg = FMA(KP634393284, Td7, KP773010453 * Td8);
			 Tda = TaT - Tb8;
			 Tdb = Tbj - Tbc;
			 Tdc = FMA(KP773010453, Tda, KP634393284 * Tdb);
			 Tdh = FNMS(KP634393284, Tda, KP773010453 * Tdb);
		    }
		    Tdd = Td9 + Tdc;
		    Tdx = Tdg - Tdh;
		    Tdi = Tdg + Tdh;
		    Tds = Tdc - Td9;
	       }
	       {
		    E Tcs, TcM, Ta5, Tct;
		    Tcs = Tbm + Tcr;
		    TcM = TcI + TcL;
		    Ta5 = W[0];
		    Tct = W[1];
		    rio[WS(ios, 1)] = FNMS(Tct, TcM, Ta5 * Tcs);
		    iio[-WS(ios, 62)] = FMA(Ta5, TcM, Tct * Tcs);
	       }
	       {
		    E Tdu, Tdy, Tdr, Tdv;
		    Tdu = Tds + Tdt;
		    Tdy = Tdw + Tdx;
		    Tdr = W[16];
		    Tdv = W[17];
		    rio[WS(ios, 9)] = FNMS(Tdv, Tdy, Tdr * Tdu);
		    iio[-WS(ios, 54)] = FMA(Tdr, Tdy, Tdv * Tdu);
	       }
	       {
		    E TdA, TdC, Tdz, TdB;
		    TdA = Tdw - Tdx;
		    TdC = Tdt - Tds;
		    Tdz = W[80];
		    TdB = W[81];
		    iio[-WS(ios, 22)] = FMA(Tdz, TdA, TdB * TdC);
		    rio[WS(ios, 41)] = FNMS(TdB, TdA, Tdz * TdC);
	       }
	       {
		    E TcO, TcQ, TcN, TcP;
		    TcO = TcI - TcL;
		    TcQ = Tcr - Tbm;
		    TcN = W[64];
		    TcP = W[65];
		    iio[-WS(ios, 30)] = FMA(TcN, TcO, TcP * TcQ);
		    rio[WS(ios, 33)] = FNMS(TcP, TcO, TcN * TcQ);
	       }
	       {
		    E TcU, TcY, TcR, TcV;
		    TcU = TcS + TcT;
		    TcY = TcW + TcX;
		    TcR = W[96];
		    TcV = W[97];
		    iio[-WS(ios, 14)] = FMA(TcR, TcU, TcV * TcY);
		    rio[WS(ios, 49)] = FNMS(TcV, TcU, TcR * TcY);
	       }
	       {
		    E Tde, Tdm, Td3, Tdf;
		    Tde = Td6 + Tdd;
		    Tdm = Tdi + Tdl;
		    Td3 = W[112];
		    Tdf = W[113];
		    iio[-WS(ios, 6)] = FMA(Td3, Tde, Tdf * Tdm);
		    rio[WS(ios, 57)] = FNMS(Tdf, Tde, Td3 * Tdm);
	       }
	       {
		    E Tdo, Tdq, Tdn, Tdp;
		    Tdo = Tdl - Tdi;
		    Tdq = Td6 - Tdd;
		    Tdn = W[48];
		    Tdp = W[49];
		    rio[WS(ios, 25)] = FNMS(Tdp, Tdq, Tdn * Tdo);
		    iio[-WS(ios, 38)] = FMA(Tdn, Tdq, Tdp * Tdo);
	       }
	       {
		    E Td0, Td2, TcZ, Td1;
		    Td0 = TcX - TcW;
		    Td2 = TcS - TcT;
		    TcZ = W[32];
		    Td1 = W[33];
		    rio[WS(ios, 17)] = FNMS(Td1, Td2, TcZ * Td0);
		    iio[-WS(ios, 46)] = FMA(TcZ, Td2, Td1 * Td0);
	       }
	  }
	  {
	       E Tfy, Thd, TgC, TgT, Tgr, Th8, Tgz, TgY, Tgb, TgD, Tgg, Tgy, TgQ, Thc, Th1;
	       E Th9;
	       {
		    E Tfi, TgS, Tfx, TgR, Tfp, Tfw;
		    Tfi = Tfa + Tfh;
		    TgS = Tgp - Tgm;
		    Tfp = FNMS(KP195090322, Tfo, KP980785280 * Tfl);
		    Tfw = FMA(KP980785280, Tfs, KP195090322 * Tfv);
		    Tfx = Tfp + Tfw;
		    TgR = Tfw - Tfp;
		    Tfy = Tfi + Tfx;
		    Thd = TgS - TgR;
		    TgC = Tfi - Tfx;
		    TgT = TgR + TgS;
	       }
	       {
		    E Tgq, TgW, Tgj, TgX, Tgh, Tgi;
		    Tgq = Tgm + Tgp;
		    TgW = Tfa - Tfh;
		    Tgh = FMA(KP195090322, Tfl, KP980785280 * Tfo);
		    Tgi = FNMS(KP195090322, Tfs, KP980785280 * Tfv);
		    Tgj = Tgh + Tgi;
		    TgX = Tgh - Tgi;
		    Tgr = Tgj + Tgq;
		    Th8 = TgW - TgX;
		    Tgz = Tgq - Tgj;
		    TgY = TgW + TgX;
	       }
	       {
		    E TfR, Tge, Tga, Tgf;
		    {
			 E TfJ, TfQ, Tg2, Tg9;
			 TfJ = TfB + TfI;
			 TfQ = TfM + TfP;
			 TfR = FNMS(KP098017140, TfQ, KP995184726 * TfJ);
			 Tge = FMA(KP098017140, TfJ, KP995184726 * TfQ);
			 Tg2 = TfU + Tg1;
			 Tg9 = Tg5 + Tg8;
			 Tga = FMA(KP995184726, Tg2, KP098017140 * Tg9);
			 Tgf = FNMS(KP098017140, Tg2, KP995184726 * Tg9);
		    }
		    Tgb = TfR + Tga;
		    TgD = Tge - Tgf;
		    Tgg = Tge + Tgf;
		    Tgy = Tga - TfR;
	       }
	       {
		    E TgM, TgZ, TgP, Th0;
		    {
			 E TgK, TgL, TgN, TgO;
			 TgK = TfP - TfM;
			 TgL = TfB - TfI;
			 TgM = FNMS(KP634393284, TgL, KP773010453 * TgK);
			 TgZ = FMA(KP773010453, TgL, KP634393284 * TgK);
			 TgN = TfU - Tg1;
			 TgO = Tg8 - Tg5;
			 TgP = FMA(KP634393284, TgN, KP773010453 * TgO);
			 Th0 = FNMS(KP634393284, TgO, KP773010453 * TgN);
		    }
		    TgQ = TgM + TgP;
		    Thc = TgZ - Th0;
		    Th1 = TgZ + Th0;
		    Th9 = TgP - TgM;
	       }
	       {
		    E Tgc, Tgs, Tf7, Tgd;
		    Tgc = Tfy + Tgb;
		    Tgs = Tgg + Tgr;
		    Tf7 = W[124];
		    Tgd = W[125];
		    iio[0] = FMA(Tf7, Tgc, Tgd * Tgs);
		    rio[WS(ios, 63)] = FNMS(Tgd, Tgc, Tf7 * Tgs);
	       }
	       {
		    E Tha, The, Th7, Thb;
		    Tha = Th8 + Th9;
		    The = Thc + Thd;
		    Th7 = W[108];
		    Thb = W[109];
		    iio[-WS(ios, 8)] = FMA(Th7, Tha, Thb * The);
		    rio[WS(ios, 55)] = FNMS(Thb, Tha, Th7 * The);
	       }
	       {
		    E Thg, Thi, Thf, Thh;
		    Thg = Thd - Thc;
		    Thi = Th8 - Th9;
		    Thf = W[44];
		    Thh = W[45];
		    rio[WS(ios, 23)] = FNMS(Thh, Thi, Thf * Thg);
		    iio[-WS(ios, 40)] = FMA(Thf, Thi, Thh * Thg);
	       }
	       {
		    E Tgu, Tgw, Tgt, Tgv;
		    Tgu = Tgr - Tgg;
		    Tgw = Tfy - Tgb;
		    Tgt = W[60];
		    Tgv = W[61];
		    rio[WS(ios, 31)] = FNMS(Tgv, Tgw, Tgt * Tgu);
		    iio[-WS(ios, 32)] = FMA(Tgt, Tgw, Tgv * Tgu);
	       }
	       {
		    E TgA, TgE, Tgx, TgB;
		    TgA = Tgy + Tgz;
		    TgE = TgC + TgD;
		    Tgx = W[28];
		    TgB = W[29];
		    rio[WS(ios, 15)] = FNMS(TgB, TgE, Tgx * TgA);
		    iio[-WS(ios, 48)] = FMA(Tgx, TgE, TgB * TgA);
	       }
	       {
		    E TgU, Th2, TgJ, TgV;
		    TgU = TgQ + TgT;
		    Th2 = TgY + Th1;
		    TgJ = W[12];
		    TgV = W[13];
		    rio[WS(ios, 7)] = FNMS(TgV, Th2, TgJ * TgU);
		    iio[-WS(ios, 56)] = FMA(TgJ, Th2, TgV * TgU);
	       }
	       {
		    E Th4, Th6, Th3, Th5;
		    Th4 = TgY - Th1;
		    Th6 = TgT - TgQ;
		    Th3 = W[76];
		    Th5 = W[77];
		    iio[-WS(ios, 24)] = FMA(Th3, Th4, Th5 * Th6);
		    rio[WS(ios, 39)] = FNMS(Th5, Th4, Th3 * Th6);
	       }
	       {
		    E TgG, TgI, TgF, TgH;
		    TgG = TgC - TgD;
		    TgI = Tgz - Tgy;
		    TgF = W[92];
		    TgH = W[93];
		    iio[-WS(ios, 16)] = FMA(TgF, TgG, TgH * TgI);
		    rio[WS(ios, 47)] = FNMS(TgH, TgG, TgF * TgI);
	       }
	  }
	  {
	       E ThJ, TiG, Ti7, Tig, ThS, TiD, Ti2, Tiv, Thy, Ti6, ThV, Ti3, Tin, TiH, Tis;
	       E TiC;
	       {
		    E ThI, Tie, ThF, Tif, ThB, ThE;
		    ThI = ThG + ThH;
		    Tie = ThM - ThN;
		    ThB = FNMS(KP555570233, ThA, KP831469612 * Thz);
		    ThE = FNMS(KP555570233, ThD, KP831469612 * ThC);
		    ThF = ThB + ThE;
		    Tif = ThE - ThB;
		    ThJ = ThF + ThI;
		    TiG = Tie - Tif;
		    Ti7 = ThI - ThF;
		    Tig = Tie + Tif;
	       }
	       {
		    E ThO, Tiu, ThR, Tit, ThP, ThQ;
		    ThO = ThM + ThN;
		    Tiu = ThH - ThG;
		    ThP = FMA(KP831469612, ThA, KP555570233 * Thz);
		    ThQ = FMA(KP831469612, ThD, KP555570233 * ThC);
		    ThR = ThP - ThQ;
		    Tit = ThP + ThQ;
		    ThS = ThO + ThR;
		    TiD = Tiu - Tit;
		    Ti2 = ThO - ThR;
		    Tiv = Tit + Tiu;
	       }
	       {
		    E Thq, ThT, Thx, ThU;
		    {
			 E Thm, Thp, Tht, Thw;
			 Thm = Thk + Thl;
			 Thp = Thn + Tho;
			 Thq = FNMS(KP290284677, Thp, KP956940335 * Thm);
			 ThT = FMA(KP956940335, Thp, KP290284677 * Thm);
			 Tht = Thr - Ths;
			 Thw = Thu + Thv;
			 Thx = FMA(KP290284677, Tht, KP956940335 * Thw);
			 ThU = FNMS(KP290284677, Thw, KP956940335 * Tht);
		    }
		    Thy = Thq + Thx;
		    Ti6 = ThT - ThU;
		    ThV = ThT + ThU;
		    Ti3 = Thx - Thq;
	       }
	       {
		    E Tij, Tiq, Tim, Tir;
		    {
			 E Tih, Tii, Tik, Til;
			 Tih = Thn - Tho;
			 Tii = Thl - Thk;
			 Tij = FNMS(KP471396736, Tii, KP881921264 * Tih);
			 Tiq = FMA(KP471396736, Tih, KP881921264 * Tii);
			 Tik = Thv - Thu;
			 Til = Ths + Thr;
			 Tim = FNMS(KP881921264, Til, KP471396736 * Tik);
			 Tir = FMA(KP471396736, Til, KP881921264 * Tik);
		    }
		    Tin = Tij + Tim;
		    TiH = Tiq - Tir;
		    Tis = Tiq + Tir;
		    TiC = Tim - Tij;
	       }
	       {
		    E ThK, ThW, Thj, ThL;
		    ThK = Thy + ThJ;
		    ThW = ThS + ThV;
		    Thj = W[4];
		    ThL = W[5];
		    rio[WS(ios, 3)] = FNMS(ThL, ThW, Thj * ThK);
		    iio[-WS(ios, 60)] = FMA(Thj, ThW, ThL * ThK);
	       }
	       {
		    E TiE, TiI, TiB, TiF;
		    TiE = TiC + TiD;
		    TiI = TiG + TiH;
		    TiB = W[20];
		    TiF = W[21];
		    rio[WS(ios, 11)] = FNMS(TiF, TiI, TiB * TiE);
		    iio[-WS(ios, 52)] = FMA(TiB, TiI, TiF * TiE);
	       }
	       {
		    E TiK, TiM, TiJ, TiL;
		    TiK = TiG - TiH;
		    TiM = TiD - TiC;
		    TiJ = W[84];
		    TiL = W[85];
		    iio[-WS(ios, 20)] = FMA(TiJ, TiK, TiL * TiM);
		    rio[WS(ios, 43)] = FNMS(TiL, TiK, TiJ * TiM);
	       }
	       {
		    E ThY, Ti0, ThX, ThZ;
		    ThY = ThS - ThV;
		    Ti0 = ThJ - Thy;
		    ThX = W[68];
		    ThZ = W[69];
		    iio[-WS(ios, 28)] = FMA(ThX, ThY, ThZ * Ti0);
		    rio[WS(ios, 35)] = FNMS(ThZ, ThY, ThX * Ti0);
	       }
	       {
		    E Ti4, Ti8, Ti1, Ti5;
		    Ti4 = Ti2 + Ti3;
		    Ti8 = Ti6 + Ti7;
		    Ti1 = W[100];
		    Ti5 = W[101];
		    iio[-WS(ios, 12)] = FMA(Ti1, Ti4, Ti5 * Ti8);
		    rio[WS(ios, 51)] = FNMS(Ti5, Ti4, Ti1 * Ti8);
	       }
	       {
		    E Tio, Tiw, Tid, Tip;
		    Tio = Tig + Tin;
		    Tiw = Tis + Tiv;
		    Tid = W[116];
		    Tip = W[117];
		    iio[-WS(ios, 4)] = FMA(Tid, Tio, Tip * Tiw);
		    rio[WS(ios, 59)] = FNMS(Tip, Tio, Tid * Tiw);
	       }
	       {
		    E Tiy, TiA, Tix, Tiz;
		    Tiy = Tiv - Tis;
		    TiA = Tig - Tin;
		    Tix = W[52];
		    Tiz = W[53];
		    rio[WS(ios, 27)] = FNMS(Tiz, TiA, Tix * Tiy);
		    iio[-WS(ios, 36)] = FMA(Tix, TiA, Tiz * Tiy);
	       }
	       {
		    E Tia, Tic, Ti9, Tib;
		    Tia = Ti7 - Ti6;
		    Tic = Ti2 - Ti3;
		    Ti9 = W[36];
		    Tib = W[37];
		    rio[WS(ios, 19)] = FNMS(Tib, Tic, Ti9 * Tia);
		    iio[-WS(ios, 44)] = FMA(Ti9, Tic, Tib * Tia);
	       }
	  }
     }
     return W;
}

static const tw_instr twinstr[] = {
     {TW_FULL, 0, 64},
     {TW_NEXT, 1, 0}
};

static const hc2hc_desc desc = { 64, "hb_64", twinstr, {808, 270, 230, 0}, &GENUS, 0, 0, 0 };

void X(codelet_hb_64) (planner *p) {
     X(khc2hc_dif_register) (p, hb_64, &desc);
}