mr2hcii_64.c

fftw-3.0.1

		    T1C = I[WS(is, 13)];
		    T1D = I[WS(is, 45)];
		    T1E = KP707106781 * (T1C - T1D);
		    T1H = KP707106781 * (T1C + T1D);
		    T1G = T1E - T1F;
		    T4m = T1E + T1F;
		    T1J = T1H + T1I;
		    T4l = T1I - T1H;
	       }
	       {
		    E T1B, T1K, T4w, T4x;
		    T1B = FNMS(KP195090322, T1A, KP980785280 * T1x);
		    T1K = FMA(KP980785280, T1G, KP195090322 * T1J);
		    T1L = T1B + T1K;
		    T3j = T1K - T1B;
		    T4w = FMA(KP831469612, T4p, KP555570233 * T4o);
		    T4x = FMA(KP831469612, T4m, KP555570233 * T4l);
		    T4y = T4w + T4x;
		    T5t = T4w - T4x;
	       }
	       {
		    E T1N, T1O, T4n, T4q;
		    T1N = FNMS(KP980785280, T1J, KP195090322 * T1G);
		    T1O = FMA(KP195090322, T1x, KP980785280 * T1A);
		    T1P = T1N - T1O;
		    T3h = T1O + T1N;
		    T4n = FNMS(KP555570233, T4m, KP831469612 * T4l);
		    T4q = FNMS(KP555570233, T4p, KP831469612 * T4o);
		    T4r = T4n - T4q;
		    T5v = T4q + T4n;
	       }
	  }
	  {
	       E Tx, T2N, T69, T6f, T1e, T6e, T2X, T30, T1Y, T2L, T2Q, T62, T2U, T31, T2H;
	       E T2K, Tw, T63;
	       Tw = Tm + Tv;
	       Tx = Td + Tw;
	       T2N = Td - Tw;
	       T63 = T35 - T34;
	       T69 = T63 - T68;
	       T6f = T63 + T68;
	       {
		    E TS, T1d, T2V, T2W;
		    TS = FNMS(KP098017140, TR, KP995184726 * TK);
		    T1d = FMA(KP995184726, T15, KP098017140 * T1c);
		    T1e = TS + T1d;
		    T6e = T1d - TS;
		    T2V = T2b - T2u;
		    T2W = T2y + T2F;
		    T2X = FNMS(KP671558954, T2W, KP740951125 * T2V);
		    T30 = FMA(KP671558954, T2V, KP740951125 * T2W);
	       }
	       {
		    E T1M, T1X, T2O, T2P;
		    T1M = T1s + T1L;
		    T1X = T1P - T1W;
		    T1Y = FMA(KP998795456, T1M, KP049067674 * T1X);
		    T2L = FNMS(KP049067674, T1M, KP998795456 * T1X);
		    T2O = FMA(KP098017140, TK, KP995184726 * TR);
		    T2P = FNMS(KP995184726, T1c, KP098017140 * T15);
		    T2Q = T2O + T2P;
		    T62 = T2P - T2O;
	       }
	       {
		    E T2S, T2T, T2v, T2G;
		    T2S = T1s - T1L;
		    T2T = T1P + T1W;
		    T2U = FMA(KP740951125, T2S, KP671558954 * T2T);
		    T31 = FNMS(KP671558954, T2S, KP740951125 * T2T);
		    T2v = T2b + T2u;
		    T2G = T2y - T2F;
		    T2H = FNMS(KP049067674, T2G, KP998795456 * T2v);
		    T2K = FMA(KP049067674, T2v, KP998795456 * T2G);
	       }
	       {
		    E T1f, T2I, T6b, T6c;
		    T1f = Tx + T1e;
		    T2I = T1Y + T2H;
		    ro[WS(ros, 31)] = T1f - T2I;
		    ro[0] = T1f + T2I;
		    T6b = T2L + T2K;
		    T6c = T62 + T69;
		    io[WS(ios, 31)] = T6b - T6c;
		    io[0] = T6b + T6c;
	       }
	       {
		    E T2J, T2M, T61, T6a;
		    T2J = Tx - T1e;
		    T2M = T2K - T2L;
		    ro[WS(ros, 16)] = T2J - T2M;
		    ro[WS(ros, 15)] = T2J + T2M;
		    T61 = T2H - T1Y;
		    T6a = T62 - T69;
		    io[WS(ios, 16)] = T61 - T6a;
		    io[WS(ios, 15)] = T61 + T6a;
	       }
	       {
		    E T2R, T2Y, T6h, T6i;
		    T2R = T2N + T2Q;
		    T2Y = T2U + T2X;
		    ro[WS(ros, 24)] = T2R - T2Y;
		    ro[WS(ros, 7)] = T2R + T2Y;
		    T6h = T31 + T30;
		    T6i = T6e + T6f;
		    io[WS(ios, 24)] = T6h - T6i;
		    io[WS(ios, 7)] = T6h + T6i;
	       }
	       {
		    E T2Z, T32, T6d, T6g;
		    T2Z = T2N - T2Q;
		    T32 = T30 - T31;
		    ro[WS(ros, 23)] = T2Z - T32;
		    ro[WS(ros, 8)] = T2Z + T32;
		    T6d = T2X - T2U;
		    T6g = T6e - T6f;
		    io[WS(ios, 23)] = T6d - T6g;
		    io[WS(ios, 8)] = T6d + T6g;
	       }
	  }
	  {
	       E T5j, T5L, T6R, T6X, T5q, T6W, T5V, T5Y, T5y, T5J, T5O, T6O, T5S, T5Z, T5F;
	       E T5I, T5i, T6P;
	       T5i = T5g - T5h;
	       T5j = T5f - T5i;
	       T5L = T5f + T5i;
	       T6P = T3Z + T3W;
	       T6R = T6P - T6Q;
	       T6X = T6P + T6Q;
	       {
		    E T5m, T5p, T5T, T5U;
		    T5m = FMA(KP290284677, T5k, KP956940335 * T5l);
		    T5p = FNMS(KP290284677, T5o, KP956940335 * T5n);
		    T5q = T5m - T5p;
		    T6W = T5p + T5m;
		    T5T = T5z + T5A;
		    T5U = T5C + T5D;
		    T5V = FNMS(KP146730474, T5U, KP989176509 * T5T);
		    T5Y = FMA(KP146730474, T5T, KP989176509 * T5U);
	       }
	       {
		    E T5u, T5x, T5M, T5N;
		    T5u = T5s - T5t;
		    T5x = T5v - T5w;
		    T5y = FMA(KP803207531, T5u, KP595699304 * T5x);
		    T5J = FNMS(KP595699304, T5u, KP803207531 * T5x);
		    T5M = FMA(KP956940335, T5o, KP290284677 * T5n);
		    T5N = FNMS(KP290284677, T5l, KP956940335 * T5k);
		    T5O = T5M + T5N;
		    T6O = T5N - T5M;
	       }
	       {
		    E T5Q, T5R, T5B, T5E;
		    T5Q = T5s + T5t;
		    T5R = T5v + T5w;
		    T5S = FMA(KP989176509, T5Q, KP146730474 * T5R);
		    T5Z = FNMS(KP146730474, T5Q, KP989176509 * T5R);
		    T5B = T5z - T5A;
		    T5E = T5C - T5D;
		    T5F = FNMS(KP595699304, T5E, KP803207531 * T5B);
		    T5I = FMA(KP595699304, T5B, KP803207531 * T5E);
	       }
	       {
		    E T5r, T5G, T6T, T6U;
		    T5r = T5j + T5q;
		    T5G = T5y + T5F;
		    ro[WS(ros, 25)] = T5r - T5G;
		    ro[WS(ros, 6)] = T5r + T5G;
		    T6T = T5J + T5I;
		    T6U = T6O + T6R;
		    io[WS(ios, 25)] = T6T - T6U;
		    io[WS(ios, 6)] = T6T + T6U;
	       }
	       {
		    E T5H, T5K, T6N, T6S;
		    T5H = T5j - T5q;
		    T5K = T5I - T5J;
		    ro[WS(ros, 22)] = T5H - T5K;
		    ro[WS(ros, 9)] = T5H + T5K;
		    T6N = T5F - T5y;
		    T6S = T6O - T6R;
		    io[WS(ios, 22)] = T6N - T6S;
		    io[WS(ios, 9)] = T6N + T6S;
	       }
	       {
		    E T5P, T5W, T6Z, T70;
		    T5P = T5L + T5O;
		    T5W = T5S + T5V;
		    ro[WS(ros, 30)] = T5P - T5W;
		    ro[WS(ros, 1)] = T5P + T5W;
		    T6Z = T5Z + T5Y;
		    T70 = T6W + T6X;
		    io[WS(ios, 30)] = T6Z - T70;
		    io[WS(ios, 1)] = T6Z + T70;
	       }
	       {
		    E T5X, T60, T6V, T6Y;
		    T5X = T5L - T5O;
		    T60 = T5Y - T5Z;
		    ro[WS(ros, 17)] = T5X - T60;
		    ro[WS(ros, 14)] = T5X + T60;
		    T6V = T5V - T5S;
		    T6Y = T6W - T6X;
		    io[WS(ios, 17)] = T6V - T6Y;
		    io[WS(ios, 14)] = T6V + T6Y;
	       }
	  }
	  {
	       E T37, T3z, T6n, T6t, T3e, T6s, T3J, T3M, T3m, T3x, T3C, T6k, T3G, T3N, T3t;
	       E T3w, T36, T6l;
	       T36 = T34 + T35;
	       T37 = T33 - T36;
	       T3z = T33 + T36;
	       T6l = Tv - Tm;
	       T6n = T6l - T6m;
	       T6t = T6l + T6m;
	       {
		    E T3a, T3d, T3H, T3I;
		    T3a = FMA(KP634393284, T38, KP773010453 * T39);
		    T3d = FNMS(KP634393284, T3c, KP773010453 * T3b);
		    T3e = T3a - T3d;
		    T6s = T3d + T3a;
		    T3H = T3n + T3o;
		    T3I = T3q + T3r;
		    T3J = FNMS(KP336889853, T3I, KP941544065 * T3H);
		    T3M = FMA(KP336889853, T3H, KP941544065 * T3I);
	       }
	       {
		    E T3i, T3l, T3A, T3B;
		    T3i = T3g - T3h;
		    T3l = T3j - T3k;
		    T3m = FMA(KP903989293, T3i, KP427555093 * T3l);
		    T3x = FNMS(KP427555093, T3i, KP903989293 * T3l);
		    T3A = FMA(KP773010453, T3c, KP634393284 * T3b);
		    T3B = FNMS(KP634393284, T39, KP773010453 * T38);
		    T3C = T3A + T3B;
		    T6k = T3B - T3A;
	       }
	       {
		    E T3E, T3F, T3p, T3s;
		    T3E = T3g + T3h;
		    T3F = T3j + T3k;
		    T3G = FMA(KP941544065, T3E, KP336889853 * T3F);
		    T3N = FNMS(KP336889853, T3E, KP941544065 * T3F);
		    T3p = T3n - T3o;
		    T3s = T3q - T3r;
		    T3t = FNMS(KP427555093, T3s, KP903989293 * T3p);
		    T3w = FMA(KP427555093, T3p, KP903989293 * T3s);
	       }
	       {
		    E T3f, T3u, T6p, T6q;
		    T3f = T37 + T3e;
		    T3u = T3m + T3t;
		    ro[WS(ros, 27)] = T3f - T3u;
		    ro[WS(ros, 4)] = T3f + T3u;
		    T6p = T3x + T3w;
		    T6q = T6k + T6n;
		    io[WS(ios, 27)] = T6p - T6q;
		    io[WS(ios, 4)] = T6p + T6q;
	       }
	       {
		    E T3v, T3y, T6j, T6o;
		    T3v = T37 - T3e;
		    T3y = T3w - T3x;
		    ro[WS(ros, 20)] = T3v - T3y;
		    ro[WS(ros, 11)] = T3v + T3y;
		    T6j = T3t - T3m;
		    T6o = T6k - T6n;
		    io[WS(ios, 20)] = T6j - T6o;
		    io[WS(ios, 11)] = T6j + T6o;
	       }
	       {
		    E T3D, T3K, T6v, T6w;
		    T3D = T3z + T3C;
		    T3K = T3G + T3J;
		    ro[WS(ros, 28)] = T3D - T3K;
		    ro[WS(ros, 3)] = T3D + T3K;
		    T6v = T3N + T3M;
		    T6w = T6s + T6t;
		    io[WS(ios, 28)] = T6v - T6w;
		    io[WS(ios, 3)] = T6v + T6w;
	       }
	       {
		    E T3L, T3O, T6r, T6u;
		    T3L = T3z - T3C;
		    T3O = T3M - T3N;
		    ro[WS(ros, 19)] = T3L - T3O;
		    ro[WS(ros, 12)] = T3L + T3O;
		    T6r = T3J - T3G;
		    T6u = T6s - T6t;
		    io[WS(ios, 19)] = T6r - T6u;
		    io[WS(ios, 12)] = T6r + T6u;
	       }
	  }
	  {
	       E T41, T4Z, T6D, T6J, T4g, T6I, T59, T5d, T4A, T4X, T52, T6y, T56, T5c, T4T;
	       E T4W, T40, T6C;
	       T40 = T3W - T3Z;
	       T41 = T3T + T40;
	       T4Z = T3T - T40;
	       T6C = T5g + T5h;
	       T6D = T6B - T6C;
	       T6J = T6C + T6B;
	       {
		    E T48, T4f, T57, T58;
		    T48 = FMA(KP881921264, T44, KP471396736 * T47);
		    T4f = FMA(KP881921264, T4b, KP471396736 * T4e);
		    T4g = T48 - T4f;
		    T6I = T48 + T4f;
		    T57 = T4K + T4H;
		    T58 = T4R + T4O;
		    T59 = FMA(KP514102744, T57, KP857728610 * T58);
		    T5d = FNMS(KP857728610, T57, KP514102744 * T58);
	       }
	       {
		    E T4s, T4z, T50, T51;
		    T4s = T4k + T4r;
		    T4z = T4v - T4y;
		    T4A = FMA(KP970031253, T4s, KP242980179 * T4z);
		    T4X = FNMS(KP242980179, T4s, KP970031253 * T4z);
		    T50 = FNMS(KP471396736, T4b, KP881921264 * T4e);
		    T51 = FNMS(KP471396736, T44, KP881921264 * T47);
		    T52 = T50 - T51;
		    T6y = T51 + T50;
	       }
	       {
		    E T54, T55, T4L, T4S;
		    T54 = T4k - T4r;
		    T55 = T4y + T4v;
		    T56 = FMA(KP514102744, T54, KP857728610 * T55);
		    T5c = FNMS(KP514102744, T55, KP857728610 * T54);
		    T4L = T4H - T4K;
		    T4S = T4O - T4R;
		    T4T = FNMS(KP242980179, T4S, KP970031253 * T4L);
		    T4W = FMA(KP242980179, T4L, KP970031253 * T4S);
	       }
	       {
		    E T4h, T4U, T6F, T6G;
		    T4h = T41 + T4g;
		    T4U = T4A + T4T;
		    ro[WS(ros, 29)] = T4h - T4U;
		    ro[WS(ros, 2)] = T4h + T4U;
		    T6F = T4X + T4W;
		    T6G = T6y + T6D;
		    io[WS(ios, 29)] = T6F - T6G;
		    io[WS(ios, 2)] = T6F + T6G;
	       }
	       {
		    E T4V, T4Y, T6x, T6E;
		    T4V = T41 - T4g;
		    T4Y = T4W - T4X;
		    ro[WS(ros, 18)] = T4V - T4Y;
		    ro[WS(ros, 13)] = T4V + T4Y;
		    T6x = T4T - T4A;
		    T6E = T6y - T6D;
		    io[WS(ios, 18)] = T6x - T6E;
		    io[WS(ios, 13)] = T6x + T6E;
	       }
	       {
		    E T53, T5a, T6L, T6M;
		    T53 = T4Z - T52;
		    T5a = T56 - T59;
		    ro[WS(ros, 21)] = T53 - T5a;
		    ro[WS(ros, 10)] = T53 + T5a;
		    T6L = T5d - T5c;
		    T6M = T6J - T6I;
		    io[WS(ios, 21)] = T6L - T6M;
		    io[WS(ios, 10)] = T6L + T6M;
	       }
	       {
		    E T5b, T5e, T6H, T6K;
		    T5b = T4Z + T52;
		    T5e = T5c + T5d;
		    ro[WS(ros, 26)] = T5b - T5e;
		    ro[WS(ros, 5)] = T5b + T5e;
		    T6H = T56 + T59;
		    T6K = T6I + T6J;
		    io[WS(ios, 5)] = -(T6H + T6K);
		    io[WS(ios, 26)] = T6K - T6H;
	       }
	  }
     }
}

static void mr2hcII_64(const R *I, R *ro, R *io, stride is, stride ros, stride ios, int v, int ivs, int ovs)
{
     int i;
     for (i = v; i > 0; --i) {
	  mr2hcII_64_0(I, ro, io, is, ros, ios);
	  I += ivs;
	  ro += ovs;
	  io += ovs;
     }
}

static const kr2hc_desc desc = { 64, "mr2hcII_64", {342, 114, 92, 0}, &GENUS, 0, 0, 0, 0, 0 };

void X(codelet_mr2hcII_64) (planner *p) {
     X(kr2hcII_register) (p, mr2hcII_64, &desc);
}