hb2_25.c

快速fft变换

					     T2I = FMA(KP559016994, T2H, T2G);
					     T6j = FNMS(KP559016994, T2H, T2G);
					     T6m = FNMS(KP559016994, T2O, T2N);
					     T2P = FMA(KP559016994, T2O, T2N);
					     TL = TB + TK;
					     T1Z = TB - TK;
					     {
						  E T6l, T7y, T6o, T7z;
						  T6l = FMA(KP951056516, T6k, T6j);
						  T7y = FNMS(KP951056516, T6k, T6j);
						  T6o = FMA(KP951056516, T6n, T6m);
						  T7z = FNMS(KP951056516, T6n, T6m);
						  T10 = Ts - TL;
						  TM = Ts + TL;
						  T2M = FNMS(KP951056516, T2L, T2I);
						  T4W = FMA(KP951056516, T2L, T2I);
						  T7K = FMA(KP939062505, T7y, T7z);
						  T7A = FNMS(KP939062505, T7z, T7y);
						  T6p = FNMS(KP549754652, T6o, T6l);
						  T6w = FMA(KP549754652, T6l, T6o);
					     }
					}
				   }
				   {
					E T34, T31, T4Y, T60, T3s, T3z, T5X;
					{
					     E T2X, T6c, T6f, T4V, T2T;
					     T2X = FNMS(KP559016994, T2W, T2V);
					     T6c = FMA(KP559016994, T2W, T2V);
					     T6f = FMA(KP559016994, T33, T32);
					     T34 = FNMS(KP559016994, T33, T32);
					     T4V = FNMS(KP951056516, T2S, T2P);
					     T2T = FMA(KP951056516, T2S, T2P);
					     {
						  E T7w, T6e, T7v, T6h;
						  T7w = FMA(KP951056516, T6d, T6c);
						  T6e = FNMS(KP951056516, T6d, T6c);
						  T7v = FMA(KP951056516, T6g, T6f);
						  T6h = FNMS(KP951056516, T6g, T6f);
						  T4X = FNMS(KP634619297, T4W, T4V);
						  T56 = FMA(KP634619297, T4V, T4W);
						  T3K = FMA(KP256756360, T2M, T2T);
						  T2U = FNMS(KP256756360, T2T, T2M);
						  T7x = FMA(KP126329378, T7w, T7v);
						  T7J = FNMS(KP126329378, T7v, T7w);
						  T6v = FNMS(KP470564281, T6e, T6h);
						  T6i = FMA(KP470564281, T6h, T6e);
						  T31 = FMA(KP951056516, T30, T2X);
						  T4Y = FNMS(KP951056516, T30, T2X);
					     }
					     T60 = FMA(KP559016994, T3r, T3q);
					     T3s = FNMS(KP559016994, T3r, T3q);
					     T3z = FNMS(KP559016994, T3y, T3x);
					     T5X = FMA(KP559016994, T3y, T3x);
					}
					{
					     E T5Z, T7r, T4Z, T38;
					     T4Z = FNMS(KP951056516, T37, T34);
					     T38 = FMA(KP951056516, T37, T34);
					     {
						  E T4O, T3w, T4P, T3D;
						  T4O = FMA(KP951056516, T3v, T3s);
						  T3w = FNMS(KP951056516, T3v, T3s);
						  T4P = FMA(KP951056516, T3C, T3z);
						  T3D = FNMS(KP951056516, T3C, T3z);
						  T50 = FNMS(KP827271945, T4Z, T4Y);
						  T57 = FMA(KP827271945, T4Y, T4Z);
						  T3L = FMA(KP634619297, T31, T38);
						  T39 = FNMS(KP634619297, T38, T31);
						  T4Q = FNMS(KP126329378, T4P, T4O);
						  T59 = FMA(KP126329378, T4O, T4P);
						  T3O = FNMS(KP939062505, T3w, T3D);
						  T3E = FMA(KP939062505, T3D, T3w);
						  T5Z = FMA(KP951056516, T5Y, T5X);
						  T7r = FNMS(KP951056516, T5Y, T5X);
					     }
					     {
						  E T3d, T3k, T64, T7s, T62;
						  T67 = FMA(KP559016994, T3c, T3b);
						  T3d = FNMS(KP559016994, T3c, T3b);
						  T3k = FNMS(KP559016994, T3j, T3i);
						  T64 = FMA(KP559016994, T3j, T3i);
						  T7s = FNMS(KP951056516, T61, T60);
						  T62 = FMA(KP951056516, T61, T60);
						  {
						       E T4S, T3h, T4R, T3o;
						       T4S = FMA(KP951056516, T3g, T3d);
						       T3h = FNMS(KP951056516, T3g, T3d);
						       T4R = FMA(KP951056516, T3n, T3k);
						       T3o = FNMS(KP951056516, T3n, T3k);
						       T7t = FNMS(KP827271945, T7s, T7r);
						       T7H = FMA(KP827271945, T7r, T7s);
						       T6y = FNMS(KP062914667, T5Z, T62);
						       T63 = FMA(KP062914667, T62, T5Z);
						       T4T = FNMS(KP470564281, T4S, T4R);
						       T5a = FMA(KP470564281, T4R, T4S);
						       T3N = FNMS(KP549754652, T3h, T3o);
						       T3p = FMA(KP549754652, T3o, T3h);
						       T66 = FNMS(KP951056516, T65, T64);
						       T7o = FMA(KP951056516, T65, T64);
						  }
					     }
					}
				   }
			      }
			 }
		    }
		    {
			 E T7q, T7G, T6J, T6I, T6q, T6b, T6B, T73, T6Q, T78, T6z, T6a;
			 cr[0] = T9 + TM;
			 {
			      E T1U, T2l, T1X, T2g, T1E, TZ, T2m, T20, T2v, T2n;
			      {
				   E T1W, T7p, T69, T1V;
				   T1W = T1S - T1T;
				   T1U = T1S + T1T;
				   T7p = FNMS(KP951056516, T68, T67);
				   T69 = FMA(KP951056516, T68, T67);
				   T1V = FNMS(KP250000000, T1U, T1R);
				   T7q = FMA(KP062914667, T7p, T7o);
				   T7G = FNMS(KP062914667, T7o, T7p);
				   T6z = FNMS(KP634619297, T66, T69);
				   T6a = FMA(KP634619297, T69, T66);
				   T2l = FNMS(KP559016994, T1W, T1V);
				   T1X = FMA(KP559016994, T1W, T1V);
				   T2g = FNMS(KP618033988, T1k, T1D);
				   T1E = FMA(KP618033988, T1D, T1k);
				   TZ = FNMS(KP250000000, TM, T9);
				   T2m = FNMS(KP618033988, T1Y, T1Z);
				   T20 = FMA(KP618033988, T1Z, T1Y);
			      }
			      ci[0] = T1R + T1U;
			      T2v = FMA(KP951056516, T2m, T2l);
			      T2n = FNMS(KP951056516, T2m, T2l);
			      {
				   E T2b, T21, T2f, T11;
				   T2b = FNMS(KP951056516, T20, T1X);
				   T21 = FMA(KP951056516, T20, T1X);
				   T2f = FNMS(KP559016994, T10, TZ);
				   T11 = FMA(KP559016994, T10, TZ);
				   {
					E T2h, T2r, T27, T1F;
					T2h = FMA(KP951056516, T2g, T2f);
					T2r = FNMS(KP951056516, T2g, T2f);
					T27 = FMA(KP951056516, T1E, T11);
					T1F = FNMS(KP951056516, T1E, T11);
					{
					     E T2o, T2i, T2w, T2s;
					     T2o = T2k * T2h;
					     T2i = T2e * T2h;
					     T2w = T2u * T2r;
					     T2s = T2q * T2r;
					     {
						  E T2c, T28, T22, T1G;
						  T2c = T2a * T27;
						  T28 = T26 * T27;
						  T22 = T1I * T1F;
						  T1G = TY * T1F;
						  ci[WS(rs, 15)] = FMA(T2q, T2v, T2w);
						  cr[WS(rs, 15)] = FNMS(T2u, T2v, T2s);
						  ci[WS(rs, 20)] = FMA(T26, T2b, T2c);
						  cr[WS(rs, 20)] = FNMS(T2a, T2b, T28);
						  ci[WS(rs, 5)] = FMA(TY, T21, T22);
						  cr[WS(rs, 5)] = FNMS(T1I, T21, T1G);
						  cr[WS(rs, 10)] = FNMS(T2k, T2n, T2i);
						  ci[WS(rs, 10)] = FMA(T2e, T2n, T2o);
					     }
					}
				   }
			      }
			 }
			 {
			      E T6x, T6A, T6O, T6P;
			      T6x = FMA(KP968479752, T6w, T6v);
			      T6J = FNMS(KP968479752, T6w, T6v);
			      T6I = FMA(KP845997307, T6z, T6y);
			      T6A = FNMS(KP845997307, T6z, T6y);
			      T6O = FNMS(KP968479752, T6p, T6i);
			      T6q = FMA(KP968479752, T6p, T6i);
			      T6b = FMA(KP845997307, T6a, T63);
			      T6P = FNMS(KP845997307, T6a, T63);
			      T6B = FNMS(KP681693190, T6A, T6x);
			      T73 = FMA(KP560319534, T6x, T6A);
			      T6Q = FMA(KP681693190, T6P, T6O);
			      T78 = FNMS(KP560319534, T6O, T6P);
			 }
			 {
			      E T7U, T8f, T7B, T7u, T82, T8k, T7Y, T7M;
			      {
				   E T7L, T7I, T80, T81;
				   {
					E T7S, T6r, T6t, T6K, T6M, T7T, T6s, T7j;
					T7S = FNMS(KP734762448, T7K, T7J);
					T7L = FMA(KP734762448, T7K, T7J);
					T6r = FMA(KP906616052, T6q, T6b);
					T6t = FNMS(KP906616052, T6q, T6b);
					T6K = FNMS(KP906616052, T6J, T6I);
					T6M = FMA(KP906616052, T6J, T6I);
					T7I = FMA(KP772036680, T7H, T7G);
					T7T = FNMS(KP772036680, T7H, T7G);
					T6s = FNMS(KP249506682, T6r, T5W);
					T7j = FMA(KP998026728, T6r, T5W);
					{
					     E T6L, T7l, T72, T6u;
					     T6L = FNMS(KP249506682, T6K, T6H);
					     T7l = FMA(KP998026728, T6K, T6H);
					     T72 = FMA(KP557913902, T6t, T6s);
					     T6u = FNMS(KP557913902, T6t, T6s);
					     {
						  E T7k, T6N, T77, T7m;
						  T7k = T4l * T7j;
						  T6N = FNMS(KP557913902, T6M, T6L);
						  T77 = FMA(KP557913902, T6M, T6L);
						  T7m = T4l * T7l;
						  {
						       E T74, T7d, T6V, T6C;
						       T74 = FNMS(KP949179823, T73, T72);
						       T7d = FMA(KP949179823, T73, T72);
						       T6V = FMA(KP860541664, T6B, T6u);
						       T6C = FNMS(KP860541664, T6B, T6u);
						       cr[WS(rs, 2)] = FNMS(T4n, T7l, T7k);
						       {
							    E T7h, T79, T6R, T6Z;
							    T7h = FNMS(KP949179823, T78, T77);
							    T79 = FMA(KP949179823, T78, T77);
							    T6R = FNMS(KP860541664, T6Q, T6N);
							    T6Z = FMA(KP860541664, T6Q, T6N);
							    ci[WS(rs, 2)] = FMA(T4n, T7j, T7m);
							    {
								 E T75, T7e, T6W, T6D;
								 T75 = T71 * T74;
								 T7e = T7c * T7d;
								 T6W = T6U * T6V;
								 T6D = T5T * T6C;
								 {
								      E T7a, T7i, T70, T6S;
								      T7a = T71 * T79;
								      T7i = T7c * T7h;
								      T70 = T6U * T6Z;
								      T6S = T5T * T6R;
								      cr[WS(rs, 12)] = FNMS(T76, T79, T75);
								      cr[WS(rs, 17)] = FNMS(T7g, T7h, T7e);
								      cr[WS(rs, 22)] = FNMS(T6Y, T6Z, T6W);
								      cr[WS(rs, 7)] = FNMS(T6E, T6R, T6D);
								      ci[WS(rs, 12)] = FMA(T76, T74, T7a);
								      ci[WS(rs, 17)] = FMA(T7g, T7d, T7i);
								      ci[WS(rs, 22)] = FMA(T6Y, T6V, T70);
								      ci[WS(rs, 7)] = FMA(T6E, T6C, T6S);
								      T7U = FNMS(KP621716863, T7T, T7S);
								      T8f = FMA(KP614372930, T7S, T7T);
								 }
							    }
						       }
						  }
					     }
					}
				   }
				   T80 = FNMS(KP734762448, T7A, T7x);
				   T7B = FMA(KP734762448, T7A, T7x);
				   T7u = FMA(KP772036680, T7t, T7q);
				   T81 = FNMS(KP772036680, T7t, T7q);
				   T82 = FNMS(KP621716863, T81, T80);
				   T8k = FMA(KP614372930, T80, T81);
				   T7Y = FNMS(KP994076283, T7L, T7I);
				   T7M = FMA(KP994076283, T7L, T7I);
			      }
			      {
				   E T5y, T5c, T51, T4U, T5f, T5E, T5o, T5i, T5k;
				   {
					E T5h, T5g, T5m, T5n, T58, T5b;
					T5h = FMA(KP912575812, T57, T56);
					T58 = FNMS(KP912575812, T57, T56);
					T5b = FNMS(KP912018591, T5a, T59);
					T5g = FMA(KP912018591, T5a, T59);
					{
					     E T7X, T7N, T7C, T7Q;
					     T7X = FNMS(KP249506682, T7M, T7F);
					     T7N = FMA(KP998026728, T7M, T7F);
					     T7C = FMA(KP994076283, T7B, T7u);
					     T7Q = FNMS(KP994076283, T7B, T7u);
					     T5y = FMA(KP525970792, T58, T5b);
					     T5c = FNMS(KP726211448, T5b, T58);
					     {
						  E T7Z, T8j, T7P, T7D;
						  T7Z = FNMS(KP557913902, T7Y, T7X);
						  T8j = FMA(KP557913902, T7Y, T7X);
						  T7P = FNMS(KP249506682, T7C, T7n);
						  T7D = FMA(KP998026728, T7C, T7n);
						  {
						       E T8b, T83, T8t, T8l;
						       T8b = FMA(KP943557151, T82, T7Z);
						       T83 = FNMS(KP943557151, T82, T7Z);
						       T8t = FMA(KP949179823, T8k, T8j);
						       T8l = FNMS(KP949179823, T8k, T8j);
						       {
							    E T8e, T7R, T7O, T7E;
							    T8e = FMA(KP557913902, T7Q, T7P);
							    T7R = FNMS(KP557913902, T7Q, T7P);
							    T7O = TR * T7D;
							    T7E = TO * T7D;
							    {
								 E T8g, T8p, T7V, T87;
								 T8g = FMA(KP949179823, T8f, T8e);
								 T8p = FNMS(KP949179823, T8f, T8e);
								 T7V = FMA(KP943557151, T7U, T7R);
								 T87 = FNMS(KP943557151, T7U, T7R);
								 ci[WS(rs, 3)] = FMA(TO, T7N, T7O);
								 cr[WS(rs, 3)] = FNMS(TR, T7N, T7E);
								 {
								      E T8m, T8h, T8u, T8q;
								      T8m = T8i * T8g;
								      T8h = T8d * T8g;
								      T8u = T8s * T8p;
								      T8q = T8o * T8p;
								      {
									   E T84, T7W, T8c, T88;
									   T84 = T4B * T7V;
									   T7W = T4z * T7V;
									   T8c = T8a * T87;
									   T88 = T86 * T87;
									   ci[WS(rs, 13)] = FMA(T8d, T8l, T8m);
									   cr[WS(rs, 13)] = FNMS(T8i, T8l, T8h);
									   ci[WS(rs, 18)] = FMA(T8o, T8t, T8u);
									   cr[WS(rs, 18)] = FNMS(T8s, T8t, T8q);
									   ci[WS(rs, 8)] = FMA(T4z, T83, T84);
									   cr[WS(rs, 8)] = FNMS(T4B, T83, T7W);
									   ci[WS(rs, 23)] = FMA(T86, T8b, T8c);
									   cr[WS(rs, 23)] = FNMS(T8a, T8b, T88);
								      }
								 }
							    }
						       }
						  }
					     }
					}
					T51 = FMA(KP912575812, T50, T4X);
					T5m = FNMS(KP912575812, T50, T4X);
					T5n = FMA(KP912018591, T4T, T4Q);
					T4U = FNMS(KP912018591, T4T, T4Q);
					T41 = FMA(KP951056516, T40, T3X);
					T5f = FNMS(KP951056516, T40, T3X);
					T5E = FMA(KP525970792, T5m, T5n);
					T5o = FNMS(KP726211448, T5n, T5m);
					T5i = FMA(KP851038619, T5h, T5g);
					T5k = FNMS(KP851038619, T5h, T5g);
				   }
				   {
					E T42, T43, T48, T49, T3M, T3P;
					T3M = FMA(KP871714437, T3L, T3K);
					T42 = FNMS(KP871714437, T3L, T3K);
					T43 = FMA(KP831864738, T3O, T3N);
					T3P = FNMS(KP831864738, T3O, T3N);
					{
					     E T5R, T5j, T54, T52;
					     T5R = FMA(KP992114701, T5i, T5f);
					     T5j = FNMS(KP248028675, T5i, T5f);
					     T54 = FNMS(KP851038619, T51, T4U);
					     T52 = FMA(KP851038619, T51, T4U);
					     T3Q = FNMS(KP559154169, T3P, T3M);
					     T4q = FMA(KP683113946, T3M, T3P);
					     {
						  E T5D, T5l, T5P, T53;
						  T5D = FMA(KP554608978, T5k, T5j);
						  T5l = FNMS(KP554608978, T5k, T5j);
						  T5P = FNMS(KP992114701, T52, T4N);
						  T53 = FMA(KP248028675, T52, T4N);
						  {
						       E T5p, T5t, T5F, T5N;
						       T5p = FNMS(KP803003575, T5o, T5l);
						       T5t = FMA(KP803003575, T5o, T5l);
						       T5F = FNMS(KP943557151, T5E, T5D);
						       T5N = FMA(KP943557151, T5E, T5D);
						       {
							    E T55, T5x, T5S, T5Q;
							    T55 = FMA(KP554608978, T54, T53);
							    T5x = FNMS(KP554608978, T54, T53);
							    T5S = TW * T5P;
							    T5Q = TS * T5P;
							    {
								 E T5J, T5z, T5r, T5d;
								 T5J = FMA(KP943557151, T5y, T5x);
								 T5z = FNMS(KP943557151, T5y, T5x);
								 T5r = FMA(KP803003575, T5c, T55);
								 T5d = FNMS(KP803003575, T5c, T55);
								 ci[WS(rs, 4)] = FMA(TS, T5R, T5S);
								 cr[WS(rs, 4)] = FNMS(TW, T5R, T5Q);
								 {
								      E T5G, T5A, T5O, T5K;
								      T5G = T5C * T5z;
								      T5A = T5w * T5z;
								      T5O = T5M * T5J;
								      T5K = T5I * T5J;
								      {
									   E T5q, T5e, T5u, T5s;
									   T5q = TX * T5d;
									   T5e = TT * T5d;
									   T5u = T25 * T5r;
									   T5s = T23 * T5r;
									   ci[WS(rs, 14)] = FMA(T5w, T5F, T5G);
									   cr[WS(rs, 14)] = FNMS(T5C, T5F, T5A);
									   ci[WS(rs, 19)] = FMA(T5I, T5N, T5O);
									   cr[WS(rs, 19)] = FNMS(T5M, T5N, T5K);
									   ci[WS(rs, 9)] = FMA(TT, T5p, T5q);
									   cr[WS(rs, 9)] = FNMS(TX, T5p, T5e);
									   ci[WS(rs, 24)] = FMA(T23, T5t, T5u);
									   cr[WS(rs, 24)] = FNMS(T25, T5t, T5s);
								      }
								 }
							    }
						       }
						  }
					     }
					}
					T48 = FNMS(KP871714437, T39, T2U);
					T3a = FMA(KP871714437, T39, T2U);
					T3F = FMA(KP831864738, T3E, T3p);
					T49 = FNMS(KP831864738, T3E, T3p);
					T4a = FMA(KP559154169, T49, T48);
					T4w = FNMS(KP683113946, T48, T49);
					T46 = FMA(KP904730450, T43, T42);