hf2_25.c

快速fft变换

						  T3u = T7 * T1y;
					     }
					     T6y = FNMS(KP559016994, T57, T56);
					     T58 = FMA(KP559016994, T57, T56);
					     T6B = FNMS(KP559016994, T50, T4Z);
					     T51 = FMA(KP559016994, T50, T4Z);
					     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 = cr[WS(rs, 14)];
						  T1P = ci[WS(rs, 14)];
						  T1T = cr[WS(rs, 19)];
						  T1X = ci[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 = cr[WS(rs, 9)];
						  T1C = ci[WS(rs, 9)];
						  T5i = T1Y - T1Q;
						  T1Z = T1Q + T1Y;
						  T3F = T3C + T3E;
						  T5p = T3E - T3C;
						  T1B = T8 * T1A;
						  T3w = T8 * T1C;
						  T1E = cr[WS(rs, 24)];
						  T1G = ci[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, T5m, T5f, T5l, T5e, T3J, T1I, T5h;
					     T23 = cr[WS(rs, 2)];
					     T25 = ci[WS(rs, 2)];
					     T1I = T1D + T1H;
					     T5h = T1H - T1D;
					     {
						  E T3A, T5o, T20, T3G;
						  T3A = T3x + T3z;
						  T5o = T3z - T3x;
						  T24 = T19 * T23;
						  T6v = FNMS(KP618033988, T5h, T5i);
						  T5j = FMA(KP618033988, T5i, T5h);
						  T20 = T1I + T1Z;
						  T5m = T1I - T1Z;
						  T6s = FNMS(KP618033988, T5o, T5p);
						  T5q = FMA(KP618033988, T5p, T5o);
						  T3G = T3A + T3F;
						  T5f = T3F - T3A;
						  T21 = T1z + T20;
						  T5l = FNMS(KP250000000, T20, T1z);
						  T3H = T3v + T3G;
						  T5e = FNMS(KP250000000, T3G, T3v);
						  T3J = T19 * T25;
					     }
					     T6r = FNMS(KP559016994, T5m, T5l);
					     T5n = FMA(KP559016994, T5m, T5l);
					     T6u = FMA(KP559016994, T5f, T5e);
					     T5g = FNMS(KP559016994, T5f, T5e);
					     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 = cr[WS(rs, 12)];
						  T2q = ci[WS(rs, 12)];
						  T2u = cr[WS(rs, 17)];
						  T2y = ci[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 = cr[WS(rs, 7)];
						  T2b = ci[WS(rs, 7)];
						  T4N = T2z - T2r;
						  T2A = T2r + T2z;
						  T3U = T3R + T3T;
						  T4U = T3R - T3T;
						  T29 = T27 * T28;
						  T3L = T27 * T2b;
						  T2f = cr[WS(rs, 22)];
						  T2j = ci[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 T8k, T6d, T6g, T8r, T6f, T8l, T6c, T8q, T69, T7r, T5Y, T8g, T8i, T66, T68;
			 E T5X, T8d, T8h;
			 {
			      E T4O, T4V, T22, T4S, T4L, T3b, T4e, T4c, T3I;
			      T8k = T3t + T3H;
			      T3I = T3t - T3H;
			      {
				   E T2l, T4M, T3P, T4T;
				   T2l = T2c + T2k;
				   T4M = T2k - T2c;
				   T3P = T3M + T3O;
				   T4T = T3O - T3M;
				   T4O = FMA(KP618033988, T4N, T4M);
				   T6d = FNMS(KP618033988, T4M, T4N);
				   {
					E T4R, T2B, T4K, T3V;
					T4R = T2A - T2l;
					T2B = T2l + T2A;
					T4V = FNMS(KP618033988, T4U, T4T);
					T6g = FMA(KP618033988, T4T, T4U);
					T4K = T3U - T3P;
					T3V = T3P + T3U;
					{
					     E T4Q, T2C, T4J, T3W, T4b;
					     T4Q = FNMS(KP250000000, T2B, T26);
					     T2C = T26 + T2B;
					     T4J = FNMS(KP250000000, T3V, T3K);
					     T3W = T3K + T3V;
					     T8r = T21 - T1v;
					     T22 = T1v + T21;
					     T4S = FNMS(KP559016994, T4R, T4Q);
					     T6f = FMA(KP559016994, T4R, T4Q);
					     T4b = T3W - T4a;
					     T8l = T3W + T4a;
					     T6c = FMA(KP559016994, T4K, T4J);
					     T4L = FNMS(KP559016994, T4K, T4J);
					     T8q = T2C - T3a;
					     T3b = T2C + T3a;
					     T4e = FNMS(KP618033988, T3I, T4b);
					     T4c = FMA(KP618033988, T4b, T3I);
					}
				   }
			      }
			      {
				   E T5H, T4t, T7V, T87, T5Q, T5P, T5D, T8e, T5A, T8f, T5K, T60, T8c, T8a, T5u;
				   E T5w, T5U, T64, T5N, T61;
				   {
					E T3e, T3d, T4h, T3c, T7T;
					T4h = FMA(KP559016994, T4g, T4f);
					T69 = FNMS(KP559016994, T4g, T4f);
					T3c = T22 + T3b;
					T3e = T22 - T3b;
					T7r = FNMS(KP559016994, T7q, T7p);
					T7T = FMA(KP559016994, T7q, T7p);
					T5H = FMA(KP951056516, T4s, T4h);
					T4t = FNMS(KP951056516, T4s, T4h);
					cr[0] = TN + T3c;
					T3d = FNMS(KP250000000, T3c, TN);
					T7V = FNMS(KP951056516, T7U, T7T);
					T87 = FMA(KP951056516, T7U, T7T);
					{
					     E T5S, T5T, T5L, T4I, T5B, T5M, T55, T5J, T5s, T5z, T4X, T5C, T5I, T5c;
					     {
						  E T5k, T5r, T4P, T4W;
						  {
						       E T4A, T4d, T3f, T4H;
						       T4A = FMA(KP951056516, T4z, T4w);
						       T5S = FNMS(KP951056516, T4z, T4w);
						       T4d = FNMS(KP559016994, T3e, T3d);
						       T3f = FMA(KP559016994, T3e, T3d);
						       T5T = FNMS(KP951056516, T4G, T4D);
						       T4H = FMA(KP951056516, T4G, T4D);
						       T5k = FNMS(KP951056516, T5j, T5g);
						       T5L = FMA(KP951056516, T5j, T5g);
						       cr[WS(rs, 5)] = FMA(KP951056516, T4c, T3f);
						       ci[WS(rs, 4)] = FNMS(KP951056516, T4c, T3f);
						       ci[WS(rs, 9)] = FMA(KP951056516, T4e, T4d);
						       cr[WS(rs, 10)] = FNMS(KP951056516, T4e, T4d);
						       T4I = FNMS(KP126329378, T4H, T4A);
						       T5B = FMA(KP126329378, T4A, T4H);
						       T5M = FNMS(KP951056516, T5q, T5n);
						       T5r = FMA(KP951056516, T5q, T5n);
						  }
						  T4P = FNMS(KP951056516, T4O, T4L);
						  T5Q = FMA(KP951056516, T4O, T4L);
						  T5P = FNMS(KP951056516, T4V, T4S);
						  T4W = FMA(KP951056516, T4V, T4S);
						  T55 = FNMS(KP951056516, T54, T51);
						  T5J = FMA(KP951056516, T54, T51);
						  T5s = FMA(KP827271945, T5r, T5k);
						  T5z = FNMS(KP827271945, T5k, T5r);
						  T4X = FNMS(KP470564281, T4W, T4P);
						  T5C = FMA(KP470564281, T4P, T4W);
						  T5I = FMA(KP951056516, T5b, T58);
						  T5c = FNMS(KP951056516, T5b, T58);
					     }
					     {
						  E T88, T4Y, T5d, T5y, T89, T5t;
						  T5D = FNMS(KP912018591, T5C, T5B);
						  T88 = FMA(KP912018591, T5C, T5B);
						  T8e = FMA(KP912018591, T4X, T4I);
						  T4Y = FNMS(KP912018591, T4X, T4I);
						  T5d = FMA(KP634619297, T5c, T55);
						  T5y = FNMS(KP634619297, T55, T5c);
						  T5A = FMA(KP912575812, T5z, T5y);
						  T89 = FNMS(KP912575812, T5z, T5y);
						  T8f = FMA(KP912575812, T5s, T5d);
						  T5t = FNMS(KP912575812, T5s, T5d);
						  T5K = FMA(KP256756360, T5J, T5I);
						  T60 = FNMS(KP256756360, T5I, T5J);
						  T8c = FNMS(KP851038619, T89, T88);
						  T8a = FMA(KP851038619, T89, T88);
						  T5u = FNMS(KP851038619, T5t, T4Y);
						  T5w = FMA(KP851038619, T5t, T4Y);
					     }
					     T5U = FMA(KP939062505, T5T, T5S);
					     T64 = FNMS(KP939062505, T5S, T5T);
					     T5N = FMA(KP634619297, T5M, T5L);
					     T61 = FNMS(KP634619297, T5L, T5M);
					}
				   }
				   {
					E T62, T7W, T83, T5O, T5R, T63;
					cr[WS(rs, 4)] = FNMS(KP992114701, T5u, T4t);
					T62 = FMA(KP871714437, T61, T60);
					T7W = FNMS(KP871714437, T61, T60);
					T83 = FNMS(KP871714437, T5N, T5K);
					T5O = FMA(KP871714437, T5N, T5K);
					T5R = FMA(KP549754652, T5Q, T5P);
					T63 = FNMS(KP549754652, T5P, T5Q);
					ci[WS(rs, 20)] = FNMS(KP992114701, T8a, T87);
					{
					     E T65, T5W, T84, T86, T81, T85, T8b;
					     {
						  E T5E, T5G, T82, T80, T7Y, T5v, T7X, T5V, T5F, T5x, T7Z;
						  T5E = FNMS(KP726211448, T5D, T5A);
						  T5G = FMA(KP525970792, T5A, T5D);
						  T65 = FNMS(KP831864738, T64, T63);
						  T7X = FMA(KP831864738, T64, T63);
						  T82 = FNMS(KP831864738, T5U, T5R);
						  T5V = FMA(KP831864738, T5U, T5R);
						  T80 = FNMS(KP904730450, T7X, T7W);
						  T7Y = FMA(KP904730450, T7X, T7W);
						  T5Y = FNMS(KP904730450, T5V, T5O);
						  T5W = FMA(KP904730450, T5V, T5O);
						  T5v = FMA(KP248028675, T5u, T4t);
						  ci[WS(rs, 23)] = FMA(KP968583161, T7Y, T7V);
						  cr[WS(rs, 1)] = FMA(KP968583161, T5W, T5H);
						  T84 = FNMS(KP683113946, T83, T82);
						  T86 = FMA(KP559154169, T82, T83);
						  T5F = FNMS(KP554608978, T5w, T5v);
						  T5x = FMA(KP554608978, T5w, T5v);
						  T7Z = FNMS(KP242145790, T7Y, T7V);
						  ci[WS(rs, 10)] = FNMS(KP943557151, T5G, T5F);
						  ci[WS(rs, 5)] = FMA(KP943557151, T5G, T5F);
						  ci[0] = FMA(KP803003575, T5E, T5x);
						  cr[WS(rs, 9)] = FNMS(KP803003575, T5E, T5x);
						  T81 = FNMS(KP541454447, T80, T7Z);
						  T85 = FMA(KP541454447, T80, T7Z);
					     }
					     T8g = FNMS(KP525970792, T8f, T8e);
					     T8i = FMA(KP726211448, T8e, T8f);
					     ci[WS(rs, 13)] = FMA(KP833417178, T84, T81);
					     cr[WS(rs, 16)] = FMS(KP833417178, T84, T81);
					     cr[WS(rs, 21)] = -(FMA(KP921177326, T86, T85));
					     ci[WS(rs, 18)] = FNMS(KP921177326, T86, T85);
					     T8b = FMA(KP248028675, T8a, T87);
					     T66 = FMA(KP559154169, T65, T62);
					     T68 = FNMS(KP683113946, T62, T65);
					     T5X = FNMS(KP242145790, T5W, T5H);
					     T8d = FNMS(KP554608978, T8c, T8b);
					     T8h = FMA(KP554608978, T8c, T8b);
					}
				   }
			      }
			 }
			 {
			      E T8s, T8u, T5Z, T67;
			      cr[WS(rs, 24)] = -(FMA(KP803003575, T8i, T8h));
			      ci[WS(rs, 15)] = FNMS(KP803003575, T8i, T8h);
			      cr[WS(rs, 19)] = FMS(KP943557151, T8g, T8d);
			      cr[WS(rs, 14)] = -(FMA(KP943557151, T8g, T8d));
			      T5Z = FMA(KP541454447, T5Y, T5X);
			      T67 = FNMS(KP541454447, T5Y, T5X);
			      cr[WS(rs, 11)] = FNMS(KP833417178, T68, T67);
			      ci[WS(rs, 8)] = FMA(KP833417178, T68, T67);
			      cr[WS(rs, 6)] = FMA(KP921177326, T66, T5Z);
			      ci[WS(rs, 3)] = FNMS(KP921177326, T66, T5Z);
			      T8s = FMA(KP618033988, T8r, T8q);
			      T8u = FNMS(KP618033988, T8q, T8r);
			      {
				   E T6X, T6T, T6b, T7H, T7v, T6Y, T72, T71, T6P, T7O, T6M, T7P, T7K, T6G, T6I;
				   E T6W, T7f, T7d, T76;
				   {
					E T74, T75, T6i, T6N, T6L, T6E, T6U, T6l, T6o, T6V, T6t, T6w;
					{
					     E T6e, T8o, T8n, T6h, T8m;
					     T6X = FNMS(KP951056516, T6d, T6c);
					     T6e = FMA(KP951056516, T6d, T6c);
					     T8o = T8k - T8l;
					     T8m = T8k + T8l;
					     T6T = FNMS(KP951056516, T6a, T69);
					     T6b = FMA(KP951056516, T6a, T69);
					     T7H = FNMS(KP951056516, T7u, T7r);
					     T7v = FMA(KP951056516, T7u, T7r);
					     ci[WS(rs, 24)] = T8m + T8j;
					     T8n = FNMS(KP250000000, T8m, T8j);
					     T6h = FMA(KP951056516, T6g, T6f);
					     T6Y = FNMS(KP951056516, T6g, T6f);
					     {
						  E T6A, T6D, T8t, T8p;
						  T74 = FMA(KP951056516, T6z, T6y);
						  T6A = FNMS(KP951056516, T6z, T6y);
						  T6D = FMA(KP951056516, T6C, T6B);
						  T75 = FNMS(KP951056516, T6C, T6B);
						  T8t = FMA(KP559016994, T8o, T8n);
						  T8p = FNMS(KP559016994, T8o, T8n);
						  T6i = FMA(KP062914667, T6h, T6e);
						  T6N = FNMS(KP062914667, T6e, T6h);
						  ci[WS(rs, 14)] = FMA(KP951056516, T8s, T8p);
						  cr[WS(rs, 15)] = FMS(KP951056516, T8s, T8p);
						  ci[WS(rs, 19)] = FMA(KP951056516, T8u, T8t);
						  cr[WS(rs, 20)] = FMS(KP951056516, T8u, T8t);
						  T6L = FNMS(KP939062505, T6A, T6D);
						  T6E = FMA(KP939062505, T6D, T6A);
					     }
					}
					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);
					T6w = FMA(KP951056516, T6v, T6u);
					T71 = FNMS(KP951056516, T6v, T6u);
					{
					     E T6q, T6F, T6O, T6p;
					     T6O = FMA(KP827271945, T6l, T6o);
					     T6p = FNMS(KP827271945, T6o, T6l);
					     {
						  E T6K, T6x, T7I, T7J;
						  T6K = FMA(KP126329378, T6t, T6w);
						  T6x = FNMS(KP126329378, T6w, T6t);
						  T7I = FMA(KP772036680, T6O, T6N);
						  T6P = FNMS(KP772036680, T6O, T6N);
						  T6q = FMA(KP772036680, T6p, T6i);
						  T7O = FNMS(KP772036680, T6p, T6i);
						  T7J = FNMS(KP734762448, T6L, T6K);
						  T6M = FMA(KP734762448, T6L, T6K);
						  T6F = FNMS(KP734762448, T6E, T6x);
						  T7P = FMA(KP734762448, T6E, T6x);
						  T7K = FMA(KP994076283, T7J, T7I);
						  T7M = FNMS(KP994076283, T7J, T7I);
					     }
					     T6G = FNMS(KP994076283, T6F, T6q);
					     T6I = FMA(KP994076283, T6F, T6q);
					}
					T6W = FMA(KP062914667, T6V, T6U);
					T7f = FNMS(KP062914667, T6U, T6V);
					T7d = FNMS(KP549754652, T74, T75);
					T76 = FMA(KP549754652, T75, T74);
				   }
				   {
					E T7h, T7C, T7e, T7D, T7y, T7A, T78, T7a;
					{
					     E T70, T77, T7g, T6Z;
					     cr[WS(rs, 3)] = FMA(KP998026728, T6G, T6b);
					     T7g = FNMS(KP634619297, T6X, T6Y);
					     T6Z = FMA(KP634619297, T6Y, T6X);
					     {
						  E T7c, T73, T7w, T7x;
						  T7c = FMA(KP470564281, T71, T72);
						  T73 = FNMS(KP470564281, T72, T71);
						  T7w = FMA(KP845997307, T7g, T7f);
						  T7h = FNMS(KP845997307, T7g, T7f);
						  T70 = FMA(KP845997307, T6Z, T6W);
						  T7C = FNMS(KP845997307, T6Z, T6W);
						  T7x = FNMS(KP968479752, T7d, T7c);
						  T7e = FMA(KP968479752, T7d, T7c);
						  T77 = FMA(KP968479752, T76, T73);
						  T7D = FNMS(KP968479752, T76, T73);
						  T7y = FMA(KP906616052, T7x, T7w);
						  T7A = FNMS(KP906616052, T7x, T7w);