hf_25.c

快速fft变换

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