q1_6.c

fftw-3.0.1

	       }
	       T4G = T4C + T4F;
	       T5l = KP866025403 * (T4F - T4C);
	       T4L = T4J + T4K;
	       T4W = KP866025403 * (T4K - T4J);
	  }
	  rio[0] = T3 + Ta;
	  iio[0] = TM + TN;
	  rio[WS(is, 1)] = TX + T14;
	  iio[WS(is, 1)] = T1G + T1H;
	  rio[WS(is, 3)] = T2L + T2S;
	  rio[WS(is, 2)] = T1R + T1Y;
	  iio[WS(is, 2)] = T2A + T2B;
	  iio[WS(is, 3)] = T3u + T3v;
	  iio[WS(is, 4)] = T4o + T4p;
	  iio[WS(is, 5)] = T5i + T5j;
	  rio[WS(is, 5)] = T4z + T4G;
	  rio[WS(is, 4)] = T3F + T3M;
	  {
	       E T1w, T1y, T1v, T1x;
	       T1w = T16 + T19;
	       T1y = T1n + T1o;
	       T1v = W[4];
	       T1x = W[5];
	       rio[WS(vs, 3) + WS(is, 1)] = FMA(T1v, T1w, T1x * T1y);
	       iio[WS(vs, 3) + WS(is, 1)] = FNMS(T1x, T1w, T1v * T1y);
	  }
	  {
	       E T58, T5a, T57, T59;
	       T58 = T4I + T4L;
	       T5a = T4Z + T50;
	       T57 = W[4];
	       T59 = W[5];
	       rio[WS(vs, 3) + WS(is, 5)] = FMA(T57, T58, T59 * T5a);
	       iio[WS(vs, 3) + WS(is, 5)] = FNMS(T59, T58, T57 * T5a);
	  }
	  {
	       E TC, TE, TB, TD;
	       TC = Tc + Tf;
	       TE = Tt + Tu;
	       TB = W[4];
	       TD = W[5];
	       rio[WS(vs, 3)] = FMA(TB, TC, TD * TE);
	       iio[WS(vs, 3)] = FNMS(TD, TC, TB * TE);
	  }
	  {
	       E T4e, T4g, T4d, T4f;
	       T4e = T3O + T3R;
	       T4g = T45 + T46;
	       T4d = W[4];
	       T4f = W[5];
	       rio[WS(vs, 3) + WS(is, 4)] = FMA(T4d, T4e, T4f * T4g);
	       iio[WS(vs, 3) + WS(is, 4)] = FNMS(T4f, T4e, T4d * T4g);
	  }
	  {
	       E T3k, T3m, T3j, T3l;
	       T3k = T2U + T2X;
	       T3m = T3b + T3c;
	       T3j = W[4];
	       T3l = W[5];
	       rio[WS(vs, 3) + WS(is, 3)] = FMA(T3j, T3k, T3l * T3m);
	       iio[WS(vs, 3) + WS(is, 3)] = FNMS(T3l, T3k, T3j * T3m);
	  }
	  {
	       E T2q, T2s, T2p, T2r;
	       T2q = T20 + T23;
	       T2s = T2h + T2i;
	       T2p = W[4];
	       T2r = W[5];
	       rio[WS(vs, 3) + WS(is, 2)] = FMA(T2p, T2q, T2r * T2s);
	       iio[WS(vs, 3) + WS(is, 2)] = FNMS(T2r, T2q, T2p * T2s);
	  }
	  {
	       E T5g, T5o, T5m, T5q, T5c, T5k;
	       T5c = FNMS(KP500000000, T4G, T4z);
	       T5g = T5c - T5f;
	       T5o = T5c + T5f;
	       T5k = FNMS(KP500000000, T5j, T5i);
	       T5m = T5k - T5l;
	       T5q = T5l + T5k;
	       {
		    E T5b, T5h, T5n, T5p;
		    T5b = W[2];
		    T5h = W[3];
		    rio[WS(vs, 2) + WS(is, 5)] = FMA(T5b, T5g, T5h * T5m);
		    iio[WS(vs, 2) + WS(is, 5)] = FNMS(T5h, T5g, T5b * T5m);
		    T5n = W[6];
		    T5p = W[7];
		    rio[WS(vs, 4) + WS(is, 5)] = FMA(T5n, T5o, T5p * T5q);
		    iio[WS(vs, 4) + WS(is, 5)] = FNMS(T5p, T5o, T5n * T5q);
	       }
	  }
	  {
	       E To, Ty, Tw, TA, Tg, Tv;
	       Tg = FNMS(KP500000000, Tf, Tc);
	       To = Tg + Tn;
	       Ty = Tg - Tn;
	       Tv = FNMS(KP500000000, Tu, Tt);
	       Tw = Tq + Tv;
	       TA = Tv - Tq;
	       {
		    E Tb, Tp, Tx, Tz;
		    Tb = W[0];
		    Tp = W[1];
		    rio[WS(vs, 1)] = FMA(Tb, To, Tp * Tw);
		    iio[WS(vs, 1)] = FNMS(Tp, To, Tb * Tw);
		    Tx = W[8];
		    Tz = W[9];
		    rio[WS(vs, 5)] = FMA(Tx, Ty, Tz * TA);
		    iio[WS(vs, 5)] = FNMS(Tz, Ty, Tx * TA);
	       }
	  }
	  {
	       E T36, T3g, T3e, T3i, T2Y, T3d;
	       T2Y = FNMS(KP500000000, T2X, T2U);
	       T36 = T2Y + T35;
	       T3g = T2Y - T35;
	       T3d = FNMS(KP500000000, T3c, T3b);
	       T3e = T38 + T3d;
	       T3i = T3d - T38;
	       {
		    E T2T, T37, T3f, T3h;
		    T2T = W[0];
		    T37 = W[1];
		    rio[WS(vs, 1) + WS(is, 3)] = FMA(T2T, T36, T37 * T3e);
		    iio[WS(vs, 1) + WS(is, 3)] = FNMS(T37, T36, T2T * T3e);
		    T3f = W[8];
		    T3h = W[9];
		    rio[WS(vs, 5) + WS(is, 3)] = FMA(T3f, T3g, T3h * T3i);
		    iio[WS(vs, 5) + WS(is, 3)] = FNMS(T3h, T3g, T3f * T3i);
	       }
	  }
	  {
	       E T2y, T2G, T2E, T2I, T2u, T2C;
	       T2u = FNMS(KP500000000, T1Y, T1R);
	       T2y = T2u - T2x;
	       T2G = T2u + T2x;
	       T2C = FNMS(KP500000000, T2B, T2A);
	       T2E = T2C - T2D;
	       T2I = T2D + T2C;
	       {
		    E T2t, T2z, T2F, T2H;
		    T2t = W[2];
		    T2z = W[3];
		    rio[WS(vs, 2) + WS(is, 2)] = FMA(T2t, T2y, T2z * T2E);
		    iio[WS(vs, 2) + WS(is, 2)] = FNMS(T2z, T2y, T2t * T2E);
		    T2F = W[6];
		    T2H = W[7];
		    rio[WS(vs, 4) + WS(is, 2)] = FMA(T2F, T2G, T2H * T2I);
		    iio[WS(vs, 4) + WS(is, 2)] = FNMS(T2H, T2G, T2F * T2I);
	       }
	  }
	  {
	       E T3s, T3A, T3y, T3C, T3o, T3w;
	       T3o = FNMS(KP500000000, T2S, T2L);
	       T3s = T3o - T3r;
	       T3A = T3o + T3r;
	       T3w = FNMS(KP500000000, T3v, T3u);
	       T3y = T3w - T3x;
	       T3C = T3x + T3w;
	       {
		    E T3n, T3t, T3z, T3B;
		    T3n = W[2];
		    T3t = W[3];
		    rio[WS(vs, 2) + WS(is, 3)] = FMA(T3n, T3s, T3t * T3y);
		    iio[WS(vs, 2) + WS(is, 3)] = FNMS(T3t, T3s, T3n * T3y);
		    T3z = W[6];
		    T3B = W[7];
		    rio[WS(vs, 4) + WS(is, 3)] = FMA(T3z, T3A, T3B * T3C);
		    iio[WS(vs, 4) + WS(is, 3)] = FNMS(T3B, T3A, T3z * T3C);
	       }
	  }
	  {
	       E T1E, T1M, T1K, T1O, T1A, T1I;
	       T1A = FNMS(KP500000000, T14, TX);
	       T1E = T1A - T1D;
	       T1M = T1A + T1D;
	       T1I = FNMS(KP500000000, T1H, T1G);
	       T1K = T1I - T1J;
	       T1O = T1J + T1I;
	       {
		    E T1z, T1F, T1L, T1N;
		    T1z = W[2];
		    T1F = W[3];
		    rio[WS(vs, 2) + WS(is, 1)] = FMA(T1z, T1E, T1F * T1K);
		    iio[WS(vs, 2) + WS(is, 1)] = FNMS(T1F, T1E, T1z * T1K);
		    T1L = W[6];
		    T1N = W[7];
		    rio[WS(vs, 4) + WS(is, 1)] = FMA(T1L, T1M, T1N * T1O);
		    iio[WS(vs, 4) + WS(is, 1)] = FNMS(T1N, T1M, T1L * T1O);
	       }
	  }
	  {
	       E T4m, T4u, T4s, T4w, T4i, T4q;
	       T4i = FNMS(KP500000000, T3M, T3F);
	       T4m = T4i - T4l;
	       T4u = T4i + T4l;
	       T4q = FNMS(KP500000000, T4p, T4o);
	       T4s = T4q - T4r;
	       T4w = T4r + T4q;
	       {
		    E T4h, T4n, T4t, T4v;
		    T4h = W[2];
		    T4n = W[3];
		    rio[WS(vs, 2) + WS(is, 4)] = FMA(T4h, T4m, T4n * T4s);
		    iio[WS(vs, 2) + WS(is, 4)] = FNMS(T4n, T4m, T4h * T4s);
		    T4t = W[6];
		    T4v = W[7];
		    rio[WS(vs, 4) + WS(is, 4)] = FMA(T4t, T4u, T4v * T4w);
		    iio[WS(vs, 4) + WS(is, 4)] = FNMS(T4v, T4u, T4t * T4w);
	       }
	  }
	  {
	       E TK, TS, TQ, TU, TG, TO;
	       TG = FNMS(KP500000000, Ta, T3);
	       TK = TG - TJ;
	       TS = TG + TJ;
	       TO = FNMS(KP500000000, TN, TM);
	       TQ = TO - TP;
	       TU = TP + TO;
	       {
		    E TF, TL, TR, TT;
		    TF = W[2];
		    TL = W[3];
		    rio[WS(vs, 2)] = FMA(TF, TK, TL * TQ);
		    iio[WS(vs, 2)] = FNMS(TL, TK, TF * TQ);
		    TR = W[6];
		    TT = W[7];
		    rio[WS(vs, 4)] = FMA(TR, TS, TT * TU);
		    iio[WS(vs, 4)] = FNMS(TT, TS, TR * TU);
	       }
	  }
	  {
	       E T2c, T2m, T2k, T2o, T24, T2j;
	       T24 = FNMS(KP500000000, T23, T20);
	       T2c = T24 + T2b;
	       T2m = T24 - T2b;
	       T2j = FNMS(KP500000000, T2i, T2h);
	       T2k = T2e + T2j;
	       T2o = T2j - T2e;
	       {
		    E T1Z, T2d, T2l, T2n;
		    T1Z = W[0];
		    T2d = W[1];
		    rio[WS(vs, 1) + WS(is, 2)] = FMA(T1Z, T2c, T2d * T2k);
		    iio[WS(vs, 1) + WS(is, 2)] = FNMS(T2d, T2c, T1Z * T2k);
		    T2l = W[8];
		    T2n = W[9];
		    rio[WS(vs, 5) + WS(is, 2)] = FMA(T2l, T2m, T2n * T2o);
		    iio[WS(vs, 5) + WS(is, 2)] = FNMS(T2n, T2m, T2l * T2o);
	       }
	  }
	  {
	       E T40, T4a, T48, T4c, T3S, T47;
	       T3S = FNMS(KP500000000, T3R, T3O);
	       T40 = T3S + T3Z;
	       T4a = T3S - T3Z;
	       T47 = FNMS(KP500000000, T46, T45);
	       T48 = T42 + T47;
	       T4c = T47 - T42;
	       {
		    E T3N, T41, T49, T4b;
		    T3N = W[0];
		    T41 = W[1];
		    rio[WS(vs, 1) + WS(is, 4)] = FMA(T3N, T40, T41 * T48);
		    iio[WS(vs, 1) + WS(is, 4)] = FNMS(T41, T40, T3N * T48);
		    T49 = W[8];
		    T4b = W[9];
		    rio[WS(vs, 5) + WS(is, 4)] = FMA(T49, T4a, T4b * T4c);
		    iio[WS(vs, 5) + WS(is, 4)] = FNMS(T4b, T4a, T49 * T4c);
	       }
	  }
	  {
	       E T1i, T1s, T1q, T1u, T1a, T1p;
	       T1a = FNMS(KP500000000, T19, T16);
	       T1i = T1a + T1h;
	       T1s = T1a - T1h;
	       T1p = FNMS(KP500000000, T1o, T1n);
	       T1q = T1k + T1p;
	       T1u = T1p - T1k;
	       {
		    E T15, T1j, T1r, T1t;
		    T15 = W[0];
		    T1j = W[1];
		    rio[WS(vs, 1) + WS(is, 1)] = FMA(T15, T1i, T1j * T1q);
		    iio[WS(vs, 1) + WS(is, 1)] = FNMS(T1j, T1i, T15 * T1q);
		    T1r = W[8];
		    T1t = W[9];
		    rio[WS(vs, 5) + WS(is, 1)] = FMA(T1r, T1s, T1t * T1u);
		    iio[WS(vs, 5) + WS(is, 1)] = FNMS(T1t, T1s, T1r * T1u);
	       }
	  }
	  {
	       E T4U, T54, T52, T56, T4M, T51;
	       T4M = FNMS(KP500000000, T4L, T4I);
	       T4U = T4M + T4T;
	       T54 = T4M - T4T;
	       T51 = FNMS(KP500000000, T50, T4Z);
	       T52 = T4W + T51;
	       T56 = T51 - T4W;
	       {
		    E T4H, T4V, T53, T55;
		    T4H = W[0];
		    T4V = W[1];
		    rio[WS(vs, 1) + WS(is, 5)] = FMA(T4H, T4U, T4V * T52);
		    iio[WS(vs, 1) + WS(is, 5)] = FNMS(T4V, T4U, T4H * T52);
		    T53 = W[8];
		    T55 = W[9];
		    rio[WS(vs, 5) + WS(is, 5)] = FMA(T53, T54, T55 * T56);
		    iio[WS(vs, 5) + WS(is, 5)] = FNMS(T55, T54, T53 * T56);
	       }
	  }
     }
     return W;
}

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

static const ct_desc desc = { 6, "q1_6", twinstr, {192, 84, 84, 0}, &GENUS, 0, 0, 0 };

void X(codelet_q1_6) (planner *p) {
     X(kdft_difsq_register) (p, q1_6, &desc);
}