hf_32.c

快速fft变换

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