hb_25.c

快速fft变换

					     {
						  E T4p, T2W, T4o, T33;
						  T4p = FMA(KP951056516, T2V, T2S);
						  T2W = FNMS(KP951056516, T2V, T2S);
						  T4o = FMA(KP951056516, T32, T2Z);
						  T33 = FNMS(KP951056516, T32, T2Z);
						  T70 = FNMS(KP827271945, T6Z, T6Y);
						  T7f = FMA(KP827271945, T6Y, T6Z);
						  T66 = FNMS(KP062914667, T5x, T5A);
						  T5B = FMA(KP062914667, T5A, T5x);
						  T4q = FNMS(KP470564281, T4p, T4o);
						  T4H = FMA(KP470564281, T4o, T4p);
						  T3s = FNMS(KP549754652, T2W, T33);
						  T34 = FMA(KP549754652, T33, T2W);
						  T5E = FNMS(KP951056516, T5D, T5C);
						  T6V = FMA(KP951056516, T5D, T5C);
					     }
					}
				   }
			      }
			 }
		    }
	       }
	       {
		    E T6X, T7e, T6A, T6F, T6C, T6G, T6B;
		    cr[0] = T9 + TM;
		    {
			 E T67, T5I, T25, T22, T1X, T26, T21;
			 {
			      E T1I, T23, T1L, T1Z, T1t, TO, T24, T1O;
			      {
				   E T1K, T6W, T5H, T1J;
				   T1K = T1G - T1H;
				   T1I = T1G + T1H;
				   T6W = FNMS(KP951056516, T5G, T5F);
				   T5H = FMA(KP951056516, T5G, T5F);
				   T1J = FNMS(KP250000000, T1I, T1F);
				   T6X = FMA(KP062914667, T6W, T6V);
				   T7e = FNMS(KP062914667, T6V, T6W);
				   T67 = FNMS(KP634619297, T5E, T5H);
				   T5I = FMA(KP634619297, T5H, T5E);
				   T23 = FNMS(KP559016994, T1K, T1J);
				   T1L = FMA(KP559016994, T1K, T1J);
				   T1Z = FNMS(KP618033988, T19, T1s);
				   T1t = FMA(KP618033988, T1s, T19);
				   TO = FNMS(KP250000000, TM, T9);
				   T24 = FNMS(KP618033988, T1M, T1N);
				   T1O = FMA(KP618033988, T1N, T1M);
			      }
			      {
				   E T2b, T2a, T1Y, TQ, T27;
				   ci[0] = T1F + T1I;
				   T2b = FMA(KP951056516, T24, T23);
				   T25 = FNMS(KP951056516, T24, T23);
				   T2a = W[29];
				   T1Y = FNMS(KP559016994, TP, TO);
				   TQ = FMA(KP559016994, TP, TO);
				   T27 = W[28];
				   {
					E T1V, T1P, T20, T1S, T1w, T1v, TN, T1Q;
					T1V = FNMS(KP951056516, T1O, T1L);
					T1P = FMA(KP951056516, T1O, T1L);
					{
					     E T28, T1u, T29, T2c;
					     T20 = FMA(KP951056516, T1Z, T1Y);
					     T28 = FNMS(KP951056516, T1Z, T1Y);
					     T1S = FMA(KP951056516, T1t, TQ);
					     T1u = FNMS(KP951056516, T1t, TQ);
					     T1w = W[9];
					     T29 = T27 * T28;
					     T2c = T2a * T28;
					     TN = W[8];
					     T1Q = T1w * T1u;
					     cr[WS(rs, 15)] = FNMS(T2a, T2b, T29);
					     ci[WS(rs, 15)] = FMA(T27, T2b, T2c);
					     T1v = TN * T1u;
					}
					ci[WS(rs, 5)] = FMA(TN, T1P, T1Q);
					{
					     E T1U, T1R, T1W, T1T;
					     T1U = W[39];
					     cr[WS(rs, 5)] = FNMS(T1w, T1P, T1v);
					     T1R = W[38];
					     T1W = T1U * T1S;
					     T22 = W[19];
					     T1T = T1R * T1S;
					     T1X = W[18];
					     ci[WS(rs, 20)] = FMA(T1R, T1V, T1W);
					     T26 = T22 * T20;
					     cr[WS(rs, 20)] = FNMS(T1U, T1V, T1T);
					     T21 = T1X * T20;
					}
				   }
			      }
			 }
			 {
			      E T6h, T6g, T5Y, T5J, T6z, T69, T6o, T6E;
			      {
				   E T6m, T6n, T65, T68;
				   T65 = FMA(KP968479752, T64, T63);
				   T6h = FNMS(KP968479752, T64, T63);
				   ci[WS(rs, 10)] = FMA(T1X, T25, T26);
				   T68 = FNMS(KP845997307, T67, T66);
				   T6g = FMA(KP845997307, T67, T66);
				   cr[WS(rs, 10)] = FNMS(T22, T25, T21);
				   T6m = FNMS(KP968479752, T5X, T5Q);
				   T5Y = FMA(KP968479752, T5X, T5Q);
				   T5J = FMA(KP845997307, T5I, T5B);
				   T6n = FNMS(KP845997307, T5I, T5B);
				   T6z = FMA(KP560319534, T65, T68);
				   T69 = FNMS(KP681693190, T68, T65);
				   T6o = FMA(KP681693190, T6n, T6m);
				   T6E = FNMS(KP560319534, T6m, T6n);
			      }
			      {
				   E T62, T6l, T6I, T6L, T6H, T6K;
				   {
					E T6Q, T6O, T6y, T6D, T6S;
					{
					     E T6N, T5Z, T61, T6i, T6k;
					     T6N = W[2];
					     T5Z = FMA(KP906616052, T5Y, T5J);
					     T61 = FNMS(KP906616052, T5Y, T5J);
					     T6i = FNMS(KP906616052, T6h, T6g);
					     T6k = FMA(KP906616052, T6h, T6g);
					     T6Q = W[3];
					     {
						  E T60, T6j, T6R, T6P;
						  T60 = FNMS(KP249506682, T5Z, T5u);
						  T6O = FMA(KP998026728, T5Z, T5u);
						  T6j = FNMS(KP249506682, T6i, T6f);
						  T6R = FMA(KP998026728, T6i, T6f);
						  T6y = FMA(KP557913902, T61, T60);
						  T62 = FNMS(KP557913902, T61, T60);
						  T6P = T6N * T6O;
						  T6l = FNMS(KP557913902, T6k, T6j);
						  T6D = FMA(KP557913902, T6k, T6j);
						  T6S = T6N * T6R;
						  cr[WS(rs, 2)] = FNMS(T6Q, T6R, T6P);
					     }
					}
					T6A = FNMS(KP949179823, T6z, T6y);
					T6I = FMA(KP949179823, T6z, T6y);
					T6L = FNMS(KP949179823, T6E, T6D);
					T6F = FMA(KP949179823, T6E, T6D);
					ci[WS(rs, 2)] = FMA(T6Q, T6O, T6S);
					T6H = W[32];
					T6K = W[33];
				   }
				   {
					E T6a, T6s, T6v, T6p, T6c, T6q, T6b, T6M, T6J, T5r;
					T6a = FNMS(KP860541664, T69, T62);
					T6s = FMA(KP860541664, T69, T62);
					T6v = FMA(KP860541664, T6o, T6l);
					T6p = FNMS(KP860541664, T6o, T6l);
					T6M = T6H * T6L;
					T6J = T6H * T6I;
					T5r = W[12];
					T6c = W[13];
					ci[WS(rs, 17)] = FMA(T6K, T6I, T6M);
					cr[WS(rs, 17)] = FNMS(T6K, T6L, T6J);
					T6q = T5r * T6p;
					T6b = T5r * T6a;
					{
					     E T6r, T6u, T6w, T6t, T6x;
					     ci[WS(rs, 7)] = FMA(T6c, T6a, T6q);
					     cr[WS(rs, 7)] = FNMS(T6c, T6p, T6b);
					     T6r = W[42];
					     T6u = W[43];
					     T6w = T6r * T6v;
					     T6t = T6r * T6s;
					     T6x = W[22];
					     T6C = W[23];
					     ci[WS(rs, 22)] = FMA(T6u, T6s, T6w);
					     cr[WS(rs, 22)] = FNMS(T6u, T6v, T6t);
					     T6G = T6x * T6F;
					     T6B = T6x * T6A;
					}
				   }
			      }
			 }
		    }
		    {
			 E T7u, T7D, T7n, T7w, T7v, T7E;
			 {
			      E T78, T7t, T7N, T71, T7C, T7S, T7y, T7k;
			      {
				   E T7j, T7g, T7A, T7B, T7r, T7s;
				   T7r = FNMS(KP734762448, T7i, T7h);
				   T7j = FMA(KP734762448, T7i, T7h);
				   T7g = FMA(KP772036680, T7f, T7e);
				   T7s = FNMS(KP772036680, T7f, T7e);
				   ci[WS(rs, 12)] = FMA(T6C, T6A, T6G);
				   cr[WS(rs, 12)] = FNMS(T6C, T6F, T6B);
				   T7A = FNMS(KP734762448, T77, T74);
				   T78 = FMA(KP734762448, T77, T74);
				   T7t = FNMS(KP621716863, T7s, T7r);
				   T7N = FMA(KP614372930, T7r, T7s);
				   T71 = FMA(KP772036680, T70, T6X);
				   T7B = FNMS(KP772036680, T70, T6X);
				   T7C = FNMS(KP621716863, T7B, T7A);
				   T7S = FMA(KP614372930, T7A, T7B);
				   T7y = FNMS(KP994076283, T7j, T7g);
				   T7k = FMA(KP994076283, T7j, T7g);
			      }
			      {
				   E T7c, T6T, T7x, T7l, T79, T7p;
				   T7c = W[5];
				   T6T = W[4];
				   T7x = FNMS(KP249506682, T7k, T7d);
				   T7l = FMA(KP998026728, T7k, T7d);
				   T79 = FMA(KP994076283, T78, T71);
				   T7p = FNMS(KP994076283, T78, T71);
				   {
					E T7z, T7Y, T7Z, T7T, T7q, T7O, T7X, T7L, T7Q, T7P, T7U;
					{
					     E T7V, T80, T7b, T7m, T7W;
					     {
						  E T7R, T7o, T7a, T7M;
						  T7V = W[34];
						  T7R = FMA(KP557913902, T7y, T7x);
						  T7z = FNMS(KP557913902, T7y, T7x);
						  T7Y = W[35];
						  T7o = FNMS(KP249506682, T79, T6U);
						  T7a = FMA(KP998026728, T79, T6U);
						  T7Z = FMA(KP949179823, T7S, T7R);
						  T7T = FNMS(KP949179823, T7S, T7R);
						  T7M = FMA(KP557913902, T7p, T7o);
						  T7q = FNMS(KP557913902, T7p, T7o);
						  T7b = T6T * T7a;
						  T7m = T7c * T7a;
						  T7W = FNMS(KP949179823, T7N, T7M);
						  T7O = FMA(KP949179823, T7N, T7M);
					     }
					     cr[WS(rs, 3)] = FNMS(T7c, T7l, T7b);
					     ci[WS(rs, 3)] = FMA(T6T, T7l, T7m);
					     T80 = T7Y * T7W;
					     T7X = T7V * T7W;
					     T7L = W[24];
					     T7Q = W[25];
					     ci[WS(rs, 18)] = FMA(T7V, T7Z, T80);
					}
					cr[WS(rs, 18)] = FNMS(T7Y, T7Z, T7X);
					T7P = T7L * T7O;
					T7U = T7Q * T7O;
					{
					     E T7J, T7F, T7I, T7H, T7K, T7G;
					     T7u = FMA(KP943557151, T7t, T7q);
					     T7G = FNMS(KP943557151, T7t, T7q);
					     cr[WS(rs, 13)] = FNMS(T7Q, T7T, T7P);
					     ci[WS(rs, 13)] = FMA(T7L, T7T, T7U);
					     T7J = FMA(KP943557151, T7C, T7z);
					     T7D = FNMS(KP943557151, T7C, T7z);
					     T7F = W[44];
					     T7I = W[45];
					     T7n = W[14];
					     T7H = T7F * T7G;
					     T7K = T7I * T7G;
					     T7w = W[15];
					     T7v = T7n * T7u;
					     cr[WS(rs, 23)] = FNMS(T7I, T7J, T7H);
					     ci[WS(rs, 23)] = FMA(T7F, T7J, T7K);
					}
				   }
			      }
			 }
			 T7E = T7w * T7u;
			 cr[WS(rs, 8)] = FNMS(T7w, T7D, T7v);
			 {
			      E T3F, T4K, T4X, T4j, T4M, T4L, T4Y;
			      {
				   E T4P, T4O, T4y, T4r, T4J, T57, T4N, T5c, T4W;
				   {
					E T4U, T4V, T4F, T4I;
					T4F = FNMS(KP912575812, T4E, T4D);
					T4P = FMA(KP912575812, T4E, T4D);
					T4O = FMA(KP912018591, T4H, T4G);
					T4I = FNMS(KP912018591, T4H, T4G);
					ci[WS(rs, 8)] = FMA(T7n, T7D, T7E);
					T4y = FMA(KP912575812, T4x, T4u);
					T4U = FNMS(KP912575812, T4x, T4u);
					T4V = FMA(KP912018591, T4q, T4n);
					T4r = FNMS(KP912018591, T4q, T4n);
					T4J = FNMS(KP726211448, T4I, T4F);
					T57 = FMA(KP525970792, T4F, T4I);
					T3F = FMA(KP951056516, T3E, T3B);
					T4N = FNMS(KP951056516, T3E, T3B);
					T5c = FMA(KP525970792, T4U, T4V);
					T4W = FNMS(KP726211448, T4V, T4U);
				   }
				   {
					E T5o, T4S, T4B, T5l, T5p, T4R, T4A, T5m, T4Q, T4z;
					T5o = W[7];
					T4Q = FMA(KP851038619, T4P, T4O);
					T4S = FNMS(KP851038619, T4P, T4O);
					T4z = FMA(KP851038619, T4y, T4r);
					T4B = FNMS(KP851038619, T4y, T4r);
					T5l = W[6];
					T5p = FMA(KP992114701, T4Q, T4N);
					T4R = FNMS(KP248028675, T4Q, T4N);
					T4A = FMA(KP248028675, T4z, T4k);
					T5m = FNMS(KP992114701, T4z, T4k);
					{
					     E T4T, T4C, T5d, T58, T55, T5a, T59, T5e;
					     {
						  E T5f, T5j, T5i, T5h, T5k, T5g;
						  T5f = W[36];
						  {
						       E T5b, T56, T5n, T5q;
						       T4T = FNMS(KP554608978, T4S, T4R);
						       T5b = FMA(KP554608978, T4S, T4R);
						       T56 = FNMS(KP554608978, T4B, T4A);
						       T4C = FMA(KP554608978, T4B, T4A);
						       T5n = T5l * T5m;
						       T5q = T5o * T5m;
						       T5j = FMA(KP943557151, T5c, T5b);
						       T5d = FNMS(KP943557151, T5c, T5b);
						       T5g = FMA(KP943557151, T57, T56);
						       T58 = FNMS(KP943557151, T57, T56);
						       cr[WS(rs, 4)] = FNMS(T5o, T5p, T5n);
						       ci[WS(rs, 4)] = FMA(T5l, T5p, T5q);
						  }
						  T5i = W[37];
						  T5h = T5f * T5g;
						  T55 = W[26];
						  T5k = T5i * T5g;
						  T5a = W[27];
						  cr[WS(rs, 19)] = FNMS(T5i, T5j, T5h);
						  T59 = T55 * T58;
						  ci[WS(rs, 19)] = FMA(T5f, T5j, T5k);
					     }
					     T5e = T5a * T58;
					     {
						  E T53, T4Z, T52, T51, T54, T50;
						  cr[WS(rs, 14)] = FNMS(T5a, T5d, T59);
						  T4K = FNMS(KP803003575, T4J, T4C);
						  T50 = FMA(KP803003575, T4J, T4C);
						  ci[WS(rs, 14)] = FMA(T55, T5d, T5e);
						  T4X = FNMS(KP803003575, T4W, T4T);
						  T53 = FMA(KP803003575, T4W, T4T);
						  T4Z = W[46];
						  T52 = W[47];
						  T4j = W[16];
						  T51 = T4Z * T50;
						  T54 = T52 * T50;
						  T4M = W[17];
						  T4L = T4j * T4K;
						  cr[WS(rs, 24)] = FNMS(T52, T53, T51);
						  ci[WS(rs, 24)] = FMA(T4Z, T53, T54);
					     }
					}
				   }
			      }
			      T4Y = T4M * T4K;
			      cr[WS(rs, 9)] = FNMS(T4M, T4X, T4L);
			      {
				   E T3G, T3H, T2P, T3k, T3Z, T3v, T3O, T44;
				   {
					E T3M, T3N, T3r, T3u;
					T3G = FNMS(KP871714437, T3q, T3p);
					T3r = FMA(KP871714437, T3q, T3p);
					T3u = FNMS(KP831864738, T3t, T3s);
					T3H = FMA(KP831864738, T3t, T3s);
					ci[WS(rs, 9)] = FMA(T4j, T4X, T4Y);
					T3M = FNMS(KP871714437, T2O, T2z);
					T2P = FMA(KP871714437, T2O, T2z);
					T3k = FMA(KP831864738, T3j, T34);
					T3N = FNMS(KP831864738, T3j, T34);
					T3Z = FMA(KP683113946, T3r, T3u);
					T3v = FNMS(KP559154169, T3u, T3r);
					T3O = FMA(KP559154169, T3N, T3M);
					T44 = FNMS(KP683113946, T3M, T3N);
				   }
				   {
					E T4g, T3K, T3n, T4d, T3J, T4h, T4e, T3m, T3I, T3l;
					T4g = W[1];
					T3K = FMA(KP904730450, T3H, T3G);
					T3I = FNMS(KP904730450, T3H, T3G);
					T3n = FNMS(KP904730450, T3k, T2P);
					T3l = FMA(KP904730450, T3k, T2P);
					T4d = W[0];
					T3J = FNMS(KP242145790, T3I, T3F);
					T4h = FMA(KP968583161, T3I, T3F);
					T4e = FMA(KP968583161, T3l, T2k);
					T3m = FNMS(KP242145790, T3l, T2k);
					{
					     E T3L, T3o, T45, T40, T3X, T42, T41, T46;
					     {
						  E T47, T4b, T4a, T49, T4c, T48;
						  T47 = W[30];
						  {
						       E T43, T3Y, T4f, T4i;
						       T43 = FNMS(KP541454447, T3K, T3J);
						       T3L = FMA(KP541454447, T3K, T3J);
						       T3o = FMA(KP541454447, T3n, T3m);
						       T3Y = FNMS(KP541454447, T3n, T3m);
						       T4f = T4d * T4e;
						       T4i = T4g * T4e;
						       T45 = FNMS(KP833417178, T44, T43);
						       T4b = FMA(KP833417178, T44, T43);
						       T40 = FNMS(KP833417178, T3Z, T3Y);
						       T48 = FMA(KP833417178, T3Z, T3Y);
						       cr[WS(rs, 1)] = FNMS(T4g, T4h, T4f);
						       ci[WS(rs, 1)] = FMA(T4d, T4h, T4i);
						  }
						  T4a = W[31];