hf_64.c

fftw-3.0.1

		    Tbs = T96 + T99;
		    Tbt = FMA(KP923879532, Tbr, KP382683432 * Tbs);
		    Tbz = FNMS(KP382683432, Tbr, KP923879532 * Tbs);
	       }
	       {
		    E TdY, TdZ, Tbo, Tbp;
		    TdY = TdW - TdX;
		    TdZ = T5u - T5F;
		    Te0 = TdY - TdZ;
		    Tee = TdZ + TdY;
		    Tbo = T8S - T8T;
		    Tbp = T8V + T8Y;
		    Tbq = FNMS(KP382683432, Tbp, KP923879532 * Tbo);
		    TbA = FMA(KP382683432, Tbo, KP923879532 * Tbp);
	       }
	  }
	  {
	       E T3X, TdB, T8a, T8d, T4v, Tdx, T80, T85, T48, TdC, T8b, T8g, T4k, Tdw, T7X;
	       E T84;
	       {
		    E T3R, T88, T3W, T89;
		    {
			 E T3O, T3Q, T3N, T3P;
			 T3O = rio[WS(ios, 5)];
			 T3Q = iio[-WS(ios, 58)];
			 T3N = W[8];
			 T3P = W[9];
			 T3R = FMA(T3N, T3O, T3P * T3Q);
			 T88 = FNMS(T3P, T3O, T3N * T3Q);
		    }
		    {
			 E T3T, T3V, T3S, T3U;
			 T3T = rio[WS(ios, 37)];
			 T3V = iio[-WS(ios, 26)];
			 T3S = W[72];
			 T3U = W[73];
			 T3W = FMA(T3S, T3T, T3U * T3V);
			 T89 = FNMS(T3U, T3T, T3S * T3V);
		    }
		    T3X = T3R + T3W;
		    TdB = T88 + T89;
		    T8a = T88 - T89;
		    T8d = T3R - T3W;
	       }
	       {
		    E T4p, T7Y, T4u, T7Z;
		    {
			 E T4m, T4o, T4l, T4n;
			 T4m = rio[WS(ios, 13)];
			 T4o = iio[-WS(ios, 50)];
			 T4l = W[24];
			 T4n = W[25];
			 T4p = FMA(T4l, T4m, T4n * T4o);
			 T7Y = FNMS(T4n, T4m, T4l * T4o);
		    }
		    {
			 E T4r, T4t, T4q, T4s;
			 T4r = rio[WS(ios, 45)];
			 T4t = iio[-WS(ios, 18)];
			 T4q = W[88];
			 T4s = W[89];
			 T4u = FMA(T4q, T4r, T4s * T4t);
			 T7Z = FNMS(T4s, T4r, T4q * T4t);
		    }
		    T4v = T4p + T4u;
		    Tdx = T7Y + T7Z;
		    T80 = T7Y - T7Z;
		    T85 = T4p - T4u;
	       }
	       {
		    E T42, T8e, T47, T8f;
		    {
			 E T3Z, T41, T3Y, T40;
			 T3Z = rio[WS(ios, 21)];
			 T41 = iio[-WS(ios, 42)];
			 T3Y = W[40];
			 T40 = W[41];
			 T42 = FMA(T3Y, T3Z, T40 * T41);
			 T8e = FNMS(T40, T3Z, T3Y * T41);
		    }
		    {
			 E T44, T46, T43, T45;
			 T44 = rio[WS(ios, 53)];
			 T46 = iio[-WS(ios, 10)];
			 T43 = W[104];
			 T45 = W[105];
			 T47 = FMA(T43, T44, T45 * T46);
			 T8f = FNMS(T45, T44, T43 * T46);
		    }
		    T48 = T42 + T47;
		    TdC = T8e + T8f;
		    T8b = T42 - T47;
		    T8g = T8e - T8f;
	       }
	       {
		    E T4e, T82, T4j, T83;
		    {
			 E T4b, T4d, T4a, T4c;
			 T4b = rio[WS(ios, 61)];
			 T4d = iio[-WS(ios, 2)];
			 T4a = W[120];
			 T4c = W[121];
			 T4e = FMA(T4a, T4b, T4c * T4d);
			 T82 = FNMS(T4c, T4b, T4a * T4d);
		    }
		    {
			 E T4g, T4i, T4f, T4h;
			 T4g = rio[WS(ios, 29)];
			 T4i = iio[-WS(ios, 34)];
			 T4f = W[56];
			 T4h = W[57];
			 T4j = FMA(T4f, T4g, T4h * T4i);
			 T83 = FNMS(T4h, T4g, T4f * T4i);
		    }
		    T4k = T4e + T4j;
		    Tdw = T82 + T83;
		    T7X = T4e - T4j;
		    T84 = T82 - T83;
	       }
	       {
		    E T49, T4w, TdA, TdD;
		    T49 = T3X + T48;
		    T4w = T4k + T4v;
		    T4x = T49 + T4w;
		    TfJ = T4w - T49;
		    TdA = T3X - T48;
		    TdD = TdB - TdC;
		    TdE = TdA + TdD;
		    TdM = TdD - TdA;
	       }
	       {
		    E TfM, TfN, T81, T86;
		    TfM = TdB + TdC;
		    TfN = Tdw + Tdx;
		    TfO = TfM - TfN;
		    Tgu = TfM + TfN;
		    T81 = T7X - T80;
		    T86 = T84 + T85;
		    T87 = FNMS(KP923879532, T86, KP382683432 * T81);
		    T8v = FMA(KP382683432, T86, KP923879532 * T81);
	       }
	       {
		    E T8c, T8h, Tb8, Tb9;
		    T8c = T8a + T8b;
		    T8h = T8d - T8g;
		    T8i = FMA(KP923879532, T8c, KP382683432 * T8h);
		    T8u = FNMS(KP923879532, T8h, KP382683432 * T8c);
		    Tb8 = T8a - T8b;
		    Tb9 = T8d + T8g;
		    Tba = FMA(KP382683432, Tb8, KP923879532 * Tb9);
		    Tbg = FNMS(KP382683432, Tb9, KP923879532 * Tb8);
	       }
	       {
		    E Tdv, Tdy, Tb5, Tb6;
		    Tdv = T4k - T4v;
		    Tdy = Tdw - Tdx;
		    Tdz = Tdv - Tdy;
		    TdN = Tdv + Tdy;
		    Tb5 = T7X + T80;
		    Tb6 = T84 - T85;
		    Tb7 = FNMS(KP382683432, Tb6, KP923879532 * Tb5);
		    Tbh = FMA(KP923879532, Tb6, KP382683432 * Tb5);
	       }
	  }
	  {
	       E T1t, Tgn, TgK, TgL, TgV, Th1, T30, Th0, T66, TgX, Tgw, TgE, TgB, TgF, Tgq;
	       E TgM;
	       {
		    E TH, T1s, TgI, TgJ;
		    TH = Tj + TG;
		    T1s = T14 + T1r;
		    T1t = TH + T1s;
		    Tgn = TH - T1s;
		    TgI = Tgt + Tgu;
		    TgJ = Tgy + Tgz;
		    TgK = TgI - TgJ;
		    TgL = TgI + TgJ;
	       }
	       {
		    E TgN, TgU, T2e, T2Z;
		    TgN = Tfq + Tfr;
		    TgU = TgO + TgT;
		    TgV = TgN + TgU;
		    Th1 = TgU - TgN;
		    T2e = T1Q + T2d;
		    T2Z = T2B + T2Y;
		    T30 = T2e + T2Z;
		    Th0 = T2Z - T2e;
	       }
	       {
		    E T4y, T65, Tgs, Tgv;
		    T4y = T3M + T4x;
		    T65 = T5j + T64;
		    T66 = T4y + T65;
		    TgX = T65 - T4y;
		    Tgs = T3M - T4x;
		    Tgv = Tgt - Tgu;
		    Tgw = Tgs + Tgv;
		    TgE = Tgv - Tgs;
	       }
	       {
		    E Tgx, TgA, Tgo, Tgp;
		    Tgx = T5j - T64;
		    TgA = Tgy - Tgz;
		    TgB = Tgx - TgA;
		    TgF = Tgx + TgA;
		    Tgo = Tfu + Tfv;
		    Tgp = TfA + TfB;
		    Tgq = Tgo - Tgp;
		    TgM = Tgo + Tgp;
	       }
	       {
		    E T31, TgW, TgH, TgY;
		    T31 = T1t + T30;
		    iio[-WS(ios, 32)] = T31 - T66;
		    rio[0] = T31 + T66;
		    TgW = TgM + TgV;
		    rio[WS(ios, 32)] = TgL - TgW;
		    iio[0] = TgL + TgW;
		    TgH = T1t - T30;
		    iio[-WS(ios, 48)] = TgH - TgK;
		    rio[WS(ios, 16)] = TgH + TgK;
		    TgY = TgV - TgM;
		    rio[WS(ios, 48)] = TgX - TgY;
		    iio[-WS(ios, 16)] = TgX + TgY;
	       }
	       {
		    E Tgr, TgC, TgZ, Th2;
		    Tgr = Tgn + Tgq;
		    TgC = KP707106781 * (Tgw + TgB);
		    iio[-WS(ios, 40)] = Tgr - TgC;
		    rio[WS(ios, 8)] = Tgr + TgC;
		    TgZ = KP707106781 * (TgE + TgF);
		    Th2 = Th0 + Th1;
		    rio[WS(ios, 40)] = TgZ - Th2;
		    iio[-WS(ios, 8)] = TgZ + Th2;
	       }
	       {
		    E TgD, TgG, Th3, Th4;
		    TgD = Tgn - Tgq;
		    TgG = KP707106781 * (TgE - TgF);
		    iio[-WS(ios, 56)] = TgD - TgG;
		    rio[WS(ios, 24)] = TgD + TgG;
		    Th3 = KP707106781 * (TgB - Tgw);
		    Th4 = Th1 - Th0;
		    rio[WS(ios, 56)] = Th3 - Th4;
		    iio[-WS(ios, 24)] = Th3 + Th4;
	       }
	  }
	  {
	       E Tft, Tg7, Tgh, Tgl, Th9, Thf, TfE, Th6, TfQ, Tg4, Tga, The, Tge, Tgk, Tg1;
	       E Tg5;
	       {
		    E Tfp, Tfs, Tgf, Tgg;
		    Tfp = Tj - TG;
		    Tfs = Tfq - Tfr;
		    Tft = Tfp - Tfs;
		    Tg7 = Tfp + Tfs;
		    Tgf = TfR + TfU;
		    Tgg = TfY + TfZ;
		    Tgh = FNMS(KP382683432, Tgg, KP923879532 * Tgf);
		    Tgl = FMA(KP923879532, Tgg, KP382683432 * Tgf);
	       }
	       {
		    E Th7, Th8, Tfy, TfD;
		    Th7 = T1r - T14;
		    Th8 = TgT - TgO;
		    Th9 = Th7 + Th8;
		    Thf = Th8 - Th7;
		    Tfy = Tfw - Tfx;
		    TfD = Tfz + TfC;
		    TfE = KP707106781 * (Tfy - TfD);
		    Th6 = KP707106781 * (Tfy + TfD);
	       }
	       {
		    E TfK, TfP, Tg8, Tg9;
		    TfK = TfI - TfJ;
		    TfP = TfL - TfO;
		    TfQ = FMA(KP923879532, TfK, KP382683432 * TfP);
		    Tg4 = FNMS(KP923879532, TfP, KP382683432 * TfK);
		    Tg8 = Tfx + Tfw;
		    Tg9 = Tfz - TfC;
		    Tga = KP707106781 * (Tg8 + Tg9);
		    The = KP707106781 * (Tg9 - Tg8);
	       }
	       {
		    E Tgc, Tgd, TfV, Tg0;
		    Tgc = TfI + TfJ;
		    Tgd = TfL + TfO;
		    Tge = FMA(KP382683432, Tgc, KP923879532 * Tgd);
		    Tgk = FNMS(KP382683432, Tgd, KP923879532 * Tgc);
		    TfV = TfR - TfU;
		    Tg0 = TfY - TfZ;
		    Tg1 = FNMS(KP923879532, Tg0, KP382683432 * TfV);
		    Tg5 = FMA(KP382683432, Tg0, KP923879532 * TfV);
	       }
	       {
		    E TfF, Tg2, Thd, Thg;
		    TfF = Tft + TfE;
		    Tg2 = TfQ + Tg1;
		    iio[-WS(ios, 44)] = TfF - Tg2;
		    rio[WS(ios, 12)] = TfF + Tg2;
		    Thd = Tg4 + Tg5;
		    Thg = The + Thf;
		    rio[WS(ios, 44)] = Thd - Thg;
		    iio[-WS(ios, 12)] = Thd + Thg;
	       }
	       {
		    E Tg3, Tg6, Thh, Thi;
		    Tg3 = Tft - TfE;
		    Tg6 = Tg4 - Tg5;
		    iio[-WS(ios, 60)] = Tg3 - Tg6;
		    rio[WS(ios, 28)] = Tg3 + Tg6;
		    Thh = Tg1 - TfQ;
		    Thi = Thf - The;
		    rio[WS(ios, 60)] = Thh - Thi;
		    iio[-WS(ios, 28)] = Thh + Thi;
	       }
	       {
		    E Tgb, Tgi, Th5, Tha;
		    Tgb = Tg7 + Tga;
		    Tgi = Tge + Tgh;
		    iio[-WS(ios, 36)] = Tgb - Tgi;
		    rio[WS(ios, 4)] = Tgb + Tgi;
		    Th5 = Tgk + Tgl;
		    Tha = Th6 + Th9;
		    rio[WS(ios, 36)] = Th5 - Tha;
		    iio[-WS(ios, 4)] = Th5 + Tha;
	       }
	       {
		    E Tgj, Tgm, Thb, Thc;
		    Tgj = Tg7 - Tga;
		    Tgm = Tgk - Tgl;
		    iio[-WS(ios, 52)] = Tgj - Tgm;
		    rio[WS(ios, 20)] = Tgj + Tgm;
		    Thb = Tgh - Tge;
		    Thc = Th9 - Th6;
		    rio[WS(ios, 52)] = Thb - Thc;
		    iio[-WS(ios, 20)] = Thb + Thc;
	       }
	  }
	  {
	       E Td1, Ten, Tdo, ThA, ThD, ThJ, Teq, ThI, Teh, TeB, Tel, Tex, TdQ, TeA, Tek;
	       E Teu;
	       {
		    E TcP, Td0, Teo, Tep;
		    TcP = TcL - TcO;
		    Td0 = KP707106781 * (TcU - TcZ);
		    Td1 = TcP - Td0;
		    Ten = TcP + Td0;
		    {
			 E Tdc, Tdn, ThB, ThC;
			 Tdc = FNMS(KP923879532, Tdb, KP382683432 * Td6);
			 Tdn = FMA(KP382683432, Tdh, KP923879532 * Tdm);
			 Tdo = Tdc - Tdn;
			 ThA = Tdc + Tdn;
			 ThB = KP707106781 * (TeF - TeE);
			 ThC = Thn - Thm;
			 ThD = ThB + ThC;
			 ThJ = ThC - ThB;
		    }
		    Teo = FMA(KP923879532, Td6, KP382683432 * Tdb);
		    Tep = FNMS(KP923879532, Tdh, KP382683432 * Tdm);
		    Teq = Teo + Tep;
		    ThI = Tep - Teo;
		    {
			 E Te7, Tev, Teg, Tew, Te6, Tef;
			 Te6 = KP707106781 * (Te0 - Te5);
			 Te7 = TdV - Te6;
			 Tev = TdV + Te6;
			 Tef = KP707106781 * (Ted - Tee);
			 Teg = Tec - Tef;
			 Tew = Tec + Tef;
			 Teh = FNMS(KP980785280, Teg, KP195090322 * Te7);
			 TeB = FMA(KP831469612, Tew, KP555570233 * Tev);
			 Tel = FMA(KP195090322, Teg, KP980785280 * Te7);
			 Tex = FNMS(KP555570233, Tew, KP831469612 * Tev);
		    }
		    {
			 E TdG, Tes, TdP, Tet, TdF, TdO;
			 TdF = KP707106781 * (Tdz - TdE);
			 TdG = Tdu - TdF;
			 Tes = Tdu + TdF;
			 TdO = KP707106781 * (TdM - TdN);
			 TdP = TdL - TdO;
			 Tet = TdL + TdO;
			 TdQ = FMA(KP980785280, TdG, KP195090322 * TdP);
			 TeA = FNMS(KP555570233, Tet, KP831469612 * Tes);
			 Tek = FNMS(KP980785280, TdP, KP195090322 * TdG);
			 Teu = FMA(KP555570233, Tes, KP831469612 * Tet);
		    }
	       }
	       {
		    E Tdp, Tei, ThH, ThK;
		    Tdp = Td1 + Tdo;
		    Tei = TdQ + Teh;
		    iio[-WS(ios, 46)] = Tdp - Tei;
		    rio[WS(ios, 14)] = Tdp + Tei;
		    ThH = Tek + Tel;
		    ThK = ThI + ThJ;
		    rio[WS(ios, 46)] = ThH - ThK;
		    iio[-WS(ios, 14)] = ThH + ThK;
	       }
	       {
		    E Tej, Tem, ThL, ThM;
		    Tej = Td1 - Tdo;
		    Tem = Tek - Tel;
		    iio[-WS(ios, 62)] = Tej - Tem;
		    rio[WS(ios, 30)] = Tej + Tem;
		    ThL = Teh - TdQ;
		    ThM = ThJ - ThI;
		    rio[WS(ios, 62)] = ThL - ThM;
		    iio[-WS(ios, 30)] = ThL + ThM;
	       }
	       {
		    E Ter, Tey, Thz, ThE;
		    Ter = Ten + Teq;
		    Tey = Teu + Tex;
		    iio[-WS(ios, 38)] = Ter - Tey;
		    rio[WS(ios, 6)] = Ter + Tey;
		    Thz = TeA + TeB;
		    ThE = ThA + ThD;
		    rio[WS(ios, 38)] = Thz - ThE;
		    iio[-WS(ios, 6)] = Thz + ThE;
	       }
	       {
		    E Tez, TeC, ThF, ThG;
		    Tez = Ten - Teq;
		    TeC = TeA - TeB;
		    iio[-WS(ios, 54)] = Tez - TeC;
		    rio[WS(ios, 22)] = Tez + TeC;
		    ThF = Tex - Teu;
		    ThG = ThD - ThA;
		    rio[WS(ios, 54)] = ThF - ThG;
		    iio[-WS(ios, 22)] = ThF + ThG;
	       }
	  }
	  {
	       E TeH, Tf9, TeO, Thk, Thp, Thv, Tfc, Thu, Tf3, Tfn, Tf7, Tfj, TeW, Tfm, Tf6;
	       E Tfg;
	       {
		    E TeD, TeG, Tfa, Tfb;
		    TeD = TcL + TcO;
		    TeG = KP707106781 * (TeE + TeF);
		    TeH = TeD - TeG;
		    Tf9 = TeD + TeG;
		    {
			 E TeK, TeN, Thl, Tho;
			 TeK = FNMS(KP382683432, TeJ, KP923879532 * TeI);
			 TeN = FMA(KP923879532, TeL, KP382683432 * TeM);
			 TeO = TeK - TeN;
			 Thk = TeK + TeN;
			 Thl = KP707106781 * (TcU + TcZ);
			 Tho = Thm + Thn;
			 Thp = Thl + Tho;
			 Thv = Tho - Thl;
		    }
		    Tfa = FMA(KP382683432, TeI, KP923879532 * TeJ);
		    Tfb = FNMS(KP382683432, TeL, KP923879532 * TeM);
		    Tfc = Tfa + Tfb;
		    Thu = Tfb - Tfa;
		    {
			 E TeZ, Tfh, Tf2, Tfi, TeY, Tf1;
			 TeY = KP707106781 * (Tee + Ted);
			 TeZ = TeX - TeY;
			 Tfh = TeX + TeY;
			 Tf1 = KP707106781 * (Te0 + Te5);
			 Tf2 = Tf0 - Tf1;
			 Tfi = Tf0 + Tf1;
			 Tf3 = FNMS(KP831469612, Tf2, KP555570233 * TeZ);
			 Tfn = FMA(KP195090322, Tfh, KP980785280 * Tfi);
			 Tf7 = FMA(KP831469612, TeZ, KP555570233 * Tf2);
			 Tfj = FNMS(KP195090322, Tfi, KP980785280 * Tfh);
		    }
		    {
			 E TeS, Tfe, TeV, Tff, TeR, TeU;
			 TeR = KP707106781 * (TdE + Tdz);
			 TeS = TeQ - TeR;
			 Tfe = TeQ + TeR;
			 TeU = KP707106781 * (TdM + TdN);
			 TeV = TeT - TeU;
			 Tff = TeT + TeU;
			 TeW = FMA(KP555570233, TeS, KP831469612 * TeV);
			 Tfm = FNMS(KP195090322, Tfe, KP980785280 * Tff);
			 Tf6 = FNMS(KP831469612, TeS, KP555570233 * TeV);
			 Tfg = FMA(KP980785280, Tfe, KP195090322 * Tff);
		    }
	       }
	       {
		    E TeP, Tf4, Tht, Thw;
		    TeP = TeH + TeO;
		    Tf4 = TeW + Tf3;
		    iio[-WS(ios, 42)] = TeP - Tf4;
		    rio[WS(ios, 10)] = TeP + Tf4;
		    Tht = Tf6 + Tf7;
		    Thw = Thu + Thv;
		    rio[WS(ios, 42)] = Tht - Thw;
		    iio[-WS(ios, 10)] = Tht + Thw;
	       }
	       {
		    E Tf5, Tf8, Thx, Thy;