hf_32.c

最新的FFT程序

					T3p = W[21];
				   }
				   {
					E T3a, T3d, T3c, T58, T3b, T5f, T3o, T39;
					T3a = rio[WS(ios, 19)];
					T3d = iio[-WS(ios, 12)];
					T5f = T3m * T3q;
					T3o = T3m * T3n;
					T39 = W[36];
					T3c = W[37];
					T5g = FNMS(T3p, T3n, T5f);
					T3r = FMA(T3p, T3q, T3o);
					T58 = T39 * T3d;
					T3b = T39 * T3a;
					T3h = rio[WS(ios, 27)];
					T3k = iio[-WS(ios, 4)];
					T59 = FNMS(T3c, T3a, T58);
					T3e = FMA(T3c, T3d, T3b);
					T3g = W[52];
					T3j = W[53];
				   }
			      }
			      T5a = T57 - T59;
			      T7c = T57 + T59;
			      T3f = T38 + T3e;
			      T55 = T38 - T3e;
			      T5d = T3g * T3k;
			      T3i = T3g * T3h;
			      T5e = FNMS(T3j, T3h, T5d);
			      T3l = FMA(T3j, T3k, T3i);
			 }
			 {
			      E T5m, T2H, T52, T30, T2Q, T2T, T2S, T5o, T2N, T4Z, T2R;
			      {
				   E T2W, T2Z, T2Y, T51, T2X;
				   {
					E T2D, T2G, T2C, T2F, T5l, T2E, T2V;
					{
					     E T5h, T7d, T5c, T3s;
					     T2D = rio[WS(ios, 31)];
					     T5s = T5a - T55;
					     T5b = T55 + T5a;
					     T5h = T5e - T5g;
					     T7d = T5e + T5g;
					     T5c = T3l - T3r;
					     T3s = T3l + T3r;
					     T2G = iio[0];
					     T7e = T7c - T7d;
					     T7O = T7c + T7d;
					     T5i = T5c - T5h;
					     T5t = T5c + T5h;
					     T3t = T3f + T3s;
					     T76 = T3s - T3f;
					     T2C = W[60];
					}
					T2F = W[61];
					T2W = rio[WS(ios, 23)];
					T2Z = iio[-WS(ios, 8)];
					T5l = T2C * T2G;
					T2E = T2C * T2D;
					T2V = W[44];
					T2Y = W[45];
					T5m = FNMS(T2F, T2D, T5l);
					T2H = FMA(T2F, T2G, T2E);
					T51 = T2V * T2Z;
					T2X = T2V * T2W;
				   }
				   {
					E T2J, T2M, T2I, T2L, T5n, T2K, T2P;
					T2J = rio[WS(ios, 15)];
					T2M = iio[-WS(ios, 16)];
					T52 = FNMS(T2Y, T2W, T51);
					T30 = FMA(T2Y, T2Z, T2X);
					T2I = W[28];
					T2L = W[29];
					T2Q = rio[WS(ios, 7)];
					T2T = iio[-WS(ios, 24)];
					T5n = T2I * T2M;
					T2K = T2I * T2J;
					T2P = W[12];
					T2S = W[13];
					T5o = FNMS(T2L, T2J, T5n);
					T2N = FMA(T2L, T2M, T2K);
					T4Z = T2P * T2T;
					T2R = T2P * T2Q;
				   }
			      }
			      {
				   E T5p, T77, T2O, T4Y, T50, T2U;
				   T5p = T5m - T5o;
				   T77 = T5m + T5o;
				   T2O = T2H + T2N;
				   T4Y = T2H - T2N;
				   T50 = FNMS(T2S, T2Q, T4Z);
				   T2U = FMA(T2S, T2T, T2R);
				   {
					E T78, T53, T5q, T31;
					T78 = T50 + T52;
					T53 = T50 - T52;
					T5q = T30 - T2U;
					T31 = T2U + T30;
					T54 = T4Y + T53;
					T6c = T4Y - T53;
					T79 = T77 - T78;
					T7N = T77 + T78;
					T32 = T2O + T31;
					T7b = T2O - T31;
					T6f = T5q - T5p;
					T5r = T5p + T5q;
				   }
			      }
			 }
		    }
	       }
	       {
		    E T6E, T8j, T6H, T8k, T7S, T8i, T8h, T7V;
		    {
			 E T7P, T7C, TO, T80, T7Z, T8e, T89, T8d, T1H, T7M, T8b, T3v, T7T, T7L;
			 {
			      E T7K, T82, T88, T7H, T2B, T3u;
			      {
				   E Tm, TN, T7X, T7Y;
				   T6E = T8 - Tl;
				   Tm = T8 + Tl;
				   TN = Tz + TM;
				   T8j = TM - Tz;
				   T7K = T7I - T7J;
				   T7X = T7I + T7J;
				   T7Y = T7N + T7O;
				   T7P = T7N - T7O;
				   T6H = T6F - T6G;
				   T82 = T6F + T6G;
				   T7C = Tm - TN;
				   TO = Tm + TN;
				   T80 = T7X + T7Y;
				   T7Z = T7X - T7Y;
				   T88 = T83 + T87;
				   T8k = T87 - T83;
			      }
			      T7H = T29 - T2A;
			      T2B = T29 + T2A;
			      T8e = T88 - T82;
			      T89 = T82 + T88;
			      T8d = T1G - T1f;
			      T1H = T1f + T1G;
			      T3u = T32 + T3t;
			      T7M = T32 - T3t;
			      T8b = T3u - T2B;
			      T3v = T2B + T3u;
			      T7T = T7K - T7H;
			      T7L = T7H + T7K;
			 }
			 {
			      E T7W, T7U, T7Q, T7F, T8a, T8c, T81, T1I;
			      T7W = TO - T1H;
			      T1I = TO + T1H;
			      T7U = T7M + T7P;
			      T7Q = T7M - T7P;
			      T81 = T7D + T7E;
			      T7F = T7D - T7E;
			      iio[-WS(ios, 16)] = T1I - T3v;
			      rio[0] = T1I + T3v;
			      T8a = T81 + T89;
			      T8c = T89 - T81;
			      {
				   E T7G, T7R, T8f, T8g;
				   T7S = T7C - T7F;
				   T7G = T7C + T7F;
				   iio[0] = T80 + T8a;
				   rio[WS(ios, 16)] = T80 - T8a;
				   rio[WS(ios, 8)] = T7W + T7Z;
				   iio[-WS(ios, 24)] = T7W - T7Z;
				   iio[-WS(ios, 8)] = T8b + T8c;
				   rio[WS(ios, 24)] = T8b - T8c;
				   T7R = T7L + T7Q;
				   T8i = T7Q - T7L;
				   T8h = T8e - T8d;
				   T8f = T8d + T8e;
				   T8g = T7T + T7U;
				   T7V = T7T - T7U;
				   rio[WS(ios, 4)] = FMA(KP707106781, T7R, T7G);
				   iio[-WS(ios, 20)] = FNMS(KP707106781, T7R, T7G);
				   iio[-WS(ios, 4)] = FMA(KP707106781, T8g, T8f);
				   rio[WS(ios, 20)] = FMS(KP707106781, T8g, T8f);
			      }
			 }
		    }
		    {
			 E T7f, T7a, T7m, T6I, T7A, T7w, T7s, T8r, T8l, T8m, T6T, T7r, T7j, T75, T8s;
			 E T7p;
			 {
			      E T7n, T6N, T6S, T7o, T7u, T7v, T6Z, T74;
			      rio[WS(ios, 12)] = FMA(KP707106781, T7V, T7S);
			      iio[-WS(ios, 28)] = FNMS(KP707106781, T7V, T7S);
			      iio[-WS(ios, 12)] = FMA(KP707106781, T8i, T8h);
			      rio[WS(ios, 28)] = FMS(KP707106781, T8i, T8h);
			      T7f = T7b - T7e;
			      T7u = T7b + T7e;
			      T7v = T79 + T76;
			      T7a = T76 - T79;
			      T7m = T6E + T6H;
			      T6I = T6E - T6H;
			      T7A = FMA(KP414213562, T7u, T7v);
			      T7w = FNMS(KP414213562, T7v, T7u);
			      T7n = T6M + T6L;
			      T6N = T6L - T6M;
			      T6S = T6O + T6R;
			      T7o = T6O - T6R;
			      T7s = T6X + T6Y;
			      T6Z = T6X - T6Y;
			      T8r = T8k - T8j;
			      T8l = T8j + T8k;
			      T8m = T6N + T6S;
			      T6T = T6N - T6S;
			      T74 = T70 - T73;
			      T7r = T70 + T73;
			      T7j = FNMS(KP414213562, T6Z, T74);
			      T75 = FMA(KP414213562, T74, T6Z);
			      T8s = T7o - T7n;
			      T7p = T7n + T7o;
			 }
			 {
			      E T7z, T7t, T8v, T8w;
			      {
				   E T7i, T6U, T7k, T7g, T8t, T8u, T7l, T7h;
				   T7i = FNMS(KP707106781, T6T, T6I);
				   T6U = FMA(KP707106781, T6T, T6I);
				   T7z = FNMS(KP414213562, T7r, T7s);
				   T7t = FMA(KP414213562, T7s, T7r);
				   T7k = FNMS(KP414213562, T7a, T7f);
				   T7g = FMA(KP414213562, T7f, T7a);
				   T8t = FMA(KP707106781, T8s, T8r);
				   T8v = FNMS(KP707106781, T8s, T8r);
				   T8u = T7k - T7j;
				   T7l = T7j + T7k;
				   T8w = T7g - T75;
				   T7h = T75 + T7g;
				   iio[-WS(ios, 6)] = FMA(KP923879532, T8u, T8t);
				   rio[WS(ios, 22)] = FMS(KP923879532, T8u, T8t);
				   rio[WS(ios, 6)] = FMA(KP923879532, T7h, T6U);
				   iio[-WS(ios, 22)] = FNMS(KP923879532, T7h, T6U);
				   iio[-WS(ios, 30)] = FMA(KP923879532, T7l, T7i);
				   rio[WS(ios, 14)] = FNMS(KP923879532, T7l, T7i);
			      }
			      {
				   E T7y, T8n, T8o, T7B, T7q, T7x;
				   T7y = FNMS(KP707106781, T7p, T7m);
				   T7q = FMA(KP707106781, T7p, T7m);
				   T7x = T7t + T7w;
				   T8q = T7w - T7t;
				   T8p = FNMS(KP707106781, T8m, T8l);
				   T8n = FMA(KP707106781, T8m, T8l);
				   iio[-WS(ios, 14)] = FMA(KP923879532, T8w, T8v);
				   rio[WS(ios, 30)] = FMS(KP923879532, T8w, T8v);
				   rio[WS(ios, 2)] = FMA(KP923879532, T7x, T7q);
				   iio[-WS(ios, 18)] = FNMS(KP923879532, T7x, T7q);
				   T8o = T7z + T7A;
				   T7B = T7z - T7A;
				   iio[-WS(ios, 2)] = FMA(KP923879532, T8o, T8n);
				   rio[WS(ios, 18)] = FMS(KP923879532, T8o, T8n);
				   rio[WS(ios, 10)] = FMA(KP923879532, T7B, T7y);
				   iio[-WS(ios, 26)] = FNMS(KP923879532, T7B, T7y);
			      }
			 }
		    }
	       }
	       {
		    E T5S, T8O, T8N, T5V, T6g, T6d, T69, T66, T8K, T8J;
		    {
			 E T5C, T3S, T8I, T4n, T8H, T8B, T8C, 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);
			      iio[-WS(ios, 10)] = FMA(KP923879532, T8q, T8p);
			      rio[WS(ios, 26)] = FMS(KP923879532, T8q, T8p);
			      T5C = FNMS(KP707106781, T3R, T3C);
			      T3S = FMA(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;
			      T8I = T4m - T47;
			      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 = T4C - T4J;
				   T4K = T4C + T4J;
				   T8C = T5D + T5E;
				   T5F = T5D - T5E;
				   T5K = FNMS(KP707106781, T5j, T54);
				   T5k = FMA(KP707106781, T5j, T54);
				   T5L = FNMS(KP707106781, T5u, T5r);
				   T5v = FMA(KP707106781, T5u, T5r);
				   T4L = FMA(KP707106781, T4K, T4v);
				   T5H = FNMS(KP707106781, T4K, T4v);
				   T5M = FMA(KP668178637, T5L, T5K);
				   T5Q = FNMS(KP668178637, T5K, T5L);
				   T5A = FMA(KP198912367, T5k, T5v);
				   T5w = FNMS(KP198912367, T5v, T5k);
				   T4V = T4T + T4U;
				   T66 = T4U - T4T;
			      }
			 }
			 {
			      E T5y, T4o, T8D, T8F, T5I, T4W;
			      T5y = FNMS(KP923879532, T4n, T3S);
			      T4o = FMA(KP923879532, T4n, T3S);
			      T8D = FMA(KP923879532, T8C, T8B);
			      T8F = FNMS(KP923879532, T8C, T8B);
			      T5I = FNMS(KP707106781, T4V, T4S);
			      T4W = FMA(KP707106781, T4V, T4S);
			      {
				   E T5O, T8L, T8M, T5R;
				   {
					E T5G, T5P, T5z, T4X, T5N, T5J;
					T5O = FMA(KP923879532, T5F, T5C);
					T5G = FNMS(KP923879532, T5F, T5C);
					T5J = FNMS(KP668178637, T5I, T5H);
					T5P = FMA(KP668178637, T5H, T5I);
					T5z = FNMS(KP198912367, T4L, T4W);
					T4X = FMA(KP198912367, T4W, T4L);
					T5N = T5J + T5M;
					T8K = T5M - T5J;
					T8J = FMA(KP923879532, T8I, T8H);
					T8L = FNMS(KP923879532, T8I, T8H);
					{
					     E T5B, T8E, T8G, T5x;
					     T5B = T5z - T5A;
					     T8E = T5z + T5A;
					     T8G = T5w - T4X;
					     T5x = T4X + T5w;
					     iio[-WS(ios, 29)] = FMA(KP831469612, T5N, T5G);
					     rio[WS(ios, 13)] = FNMS(KP831469612, T5N, T5G);
					     rio[WS(ios, 9)] = FMA(KP980785280, T5B, T5y);
					     iio[-WS(ios, 25)] = FNMS(KP980785280, T5B, T5y);
					     iio[-WS(ios, 1)] = FMA(KP980785280, T8E, T8D);
					     rio[WS(ios, 17)] = FMS(KP980785280, T8E, T8D);
					     iio[-WS(ios, 9)] = FMA(KP980785280, T8G, T8F);
					     rio[WS(ios, 25)] = FMS(KP980785280, T8G, T8F);
					     rio[WS(ios, 1)] = FMA(KP980785280, T5x, T4o);
					     iio[-WS(ios, 17)] = FNMS(KP980785280, T5x, T4o);
					     T8M = T5P + T5Q;
					     T5R = T5P - T5Q;
					}
				   }
				   iio[-WS(ios, 13)] = FNMS(KP831469612, T8M, T8L);
				   rio[WS(ios, 29)] = -(FMA(KP831469612, T8M, T8L));
				   rio[WS(ios, 5)] = FMA(KP831469612, T5R, T5O);
				   iio[-WS(ios, 21)] = 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);
			      iio[-WS(ios, 5)] = FMA(KP831469612, T8K, T8J);
			      rio[WS(ios, 21)] = 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;
					     rio[WS(ios, 7)] = FMA(KP980785280, T6z, T6s);
					     iio[-WS(ios, 23)] = FNMS(KP980785280, T6z, T6s);
					     rio[WS(ios, 11)] = FMA(KP831469612, T6n, T6k);
					     iio[-WS(ios, 27)] = FNMS(KP831469612, T6n, T6k);
					     iio[-WS(ios, 3)] = FMA(KP831469612, T8S, T8R);
					     rio[WS(ios, 19)] = FMS(KP831469612, T8S, T8R);
					     iio[-WS(ios, 11)] = FMA(KP831469612, T8U, T8T);
					     rio[WS(ios, 27)] = FMS(KP831469612, T8U, T8T);
					     rio[WS(ios, 3)] = FMA(KP831469612, T6j, T64);
					     iio[-WS(ios, 19)] = FNMS(KP831469612, T6j, T64);
					     T8Y = T6C - T6B;
					     T6D = T6B + T6C;
					}
				   }
				   iio[-WS(ios, 7)] = FMA(KP980785280, T8Y, T8X);
				   rio[WS(ios, 23)] = FMS(KP980785280, T8Y, T8X);
				   iio[-WS(ios, 31)] = FMA(KP980785280, T6D, T6A);
				   rio[WS(ios, 15)] = FNMS(KP980785280, T6D, T6A);
			      }
			 }