t2_32.c

fftw-3.0.1

			 E T45, T46, T48, T4b;
			 T45 = T43 - T44;
			 T46 = T20 - T27;
			 T47 = T45 + T46;
			 T5B = T45 - T46;
			 T48 = T1J - T1U;
			 T4b = T49 - T4a;
			 T4c = T48 - T4b;
			 T5C = T48 + T4b;
		    }
		    T4d = FMA(KP382683432, T47, KP923879532 * T4c);
		    T61 = FNMS(KP382683432, T5B, KP923879532 * T5C);
		    T5f = FNMS(KP923879532, T47, KP382683432 * T4c);
		    T5D = FMA(KP923879532, T5B, KP382683432 * T5C);
	       }
	       {
		    E T1E, T7e, T29, T6p, T6s, T7f;
		    T1E = T1s + T1D;
		    T7e = T6k + T6l;
		    {
			 E T1V, T28, T6q, T6r;
			 T1V = T1J + T1U;
			 T28 = T20 + T27;
			 T29 = T1V + T28;
			 T6p = T1V - T28;
			 T6q = T43 + T44;
			 T6r = T49 + T4a;
			 T6s = T6q - T6r;
			 T7f = T6q + T6r;
		    }
		    T2a = T1E + T29;
		    T6t = T6p + T6s;
		    T7O = T29 - T1E;
		    T7C = T7e + T7f;
		    T7g = T7e - T7f;
		    T6Z = T6p - T6s;
	       }
	       {
		    E T4k, T5J, T4B, T5G, T4v, T5H, T4E, T5K, T4j, T4A;
		    T4j = T2n - T2w;
		    T4k = T4i + T4j;
		    T5J = T4i - T4j;
		    T4A = T4y - T4z;
		    T4B = T4x - T4A;
		    T5G = T4x + T4A;
		    {
			 E T4p, T4u, T4C, T4D;
			 T4p = T4l - T4o;
			 T4u = T4q + T4t;
			 T4v = KP707106781 * (T4p - T4u);
			 T5H = KP707106781 * (T4u + T4p);
			 T4C = T4t - T4q;
			 T4D = T4l + T4o;
			 T4E = KP707106781 * (T4C - T4D);
			 T5K = KP707106781 * (T4C + T4D);
		    }
		    T4w = T4k - T4v;
		    T64 = T5G + T5H;
		    T65 = T5J + T5K;
		    T4F = T4B - T4E;
		    T5i = T4k + T4v;
		    T5I = T5G - T5H;
		    T5L = T5J - T5K;
		    T5j = T4B + T4E;
	       }
	       {
		    E T2y, T6B, T6y, T7j, T2R, T6z, T6E, T7k, T2x, T6x;
		    T2x = T2n + T2w;
		    T2y = T2k + T2x;
		    T6B = T2k - T2x;
		    T6x = T4y + T4z;
		    T6y = T6w - T6x;
		    T7j = T6w + T6x;
		    {
			 E T2H, T2Q, T6C, T6D;
			 T2H = T2B + T2G;
			 T2Q = T2M + T2P;
			 T2R = T2H + T2Q;
			 T6z = T2Q - T2H;
			 T6C = T4r + T4s;
			 T6D = T4m + T4n;
			 T6E = T6C - T6D;
			 T7k = T6C + T6D;
		    }
		    T2S = T2y + T2R;
		    T7l = T7j - T7k;
		    T7y = T7j + T7k;
		    T6A = T6y - T6z;
		    T6F = T6B - T6E;
		    T73 = T6B + T6E;
		    T7i = T2y - T2R;
		    T72 = T6y + T6z;
	       }
	       {
		    E T4L, T5N, T55, T5O, T4W, T5R, T52, T5Q;
		    {
			 E T4H, T4K, T53, T54;
			 T4H = T31 - T36;
			 T4K = T4I - T4J;
			 T4L = T4H - T4K;
			 T5N = T4H + T4K;
			 T53 = T4R - T4U;
			 T54 = T4P + T4O;
			 T55 = KP707106781 * (T53 - T54);
			 T5O = KP707106781 * (T54 + T53);
		    }
		    {
			 E T4Q, T4V, T50, T51;
			 T4Q = T4O - T4P;
			 T4V = T4R + T4U;
			 T4W = KP707106781 * (T4Q - T4V);
			 T5R = KP707106781 * (T4Q + T4V);
			 T50 = T4Y - T4Z;
			 T51 = T3a - T3f;
			 T52 = T50 + T51;
			 T5Q = T50 - T51;
		    }
		    T4X = T4L - T4W;
		    T67 = T5N + T5O;
		    T68 = T5Q + T5R;
		    T56 = T52 - T55;
		    T5l = T4L + T4W;
		    T5P = T5N - T5O;
		    T5S = T5Q - T5R;
		    T5m = T52 + T55;
	       }
	       {
		    E T3y, T6P, T6K, T7p, T3h, T6H, T6O, T7o, T3x, T6J;
		    T3x = T3t + T3w;
		    T3y = T3q + T3x;
		    T6P = T3x - T3q;
		    T6J = T4S + T4T;
		    T6K = T6I - T6J;
		    T7p = T6I + T6J;
		    {
			 E T37, T3g, T6M, T6N;
			 T37 = T31 + T36;
			 T3g = T3a + T3f;
			 T3h = T37 + T3g;
			 T6H = T37 - T3g;
			 T6M = T4Y + T4Z;
			 T6N = T4I + T4J;
			 T6O = T6M - T6N;
			 T7o = T6M + T6N;
		    }
		    T3z = T3h + T3y;
		    T7q = T7o - T7p;
		    T7z = T7o + T7p;
		    T6L = T6H - T6K;
		    T6Q = T6O - T6P;
		    T76 = T6O + T6P;
		    T7n = T3h - T3y;
		    T75 = T6H + T6K;
	       }
	       {
		    E T3A, T7A, T2b, T7x, T1j;
		    T3A = T2S + T3z;
		    T7A = T7y - T7z;
		    T1j = TL + T1i;
		    T2b = T1j + T2a;
		    T7x = T1j - T2a;
		    ri[WS(ios, 16)] = T2b - T3A;
		    ri[WS(ios, 8)] = T7x + T7A;
		    ri[0] = T2b + T3A;
		    ri[WS(ios, 24)] = T7x - T7A;
	       }
	       {
		    E T7B, T7L, T7K, T7M, T7J;
		    T7B = T7y + T7z;
		    T7L = T3z - T2S;
		    T7J = T7D + T7I;
		    T7K = T7C + T7J;
		    T7M = T7J - T7C;
		    ii[0] = T7B + T7K;
		    ii[WS(ios, 24)] = T7M - T7L;
		    ii[WS(ios, 16)] = T7K - T7B;
		    ii[WS(ios, 8)] = T7L + T7M;
	       }
	       {
		    E T7h, T7t, T7Q, T7S, T7s, T7R, T7w, T7N, T7d, T7P;
		    T7d = TL - T1i;
		    T7h = T7d + T7g;
		    T7t = T7d - T7g;
		    T7P = T7I - T7D;
		    T7Q = T7O + T7P;
		    T7S = T7P - T7O;
		    {
			 E T7m, T7r, T7u, T7v;
			 T7m = T7i + T7l;
			 T7r = T7n - T7q;
			 T7s = KP707106781 * (T7m + T7r);
			 T7R = KP707106781 * (T7r - T7m);
			 T7u = T7l - T7i;
			 T7v = T7n + T7q;
			 T7w = KP707106781 * (T7u - T7v);
			 T7N = KP707106781 * (T7u + T7v);
		    }
		    ri[WS(ios, 20)] = T7h - T7s;
		    ii[WS(ios, 20)] = T7Q - T7N;
		    ri[WS(ios, 4)] = T7h + T7s;
		    ii[WS(ios, 4)] = T7N + T7Q;
		    ri[WS(ios, 28)] = T7t - T7w;
		    ii[WS(ios, 28)] = T7S - T7R;
		    ri[WS(ios, 12)] = T7t + T7w;
		    ii[WS(ios, 12)] = T7R + T7S;
	       }
	       {
		    E T71, T79, T7Y, T80, T78, T7Z, T7c, T7T;
		    {
			 E T6X, T70, T7U, T7X;
			 T6X = T6f + T6i;
			 T70 = KP707106781 * (T6Y + T6Z);
			 T71 = T6X + T70;
			 T79 = T6X - T70;
			 T7U = KP707106781 * (T6o + T6t);
			 T7X = T7V + T7W;
			 T7Y = T7U + T7X;
			 T80 = T7X - T7U;
		    }
		    {
			 E T74, T77, T7a, T7b;
			 T74 = FMA(KP382683432, T72, KP923879532 * T73);
			 T77 = FNMS(KP382683432, T76, KP923879532 * T75);
			 T78 = T74 + T77;
			 T7Z = T77 - T74;
			 T7a = FNMS(KP382683432, T73, KP923879532 * T72);
			 T7b = FMA(KP923879532, T76, KP382683432 * T75);
			 T7c = T7a - T7b;
			 T7T = T7a + T7b;
		    }
		    ri[WS(ios, 18)] = T71 - T78;
		    ii[WS(ios, 18)] = T7Y - T7T;
		    ri[WS(ios, 2)] = T71 + T78;
		    ii[WS(ios, 2)] = T7T + T7Y;
		    ri[WS(ios, 26)] = T79 - T7c;
		    ii[WS(ios, 26)] = T80 - T7Z;
		    ri[WS(ios, 10)] = T79 + T7c;
		    ii[WS(ios, 10)] = T7Z + T80;
	       }
	       {
		    E T4f, T59, T8y, T8A, T58, T8z, T5c, T8v;
		    {
			 E T3R, T4e, T8w, T8x;
			 T3R = T3F - T3Q;
			 T4e = T42 - T4d;
			 T4f = T3R + T4e;
			 T59 = T3R - T4e;
			 T8w = T5f - T5e;
			 T8x = T8q - T8p;
			 T8y = T8w + T8x;
			 T8A = T8x - T8w;
		    }
		    {
			 E T4G, T57, T5a, T5b;
			 T4G = FMA(KP980785280, T4w, KP195090322 * T4F);
			 T57 = FNMS(KP980785280, T56, KP195090322 * T4X);
			 T58 = T4G + T57;
			 T8z = T57 - T4G;
			 T5a = FNMS(KP980785280, T4F, KP195090322 * T4w);
			 T5b = FMA(KP195090322, T56, KP980785280 * T4X);
			 T5c = T5a - T5b;
			 T8v = T5a + T5b;
		    }
		    ri[WS(ios, 23)] = T4f - T58;
		    ii[WS(ios, 23)] = T8y - T8v;
		    ri[WS(ios, 7)] = T4f + T58;
		    ii[WS(ios, 7)] = T8v + T8y;
		    ri[WS(ios, 31)] = T59 - T5c;
		    ii[WS(ios, 31)] = T8A - T8z;
		    ri[WS(ios, 15)] = T59 + T5c;
		    ii[WS(ios, 15)] = T8z + T8A;
	       }
	       {
		    E T5F, T5V, T8k, T8m, T5U, T8l, T5Y, T8h;
		    {
			 E T5x, T5E, T8i, T8j;
			 T5x = T5t - T5w;
			 T5E = T5A - T5D;
			 T5F = T5x + T5E;
			 T5V = T5x - T5E;
			 T8i = T61 - T60;
			 T8j = T8c - T89;
			 T8k = T8i + T8j;
			 T8m = T8j - T8i;
		    }
		    {
			 E T5M, T5T, T5W, T5X;
			 T5M = FMA(KP555570233, T5I, KP831469612 * T5L);
			 T5T = FNMS(KP831469612, T5S, KP555570233 * T5P);
			 T5U = T5M + T5T;
			 T8l = T5T - T5M;
			 T5W = FNMS(KP831469612, T5I, KP555570233 * T5L);
			 T5X = FMA(KP831469612, T5P, KP555570233 * T5S);
			 T5Y = T5W - T5X;
			 T8h = T5W + T5X;
		    }
		    ri[WS(ios, 21)] = T5F - T5U;
		    ii[WS(ios, 21)] = T8k - T8h;
		    ri[WS(ios, 5)] = T5F + T5U;
		    ii[WS(ios, 5)] = T8h + T8k;
		    ri[WS(ios, 29)] = T5V - T5Y;
		    ii[WS(ios, 29)] = T8m - T8l;
		    ri[WS(ios, 13)] = T5V + T5Y;
		    ii[WS(ios, 13)] = T8l + T8m;
	       }
	       {
		    E T6v, T6T, T84, T86, T6S, T85, T6W, T81;
		    {
			 E T6j, T6u, T82, T83;
			 T6j = T6f - T6i;
			 T6u = KP707106781 * (T6o - T6t);
			 T6v = T6j + T6u;
			 T6T = T6j - T6u;
			 T82 = KP707106781 * (T6Z - T6Y);
			 T83 = T7W - T7V;
			 T84 = T82 + T83;
			 T86 = T83 - T82;
		    }
		    {
			 E T6G, T6R, T6U, T6V;
			 T6G = FMA(KP923879532, T6A, KP382683432 * T6F);
			 T6R = FNMS(KP923879532, T6Q, KP382683432 * T6L);
			 T6S = T6G + T6R;
			 T85 = T6R - T6G;
			 T6U = FNMS(KP923879532, T6F, KP382683432 * T6A);
			 T6V = FMA(KP382683432, T6Q, KP923879532 * T6L);
			 T6W = T6U - T6V;
			 T81 = T6U + T6V;
		    }
		    ri[WS(ios, 22)] = T6v - T6S;
		    ii[WS(ios, 22)] = T84 - T81;
		    ri[WS(ios, 6)] = T6v + T6S;
		    ii[WS(ios, 6)] = T81 + T84;
		    ri[WS(ios, 30)] = T6T - T6W;
		    ii[WS(ios, 30)] = T86 - T85;
		    ri[WS(ios, 14)] = T6T + T6W;
		    ii[WS(ios, 14)] = T85 + T86;
	       }
	       {
		    E T5h, T5p, T8s, T8u, T5o, T8t, T5s, T8n;
		    {
			 E T5d, T5g, T8o, T8r;
			 T5d = T3F + T3Q;
			 T5g = T5e + T5f;
			 T5h = T5d + T5g;
			 T5p = T5d - T5g;
			 T8o = T42 + T4d;
			 T8r = T8p + T8q;
			 T8s = T8o + T8r;
			 T8u = T8r - T8o;
		    }
		    {
			 E T5k, T5n, T5q, T5r;
			 T5k = FMA(KP555570233, T5i, KP831469612 * T5j);
			 T5n = FNMS(KP555570233, T5m, KP831469612 * T5l);
			 T5o = T5k + T5n;
			 T8t = T5n - T5k;
			 T5q = FNMS(KP555570233, T5j, KP831469612 * T5i);
			 T5r = FMA(KP831469612, T5m, KP555570233 * T5l);
			 T5s = T5q - T5r;
			 T8n = T5q + T5r;
		    }
		    ri[WS(ios, 19)] = T5h - T5o;
		    ii[WS(ios, 19)] = T8s - T8n;
		    ri[WS(ios, 3)] = T5h + T5o;
		    ii[WS(ios, 3)] = T8n + T8s;
		    ri[WS(ios, 27)] = T5p - T5s;
		    ii[WS(ios, 27)] = T8u - T8t;
		    ri[WS(ios, 11)] = T5p + T5s;
		    ii[WS(ios, 11)] = T8t + T8u;
	       }
	       {
		    E T63, T6b, T8e, T8g, T6a, T8f, T6e, T87;
		    {
			 E T5Z, T62, T88, T8d;
			 T5Z = T5t + T5w;
			 T62 = T60 + T61;
			 T63 = T5Z + T62;
			 T6b = T5Z - T62;
			 T88 = T5A + T5D;
			 T8d = T89 + T8c;
			 T8e = T88 + T8d;
			 T8g = T8d - T88;
		    }
		    {
			 E T66, T69, T6c, T6d;
			 T66 = FMA(KP980785280, T64, KP195090322 * T65);
			 T69 = FNMS(KP195090322, T68, KP980785280 * T67);
			 T6a = T66 + T69;
			 T8f = T69 - T66;
			 T6c = FNMS(KP195090322, T64, KP980785280 * T65);
			 T6d = FMA(KP195090322, T67, KP980785280 * T68);
			 T6e = T6c - T6d;
			 T87 = T6c + T6d;
		    }
		    ri[WS(ios, 17)] = T63 - T6a;
		    ii[WS(ios, 17)] = T8e - T87;
		    ri[WS(ios, 1)] = T63 + T6a;
		    ii[WS(ios, 1)] = T87 + T8e;
		    ri[WS(ios, 25)] = T6b - T6e;
		    ii[WS(ios, 25)] = T8g - T8f;
		    ri[WS(ios, 9)] = T6b + T6e;
		    ii[WS(ios, 9)] = T8f + T8g;
	       }
	  }
     }
     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_NEXT, 1, 0}
};

static const ct_desc desc = { 32, "t2_32", twinstr, {376, 168, 112, 0}, &GENUS, 0, 0, 0 };

void X(codelet_t2_32) (planner *p) {
     X(kdft_dit_register) (p, t2_32, &desc);
}