hf_25.c

快速fft变换

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

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

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

void X(codelet_hf_25) (planner *p) {
     X(khc2hc_register) (p, hf_25, &desc);
}
#endif				/* HAVE_FMA */