t1_25.c

快速fft变换

		    T36 = KP559016994 * (T34 - T35);
		    T39 = FNMS(KP250000000, T38, T37);
		    T3a = T36 + T39;
		    T4S = T39 - T36;
		    T3e = T3a - T3d;
		    T5v = T4S - T4R;
		    T4k = T3d + T3a;
		    T4T = T4R + T4S;
	       }
	  }
	  {
	       E T1r, T3v, T3m, T3p, T3A, T3z, T3s, T3t, T3w, T1C, T1N, T1O;
	       {
		    E T1o, T1q, T1n, T1p;
		    T1o = ri[WS(rs, 2)];
		    T1q = ii[WS(rs, 2)];
		    T1n = W[2];
		    T1p = W[3];
		    T1r = FMA(T1n, T1o, T1p * T1q);
		    T3v = FNMS(T1p, T1o, T1n * T1q);
	       }
	       {
		    E T1w, T3k, T1M, T3o, T1B, T3l, T1H, T3n;
		    {
			 E T1t, T1v, T1s, T1u;
			 T1t = ri[WS(rs, 7)];
			 T1v = ii[WS(rs, 7)];
			 T1s = W[12];
			 T1u = W[13];
			 T1w = FMA(T1s, T1t, T1u * T1v);
			 T3k = FNMS(T1u, T1t, T1s * T1v);
		    }
		    {
			 E T1J, T1L, T1I, T1K;
			 T1J = ri[WS(rs, 17)];
			 T1L = ii[WS(rs, 17)];
			 T1I = W[32];
			 T1K = W[33];
			 T1M = FMA(T1I, T1J, T1K * T1L);
			 T3o = FNMS(T1K, T1J, T1I * T1L);
		    }
		    {
			 E T1y, T1A, T1x, T1z;
			 T1y = ri[WS(rs, 22)];
			 T1A = ii[WS(rs, 22)];
			 T1x = W[42];
			 T1z = W[43];
			 T1B = FMA(T1x, T1y, T1z * T1A);
			 T3l = FNMS(T1z, T1y, T1x * T1A);
		    }
		    {
			 E T1E, T1G, T1D, T1F;
			 T1E = ri[WS(rs, 12)];
			 T1G = ii[WS(rs, 12)];
			 T1D = W[22];
			 T1F = W[23];
			 T1H = FMA(T1D, T1E, T1F * T1G);
			 T3n = FNMS(T1F, T1E, T1D * T1G);
		    }
		    T3m = T3k - T3l;
		    T3p = T3n - T3o;
		    T3A = T1H - T1M;
		    T3z = T1w - T1B;
		    T3s = T3k + T3l;
		    T3t = T3n + T3o;
		    T3w = T3s + T3t;
		    T1C = T1w + T1B;
		    T1N = T1H + T1M;
		    T1O = T1C + T1N;
	       }
	       T1P = T1r + T1O;
	       T5Z = T3v + T3w;
	       {
		    E T3q, T50, T3j, T4Z, T3h, T3i;
		    T3q = FMA(KP951056516, T3m, KP587785252 * T3p);
		    T50 = FNMS(KP587785252, T3m, KP951056516 * T3p);
		    T3h = KP559016994 * (T1C - T1N);
		    T3i = FNMS(KP250000000, T1O, T1r);
		    T3j = T3h + T3i;
		    T4Z = T3i - T3h;
		    T3r = T3j + T3q;
		    T5z = T4Z + T50;
		    T4o = T3j - T3q;
		    T51 = T4Z - T50;
	       }
	       {
		    E T3B, T52, T3y, T53, T3u, T3x;
		    T3B = FMA(KP951056516, T3z, KP587785252 * T3A);
		    T52 = FNMS(KP587785252, T3z, KP951056516 * T3A);
		    T3u = KP559016994 * (T3s - T3t);
		    T3x = FNMS(KP250000000, T3w, T3v);
		    T3y = T3u + T3x;
		    T53 = T3x - T3u;
		    T3C = T3y - T3B;
		    T5A = T53 - T52;
		    T4n = T3B + T3y;
		    T54 = T52 + T53;
	       }
	  }
	  {
	       E T62, T64, Tp, T2k, T5T, T5U, T63, T5V;
	       {
		    E T5Y, T61, T1m, T2j;
		    T5Y = T5W - T5X;
		    T61 = T5Z - T60;
		    T62 = FMA(KP951056516, T5Y, KP587785252 * T61);
		    T64 = FNMS(KP587785252, T5Y, KP951056516 * T61);
		    Tp = T1 + To;
		    T1m = TS + T1l;
		    T2j = T1P + T2i;
		    T2k = T1m + T2j;
		    T5T = KP559016994 * (T1m - T2j);
		    T5U = FNMS(KP250000000, T2k, Tp);
	       }
	       ri[0] = Tp + T2k;
	       T63 = T5U - T5T;
	       ri[WS(rs, 10)] = T63 - T64;
	       ri[WS(rs, 15)] = T63 + T64;
	       T5V = T5T + T5U;
	       ri[WS(rs, 20)] = T5V - T62;
	       ri[WS(rs, 5)] = T5V + T62;
	  }
	  {
	       E T6i, T6j, T6c, T67, T6d, T6e, T6k, T6f;
	       {
		    E T6g, T6h, T65, T66;
		    T6g = TS - T1l;
		    T6h = T1P - T2i;
		    T6i = FMA(KP951056516, T6g, KP587785252 * T6h);
		    T6j = FNMS(KP587785252, T6g, KP951056516 * T6h);
		    T6c = T6a + T6b;
		    T65 = T5W + T5X;
		    T66 = T5Z + T60;
		    T67 = T65 + T66;
		    T6d = KP559016994 * (T65 - T66);
		    T6e = FNMS(KP250000000, T67, T6c);
	       }
	       ii[0] = T67 + T6c;
	       T6k = T6e - T6d;
	       ii[WS(rs, 10)] = T6j + T6k;
	       ii[WS(rs, 15)] = T6k - T6j;
	       T6f = T6d + T6e;
	       ii[WS(rs, 5)] = T6f - T6i;
	       ii[WS(rs, 20)] = T6i + T6f;
	  }
	  {
	       E T2v, T4f, T6u, T6G, T42, T6z, T43, T6y, T4A, T6H, T4D, T6F, T4u, T6L, T4v;
	       E T6K, T48, T6v, T4b, T6n, T2n, T6q;
	       T2n = T2l + T2m;
	       T2v = T2n + T2u;
	       T4f = T2n - T2u;
	       T6q = T6o + T6p;
	       T6u = T6q - T6t;
	       T6G = T6t + T6q;
	       {
		    E T2S, T3f, T3g, T3D, T40, T41;
		    T2S = FMA(KP968583161, T2G, KP248689887 * T2R);
		    T3f = FMA(KP535826794, T33, KP844327925 * T3e);
		    T3g = T2S + T3f;
		    T3D = FMA(KP876306680, T3r, KP481753674 * T3C);
		    T40 = FMA(KP728968627, T3O, KP684547105 * T3Z);
		    T41 = T3D + T40;
		    T42 = T3g + T41;
		    T6z = T3D - T40;
		    T43 = KP559016994 * (T3g - T41);
		    T6y = T2S - T3f;
	       }
	       {
		    E T4y, T4z, T6D, T4B, T4C, T6E;
		    T4y = FNMS(KP844327925, T4g, KP535826794 * T4h);
		    T4z = FNMS(KP637423989, T4k, KP770513242 * T4j);
		    T6D = T4y + T4z;
		    T4B = FMA(KP125333233, T4r, KP992114701 * T4q);
		    T4C = FMA(KP904827052, T4o, KP425779291 * T4n);
		    T6E = T4C + T4B;
		    T4A = T4y - T4z;
		    T6H = KP559016994 * (T6D + T6E);
		    T4D = T4B - T4C;
		    T6F = T6D - T6E;
	       }
	       {
		    E T4i, T4l, T4m, T4p, T4s, T4t;
		    T4i = FMA(KP535826794, T4g, KP844327925 * T4h);
		    T4l = FMA(KP637423989, T4j, KP770513242 * T4k);
		    T4m = T4i - T4l;
		    T4p = FNMS(KP425779291, T4o, KP904827052 * T4n);
		    T4s = FNMS(KP992114701, T4r, KP125333233 * T4q);
		    T4t = T4p + T4s;
		    T4u = T4m + T4t;
		    T6L = T4p - T4s;
		    T4v = KP559016994 * (T4m - T4t);
		    T6K = T4i + T4l;
	       }
	       {
		    E T46, T47, T6l, T49, T4a, T6m;
		    T46 = FNMS(KP248689887, T2G, KP968583161 * T2R);
		    T47 = FNMS(KP844327925, T33, KP535826794 * T3e);
		    T6l = T46 + T47;
		    T49 = FNMS(KP481753674, T3r, KP876306680 * T3C);
		    T4a = FNMS(KP684547105, T3O, KP728968627 * T3Z);
		    T6m = T49 + T4a;
		    T48 = T46 - T47;
		    T6v = KP559016994 * (T6l - T6m);
		    T4b = T49 - T4a;
		    T6n = T6l + T6m;
	       }
	       ri[WS(rs, 1)] = T2v + T42;
	       ii[WS(rs, 1)] = T6n + T6u;
	       ri[WS(rs, 4)] = T4f + T4u;
	       ii[WS(rs, 4)] = T6F + T6G;
	       {
		    E T4c, T4e, T45, T4d, T44;
		    T4c = FMA(KP951056516, T48, KP587785252 * T4b);
		    T4e = FNMS(KP587785252, T48, KP951056516 * T4b);
		    T44 = FNMS(KP250000000, T42, T2v);
		    T45 = T43 + T44;
		    T4d = T44 - T43;
		    ri[WS(rs, 21)] = T45 - T4c;
		    ri[WS(rs, 16)] = T4d + T4e;
		    ri[WS(rs, 6)] = T45 + T4c;
		    ri[WS(rs, 11)] = T4d - T4e;
	       }
	       {
		    E T6A, T6B, T6x, T6C, T6w;
		    T6A = FMA(KP951056516, T6y, KP587785252 * T6z);
		    T6B = FNMS(KP587785252, T6y, KP951056516 * T6z);
		    T6w = FNMS(KP250000000, T6n, T6u);
		    T6x = T6v + T6w;
		    T6C = T6w - T6v;
		    ii[WS(rs, 6)] = T6x - T6A;
		    ii[WS(rs, 16)] = T6C - T6B;
		    ii[WS(rs, 21)] = T6A + T6x;
		    ii[WS(rs, 11)] = T6B + T6C;
	       }
	       {
		    E T4E, T4G, T4x, T4F, T4w;
		    T4E = FMA(KP951056516, T4A, KP587785252 * T4D);
		    T4G = FNMS(KP587785252, T4A, KP951056516 * T4D);
		    T4w = FNMS(KP250000000, T4u, T4f);
		    T4x = T4v + T4w;
		    T4F = T4w - T4v;
		    ri[WS(rs, 24)] = T4x - T4E;
		    ri[WS(rs, 19)] = T4F + T4G;
		    ri[WS(rs, 9)] = T4x + T4E;
		    ri[WS(rs, 14)] = T4F - T4G;
	       }
	       {
		    E T6M, T6N, T6J, T6O, T6I;
		    T6M = FMA(KP951056516, T6K, KP587785252 * T6L);
		    T6N = FNMS(KP587785252, T6K, KP951056516 * T6L);
		    T6I = FNMS(KP250000000, T6F, T6G);
		    T6J = T6H + T6I;
		    T6O = T6I - T6H;
		    ii[WS(rs, 9)] = T6J - T6M;
		    ii[WS(rs, 19)] = T6O - T6N;
		    ii[WS(rs, 24)] = T6M + T6J;
		    ii[WS(rs, 14)] = T6N + T6O;
	       }
	  }
	  {
	       E T4J, T5r, T6U, T76, T5e, T6Z, T5f, T6Y, T5M, T77, T5P, T75, T5G, T7b, T5H;
	       E T7a, T5k, T6V, T5n, T6R, T4H, T6T;
	       T4H = T2m - T2l;
	       T4J = T4H - T4I;
	       T5r = T4H + T4I;
	       T6T = T6p - T6o;
	       T6U = T6S + T6T;
	       T76 = T6T - T6S;
	       {
		    E T4Q, T4X, T4Y, T55, T5c, T5d;
		    T4Q = FMA(KP876306680, T4M, KP481753674 * T4P);
		    T4X = FNMS(KP425779291, T4W, KP904827052 * T4T);
		    T4Y = T4Q + T4X;
		    T55 = FMA(KP535826794, T51, KP844327925 * T54);
		    T5c = FMA(KP062790519, T58, KP998026728 * T5b);
		    T5d = T55 + T5c;
		    T5e = T4Y + T5d;
		    T6Z = T55 - T5c;
		    T5f = KP559016994 * (T4Y - T5d);
		    T6Y = T4Q - T4X;
	       }
	       {
		    E T5K, T5L, T73, T5N, T5O, T74;
		    T5K = FNMS(KP684547105, T5s, KP728968627 * T5t);
		    T5L = FMA(KP125333233, T5w, KP992114701 * T5v);
		    T73 = T5K - T5L;
		    T5N = FNMS(KP998026728, T5z, KP062790519 * T5A);
		    T5O = FMA(KP770513242, T5D, KP637423989 * T5C);
		    T74 = T5N - T5O;
		    T5M = T5K + T5L;
		    T77 = KP559016994 * (T73 - T74);
		    T5P = T5N + T5O;
		    T75 = T73 + T74;
	       }
	       {
		    E T5u, T5x, T5y, T5B, T5E, T5F;
		    T5u = FMA(KP728968627, T5s, KP684547105 * T5t);
		    T5x = FNMS(KP992114701, T5w, KP125333233 * T5v);
		    T5y = T5u + T5x;
		    T5B = FMA(KP062790519, T5z, KP998026728 * T5A);
		    T5E = FNMS(KP637423989, T5D, KP770513242 * T5C);
		    T5F = T5B + T5E;
		    T5G = T5y + T5F;
		    T7b = T5B - T5E;
		    T5H = KP559016994 * (T5y - T5F);
		    T7a = T5u - T5x;
	       }
	       {
		    E T5i, T5j, T6P, T5l, T5m, T6Q;
		    T5i = FNMS(KP481753674, T4M, KP876306680 * T4P);
		    T5j = FMA(KP904827052, T4W, KP425779291 * T4T);
		    T6P = T5i - T5j;
		    T5l = FNMS(KP844327925, T51, KP535826794 * T54);
		    T5m = FNMS(KP998026728, T58, KP062790519 * T5b);
		    T6Q = T5l + T5m;
		    T5k = T5i + T5j;
		    T6V = KP559016994 * (T6P - T6Q);
		    T5n = T5l - T5m;
		    T6R = T6P + T6Q;
	       }
	       ri[WS(rs, 2)] = T4J + T5e;
	       ii[WS(rs, 2)] = T6R + T6U;
	       ri[WS(rs, 3)] = T5r + T5G;
	       ii[WS(rs, 3)] = T75 + T76;
	       {
		    E T5o, T5q, T5h, T5p, T5g;
		    T5o = FMA(KP951056516, T5k, KP587785252 * T5n);
		    T5q = FNMS(KP587785252, T5k, KP951056516 * T5n);
		    T5g = FNMS(KP250000000, T5e, T4J);
		    T5h = T5f + T5g;
		    T5p = T5g - T5f;
		    ri[WS(rs, 22)] = T5h - T5o;
		    ri[WS(rs, 17)] = T5p + T5q;
		    ri[WS(rs, 7)] = T5h + T5o;
		    ri[WS(rs, 12)] = T5p - T5q;
	       }
	       {
		    E T70, T71, T6X, T72, T6W;
		    T70 = FMA(KP951056516, T6Y, KP587785252 * T6Z);
		    T71 = FNMS(KP587785252, T6Y, KP951056516 * T6Z);
		    T6W = FNMS(KP250000000, T6R, T6U);
		    T6X = T6V + T6W;
		    T72 = T6W - T6V;
		    ii[WS(rs, 7)] = T6X - T70;
		    ii[WS(rs, 17)] = T72 - T71;
		    ii[WS(rs, 22)] = T70 + T6X;
		    ii[WS(rs, 12)] = T71 + T72;
	       }
	       {
		    E T5Q, T5S, T5J, T5R, T5I;
		    T5Q = FMA(KP951056516, T5M, KP587785252 * T5P);
		    T5S = FNMS(KP587785252, T5M, KP951056516 * T5P);
		    T5I = FNMS(KP250000000, T5G, T5r);
		    T5J = T5H + T5I;
		    T5R = T5I - T5H;
		    ri[WS(rs, 23)] = T5J - T5Q;
		    ri[WS(rs, 18)] = T5R + T5S;
		    ri[WS(rs, 8)] = T5J + T5Q;
		    ri[WS(rs, 13)] = T5R - T5S;
	       }
	       {
		    E T7c, T7d, T79, T7e, T78;
		    T7c = FMA(KP951056516, T7a, KP587785252 * T7b);
		    T7d = FNMS(KP587785252, T7a, KP951056516 * T7b);
		    T78 = FNMS(KP250000000, T75, T76);
		    T79 = T77 + T78;
		    T7e = T78 - T77;
		    ii[WS(rs, 8)] = T79 - T7c;
		    ii[WS(rs, 18)] = T7e - T7d;
		    ii[WS(rs, 23)] = T7c + T79;
		    ii[WS(rs, 13)] = T7d + T7e;
	       }
	  }
     }
}

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

static const ct_desc desc = { 25, "t1_25", twinstr, &GENUS, {260, 140, 140, 0}, 0, 0, 0 };

void X(codelet_t1_25) (planner *p) {
     X(kdft_dit_register) (p, t1_25, &desc);
}
#endif				/* HAVE_FMA */