t1_25.c

快速fft变换

			      }
			      {
				   E T1I, T1L, T1O, T3h, T1J, T1R, T1N, T1K, T1Q, T3j, T1P;
				   {
					E T1C, T1F, T1B, T1E;
					{
					     E T4g, T4n, T4f, T4m, T4p, T36;
					     T1C = ri[WS(rs, 2)];
					     T4p = T35 - T33;
					     T36 = T33 + T35;
					     {
						  E T4i, T1k, T3c, T1y;
						  T4i = T1j - T1d;
						  T1k = T1d + T1j;
						  T5Y = FNMS(KP618033988, T4p, T4q);
						  T4r = FMA(KP618033988, T4q, T4p);
						  T3c = T36 + T3b;
						  T4g = T3b - T36;
						  T61 = FNMS(KP618033988, T4i, T4j);
						  T4k = FMA(KP618033988, T4j, T4i);
						  T1y = T1k + T1x;
						  T4n = T1k - T1x;
						  T3d = T31 + T3c;
						  T4f = FNMS(KP250000000, T3c, T31);
						  T1z = T17 + T1y;
						  T4m = FNMS(KP250000000, T1y, T17);
						  T1F = ii[WS(rs, 2)];
					     }
					     T60 = FMA(KP559016994, T4g, T4f);
					     T4h = FNMS(KP559016994, T4g, T4f);
					     T5X = FNMS(KP559016994, T4n, T4m);
					     T4o = FMA(KP559016994, T4n, T4m);
					     T1B = W[2];
					}
					T1E = W[3];
					{
					     E T1V, T1Y, T21, T3m, T1W, T24, T20, T1X, T23, T3f, T1D, T1U, T1H;
					     T1V = ri[WS(rs, 12)];
					     T1Y = ii[WS(rs, 12)];
					     T3f = T1B * T1F;
					     T1D = T1B * T1C;
					     T1U = W[22];
					     T21 = ri[WS(rs, 17)];
					     T3g = FNMS(T1E, T1C, T3f);
					     T1G = FMA(T1E, T1F, T1D);
					     T3m = T1U * T1Y;
					     T1W = T1U * T1V;
					     T24 = ii[WS(rs, 17)];
					     T20 = W[32];
					     T1X = W[23];
					     T23 = W[33];
					     {
						  E T3n, T1Z, T3p, T25, T3o, T22;
						  T1I = ri[WS(rs, 7)];
						  T3o = T20 * T24;
						  T22 = T20 * T21;
						  T3n = FNMS(T1X, T1V, T3m);
						  T1Z = FMA(T1X, T1Y, T1W);
						  T3p = FNMS(T23, T21, T3o);
						  T25 = FMA(T23, T24, T22);
						  T1L = ii[WS(rs, 7)];
						  T1H = W[12];
						  T3q = T3n + T3p;
						  T4z = T3n - T3p;
						  T4G = T25 - T1Z;
						  T26 = T1Z + T25;
					     }
					     T1O = ri[WS(rs, 22)];
					     T3h = T1H * T1L;
					     T1J = T1H * T1I;
					     T1R = ii[WS(rs, 22)];
					     T1N = W[42];
					     T1K = W[13];
					     T1Q = W[43];
					}
				   }
				   T3j = T1N * T1R;
				   T1P = T1N * T1O;
				   T3i = FNMS(T1K, T1I, T3h);
				   T1M = FMA(T1K, T1L, T1J);
				   T3k = FNMS(T1Q, T1O, T3j);
				   T1S = FMA(T1Q, T1R, T1P);
			      }
			 }
		    }
	       }
	       {
		    E T6R, T5M, T4A, T5J, T4H, T6S, T5I, T4E, T5L, T4x, T3K, T3I, T2K, T74, T76;
		    E T2J;
		    {
			 E T1A, T72, T73, T2H, T28, T2I;
			 {
			      E T3e, T4D, T4w, T4C, T4v, T3H, T4y, T3l;
			      T6R = T2Z + T3d;
			      T3e = T2Z - T3d;
			      T4y = T3k - T3i;
			      T3l = T3i + T3k;
			      {
				   E T4F, T1T, T3r, T27, T3s;
				   T4F = T1S - T1M;
				   T1T = T1M + T1S;
				   T5M = FMA(KP618033988, T4y, T4z);
				   T4A = FNMS(KP618033988, T4z, T4y);
				   T3r = T3l + T3q;
				   T4D = T3q - T3l;
				   T5J = FNMS(KP618033988, T4F, T4G);
				   T4H = FMA(KP618033988, T4G, T4F);
				   T27 = T1T + T26;
				   T4w = T26 - T1T;
				   T3s = T3g + T3r;
				   T4C = FNMS(KP250000000, T3r, T3g);
				   T28 = T1G + T27;
				   T4v = FNMS(KP250000000, T27, T1G);
				   T3H = T3s - T3G;
				   T6S = T3s + T3G;
			      }
			      T5I = FMA(KP559016994, T4D, T4C);
			      T4E = FNMS(KP559016994, T4D, T4C);
			      T5L = FMA(KP559016994, T4w, T4v);
			      T4x = FNMS(KP559016994, T4w, T4v);
			      T3K = FNMS(KP618033988, T3e, T3H);
			      T3I = FMA(KP618033988, T3H, T3e);
			 }
			 T1A = T11 + T1z;
			 T72 = T11 - T1z;
			 T73 = T28 - T2G;
			 T2H = T28 + T2G;
			 T2I = T1A + T2H;
			 T2K = T1A - T2H;
			 T74 = FMA(KP618033988, T73, T72);
			 T76 = FNMS(KP618033988, T72, T73);
			 ri[0] = Tt + T2I;
			 T2J = FNMS(KP250000000, T2I, Tt);
		    }
		    {
			 E T5F, T7B, T7u, T5E, T5C, T7A, T7y, T7t, T5D, T5v;
			 {
			      E T3Z, T5d, T7p, T7d, T5m, T5l, T56, T7k, T59, T7l, T5z, T5g, T7g, T7i, T52;
			      E T50, T5x, T5q, T5A, T5j, T70, T6Z, T3N;
			      T5F = FNMS(KP559016994, T3M, T3L);
			      T3N = FMA(KP559016994, T3M, T3L);
			      {
				   E T79, T3J, T2L, T6T;
				   T79 = FMA(KP559016994, T78, T77);
				   T7B = FNMS(KP559016994, T78, T77);
				   T3J = FNMS(KP559016994, T2K, T2J);
				   T2L = FMA(KP559016994, T2K, T2J);
				   T6T = T6R + T6S;
				   T70 = T6R - T6S;
				   T3Z = FMA(KP951056516, T3Y, T3N);
				   T5d = FNMS(KP951056516, T3Y, T3N);
				   ri[WS(rs, 5)] = FMA(KP951056516, T3I, T2L);
				   ri[WS(rs, 20)] = FNMS(KP951056516, T3I, T2L);
				   ri[WS(rs, 15)] = FMA(KP951056516, T3K, T3J);
				   ri[WS(rs, 10)] = FNMS(KP951056516, T3K, T3J);
				   ii[0] = T6T + T6Y;
				   T6Z = FNMS(KP250000000, T6T, T6Y);
				   T7p = FMA(KP951056516, T7c, T79);
				   T7d = FNMS(KP951056516, T7c, T79);
			      }
			      {
				   E T5e, T54, T4e, T5f, T5o, T5p, T5i, T4B, T58, T4Y, T55, T4t, T4I, T5h;
				   {
					E T4Q, T4X, T4l, T4s;
					{
					     E T46, T71, T75, T4d;
					     T5m = FNMS(KP951056516, T45, T42);
					     T46 = FMA(KP951056516, T45, T42);
					     T71 = FMA(KP559016994, T70, T6Z);
					     T75 = FNMS(KP559016994, T70, T6Z);
					     T4d = FMA(KP951056516, T4c, T49);
					     T5l = FNMS(KP951056516, T4c, T49);
					     T5e = FMA(KP951056516, T4P, T4M);
					     T4Q = FNMS(KP951056516, T4P, T4M);
					     ii[WS(rs, 20)] = FMA(KP951056516, T74, T71);
					     ii[WS(rs, 5)] = FNMS(KP951056516, T74, T71);
					     ii[WS(rs, 15)] = FNMS(KP951056516, T76, T75);
					     ii[WS(rs, 10)] = FMA(KP951056516, T76, T75);
					     T54 = FNMS(KP256756360, T46, T4d);
					     T4e = FMA(KP256756360, T4d, T46);
					     T4X = FNMS(KP951056516, T4W, T4T);
					     T5f = FMA(KP951056516, T4W, T4T);
					}
					T5o = FNMS(KP951056516, T4k, T4h);
					T4l = FMA(KP951056516, T4k, T4h);
					T4s = FNMS(KP951056516, T4r, T4o);
					T5p = FMA(KP951056516, T4r, T4o);
					T5i = FMA(KP951056516, T4A, T4x);
					T4B = FNMS(KP951056516, T4A, T4x);
					T58 = FNMS(KP939062505, T4Q, T4X);
					T4Y = FMA(KP939062505, T4X, T4Q);
					T55 = FNMS(KP634619297, T4l, T4s);
					T4t = FMA(KP634619297, T4s, T4l);
					T4I = FMA(KP951056516, T4H, T4E);
					T5h = FNMS(KP951056516, T4H, T4E);
				   }
				   {
					E T7e, T4u, T57, T4J, T7f, T4Z;
					T7e = FNMS(KP871714437, T55, T54);
					T56 = FMA(KP871714437, T55, T54);
					T4u = FMA(KP871714437, T4t, T4e);
					T7k = FNMS(KP871714437, T4t, T4e);
					T57 = FNMS(KP549754652, T4B, T4I);
					T4J = FMA(KP549754652, T4I, T4B);
					T7f = FMA(KP831864738, T58, T57);
					T59 = FNMS(KP831864738, T58, T57);
					T4Z = FMA(KP831864738, T4Y, T4J);
					T7l = FNMS(KP831864738, T4Y, T4J);
					T5z = FMA(KP126329378, T5e, T5f);
					T5g = FNMS(KP126329378, T5f, T5e);
					T7g = FMA(KP904730450, T7f, T7e);
					T7i = FNMS(KP904730450, T7f, T7e);
					T52 = FNMS(KP904730450, T4Z, T4u);
					T50 = FMA(KP904730450, T4Z, T4u);
				   }
				   T5x = FNMS(KP827271945, T5o, T5p);
				   T5q = FMA(KP827271945, T5p, T5o);
				   T5A = FMA(KP470564281, T5h, T5i);
				   T5j = FNMS(KP470564281, T5i, T5h);
			      }
			      {
				   E T7q, T5B, T5k, T7x, T5w, T5n;
				   ri[WS(rs, 1)] = FMA(KP968583161, T50, T3Z);
				   T7q = FMA(KP912018591, T5A, T5z);
				   T5B = FNMS(KP912018591, T5A, T5z);
				   T5k = FNMS(KP912018591, T5j, T5g);
				   T7x = FMA(KP912018591, T5j, T5g);
				   T5w = FNMS(KP634619297, T5l, T5m);
				   T5n = FMA(KP634619297, T5m, T5l);
				   ii[WS(rs, 1)] = FMA(KP968583161, T7g, T7d);
				   {
					E T5y, T7w, T7s, T5s, T5u, T7o, T7m, T7n, T7j, T5t;
					{
					     E T5c, T5a, T51, T7r, T5r, T53, T5b, T7h;
					     T5c = FNMS(KP683113946, T56, T59);
					     T5a = FMA(KP559154169, T59, T56);
					     T7r = FNMS(KP912575812, T5x, T5w);
					     T5y = FMA(KP912575812, T5x, T5w);
					     T5r = FNMS(KP912575812, T5q, T5n);
					     T7w = FMA(KP912575812, T5q, T5n);
					     T7s = FMA(KP851038619, T7r, T7q);
					     T7u = FNMS(KP851038619, T7r, T7q);
					     T5s = FNMS(KP851038619, T5r, T5k);
					     T5u = FMA(KP851038619, T5r, T5k);
					     T51 = FNMS(KP242145790, T50, T3Z);
					     ii[WS(rs, 4)] = FNMS(KP992114701, T7s, T7p);
					     ri[WS(rs, 4)] = FNMS(KP992114701, T5s, T5d);
					     T7o = FNMS(KP683113946, T7k, T7l);
					     T7m = FMA(KP559154169, T7l, T7k);
					     T53 = FMA(KP541454447, T52, T51);
					     T5b = FNMS(KP541454447, T52, T51);
					     T7h = FNMS(KP242145790, T7g, T7d);
					     ri[WS(rs, 11)] = FNMS(KP833417178, T5c, T5b);
					     ri[WS(rs, 16)] = FMA(KP833417178, T5c, T5b);
					     ri[WS(rs, 21)] = FNMS(KP921177326, T5a, T53);
					     ri[WS(rs, 6)] = FMA(KP921177326, T5a, T53);
					     T7n = FNMS(KP541454447, T7i, T7h);
					     T7j = FMA(KP541454447, T7i, T7h);
					}
					T5E = FMA(KP525970792, T5y, T5B);
					T5C = FNMS(KP726211448, T5B, T5y);
					ii[WS(rs, 21)] = FMA(KP921177326, T7m, T7j);
					ii[WS(rs, 6)] = FNMS(KP921177326, T7m, T7j);
					ii[WS(rs, 11)] = FMA(KP833417178, T7o, T7n);
					ii[WS(rs, 16)] = FNMS(KP833417178, T7o, T7n);
					T5t = FMA(KP248028675, T5s, T5d);
					T7A = FNMS(KP525970792, T7w, T7x);
					T7y = FMA(KP726211448, T7x, T7w);
					T7t = FMA(KP248028675, T7s, T7p);
					T5D = FNMS(KP554608978, T5u, T5t);
					T5v = FMA(KP554608978, T5u, T5t);
				   }
			      }
			 }
			 {
			      E T5H, T6p, T7P, T7D, T6y, T6x, T6l, T7X, T6i, T7W, T6L, T6s, T7S, T7U, T6e;
			      E T6c, T6J, T6C, T6M, T6v, T7z, T7v;
			      ri[WS(rs, 14)] = FNMS(KP943557151, T5E, T5D);
			      ri[WS(rs, 19)] = FMA(KP943557151, T5E, T5D);
			      ri[WS(rs, 24)] = FMA(KP803003575, T5C, T5v);
			      ri[WS(rs, 9)] = FNMS(KP803003575, T5C, T5v);
			      T7z = FNMS(KP554608978, T7u, T7t);
			      T7v = FMA(KP554608978, T7u, T7t);
			      T5H = FMA(KP951056516, T5G, T5F);
			      T6p = FNMS(KP951056516, T5G, T5F);
			      ii[WS(rs, 14)] = FMA(KP943557151, T7A, T7z);
			      ii[WS(rs, 19)] = FNMS(KP943557151, T7A, T7z);
			      ii[WS(rs, 24)] = FMA(KP803003575, T7y, T7v);
			      ii[WS(rs, 9)] = FNMS(KP803003575, T7y, T7v);
			      {
				   E T6t, T6u, T6A, T6j, T5O, T6B, T6q, T6r, T5Z, T6h, T6a, T6k, T5V, T62;
				   {
					E T66, T69, T5K, T5N, T5R, T5U;
					T6t = FNMS(KP951056516, T5J, T5I);
					T5K = FMA(KP951056516, T5J, T5I);
					T5N = FMA(KP951056516, T5M, T5L);
					T6u = FNMS(KP951056516, T5M, T5L);
					T6A = FMA(KP951056516, T65, T64);
					T66 = FNMS(KP951056516, T65, T64);
					T7P = FNMS(KP951056516, T7C, T7B);
					T7D = FMA(KP951056516, T7C, T7B);
					T6j = FNMS(KP062914667, T5K, T5N);
					T5O = FMA(KP062914667, T5N, T5K);
					T69 = FMA(KP951056516, T68, T67);
					T6B = FNMS(KP951056516, T68, T67);
					T6q = FMA(KP951056516, T5Q, T5P);
					T5R = FNMS(KP951056516, T5Q, T5P);
					T5U = FNMS(KP951056516, T5T, T5S);
					T6r = FMA(KP951056516, T5T, T5S);
					T6y = FMA(KP951056516, T5Y, T5X);
					T5Z = FNMS(KP951056516, T5Y, T5X);
					T6h = FNMS(KP939062505, T66, T69);
					T6a = FMA(KP939062505, T69, T66);
					T6k = FMA(KP827271945, T5R, T5U);
					T5V = FNMS(KP827271945, T5U, T5R);
					T62 = FMA(KP951056516, T61, T60);
					T6x = FNMS(KP951056516, T61, T60);
				   }
				   {
					E T7Q, T5W, T6g, T63, T7R, T6b;
					T7Q = FMA(KP772036680, T6k, T6j);
					T6l = FNMS(KP772036680, T6k, T6j);
					T5W = FMA(KP772036680, T5V, T5O);
					T7X = FNMS(KP772036680, T5V, T5O);
					T6g = FMA(KP126329378, T5Z, T62);
					T63 = FNMS(KP126329378, T62, T5Z);
					T7R = FNMS(KP734762448, T6h, T6g);
					T6i = FMA(KP734762448, T6h, T6g);
					T6b = FNMS(KP734762448, T6a, T63);
					T7W = FMA(KP734762448, T6a, T63);
					T6L = FNMS(KP062914667, T6q, T6r);
					T6s = FMA(KP062914667, T6r, T6q);
					T7S = FMA(KP994076283, T7R, T7Q);
					T7U = FNMS(KP994076283, T7R, T7Q);
					T6e = FMA(KP994076283, T6b, T5W);
					T6c = FNMS(KP994076283, T6b, T5W);
				   }
				   T6J = FNMS(KP549754652, T6A, T6B);
				   T6C = FMA(KP549754652, T6B, T6A);
				   T6M = FNMS(KP634619297, T6t, T6u);
				   T6v = FMA(KP634619297, T6u, T6t);
			      }
			      {
				   E T7E, T6N, T6w, T7L, T6I, T6z;
				   ri[WS(rs, 3)] = FMA(KP998026728, T6c, T5H);
				   T7E = FMA(KP845997307, T6M, T6L);
				   T6N = FNMS(KP845997307, T6M, T6L);
				   T6w = FMA(KP845997307, T6v, T6s);
				   T7L = FNMS(KP845997307, T6v, T6s);
				   T6I = FMA(KP470564281, T6x, T6y);
				   T6z = FNMS(KP470564281, T6y, T6x);
				   ii[WS(rs, 3)] = FNMS(KP998026728, T7S, T7P);
				   {
					E T6K, T7K, T7G, T6E, T6G, T80, T7Y, T7Z, T7V, T6F;
					{
					     E T6o, T6m, T6d, T7F, T6D, T6f, T6n, T7T;
					     T6o = FMA(KP614372930, T6i, T6l);
					     T6m = FNMS(KP621716863, T6l, T6i);
					     T7F = FNMS(KP968479752, T6J, T6I);
					     T6K = FMA(KP968479752, T6J, T6I);
					     T6D = FMA(KP968479752, T6C, T6z);
					     T7K = FNMS(KP968479752, T6C, T6z);
					     T7G = FMA(KP906616052, T7F, T7E);
					     T7I = FNMS(KP906616052, T7F, T7E);
					     T6E = FMA(KP906616052, T6D, T6w);
					     T6G = FNMS(KP906616052, T6D, T6w);
					     T6d = FNMS(KP249506682, T6c, T5H);
					     ii[WS(rs, 2)] = FNMS(KP998026728, T7G, T7D);
					     ri[WS(rs, 2)] = FMA(KP998026728, T6E, T6p);
					     T80 = FNMS(KP614372930, T7W, T7X);
					     T7Y = FMA(KP621716863, T7X, T7W);
					     T6f = FNMS(KP557913902, T6e, T6d);
					     T6n = FMA(KP557913902, T6e, T6d);
					     T7T = FMA(KP249506682, T7S, T7P);
					     ri[WS(rs, 18)] = FNMS(KP949179823, T6o, T6n);
					     ri[WS(rs, 13)] = FMA(KP949179823, T6o, T6n);
					     ri[WS(rs, 8)] = FMA(KP943557151, T6m, T6f);
					     ri[WS(rs, 23)] = FNMS(KP943557151, T6m, T6f);
					     T7Z = FNMS(KP557913902, T7U, T7T);
					     T7V = FMA(KP557913902, T7U, T7T);
					}
					T6Q = FNMS(KP560319534, T6K, T6N);
					T6O = FMA(KP681693190, T6N, T6K);
					ii[WS(rs, 23)] = FMA(KP943557151, T7Y, T7V);
					ii[WS(rs, 8)] = FNMS(KP943557151, T7Y, T7V);
					ii[WS(rs, 13)] = FMA(KP949179823, T80, T7Z);
					ii[WS(rs, 18)] = FNMS(KP949179823, T80, T7Z);
					T6F = FNMS(KP249506682, T6E, T6p);
					T7O = FNMS(KP560319534, T7K, T7L);
					T7M = FMA(KP681693190, T7L, T7K);