t2_25.c

快速fft变换

						  T3t = T3h + T3s;
						  T4B = FNMS(KP250000000, T3s, T3h);
						  T3u = T7 * T1y;
					     }
					     T6y = FNMS(KP559016994, T4v, T4u);
					     T4w = FMA(KP559016994, T4v, T4u);
					     T6B = FNMS(KP559016994, T4C, T4B);
					     T4D = FMA(KP559016994, T4C, T4B);
					     T1z = FMA(Tb, T1y, T1x);
					     T3v = FNMS(Tb, T1w, T3u);
					}
					{
					     E T1Q, T3C, T1Y, T3E;
					     {
						  E T1L, T1P, T1T, T1X, T1M, T3B, T1U, T3D;
						  T1L = ri[WS(rs, 14)];
						  T1P = ii[WS(rs, 14)];
						  T1T = ri[WS(rs, 19)];
						  T1X = ii[WS(rs, 19)];
						  T1M = T1K * T1L;
						  T3B = T1K * T1P;
						  T1U = T1S * T1T;
						  T3D = T1S * T1X;
						  T1Q = FMA(T1O, T1P, T1M);
						  T3C = FNMS(T1O, T1L, T3B);
						  T1Y = FMA(T1W, T1X, T1U);
						  T3E = FNMS(T1W, T1T, T3D);
					     }
					     {
						  E T1A, T1C, T1B, T3w;
						  T1A = ri[WS(rs, 9)];
						  T1C = ii[WS(rs, 9)];
						  T4N = T1Y - T1Q;
						  T1Z = T1Q + T1Y;
						  T3F = T3C + T3E;
						  T4U = T3E - T3C;
						  T1B = T8 * T1A;
						  T3w = T8 * T1C;
						  T1E = ri[WS(rs, 24)];
						  T1G = ii[WS(rs, 24)];
						  T1D = FMA(Tc, T1C, T1B);
						  T3x = FNMS(Tc, T1A, T3w);
						  T1F = Tk * T1E;
						  T3y = Tk * T1G;
					     }
					}
					T1H = FMA(Tm, T1G, T1F);
					T3z = FNMS(Tm, T1E, T3y);
				   }
				   {
					E T2f, T2j, T2g, T3N;
					{
					     E T23, T25, T24, T4R, T4K, T4Q, T4J, T3J, T1I, T4M;
					     T23 = ri[WS(rs, 2)];
					     T25 = ii[WS(rs, 2)];
					     T1I = T1D + T1H;
					     T4M = T1H - T1D;
					     {
						  E T3A, T4T, T20, T3G;
						  T3A = T3x + T3z;
						  T4T = T3z - T3x;
						  T24 = T19 * T23;
						  T6v = FNMS(KP618033988, T4M, T4N);
						  T4O = FMA(KP618033988, T4N, T4M);
						  T20 = T1I + T1Z;
						  T4R = T1I - T1Z;
						  T6s = FNMS(KP618033988, T4T, T4U);
						  T4V = FMA(KP618033988, T4U, T4T);
						  T3G = T3A + T3F;
						  T4K = T3F - T3A;
						  T21 = T1z + T20;
						  T4Q = FNMS(KP250000000, T20, T1z);
						  T3H = T3v + T3G;
						  T4J = FNMS(KP250000000, T3G, T3v);
						  T3J = T19 * T25;
					     }
					     T6r = FNMS(KP559016994, T4R, T4Q);
					     T4S = FMA(KP559016994, T4R, T4Q);
					     T6u = FMA(KP559016994, T4K, T4J);
					     T4L = FNMS(KP559016994, T4K, T4J);
					     T26 = FMA(T1b, T25, T24);
					     T3K = FNMS(T1b, T23, T3J);
					}
					{
					     E T2r, T3R, T2z, T3T;
					     {
						  E T2n, T2q, T2u, T2y, T2o, T3Q, T2v, T3S;
						  T2n = ri[WS(rs, 12)];
						  T2q = ii[WS(rs, 12)];
						  T2u = ri[WS(rs, 17)];
						  T2y = ii[WS(rs, 17)];
						  T2o = T2m * T2n;
						  T3Q = T2m * T2q;
						  T2v = T2t * T2u;
						  T3S = T2t * T2y;
						  T2r = FMA(T2p, T2q, T2o);
						  T3R = FNMS(T2p, T2n, T3Q);
						  T2z = FMA(T2x, T2y, T2v);
						  T3T = FNMS(T2x, T2u, T3S);
					     }
					     {
						  E T28, T2b, T29, T3L;
						  T28 = ri[WS(rs, 7)];
						  T2b = ii[WS(rs, 7)];
						  T5a = T2z - T2r;
						  T2A = T2r + T2z;
						  T3U = T3R + T3T;
						  T53 = T3R - T3T;
						  T29 = T27 * T28;
						  T3L = T27 * T2b;
						  T2f = ri[WS(rs, 22)];
						  T2j = ii[WS(rs, 22)];
						  T2c = FMA(T2a, T2b, T29);
						  T3M = FNMS(T2a, T28, T3L);
						  T2g = T2e * T2f;
						  T3N = T2e * T2j;
					     }
					}
					T2k = FMA(T2i, T2j, T2g);
					T3O = FNMS(T2i, T2f, T3N);
				   }
			      }
			 }
		    }
		    {
			 E T7l, T5b, T6d, T54, T6g, T51, T6f, T7m, T6c, T58, T4e, T4c, T7A, T7y, T4d;
			 E T3f;
			 {
			      E T7w, T22, T7x, T3b, T3I, T3c, T3e, T3d;
			      T7l = T3t + T3H;
			      T3I = T3t - T3H;
			      {
				   E T2l, T59, T3P, T52;
				   T2l = T2c + T2k;
				   T59 = T2k - T2c;
				   T3P = T3M + T3O;
				   T52 = T3O - T3M;
				   T5b = FMA(KP618033988, T5a, T59);
				   T6d = FNMS(KP618033988, T59, T5a);
				   {
					E T50, T2B, T57, T3V;
					T50 = T2A - T2l;
					T2B = T2l + T2A;
					T54 = FNMS(KP618033988, T53, T52);
					T6g = FMA(KP618033988, T52, T53);
					T57 = T3U - T3P;
					T3V = T3P + T3U;
					{
					     E T4Z, T2C, T56, T3W, T4b;
					     T4Z = FNMS(KP250000000, T2B, T26);
					     T2C = T26 + T2B;
					     T56 = FNMS(KP250000000, T3V, T3K);
					     T3W = T3K + T3V;
					     T7w = T1v - T21;
					     T22 = T1v + T21;
					     T51 = FNMS(KP559016994, T50, T4Z);
					     T6f = FMA(KP559016994, T50, T4Z);
					     T4b = T3W - T4a;
					     T7m = T3W + T4a;
					     T6c = FMA(KP559016994, T57, T56);
					     T58 = FNMS(KP559016994, T57, T56);
					     T7x = T2C - T3a;
					     T3b = T2C + T3a;
					     T4e = FNMS(KP618033988, T3I, T4b);
					     T4c = FMA(KP618033988, T4b, T3I);
					}
				   }
			      }
			      T3c = T22 + T3b;
			      T3e = T22 - T3b;
			      ri[0] = TN + T3c;
			      T3d = FNMS(KP250000000, T3c, TN);
			      T7A = FNMS(KP618033988, T7w, T7x);
			      T7y = FMA(KP618033988, T7x, T7w);
			      T4d = FNMS(KP559016994, T3e, T3d);
			      T3f = FMA(KP559016994, T3e, T3d);
			 }
			 {
			      E T69, T85, T7Y, T68, T66, T84, T82, T7X, T67, T5Z;
			      {
				   E T4t, T5H, T5Q, T7T, T7H, T5P, T5M, T5L, T5A, T7O, T5D, T7P, T7K, T7M, T5u;
				   E T5w, T5K, T63, T61, T5U, T7D, T7z, T7v;
				   {
					E T7u, T7t, T4h, T7n;
					T69 = FNMS(KP559016994, T4g, T4f);
					T4h = FMA(KP559016994, T4g, T4f);
					T7u = T7l - T7m;
					T7n = T7l + T7m;
					ri[WS(rs, 5)] = FMA(KP951056516, T4c, T3f);
					ri[WS(rs, 20)] = FNMS(KP951056516, T4c, T3f);
					ri[WS(rs, 15)] = FMA(KP951056516, T4e, T4d);
					ri[WS(rs, 10)] = FNMS(KP951056516, T4e, T4d);
					ii[0] = T7n + T7s;
					T7t = FNMS(KP250000000, T7n, T7s);
					T4t = FMA(KP951056516, T4s, T4h);
					T5H = FNMS(KP951056516, T4s, T4h);
					T7D = FMA(KP559016994, T7C, T7B);
					T85 = FNMS(KP559016994, T7C, T7B);
					T7z = FNMS(KP559016994, T7u, T7t);
					T7v = FMA(KP559016994, T7u, T7t);
				   }
				   {
					E T5I, T5J, T5S, T4P, T5y, T4I, T5C, T5s, T4W, T5T, T55, T5c;
					{
					     E T4A, T4H, T5k, T5r;
					     T5Q = FNMS(KP951056516, T4z, T4w);
					     T4A = FMA(KP951056516, T4z, T4w);
					     T7T = FMA(KP951056516, T7G, T7D);
					     T7H = FNMS(KP951056516, T7G, T7D);
					     ii[WS(rs, 20)] = FMA(KP951056516, T7y, T7v);
					     ii[WS(rs, 5)] = FNMS(KP951056516, T7y, T7v);
					     ii[WS(rs, 15)] = FNMS(KP951056516, T7A, T7z);
					     ii[WS(rs, 10)] = FMA(KP951056516, T7A, T7z);
					     T4H = FMA(KP951056516, T4G, T4D);
					     T5P = FNMS(KP951056516, T4G, T4D);
					     T5I = FMA(KP951056516, T5j, T5g);
					     T5k = FNMS(KP951056516, T5j, T5g);
					     T5r = FNMS(KP951056516, T5q, T5n);
					     T5J = FMA(KP951056516, T5q, T5n);
					     T5S = FNMS(KP951056516, T4O, T4L);
					     T4P = FMA(KP951056516, T4O, T4L);
					     T5y = FNMS(KP256756360, T4A, T4H);
					     T4I = FMA(KP256756360, T4H, T4A);
					     T5C = FNMS(KP939062505, T5k, T5r);
					     T5s = FMA(KP939062505, T5r, T5k);
					     T4W = FNMS(KP951056516, T4V, T4S);
					     T5T = FMA(KP951056516, T4V, T4S);
					     T5M = FMA(KP951056516, T54, T51);
					     T55 = FNMS(KP951056516, T54, T51);
					     T5c = FMA(KP951056516, T5b, T58);
					     T5L = FNMS(KP951056516, T5b, T58);
					}
					{
					     E T4Y, T5t, T5z, T4X;
					     T5z = FNMS(KP634619297, T4P, T4W);
					     T4X = FMA(KP634619297, T4W, T4P);
					     {
						  E T5B, T5d, T7I, T7J;
						  T5B = FNMS(KP549754652, T55, T5c);
						  T5d = FMA(KP549754652, T5c, T55);
						  T7I = FNMS(KP871714437, T5z, T5y);
						  T5A = FMA(KP871714437, T5z, T5y);
						  T4Y = FMA(KP871714437, T4X, T4I);
						  T7O = FNMS(KP871714437, T4X, T4I);
						  T7J = FMA(KP831864738, T5C, T5B);
						  T5D = FNMS(KP831864738, T5C, T5B);
						  T5t = FMA(KP831864738, T5s, T5d);
						  T7P = FNMS(KP831864738, T5s, T5d);
						  T7K = FMA(KP904730450, T7J, T7I);
						  T7M = FNMS(KP904730450, T7J, T7I);
					     }
					     T5u = FMA(KP904730450, T5t, T4Y);
					     T5w = FNMS(KP904730450, T5t, T4Y);
					}
					T5K = FNMS(KP126329378, T5J, T5I);
					T63 = FMA(KP126329378, T5I, T5J);
					T61 = FNMS(KP827271945, T5S, T5T);
					T5U = FMA(KP827271945, T5T, T5S);
				   }
				   {
					E T65, T81, T62, T80, T7W, T5W, T5Y;
					{
					     E T5O, T5V, T64, T5N;
					     ri[WS(rs, 1)] = FMA(KP968583161, T5u, T4t);
					     T64 = FMA(KP470564281, T5L, T5M);
					     T5N = FNMS(KP470564281, T5M, T5L);
					     {
						  E T60, T5R, T7U, T7V;
						  T60 = FNMS(KP634619297, T5P, T5Q);
						  T5R = FMA(KP634619297, T5Q, T5P);
						  T7U = FMA(KP912018591, T64, T63);
						  T65 = FNMS(KP912018591, T64, T63);
						  T5O = FNMS(KP912018591, T5N, T5K);
						  T81 = FMA(KP912018591, T5N, T5K);
						  T7V = FNMS(KP912575812, T61, T60);
						  T62 = FMA(KP912575812, T61, T60);
						  T5V = FNMS(KP912575812, T5U, T5R);
						  T80 = FMA(KP912575812, T5U, T5R);
						  T7W = FMA(KP851038619, T7V, T7U);
						  T7Y = FNMS(KP851038619, T7V, T7U);
						  ii[WS(rs, 1)] = FMA(KP968583161, T7K, T7H);
					     }
					     T5W = FNMS(KP851038619, T5V, T5O);
					     T5Y = FMA(KP851038619, T5V, T5O);
					}
					{
					     E T5G, T5E, T7S, T7Q, T7L, T5F, T5x, T5v, T5X, T7R, T7N;
					     T5G = FNMS(KP683113946, T5A, T5D);
					     T5E = FMA(KP559154169, T5D, T5A);
					     ii[WS(rs, 4)] = FNMS(KP992114701, T7W, T7T);
					     ri[WS(rs, 4)] = FNMS(KP992114701, T5W, T5H);
					     T5v = FNMS(KP242145790, T5u, T4t);
					     T7S = FNMS(KP683113946, T7O, T7P);
					     T7Q = FMA(KP559154169, T7P, T7O);
					     T7L = FNMS(KP242145790, T7K, T7H);
					     T5F = FNMS(KP541454447, T5w, T5v);
					     T5x = FMA(KP541454447, T5w, T5v);
					     T68 = FMA(KP525970792, T62, T65);
					     T66 = FNMS(KP726211448, T65, T62);
					     ri[WS(rs, 11)] = FNMS(KP833417178, T5G, T5F);
					     ri[WS(rs, 16)] = FMA(KP833417178, T5G, T5F);
					     ri[WS(rs, 21)] = FNMS(KP921177326, T5E, T5x);
					     ri[WS(rs, 6)] = FMA(KP921177326, T5E, T5x);
					     T7R = FNMS(KP541454447, T7M, T7L);
					     T7N = FMA(KP541454447, T7M, T7L);
					     T5X = FMA(KP248028675, T5W, T5H);
					     ii[WS(rs, 11)] = FMA(KP833417178, T7S, T7R);
					     ii[WS(rs, 16)] = FNMS(KP833417178, T7S, T7R);
					     ii[WS(rs, 21)] = FMA(KP921177326, T7Q, T7N);
					     ii[WS(rs, 6)] = FNMS(KP921177326, T7Q, T7N);
					     T84 = FNMS(KP525970792, T80, T81);
					     T82 = FMA(KP726211448, T81, T80);
					     T7X = FMA(KP248028675, T7W, T7T);
					     T67 = FNMS(KP554608978, T5Y, T5X);
					     T5Z = FMA(KP554608978, T5Y, T5X);
					}
				   }
			      }
			      {
				   E T6b, T6T, T8j, T87, T72, T71, T6P, T8r, T6M, T8q, T7f, T6W, T8m, T8o, T6I;
				   E T6G, T7d, T76, T7g, T6Z, T83, T7Z;
				   ri[WS(rs, 14)] = FNMS(KP943557151, T68, T67);
				   ri[WS(rs, 19)] = FMA(KP943557151, T68, T67);
				   ri[WS(rs, 24)] = FMA(KP803003575, T66, T5Z);
				   ri[WS(rs, 9)] = FNMS(KP803003575, T66, T5Z);
				   T83 = FNMS(KP554608978, T7Y, T7X);
				   T7Z = FMA(KP554608978, T7Y, T7X);
				   T6b = FMA(KP951056516, T6a, T69);
				   T6T = FNMS(KP951056516, T6a, T69);
				   ii[WS(rs, 14)] = FMA(KP943557151, T84, T83);
				   ii[WS(rs, 19)] = FNMS(KP943557151, T84, T83);
				   ii[WS(rs, 24)] = FMA(KP803003575, T82, T7Z);
				   ii[WS(rs, 9)] = FNMS(KP803003575, T82, T7Z);
				   {
					E T6X, T6Y, T74, T6N, T6i, T75, T6U, T6V, T6t, T6L, T6E, T6O, T6p, T6w;
					{
					     E T6A, T6D, T6e, T6h, T6l, T6o;
					     T6X = FNMS(KP951056516, T6d, T6c);
					     T6e = FMA(KP951056516, T6d, T6c);
					     T6h = FMA(KP951056516, T6g, T6f);
					     T6Y = FNMS(KP951056516, T6g, T6f);
					     T74 = FMA(KP951056516, T6z, T6y);
					     T6A = FNMS(KP951056516, T6z, T6y);
					     T8j = FNMS(KP951056516, T86, T85);
					     T87 = FMA(KP951056516, T86, T85);
					     T6N = FNMS(KP062914667, T6e, T6h);
					     T6i = FMA(KP062914667, T6h, T6e);
					     T6D = FMA(KP951056516, T6C, T6B);
					     T75 = FNMS(KP951056516, T6C, T6B);
					     T6U = FMA(KP951056516, T6k, T6j);
					     T6l = FNMS(KP951056516, T6k, T6j);
					     T6o = FNMS(KP951056516, T6n, T6m);
					     T6V = FMA(KP951056516, T6n, T6m);
					     T72 = FMA(KP951056516, T6s, T6r);
					     T6t = FNMS(KP951056516, T6s, T6r);
					     T6L = FNMS(KP939062505, T6A, T6D);
					     T6E = FMA(KP939062505, T6D, T6A);
					     T6O = FMA(KP827271945, T6l, T6o);
					     T6p = FNMS(KP827271945, T6o, T6l);
					     T6w = FMA(KP951056516, T6v, T6u);
					     T71 = FNMS(KP951056516, T6v, T6u);
					}
					{
					     E T8k, T6q, T6K, T6x, T8l, T6F;
					     T8k = FMA(KP772036680, T6O, T6N);
					     T6P = FNMS(KP772036680, T6O, T6N);
					     T6q = FMA(KP772036680, T6p, T6i);
					     T8r = FNMS(KP772036680, T6p, T6i);
					     T6K = FMA(KP126329378, T6t, T6w);
					     T6x = FNMS(KP126329378, T6w, T6t);
					     T8l = FNMS(KP734762448, T6L, T6K);
					     T6M = FMA(KP734762448, T6L, T6K);
					     T6F = FNMS(KP734762448, T6E, T6x);
					     T8q = FMA(KP734762448, T6E, T6x);
					     T7f = FNMS(KP062914667, T6U, T6V);
					     T6W = FMA(KP062914667, T6V, T6U);
					     T8m = FMA(KP994076283, T8l, T8k);
					     T8o = FNMS(KP994076283, T8l, T8k);
					     T6I = FMA(KP994076283, T6F, T6q);
					     T6G = FNMS(KP994076283, T6F, T6q);
					}
					T7d = FNMS(KP549754652, T74, T75);
					T76 = FMA(KP549754652, T75, T74);
					T7g = FNMS(KP634619297, T6X, T6Y);
					T6Z = FMA(KP634619297, T6Y, T6X);
				   }
				   {
					E T88, T7h, T70, T8f, T7c, T73;
					ri[WS(rs, 3)] = FMA(KP998026728, T6G, T6b);
					T88 = FMA(KP845997307, T7g, T7f);
					T7h = FNMS(KP845997307, T7g, T7f);
					T70 = FMA(KP845997307, T6Z, T6W);
					T8f = FNMS(KP845997307, T6Z, T6W);
					T7c = FMA(KP470564281, T71, T72);
					T73 = FNMS(KP470564281, T72, T71);
					ii[WS(rs, 3)] = FNMS(KP998026728, T8m, T8j);
					{
					     E T7e, T8e, T8a, T78, T7a, T8u, T8s, T8t, T8p, T79;
					     {
						  E T6S, T6Q, T6H, T89, T77, T6J, T6R, T8n;
						  T6S = FMA(KP614372930, T6M, T6P);
						  T6Q = FNMS(KP621716863, T6P, T6M);
						  T89 = FNMS(KP968479752, T7d, T7c);
						  T7e = FMA(KP968479752, T7d, T7c);