hf2_32.c

最新的FFT程序

					T5b = T58 - T5a;
					T5e = T34 * T32;
					T33 = T31 * T32;
				   }
				   {
					E T3N, T3K, T3O, T5H, T5Q, T41, T3R, T3U, T3W;
					{
					     E T3H, T3I, T3J, T40, T3Y, T5G, T5P, T3Z;
					     {
						  E T5i, T7C, T5d, T3b, T5f, T36;
						  T3H = rio[WS(ios, 3)];
						  T5f = FMA(T31, T35, T5e);
						  T36 = FNMS(T34, T35, T33);
						  T5t = T5b - T56;
						  T5c = T56 + T5b;
						  T5i = T5f - T5h;
						  T7C = T5f + T5h;
						  T5d = T36 - T3a;
						  T3b = T36 + T3a;
						  T3I = T4 * T3H;
						  T7D = T7B - T7C;
						  T8j = T7B + T7C;
						  T5j = T5d - T5i;
						  T5u = T5d + T5i;
						  T3c = T30 + T3b;
						  T7y = T3b - T30;
						  T3J = iio[-WS(ios, 28)];
					     }
					     T40 = iio[-WS(ios, 20)];
					     T3Y = rio[WS(ios, 11)];
					     T3N = rio[WS(ios, 19)];
					     T3K = FMA(T7, T3J, T3I);
					     T5G = T4 * T3J;
					     T5P = Ti * T3Y;
					     T3Z = Te * T3Y;
					     T3O = T3M * T3N;
					     T5H = FNMS(T7, T3H, T5G);
					     T5Q = FMA(Te, T40, T5P);
					     T41 = FNMS(Ti, T40, T3Z);
					     T3R = iio[-WS(ios, 12)];
					     T3U = rio[WS(ios, 27)];
					     T3W = iio[-WS(ios, 4)];
					}
					{
					     E T3l, T5W, T5C, T3E, T3v, T3x, T3t, T5Y, T3w, T5z;
					     {
						  E T3g, T3h, T3k, T3D, T3A, T3s, T3p, T5X;
						  {
						       E T3T, T5F, T5J, T3X, T5O, T7M, T5K;
						       T3g = rio[WS(ios, 31)];
						       {
							    E T3S, T5I, T3V, T5N;
							    T3S = FMA(T3Q, T3R, T3O);
							    T5I = T3M * T3R;
							    T3V = T2 * T3U;
							    T5N = T2 * T3W;
							    T3T = T3K + T3S;
							    T5F = T3K - T3S;
							    T5J = FNMS(T3Q, T3N, T5I);
							    T3X = FMA(Tg, T3W, T3V);
							    T5O = FNMS(Tg, T3U, T5N);
							    T3h = T3f * T3g;
						       }
						       T7M = T5H + T5J;
						       T5K = T5H - T5J;
						       {
							    E T42, T5M, T7N, T5R;
							    T42 = T3X + T41;
							    T5M = T3X - T41;
							    T7N = T5O + T5Q;
							    T5R = T5O - T5Q;
							    T5L = T5F + T5K;
							    T62 = T5K - T5F;
							    T43 = T3T + T42;
							    T7G = T42 - T3T;
							    T5S = T5M - T5R;
							    T63 = T5M + T5R;
							    T7O = T7M - T7N;
							    T8o = T7M + T7N;
							    T3k = iio[0];
						       }
						  }
						  T3D = iio[-WS(ios, 8)];
						  T3A = rio[WS(ios, 23)];
						  T3s = iio[-WS(ios, 16)];
						  T3l = FMA(T3j, T3k, T3h);
						  {
						       E T5V, T5B, T3B, T3o;
						       T5V = T3f * T3k;
						       T5B = T3C * T3A;
						       T3B = T3z * T3A;
						       T3o = rio[WS(ios, 15)];
						       T5W = FNMS(T3j, T3g, T5V);
						       T5C = FMA(T3z, T3D, T5B);
						       T3E = FNMS(T3C, T3D, T3B);
						       T3p = T3n * T3o;
						       T5X = T3r * T3o;
						  }
						  T3v = rio[WS(ios, 7)];
						  T3x = iio[-WS(ios, 24)];
						  T3t = FNMS(T3r, T3s, T3p);
						  T5Y = FMA(T3n, T3s, T5X);
						  T3w = TP * T3v;
						  T5z = TP * T3x;
					     }
					     {
						  E T5y, T3u, T7H, T5Z, T3y, T5A;
						  T5y = T3l - T3t;
						  T3u = T3l + T3t;
						  T7H = T5W + T5Y;
						  T5Z = T5W - T5Y;
						  T3y = FMA(TR, T3x, T3w);
						  T5A = FNMS(TR, T3v, T5z);
						  {
						       E T60, T3F, T5D, T7I;
						       T60 = T3E - T3y;
						       T3F = T3y + T3E;
						       T5D = T5A - T5C;
						       T7I = T5A + T5C;
						       T3G = T3u + T3F;
						       T7L = T3u - T3F;
						       T5E = T5y + T5D;
						       T6M = T5y - T5D;
						       T7J = T7H - T7I;
						       T8n = T7H + T7I;
						       T6P = T60 - T5Z;
						       T61 = T5Z + T60;
						  }
					     }
					}
				   }
			      }
			 }
		    }
		    {
			 E T7e, T8T, T7h, T8U, T8s, T8S, T8R, T8v;
			 {
			      E T8p, T8c, T1i, T8A, T8z, T8O, T8J, T8N, T2l, T8m, T8L, T45, T8t, T8l;
			      {
				   E T8k, T8C, T8I, T8h, T3d, T44;
				   {
					E TI, T1h, T8x, T8y;
					T7e = Tq - TH;
					TI = Tq + TH;
					T1h = TZ + T1g;
					T8T = T1g - TZ;
					T8k = T8i - T8j;
					T8x = T8i + T8j;
					T8y = T8n + T8o;
					T8p = T8n - T8o;
					T7h = T7f - T7g;
					T8C = T7f + T7g;
					T8c = TI - T1h;
					T1i = TI + T1h;
					T8A = T8x + T8y;
					T8z = T8x - T8y;
					T8I = T8D + T8H;
					T8U = T8H - T8D;
				   }
				   T8h = T2N - T3c;
				   T3d = T2N + T3c;
				   T8O = T8I - T8C;
				   T8J = T8C + T8I;
				   T8N = T2k - T1J;
				   T2l = T1J + T2k;
				   T44 = T3G + T43;
				   T8m = T3G - T43;
				   T8L = T44 - T3d;
				   T45 = T3d + T44;
				   T8t = T8k - T8h;
				   T8l = T8h + T8k;
			      }
			      {
				   E T8w, T8u, T8q, T8f, T8K, T8M, T8B, T2m;
				   T8w = T1i - T2l;
				   T2m = T1i + T2l;
				   T8u = T8m + T8p;
				   T8q = T8m - T8p;
				   T8B = T8d + T8e;
				   T8f = T8d - T8e;
				   iio[-WS(ios, 16)] = T2m - T45;
				   rio[0] = T2m + T45;
				   T8K = T8B + T8J;
				   T8M = T8J - T8B;
				   {
					E T8g, T8r, T8P, T8Q;
					T8s = T8c - T8f;
					T8g = T8c + T8f;
					iio[0] = T8A + T8K;
					rio[WS(ios, 16)] = T8A - T8K;
					rio[WS(ios, 8)] = T8w + T8z;
					iio[-WS(ios, 24)] = T8w - T8z;
					iio[-WS(ios, 8)] = T8L + T8M;
					rio[WS(ios, 24)] = T8L - T8M;
					T8r = T8l + T8q;
					T8S = T8q - T8l;
					T8R = T8O - T8N;
					T8P = T8N + T8O;
					T8Q = T8t + T8u;
					T8v = T8t - T8u;
					rio[WS(ios, 4)] = FMA(KP707106781, T8r, T8g);
					iio[-WS(ios, 20)] = FNMS(KP707106781, T8r, T8g);
					iio[-WS(ios, 4)] = FMA(KP707106781, T8Q, T8P);
					rio[WS(ios, 20)] = FMS(KP707106781, T8Q, T8P);
				   }
			      }
			 }
			 {
			      E T7P, T7K, T7W, T7i, T8a, T86, T82, T91, T8V, T8W, T7t, T81, T7T, T7F, T92;
			      E T7Z;
			      {
				   E T7X, T7n, T7s, T7Y, T84, T85, T7z, T7E;
				   rio[WS(ios, 12)] = FMA(KP707106781, T8v, T8s);
				   iio[-WS(ios, 28)] = FNMS(KP707106781, T8v, T8s);
				   iio[-WS(ios, 12)] = FMA(KP707106781, T8S, T8R);
				   rio[WS(ios, 28)] = FMS(KP707106781, T8S, T8R);
				   T7P = T7L - T7O;
				   T84 = T7L + T7O;
				   T85 = T7J + T7G;
				   T7K = T7G - T7J;
				   T7W = T7e + T7h;
				   T7i = T7e - T7h;
				   T8a = FMA(KP414213562, T84, T85);
				   T86 = FNMS(KP414213562, T85, T84);
				   T7X = T7m + T7l;
				   T7n = T7l - T7m;
				   T7s = T7o + T7r;
				   T7Y = T7o - T7r;
				   T82 = T7x + T7y;
				   T7z = T7x - T7y;
				   T91 = T8U - T8T;
				   T8V = T8T + T8U;
				   T8W = T7n + T7s;
				   T7t = T7n - T7s;
				   T7E = T7A - T7D;
				   T81 = T7A + T7D;
				   T7T = FNMS(KP414213562, T7z, T7E);
				   T7F = FMA(KP414213562, T7E, T7z);
				   T92 = T7Y - T7X;
				   T7Z = T7X + T7Y;
			      }
			      {
				   E T89, T83, T95, T96;
				   {
					E T7S, T7u, T7U, T7Q, T93, T94, T7V, T7R;
					T7S = FNMS(KP707106781, T7t, T7i);
					T7u = FMA(KP707106781, T7t, T7i);
					T89 = FNMS(KP414213562, T81, T82);
					T83 = FMA(KP414213562, T82, T81);
					T7U = FNMS(KP414213562, T7K, T7P);
					T7Q = FMA(KP414213562, T7P, T7K);
					T93 = FMA(KP707106781, T92, T91);
					T95 = FNMS(KP707106781, T92, T91);
					T94 = T7U - T7T;
					T7V = T7T + T7U;
					T96 = T7Q - T7F;
					T7R = T7F + T7Q;
					iio[-WS(ios, 6)] = FMA(KP923879532, T94, T93);
					rio[WS(ios, 22)] = FMS(KP923879532, T94, T93);
					rio[WS(ios, 6)] = FMA(KP923879532, T7R, T7u);
					iio[-WS(ios, 22)] = FNMS(KP923879532, T7R, T7u);
					iio[-WS(ios, 30)] = FMA(KP923879532, T7V, T7S);
					rio[WS(ios, 14)] = FNMS(KP923879532, T7V, T7S);
				   }
				   {
					E T88, T8X, T8Y, T8b, T80, T87;
					T88 = FNMS(KP707106781, T7Z, T7W);
					T80 = FMA(KP707106781, T7Z, T7W);
					T87 = T83 + T86;
					T90 = T86 - T83;
					T8Z = FNMS(KP707106781, T8W, T8V);
					T8X = FMA(KP707106781, T8W, T8V);
					iio[-WS(ios, 14)] = FMA(KP923879532, T96, T95);
					rio[WS(ios, 30)] = FMS(KP923879532, T96, T95);
					rio[WS(ios, 2)] = FMA(KP923879532, T87, T80);
					iio[-WS(ios, 18)] = FNMS(KP923879532, T87, T80);
					T8Y = T89 + T8a;
					T8b = T89 - T8a;
					iio[-WS(ios, 2)] = FMA(KP923879532, T8Y, T8X);
					rio[WS(ios, 18)] = FMS(KP923879532, T8Y, T8X);
					rio[WS(ios, 10)] = FMA(KP923879532, T8b, T88);
					iio[-WS(ios, 26)] = FNMS(KP923879532, T8b, T88);
				   }
			      }
			 }
		    }
		    {
			 E T6s, T9o, T9n, T6v, T6Q, T6N, T6J, T6G, T9k, T9j;
			 {
			      E T6c, T4s, T9i, T4X, T9h, T9b, T9c, T6f, T5l, T6h, T6m, T6q, T6a, T66, T5v;
			      {
				   E T6d, T4H, T4W, T6e, T99, T9a, T4c, T4r, T5T, T64;
				   T6s = T46 - T4b;
				   T4c = T46 + T4b;
				   T4r = T4j + T4q;
				   T9o = T4q - T4j;
				   T6d = FNMS(KP414213562, T4z, T4G);
				   T4H = FMA(KP414213562, T4G, T4z);
				   iio[-WS(ios, 10)] = FMA(KP923879532, T90, T8Z);
				   rio[WS(ios, 26)] = FMS(KP923879532, T90, T8Z);
				   T6c = FNMS(KP707106781, T4r, T4c);
				   T4s = FMA(KP707106781, T4r, T4c);
				   T4W = FNMS(KP414213562, T4V, T4O);
				   T6e = FMA(KP414213562, T4O, T4V);
				   T9n = T98 + T97;
				   T99 = T97 - T98;
				   T9a = T6t + T6u;
				   T6v = T6t - T6u;
				   T9i = T4W - T4H;
				   T4X = T4H + T4W;
				   T9h = FNMS(KP707106781, T9a, T99);
				   T9b = FMA(KP707106781, T9a, T99);
				   T6Q = T5S - T5L;
				   T5T = T5L + T5S;
				   T64 = T62 + T63;
				   T6N = T63 - T62;
				   {
					E T6k, T5U, T6l, T65, T5k;
					T6J = T5c - T5j;
					T5k = T5c + T5j;
					T9c = T6d + T6e;
					T6f = T6d - T6e;
					T6k = FNMS(KP707106781, T5T, T5E);
					T5U = FMA(KP707106781, T5T, T5E);
					T6l = FNMS(KP707106781, T64, T61);
					T65 = FMA(KP707106781, T64, T61);
					T5l = FMA(KP707106781, T5k, T55);
					T6h = FNMS(KP707106781, T5k, T55);
					T6m = FMA(KP668178637, T6l, T6k);
					T6q = FNMS(KP668178637, T6k, T6l);
					T6a = FMA(KP198912367, T5U, T65);
					T66 = FNMS(KP198912367, T65, T5U);
					T5v = T5t + T5u;
					T6G = T5u - T5t;
				   }
			      }
			      {
				   E T68, T4Y, T9d, T9f, T6i, T5w;
				   T68 = FNMS(KP923879532, T4X, T4s);
				   T4Y = FMA(KP923879532, T4X, T4s);
				   T9d = FMA(KP923879532, T9c, T9b);
				   T9f = FNMS(KP923879532, T9c, T9b);
				   T6i = FNMS(KP707106781, T5v, T5s);
				   T5w = FMA(KP707106781, T5v, T5s);
				   {
					E T6o, T9l, T9m, T6r;
					{
					     E T6g, T6p, T69, T5x, T6n, T6j;
					     T6o = FMA(KP923879532, T6f, T6c);
					     T6g = FNMS(KP923879532, T6f, T6c);
					     T6j = FNMS(KP668178637, T6i, T6h);
					     T6p = FMA(KP668178637, T6h, T6i);
					     T69 = FNMS(KP198912367, T5l, T5w);
					     T5x = FMA(KP198912367, T5w, T5l);
					     T6n = T6j + T6m;
					     T9k = T6m - T6j;
					     T9j = FMA(KP923879532, T9i, T9h);
					     T9l = FNMS(KP923879532, T9i, T9h);
					     {
						  E T6b, T9e, T9g, T67;
						  T6b = T69 - T6a;
						  T9e = T69 + T6a;
						  T9g = T66 - T5x;
						  T67 = T5x + T66;
						  iio[-WS(ios, 29)] = FMA(KP831469612, T6n, T6g);
						  rio[WS(ios, 13)] = FNMS(KP831469612, T6n, T6g);
						  rio[WS(ios, 9)] = FMA(KP980785280, T6b, T68);
						  iio[-WS(ios, 25)] = FNMS(KP980785280, T6b, T68);
						  iio[-WS(ios, 1)] = FMA(KP980785280, T9e, T9d);
						  rio[WS(ios, 17)] = FMS(KP980785280, T9e, T9d);
						  iio[-WS(ios, 9)] = FMA(KP980785280, T9g, T9f);
						  rio[WS(ios, 25)] = FMS(KP980785280, T9g, T9f);
						  rio[WS(ios, 1)] = FMA(KP980785280, T67, T4Y);
						  iio[-WS(ios, 17)] = FNMS(KP980785280, T67, T4Y);
						  T9m = T6p + T6q;
						  T6r = T6p - T6q;
					     }
					}
					iio[-WS(ios, 13)] = FNMS(KP831469612, T9m, T9l);
					rio[WS(ios, 29)] = -(FMA(KP831469612, T9m, T9l));
					rio[WS(ios, 5)] = FMA(KP831469612, T6r, T6o);
					iio[-WS(ios, 21)] = FNMS(KP831469612, T6r, T6o);
				   }
			      }
			 }
			 {
			      E T6Y, T6w, T9w, T6D, T9v, T9p, T9q, T71, T6H, T74, T78, T7c, T6W, T6S;
			      {
				   E T6Z, T6z, T6C, T70;
				   T6Z = FNMS(KP414213562, T6x, T6y);
				   T6z = FMA(KP414213562, T6y, T6x);
				   iio[-WS(ios, 5)] = FMA(KP831469612, T9k, T9j);
				   rio[WS(ios, 21)] = FMS(KP831469612, T9k, T9j);
				   T6Y = FNMS(KP707106781, T6v, T6s);
				   T6w = FMA(KP707106781, T6v, T6s);
				   T6C = FNMS(KP414213562, T6B, T6A);
				   T70 = FMA(KP414213562, T6A, T6B);
				   T9w = T6z + T6C;
				   T6D = T6z - T6C;
				   T9v = FNMS(KP707106781, T9o, T9n);
				   T9p = FMA(KP707106781, T9o, T9n);
				   {
					E T77, T6O, T76, T6R;
					T9q = T70 - T6Z;
					T71 = T6Z + T70;
					T77 = FMA(KP707106781, T6N, T6M);
					T6O = FNMS(KP707106781, T6N, T6M);
					T76 = FMA(KP707106781, T6Q, T6P);
					T6R = FNMS(KP707106781, T6Q, T6P);
					T6H = FNMS(KP707106781, T6G, T6F);
					T74 = FMA(KP707106781, T6G, T6F);
					T78 = FMA(KP198912367, T77, T76);
					T7c = FNMS(KP198912367, T76, T77);
					T6W = FNMS(KP668178637, T6O, T6R);
					T6S = FMA(KP668178637, T6R, T6O);
				   }
			      }
			      {
				   E T6U, T6E, T9r, T9t, T73, T6K;
				   T6U = FNMS(KP923879532, T6D, T6w);
				   T6E = FMA(KP923879532, T6D, T6w);
				   T9r = FMA(KP923879532, T9q, T9p);
				   T9t = FNMS(KP923879532, T9q, T9p);
				   T73 = FMA(KP707106781, T6J, T6I);
				   T6K = FNMS(KP707106781, T6J, T6I);
				   {
					E T7a, T9x, T9y, T7d;
					{
					     E T72, T7b, T6V, T6L, T79, T75;
					     T7a = FMA(KP923879532, T71, T6Y);
					     T72 = FNMS(KP923879532, T71, T6Y);
					     T75 = FMA(KP198912367, T74, T73);
					     T7b = FNMS(KP198912367, T73, T74);
					     T6V = FNMS(KP668178637, T6H, T6K);
					     T6L = FMA(KP668178637, T6K, T6H);
					     T79 = T75 + T78;
					     T9A = T78 - T75;
					     T9z = FMA(KP923879532, T9w, T9v);
					     T9x = FNMS(KP923879532, T9w, T9v);
					     {
						  E T6X, T9s, T9u, T6T;
						  T6X = T6V + T6W;
						  T9s = T6V - T6W;
						  T9u = T6S - T6L;
						  T6T = T6L + T6S;
						  rio[WS(ios, 7)] = FMA(KP980785280, T79, T72);
						  iio[-WS(ios, 23)] = FNMS(KP980785280, T79, T72);
						  rio[WS(ios, 11)] = FMA(KP831469612, T6X, T6U);
						  iio[-WS(ios, 27)] = FNMS(KP831469612, T6X, T6U);
						  iio[-WS(ios, 3)] = FMA(KP831469612, T9s, T9r);
						  rio[WS(ios, 19)] = FMS(KP831469612, T9s, T9r);
						  iio[-WS(ios, 11)] = FMA(KP831469612, T9u, T9t);
						  rio[WS(ios, 27)] = FMS(KP831469612, T9u, T9t);
						  rio[WS(ios, 3)] = FMA(KP831469612, T6T, T6E);
						  iio[-WS(ios, 19)] = FNMS(KP831469612, T6T, T6E);
						  T9y = T7c - T7b;
						  T7d = T7b + T7c;
					     }
					}