hc2cf_32.c

快速fft变换

			      }
			      {
				   E T3a, T3d, T39, T3c, T58, T3b, T3g;
				   T3a = Ip[WS(rs, 9)];
				   T3d = Im[WS(rs, 9)];
				   T5g = FNMS(T3p, T3n, T5f);
				   T3r = FMA(T3p, T3q, T3o);
				   T39 = W[36];
				   T3c = W[37];
				   T3h = Ip[WS(rs, 13)];
				   T3k = Im[WS(rs, 13)];
				   T58 = T39 * T3d;
				   T3b = T39 * T3a;
				   T3g = W[52];
				   T3j = W[53];
				   T59 = FNMS(T3c, T3a, T58);
				   T3e = FMA(T3c, T3d, T3b);
				   T5d = T3g * T3k;
				   T3i = T3g * T3h;
			      }
			 }
			 {
			      E T5a, T7c, T3f, T55, T5e, T3l, T7d, T3s;
			      T5a = T57 - T59;
			      T7c = T57 + T59;
			      T3f = T38 + T3e;
			      T55 = T38 - T3e;
			      T5e = FNMS(T3j, T3h, T5d);
			      T3l = FMA(T3j, T3k, T3i);
			      T5h = T5e - T5g;
			      T7d = T5e + T5g;
			      T3s = T3l + T3r;
			      T5c = T3l - T3r;
			      T5s = T5a - T55;
			      T5b = T55 + T5a;
			      T7O = T7c + T7d;
			      T7e = T7c - T7d;
			      T76 = T3s - T3f;
			      T3t = T3f + T3s;
			 }
		    }
		    {
			 E T4y, T2f, T2o, T2r, T4A, T2l, T2n, T2q, T4E, T2p;
			 {
			      E T2u, T2x, T2t, T2w;
			      {
				   E T2b, T2e, T2d, T4x, T2c, T2a;
				   T2b = Ip[WS(rs, 2)];
				   T2e = Im[WS(rs, 2)];
				   T2a = W[8];
				   T5t = T5c + T5h;
				   T5i = T5c - T5h;
				   T2d = W[9];
				   T4x = T2a * T2e;
				   T2c = T2a * T2b;
				   T2u = Ip[WS(rs, 6)];
				   T2x = Im[WS(rs, 6)];
				   T4y = FNMS(T2d, T2b, T4x);
				   T2f = FMA(T2d, T2e, T2c);
				   T2t = W[24];
				   T2w = W[25];
			      }
			      {
				   E T2h, T2k, T2j, T4z, T2i, T4G, T2v, T2g;
				   T2h = Ip[WS(rs, 10)];
				   T2k = Im[WS(rs, 10)];
				   T4G = T2t * T2x;
				   T2v = T2t * T2u;
				   T2g = W[40];
				   T2j = W[41];
				   T4H = FNMS(T2w, T2u, T4G);
				   T2y = FMA(T2w, T2x, T2v);
				   T4z = T2g * T2k;
				   T2i = T2g * T2h;
				   T2o = Ip[WS(rs, 14)];
				   T2r = Im[WS(rs, 14)];
				   T4A = FNMS(T2j, T2h, T4z);
				   T2l = FMA(T2j, T2k, T2i);
				   T2n = W[56];
				   T2q = W[57];
			      }
			 }
			 T4B = T4y - T4A;
			 T71 = T4y + T4A;
			 T2m = T2f + T2l;
			 T4w = T2f - T2l;
			 T4E = T2n * T2r;
			 T2p = T2n * T2o;
			 T4F = FNMS(T2q, T2o, T4E);
			 T2s = FMA(T2q, T2r, T2p);
		    }
	       }
	       {
		    E T4T, T4C, T4J, T4U, T7y, T8q, T8p, T7B;
		    {
			 E T6E, T8j, T73, T6Y, T6H, T8k, T8i, T8h;
			 {
			      E T7C, TO, T80, T7Z, T8e, T89, T8d, T1H, T8b, T3v, T7T, T7L, T7U, T7Q, T2A;
			      E T7K, T7P, T7W, T1I;
			      {
				   E T7X, T7Y, T7J, T82, T88;
				   {
					E Tm, T4I, T72, T4D, T2z, TN;
					T6E = T8 - Tl;
					Tm = T8 + Tl;
					T4T = T4B - T4w;
					T4C = T4w + T4B;
					T4I = T4F - T4H;
					T72 = T4F + T4H;
					T4D = T2s - T2y;
					T2z = T2s + T2y;
					TN = Tz + TM;
					T8j = TM - Tz;
					T73 = T71 - T72;
					T7J = T71 + T72;
					T4J = T4D - T4I;
					T4U = T4D + T4I;
					T2A = T2m + T2z;
					T6Y = T2z - T2m;
					T7C = Tm - TN;
					TO = Tm + TN;
				   }
				   T7K = T7I - T7J;
				   T7X = T7I + T7J;
				   T7Y = T7N + T7O;
				   T7P = T7N - T7O;
				   T6H = T6F - T6G;
				   T82 = T6F + T6G;
				   T88 = T83 + T87;
				   T8k = T87 - T83;
				   T80 = T7X + T7Y;
				   T7Z = T7X - T7Y;
				   T8e = T88 - T82;
				   T89 = T82 + T88;
			      }
			      {
				   E T7H, T7M, T2B, T3u;
				   T7H = T29 - T2A;
				   T2B = T29 + T2A;
				   T3u = T32 + T3t;
				   T7M = T32 - T3t;
				   T8d = T1G - T1f;
				   T1H = T1f + T1G;
				   T8b = T3u - T2B;
				   T3v = T2B + T3u;
				   T7T = T7K - T7H;
				   T7L = T7H + T7K;
				   T7U = T7M + T7P;
				   T7Q = T7M - T7P;
			      }
			      T7W = TO - T1H;
			      T1I = TO + T1H;
			      {
				   E T7S, T8f, T8g, T7V;
				   {
					E T7R, T8c, T8a, T7G, T81, T7F;
					T8i = T7Q - T7L;
					T7R = T7L + T7Q;
					T81 = T7D + T7E;
					T7F = T7D - T7E;
					Rp[0] = T1I + T3v;
					Rm[WS(rs, 15)] = T1I - T3v;
					Rp[WS(rs, 8)] = T7W + T7Z;
					Rm[WS(rs, 7)] = T7W - T7Z;
					T8c = T89 - T81;
					T8a = T81 + T89;
					T7G = T7C + T7F;
					T7S = T7C - T7F;
					T8h = T8e - T8d;
					T8f = T8d + T8e;
					Ip[WS(rs, 8)] = T8b + T8c;
					Im[WS(rs, 7)] = T8b - T8c;
					Ip[0] = T80 + T8a;
					Im[WS(rs, 15)] = T80 - T8a;
					Rp[WS(rs, 4)] = FMA(KP707106781, T7R, T7G);
					Rm[WS(rs, 11)] = FNMS(KP707106781, T7R, T7G);
					T8g = T7T + T7U;
					T7V = T7T - T7U;
				   }
				   Ip[WS(rs, 4)] = FMA(KP707106781, T8g, T8f);
				   Im[WS(rs, 11)] = FMS(KP707106781, T8g, T8f);
				   Rp[WS(rs, 12)] = FMA(KP707106781, T7V, T7S);
				   Rm[WS(rs, 3)] = FNMS(KP707106781, T7V, T7S);
			      }
			 }
			 {
			      E T7f, T7m, T6I, T7a, T7A, T7w, T8r, T8l, T8m, T6T, T7j, T75, T8s, T7p, T7z;
			      E T7t;
			      {
				   E T7n, T6N, T6S, T7o, T7u, T7v;
				   T7f = T7b - T7e;
				   T7u = T7b + T7e;
				   Ip[WS(rs, 12)] = FMA(KP707106781, T8i, T8h);
				   Im[WS(rs, 3)] = FMS(KP707106781, T8i, T8h);
				   T7m = T6E + T6H;
				   T6I = T6E - T6H;
				   T7v = T79 + T76;
				   T7a = T76 - T79;
				   T7n = T6M + T6L;
				   T6N = T6L - T6M;
				   T7A = FMA(KP414213562, T7u, T7v);
				   T7w = FNMS(KP414213562, T7v, T7u);
				   T8r = T8k - T8j;
				   T8l = T8j + T8k;
				   T6S = T6O + T6R;
				   T7o = T6O - T6R;
				   {
					E T7s, T7r, T6Z, T74;
					T7s = T6X + T6Y;
					T6Z = T6X - T6Y;
					T74 = T70 - T73;
					T7r = T70 + T73;
					T8m = T6N + T6S;
					T6T = T6N - T6S;
					T7j = FNMS(KP414213562, T6Z, T74);
					T75 = FMA(KP414213562, T74, T6Z);
					T8s = T7o - T7n;
					T7p = T7n + T7o;
					T7z = FNMS(KP414213562, T7r, T7s);
					T7t = FMA(KP414213562, T7s, T7r);
				   }
			      }
			      {
				   E T7i, T6U, T8t, T8v, T7k, T7g;
				   T7i = FNMS(KP707106781, T6T, T6I);
				   T6U = FMA(KP707106781, T6T, T6I);
				   T8t = FMA(KP707106781, T8s, T8r);
				   T8v = FNMS(KP707106781, T8s, T8r);
				   T7k = FNMS(KP414213562, T7a, T7f);
				   T7g = FMA(KP414213562, T7f, T7a);
				   {
					E T7q, T7x, T8n, T8o;
					T7y = FNMS(KP707106781, T7p, T7m);
					T7q = FMA(KP707106781, T7p, T7m);
					{
					     E T7l, T8u, T8w, T7h;
					     T7l = T7j + T7k;
					     T8u = T7k - T7j;
					     T8w = T7g - T75;
					     T7h = T75 + T7g;
					     Rm[WS(rs, 1)] = FMA(KP923879532, T7l, T7i);
					     Rp[WS(rs, 14)] = FNMS(KP923879532, T7l, T7i);
					     Ip[WS(rs, 6)] = FMA(KP923879532, T8u, T8t);
					     Im[WS(rs, 9)] = FMS(KP923879532, T8u, T8t);
					     Ip[WS(rs, 14)] = FMA(KP923879532, T8w, T8v);
					     Im[WS(rs, 1)] = FMS(KP923879532, T8w, T8v);
					     Rp[WS(rs, 6)] = FMA(KP923879532, T7h, T6U);
					     Rm[WS(rs, 9)] = FNMS(KP923879532, T7h, T6U);
					     T7x = T7t + T7w;
					     T8q = T7w - T7t;
					}
					T8p = FNMS(KP707106781, T8m, T8l);
					T8n = FMA(KP707106781, T8m, T8l);
					T8o = T7z + T7A;
					T7B = T7z - T7A;
					Rp[WS(rs, 2)] = FMA(KP923879532, T7x, T7q);
					Rm[WS(rs, 13)] = FNMS(KP923879532, T7x, T7q);
					Ip[WS(rs, 2)] = FMA(KP923879532, T8o, T8n);
					Im[WS(rs, 13)] = FMS(KP923879532, T8o, T8n);
				   }
			      }
			 }
		    }
		    {
			 E T5S, T8O, T8N, T5V, T6g, T6d, T69, T66, T8K, T8J;
			 {
			      E T5C, T3S, T8I, T4n, T8H, T8B, T8C, T5F, T5k, T5K, T5u, T4K, T4V;
			      {
				   E T5D, T5E, T8z, T8A, T5j;
				   {
					E T3C, T3R, T47, T4m;
					T5S = T3w - T3B;
					T3C = T3w + T3B;
					Rp[WS(rs, 10)] = FMA(KP923879532, T7B, T7y);
					Rm[WS(rs, 5)] = FNMS(KP923879532, T7B, T7y);
					Ip[WS(rs, 10)] = FMA(KP923879532, T8q, T8p);
					Im[WS(rs, 5)] = FMS(KP923879532, T8q, T8p);
					T3R = T3J + T3Q;
					T8O = T3Q - T3J;
					T5D = FNMS(KP414213562, T3Z, T46);
					T47 = FMA(KP414213562, T46, T3Z);
					T4m = FNMS(KP414213562, T4l, T4e);
					T5E = FMA(KP414213562, T4e, T4l);
					T8N = T8y + T8x;
					T8z = T8x - T8y;
					T5C = FNMS(KP707106781, T3R, T3C);
					T3S = FMA(KP707106781, T3R, T3C);
					T8I = T4m - T47;
					T4n = T47 + T4m;
					T8A = T5T + T5U;
					T5V = T5T - T5U;
				   }
				   T6g = T5i - T5b;
				   T5j = T5b + T5i;
				   T8H = FNMS(KP707106781, T8A, T8z);
				   T8B = FMA(KP707106781, T8A, T8z);
				   T8C = T5D + T5E;
				   T5F = T5D - T5E;
				   T5k = FMA(KP707106781, T5j, T54);
				   T5K = FNMS(KP707106781, T5j, T54);
				   T5u = T5s + T5t;
				   T6d = T5t - T5s;
				   T69 = T4C - T4J;
				   T4K = T4C + T4J;
				   T4V = T4T + T4U;
				   T66 = T4U - T4T;
			      }
			      {
				   E T5M, T5Q, T5J, T5P, T8F, T8G;
				   {
					E T5y, T4o, T5A, T5w, T5z, T4X, T8D, T5L, T5v, T8E, T5B, T5x;
					T5y = FNMS(KP923879532, T4n, T3S);
					T4o = FMA(KP923879532, T4n, T3S);
					T5L = FNMS(KP707106781, T5u, T5r);
					T5v = FMA(KP707106781, T5u, T5r);
					{
					     E T5H, T4L, T5I, T4W;
					     T5H = FNMS(KP707106781, T4K, T4v);
					     T4L = FMA(KP707106781, T4K, T4v);
					     T5I = FNMS(KP707106781, T4V, T4S);
					     T4W = FMA(KP707106781, T4V, T4S);
					     T5M = FMA(KP668178637, T5L, T5K);
					     T5Q = FNMS(KP668178637, T5K, T5L);
					     T5A = FMA(KP198912367, T5k, T5v);
					     T5w = FNMS(KP198912367, T5v, T5k);
					     T5J = FNMS(KP668178637, T5I, T5H);
					     T5P = FMA(KP668178637, T5H, T5I);
					     T5z = FNMS(KP198912367, T4L, T4W);
					     T4X = FMA(KP198912367, T4W, T4L);
					}
					T8D = FMA(KP923879532, T8C, T8B);
					T8F = FNMS(KP923879532, T8C, T8B);
					T8E = T5z + T5A;
					T5B = T5z - T5A;
					T8G = T5w - T4X;
					T5x = T4X + T5w;
					Ip[WS(rs, 1)] = FMA(KP980785280, T8E, T8D);
					Im[WS(rs, 14)] = FMS(KP980785280, T8E, T8D);
					Rp[WS(rs, 1)] = FMA(KP980785280, T5x, T4o);
					Rm[WS(rs, 14)] = FNMS(KP980785280, T5x, T4o);
					Rp[WS(rs, 9)] = FMA(KP980785280, T5B, T5y);
					Rm[WS(rs, 6)] = FNMS(KP980785280, T5B, T5y);
				   }
				   {
					E T5O, T8L, T8M, T5R, T5G, T5N;
					T5O = FMA(KP923879532, T5F, T5C);
					T5G = FNMS(KP923879532, T5F, T5C);
					T5N = T5J + T5M;
					T8K = T5M - T5J;
					T8J = FMA(KP923879532, T8I, T8H);
					T8L = FNMS(KP923879532, T8I, T8H);
					Ip[WS(rs, 9)] = FMA(KP980785280, T8G, T8F);
					Im[WS(rs, 6)] = FMS(KP980785280, T8G, T8F);
					Rm[WS(rs, 2)] = FMA(KP831469612, T5N, T5G);
					Rp[WS(rs, 13)] = FNMS(KP831469612, T5N, T5G);
					T8M = T5P + T5Q;
					T5R = T5P - T5Q;
					Ip[WS(rs, 13)] = FNMS(KP831469612, T8M, T8L);
					Im[WS(rs, 2)] = -(FMA(KP831469612, T8M, T8L));
					Rp[WS(rs, 5)] = FMA(KP831469612, T5R, T5O);
					Rm[WS(rs, 10)] = FNMS(KP831469612, T5R, T5O);
				   }
			      }
			 }
			 {
			      E T6o, T5W, T8W, T63, T8V, T8P, T8Q, T6r, T67, T6u, T6y, T6C, T6m, T6i;
			      {
				   E T6p, T5Z, T62, T6q;
				   T6p = FNMS(KP414213562, T5X, T5Y);
				   T5Z = FMA(KP414213562, T5Y, T5X);
				   Ip[WS(rs, 5)] = FMA(KP831469612, T8K, T8J);
				   Im[WS(rs, 10)] = FMS(KP831469612, T8K, T8J);
				   T6o = FNMS(KP707106781, T5V, T5S);
				   T5W = FMA(KP707106781, T5V, T5S);
				   T62 = FNMS(KP414213562, T61, T60);
				   T6q = FMA(KP414213562, T60, T61);
				   T8W = T5Z + T62;
				   T63 = T5Z - T62;
				   T8V = FNMS(KP707106781, T8O, T8N);
				   T8P = FMA(KP707106781, T8O, T8N);
				   {
					E T6x, T6e, T6w, T6h;
					T8Q = T6q - T6p;
					T6r = T6p + T6q;
					T6x = FMA(KP707106781, T6d, T6c);
					T6e = FNMS(KP707106781, T6d, T6c);
					T6w = FMA(KP707106781, T6g, T6f);
					T6h = FNMS(KP707106781, T6g, T6f);
					T67 = FNMS(KP707106781, T66, T65);
					T6u = FMA(KP707106781, T66, T65);
					T6y = FMA(KP198912367, T6x, T6w);
					T6C = FNMS(KP198912367, T6w, T6x);
					T6m = FNMS(KP668178637, T6e, T6h);
					T6i = FMA(KP668178637, T6h, T6e);
				   }
			      }
			      {
				   E T6k, T64, T8R, T8T, T6t, T6a;
				   T6k = FNMS(KP923879532, T63, T5W);
				   T64 = FMA(KP923879532, T63, T5W);
				   T8R = FMA(KP923879532, T8Q, T8P);
				   T8T = FNMS(KP923879532, T8Q, T8P);
				   T6t = FMA(KP707106781, T69, T68);
				   T6a = FNMS(KP707106781, T69, T68);
				   {
					E T6A, T8X, T8Y, T6D;
					{
					     E T6s, T6B, T6l, T6b, T6z, T6v;
					     T6A = FMA(KP923879532, T6r, T6o);
					     T6s = FNMS(KP923879532, T6r, T6o);
					     T6v = FMA(KP198912367, T6u, T6t);
					     T6B = FNMS(KP198912367, T6t, T6u);
					     T6l = FNMS(KP668178637, T67, T6a);
					     T6b = FMA(KP668178637, T6a, T67);
					     T6z = T6v + T6y;
					     T90 = T6y - T6v;
					     T8Z = FMA(KP923879532, T8W, T8V);
					     T8X = FNMS(KP923879532, T8W, T8V);
					     {
						  E T6n, T8S, T8U, T6j;
						  T6n = T6l + T6m;
						  T8S = T6l - T6m;
						  T8U = T6i - T6b;
						  T6j = T6b + T6i;
						  Rp[WS(rs, 7)] = FMA(KP980785280, T6z, T6s);
						  Rm[WS(rs, 8)] = FNMS(KP980785280, T6z, T6s);
						  Rp[WS(rs, 11)] = FMA(KP831469612, T6n, T6k);
						  Rm[WS(rs, 4)] = FNMS(KP831469612, T6n, T6k);
						  Ip[WS(rs, 3)] = FMA(KP831469612, T8S, T8R);
						  Im[WS(rs, 12)] = FMS(KP831469612, T8S, T8R);
						  Ip[WS(rs, 11)] = FMA(KP831469612, T8U, T8T);
						  Im[WS(rs, 4)] = FMS(KP831469612, T8U, T8T);
						  Rp[WS(rs, 3)] = FMA(KP831469612, T6j, T64);
						  Rm[WS(rs, 12)] = FNMS(KP831469612, T6j, T64);
						  T8Y = T6C - T6B;
						  T6D = T6B + T6C;
					     }
					}
					Ip[WS(rs, 7)] = FMA(KP980785280, T8Y, T8X);
					Im[WS(rs, 8)] = FMS(KP980785280, T8Y, T8X);
					Rm[0] = FMA(KP980785280, T6D, T6A);
					Rp[WS(rs, 15)] = FNMS(KP980785280, T6D, T6A);