t2_32.c

最新的FFT程序

					     E T3K, T5G, T3S, T5I, T3V, T3W;
					     {
						  E T40, T3Y, T3N, T3R;
						  {
						       E T7v, T5m, T7w, T5r;
						       {
							    E T3H, T53, T2z, T3I, T3J, T2M, T54, T5o, T5F, T2D, T5n;
							    T3H = ri[WS(ios, 3)];
							    T7v = T50 + T52;
							    T53 = T50 - T52;
							    T2z = T2q + T2y;
							    T5m = T2q - T2y;
							    T2D = FMA(Tb, T2C, T2B);
							    T5n = T3 * T2C;
							    T3I = T4 * T3H;
							    T3J = ii[WS(ios, 3)];
							    T2M = T2D + T2L;
							    T54 = T2D - T2L;
							    T5o = FNMS(Tb, T2A, T5n);
							    T5F = T4 * T3J;
							    T55 = T53 - T54;
							    T6I = T53 + T54;
							    T2N = T2z + T2M;
							    T7A = T2z - T2M;
							    T7w = T5o + T5q;
							    T5r = T5o - T5q;
							    T3K = FMA(T7, T3J, T3I);
							    T5G = FNMS(T7, T3H, T5F);
						       }
						       T5s = T5m + T5r;
						       T6F = T5m - T5r;
						       T7x = T7v - T7w;
						       T8i = T7v + T7w;
						       T40 = ii[WS(ios, 11)];
						       T3Y = ri[WS(ios, 11)];
						  }
						  T3N = ri[WS(ios, 19)];
						  T3R = ii[WS(ios, 19)];
						  T3U = ri[WS(ios, 27)];
						  {
						       E T5P, T3Z, T3O, T5H;
						       T5P = Ti * T3Y;
						       T3Z = Te * T3Y;
						       T3O = T3M * T3N;
						       T5H = T3M * T3R;
						       T5Q = FMA(Te, T40, T5P);
						       T41 = FNMS(Ti, T40, T3Z);
						       T3S = FMA(T3Q, T3R, T3O);
						       T5I = FNMS(T3Q, T3N, T5H);
						       T3V = T2 * T3U;
						  }
						  T3W = ii[WS(ios, 27)];
					     }
					     T3T = T3K + T3S;
					     T5K = T3K - T3S;
					     T5J = T5G - T5I;
					     T7M = T5G + T5I;
					     T3X = FMA(Tg, T3W, T3V);
					     T5N = T2 * T3W;
					}
					{
					     E T37, T38, T39, T5g, T2Y, T2V;
					     {
						  E T2Q, T5M, T5R, T7N, T2O, T42, T5O, T56, T2P;
						  T2Q = ii[WS(ios, 5)];
						  T42 = T3X + T41;
						  T5M = T3X - T41;
						  T5O = FNMS(Tg, T3U, T5N);
						  T62 = T5K + T5J;
						  T5L = T5J - T5K;
						  T7J = T42 - T3T;
						  T43 = T3T + T42;
						  T5R = T5O - T5Q;
						  T7N = T5O + T5Q;
						  T2O = ri[WS(ios, 5)];
						  T37 = ri[WS(ios, 13)];
						  T63 = T5M - T5R;
						  T5S = T5M + T5R;
						  T8o = T7M + T7N;
						  T7O = T7M - T7N;
						  T56 = T2b * T2O;
						  T2P = T29 * T2O;
						  T38 = T1R * T37;
						  T39 = ii[WS(ios, 13)];
						  T57 = FMA(T29, T2Q, T56);
						  T2R = FNMS(T2b, T2Q, T2P);
					     }
					     T2U = ri[WS(ios, 21)];
					     T3a = FMA(T1U, T39, T38);
					     T5g = T1R * T39;
					     T2Y = ii[WS(ios, 21)];
					     T2V = T2T * T2U;
					     T35 = ii[WS(ios, 29)];
					     T5h = FNMS(T1U, T37, T5g);
					     T58 = T2T * T2Y;
					     T2Z = FMA(T2X, T2Y, T2V);
					     T32 = ri[WS(ios, 29)];
					}
				   }
			      }
			 }
		    }
		    {
			 E T5t, T5c, T5j, T5u, T88, T90, T8Z, T8b;
			 {
			      E T7e, T8T, T7D, T7y, T7h, T8U, T8S, T8R;
			      {
				   E T8c, T1i, T8A, T8z, T8O, T8J, T8N, T2l, T8L, T45, T8t, T8l, T8u, T8q, T3c;
				   E T8k, T8p, T8w, T2m;
				   {
					E T8x, T8y, T8j, T8C, T8I;
					{
					     E TI, T30, T7B, T5i, T7C, T5d, T3b, T1h, T59, T5b, T5e;
					     TI = Tq + TH;
					     T7e = Tq - TH;
					     T59 = FNMS(T2X, T2U, T58);
					     T5b = T2R - T2Z;
					     T30 = T2R + T2Z;
					     T5e = T34 * T32;
					     {
						  E T33, T5a, T5f, T36;
						  T33 = T31 * T32;
						  T5a = T57 - T59;
						  T7B = T57 + T59;
						  T5f = FMA(T31, T35, T5e);
						  T36 = FNMS(T34, T35, T33);
						  T5t = T5b + T5a;
						  T5c = T5a - T5b;
						  T5i = T5f - T5h;
						  T7C = T5f + T5h;
						  T5d = T36 - T3a;
						  T3b = T36 + T3a;
						  T1h = TZ + T1g;
						  T8T = T1g - TZ;
					     }
					     T7D = T7B - T7C;
					     T8j = T7B + T7C;
					     T5j = T5d + T5i;
					     T5u = T5d - T5i;
					     T3c = T30 + T3b;
					     T7y = T3b - T30;
					     T8c = TI - T1h;
					     T1i = TI + T1h;
					}
					T8k = T8i - T8j;
					T8x = T8i + T8j;
					T8y = T8n + T8o;
					T8p = T8n - T8o;
					T7h = T7f - T7g;
					T8C = T7f + T7g;
					T8I = T8D + T8H;
					T8U = T8H - T8D;
					T8A = T8x + T8y;
					T8z = T8x - T8y;
					T8O = T8I - T8C;
					T8J = T8C + T8I;
				   }
				   {
					E T8h, T8m, T3d, T44;
					T8h = T2N - T3c;
					T3d = T2N + T3c;
					T44 = T3G + T43;
					T8m = T3G - T43;
					T8N = T2k - T1J;
					T2l = T1J + T2k;
					T8L = T44 - T3d;
					T45 = T3d + T44;
					T8t = T8k - T8h;
					T8l = T8h + T8k;
					T8u = T8m + T8p;
					T8q = T8m - T8p;
				   }
				   T8w = T1i - T2l;
				   T2m = T1i + T2l;
				   {
					E T8s, T8P, T8Q, T8v;
					{
					     E T8r, T8M, T8K, T8g, T8B, T8f;
					     T8S = T8q - T8l;
					     T8r = T8l + T8q;
					     T8B = T8d + T8e;
					     T8f = T8d - T8e;
					     ri[0] = T2m + T45;
					     ri[WS(ios, 16)] = T2m - T45;
					     ri[WS(ios, 8)] = T8w + T8z;
					     ri[WS(ios, 24)] = T8w - T8z;
					     T8M = T8J - T8B;
					     T8K = T8B + T8J;
					     T8g = T8c + T8f;
					     T8s = T8c - T8f;
					     T8R = T8O - T8N;
					     T8P = T8N + T8O;
					     ii[WS(ios, 24)] = T8M - T8L;
					     ii[WS(ios, 8)] = T8L + T8M;
					     ii[WS(ios, 16)] = T8K - T8A;
					     ii[0] = T8A + T8K;
					     ri[WS(ios, 4)] = FMA(KP707106781, T8r, T8g);
					     ri[WS(ios, 20)] = FNMS(KP707106781, T8r, T8g);
					     T8Q = T8t + T8u;
					     T8v = T8t - T8u;
					}
					ii[WS(ios, 20)] = FNMS(KP707106781, T8Q, T8P);
					ii[WS(ios, 4)] = FMA(KP707106781, T8Q, T8P);
					ri[WS(ios, 12)] = FMA(KP707106781, T8v, T8s);
					ri[WS(ios, 28)] = FNMS(KP707106781, T8v, T8s);
				   }
			      }
			      {
				   E T7P, T7W, T7i, T7K, T8a, T86, T91, T8V, T8W, T7t, T7T, T7F, T92, T7Z, T89;
				   E T83;
				   {
					E T7X, T7n, T7s, T7Y, T84, T85;
					T7P = T7L - T7O;
					T84 = T7L + T7O;
					ii[WS(ios, 28)] = FNMS(KP707106781, T8S, T8R);
					ii[WS(ios, 12)] = FMA(KP707106781, T8S, T8R);
					T7W = T7e + T7h;
					T7i = T7e - T7h;
					T85 = T7I + T7J;
					T7K = T7I - T7J;
					T7X = T7m + T7l;
					T7n = T7l - T7m;
					T8a = FMA(KP414213562, T84, T85);
					T86 = FNMS(KP414213562, T85, T84);
					T91 = T8U - T8T;
					T8V = T8T + T8U;
					T7s = T7o + T7r;
					T7Y = T7o - T7r;
					{
					     E T82, T81, T7z, T7E;
					     T82 = T7x + T7y;
					     T7z = T7x - T7y;
					     T7E = T7A - T7D;
					     T81 = T7A + T7D;
					     T8W = T7n + T7s;
					     T7t = T7n - T7s;
					     T7T = FNMS(KP414213562, T7z, T7E);
					     T7F = FMA(KP414213562, T7E, T7z);
					     T92 = T7Y - T7X;
					     T7Z = T7X + T7Y;
					     T89 = FNMS(KP414213562, T81, T82);
					     T83 = FMA(KP414213562, T82, T81);
					}
				   }
				   {
					E T7S, T7u, T93, T95, T7U, T7Q;
					T7S = FNMS(KP707106781, T7t, T7i);
					T7u = FMA(KP707106781, T7t, T7i);
					T93 = FMA(KP707106781, T92, T91);
					T95 = FNMS(KP707106781, T92, T91);
					T7U = FMA(KP414213562, T7K, T7P);
					T7Q = FNMS(KP414213562, T7P, T7K);
					{
					     E T80, T87, T8X, T8Y;
					     T88 = FNMS(KP707106781, T7Z, T7W);
					     T80 = FMA(KP707106781, T7Z, T7W);
					     {
						  E T7V, T94, T96, T7R;
						  T7V = T7T + T7U;
						  T94 = T7U - T7T;
						  T96 = T7F + T7Q;
						  T7R = T7F - T7Q;
						  ri[WS(ios, 30)] = FMA(KP923879532, T7V, T7S);
						  ri[WS(ios, 14)] = FNMS(KP923879532, T7V, T7S);
						  ii[WS(ios, 22)] = FNMS(KP923879532, T94, T93);
						  ii[WS(ios, 6)] = FMA(KP923879532, T94, T93);
						  ii[WS(ios, 30)] = FMA(KP923879532, T96, T95);
						  ii[WS(ios, 14)] = FNMS(KP923879532, T96, T95);
						  ri[WS(ios, 6)] = FMA(KP923879532, T7R, T7u);
						  ri[WS(ios, 22)] = FNMS(KP923879532, T7R, T7u);
						  T87 = T83 + T86;
						  T90 = T86 - T83;
					     }
					     T8Z = FNMS(KP707106781, T8W, T8V);
					     T8X = FMA(KP707106781, T8W, T8V);
					     T8Y = T89 + T8a;
					     T8b = T89 - T8a;
					     ri[WS(ios, 2)] = FMA(KP923879532, T87, T80);
					     ri[WS(ios, 18)] = FNMS(KP923879532, T87, T80);
					     ii[WS(ios, 18)] = FNMS(KP923879532, T8Y, T8X);
					     ii[WS(ios, 2)] = FMA(KP923879532, T8Y, T8X);
					}
				   }
			      }
			 }
			 {
			      E T6s, T9o, T9n, T6v, T6N, T6Q, T6G, T6J, T9g, T9f;
			      {
				   E T6c, T4s, T9c, T4X, T9h, T9b, T9i, T6f, T5U, T6l, T64, T5k, T5v;
				   {
					E T6d, T6e, T99, T9a, T5T;
					{
					     E T4c, T4r, T4H, T4W;
					     T6s = T46 - T4b;
					     T4c = T46 + T4b;
					     ri[WS(ios, 10)] = FMA(KP923879532, T8b, T88);
					     ri[WS(ios, 26)] = FNMS(KP923879532, T8b, T88);
					     ii[WS(ios, 26)] = FNMS(KP923879532, T90, T8Z);
					     ii[WS(ios, 10)] = FMA(KP923879532, T90, T8Z);
					     T4r = T4j + T4q;
					     T9o = T4q - T4j;
					     T6d = FMA(KP414213562, T4z, T4G);
					     T4H = FNMS(KP414213562, T4G, T4z);
					     T4W = FMA(KP414213562, T4V, T4O);
					     T6e = FNMS(KP414213562, T4O, T4V);
					     T9n = T98 + T97;
					     T99 = T97 - T98;
					     T6c = FMA(KP707106781, T4r, T4c);
					     T4s = FNMS(KP707106781, T4r, T4c);
					     T9c = T4H + T4W;
					     T4X = T4H - T4W;
					     T9a = T6t + T6u;
					     T6v = T6t - T6u;
					}
					T6N = T5S - T5L;
					T5T = T5L + T5S;
					T9h = FNMS(KP707106781, T9a, T99);
					T9b = FMA(KP707106781, T9a, T99);
					T9i = T6e - T6d;
					T6f = T6d + T6e;
					T5U = FNMS(KP707106781, T5T, T5E);
					T6l = FMA(KP707106781, T5T, T5E);
					T64 = T62 + T63;
					T6Q = T62 - T63;
					T6G = T5j - T5c;
					T5k = T5c + T5j;
					T5v = T5t + T5u;
					T6J = T5t - T5u;
				   }
				   {
					E T6m, T6q, T6j, T6p, T9l, T9m;
					{
					     E T68, T4Y, T6a, T66, T69, T5x, T9j, T6k, T65, T9k, T6b, T67;
					     T68 = FNMS(KP923879532, T4X, T4s);
					     T4Y = FMA(KP923879532, T4X, T4s);
					     T6k = FMA(KP707106781, T64, T61);
					     T65 = FNMS(KP707106781, T64, T61);
					     {
						  E T6i, T5l, T6h, T5w;
						  T6i = FMA(KP707106781, T5k, T55);
						  T5l = FNMS(KP707106781, T5k, T55);
						  T6h = FMA(KP707106781, T5v, T5s);
						  T5w = FNMS(KP707106781, T5v, T5s);
						  T6m = FNMS(KP198912367, T6l, T6k);
						  T6q = FMA(KP198912367, T6k, T6l);
						  T6a = FMA(KP668178637, T5U, T65);
						  T66 = FNMS(KP668178637, T65, T5U);
						  T6j = FMA(KP198912367, T6i, T6h);
						  T6p = FNMS(KP198912367, T6h, T6i);
						  T69 = FNMS(KP668178637, T5l, T5w);
						  T5x = FMA(KP668178637, T5w, T5l);
					     }
					     T9j = FMA(KP923879532, T9i, T9h);
					     T9l = FNMS(KP923879532, T9i, T9h);
					     T9k = T6a - T69;
					     T6b = T69 + T6a;
					     T9m = T5x + T66;
					     T67 = T5x - T66;
					     ii[WS(ios, 21)] = FNMS(KP831469612, T9k, T9j);
					     ii[WS(ios, 5)] = FMA(KP831469612, T9k, T9j);
					     ri[WS(ios, 5)] = FMA(KP831469612, T67, T4Y);
					     ri[WS(ios, 21)] = FNMS(KP831469612, T67, T4Y);
					     ri[WS(ios, 29)] = FMA(KP831469612, T6b, T68);
					     ri[WS(ios, 13)] = FNMS(KP831469612, T6b, T68);
					}
					{
					     E T6o, T9d, T9e, T6r, T6g, T6n;
					     T6o = FNMS(KP923879532, T6f, T6c);
					     T6g = FMA(KP923879532, T6f, T6c);
					     T6n = T6j + T6m;
					     T9g = T6m - T6j;
					     T9f = FNMS(KP923879532, T9c, T9b);
					     T9d = FMA(KP923879532, T9c, T9b);
					     ii[WS(ios, 29)] = FMA(KP831469612, T9m, T9l);
					     ii[WS(ios, 13)] = FNMS(KP831469612, T9m, T9l);
					     ri[WS(ios, 1)] = FMA(KP980785280, T6n, T6g);
					     ri[WS(ios, 17)] = FNMS(KP980785280, T6n, T6g);
					     T9e = T6p + T6q;
					     T6r = T6p - T6q;
					     ii[WS(ios, 17)] = FNMS(KP980785280, T9e, T9d);
					     ii[WS(ios, 1)] = FMA(KP980785280, T9e, T9d);
					     ri[WS(ios, 9)] = FMA(KP980785280, T6r, T6o);
					     ri[WS(ios, 25)] = FNMS(KP980785280, 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);
					ii[WS(ios, 25)] = FNMS(KP980785280, T9g, T9f);
					ii[WS(ios, 9)] = FMA(KP980785280, T9g, T9f);
					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 = FNMS(KP198912367, T77, T76);
					     T7c = FMA(KP198912367, T76, T77);
					     T6W = FMA(KP668178637, T6O, T6R);
					     T6S = FNMS(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 = T75 + T78;
						  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;
						       ri[WS(ios, 7)] = FMA(KP980785280, T79, T72);
						       ri[WS(ios, 23)] = FNMS(KP980785280, T79, T72);
						       ri[WS(ios, 11)] = FMA(KP831469612, T6X, T6U);
						       ri[WS(ios, 27)] = FNMS(KP831469612, T6X, T6U);
						       ii[WS(ios, 19)] = FNMS(KP831469612, T9s, T9r);
						       ii[WS(ios, 3)] = FMA(KP831469612, T9s, T9r);
						       ii[WS(ios, 27)] = FNMS(KP831469612, T9u, T9t);
						       ii[WS(ios, 11)] = FMA(KP831469612, T9u, T9t);
						       ri[WS(ios, 3)] = FMA(KP831469612, T6T, T6E);
						       ri[WS(ios, 19)] = FNMS(KP831469612, T6T, T6E);
						       T9y = T7c - T7b;
						       T7d = T7b + T7c;
						  }
					     }