hc2cb2_32.c

快速fft变换

				   {
					E T1G, T1H, Tg, Th, Tj, Tk;
					Tg = Rp[WS(rs, 2)];
					Th = Rm[WS(rs, 13)];
					T1G = Ip[WS(rs, 2)];
					T6Y = T6s + T6z;
					T6A = T6s - T6z;
					T1K = Tg - Th;
					Ti = Tg + Th;
					T1H = Im[WS(rs, 13)];
					Tj = Rp[WS(rs, 10)];
					Tk = Rm[WS(rs, 5)];
					T1L = Ip[WS(rs, 10)];
					T6n = T1G - T1H;
					T1I = T1G + T1H;
					T1F = Tj - Tk;
					Tl = Tj + Tk;
					T1M = Im[WS(rs, 5)];
				   }
				   {
					E T1v, T1w, Tq, Tr;
					{
					     E Tn, T1N, T6o, To;
					     Tn = Rm[WS(rs, 1)];
					     T4w = T1I - T1F;
					     T1J = T1F + T1I;
					     T7Z = Ti - Tl;
					     Tm = Ti + Tl;
					     T1N = T1L + T1M;
					     T6o = T1L - T1M;
					     To = Rp[WS(rs, 14)];
					     T1v = Ip[WS(rs, 14)];
					     T6p = T6n + T6o;
					     T80 = T6n - T6o;
					     T4x = T1K + T1N;
					     T1O = T1K - T1N;
					     T1z = Tn - To;
					     Tp = Tn + To;
					     T1w = Im[WS(rs, 1)];
					}
					Tq = Rp[WS(rs, 6)];
					Tr = Rm[WS(rs, 9)];
					T1A = Ip[WS(rs, 6)];
					T6k = T1v - T1w;
					T1x = T1v + T1w;
					T1u = Tq - Tr;
					Ts = Tq + Tr;
					T1B = Im[WS(rs, 9)];
				   }
			      }
			      {
				   E T4z, T6m, T4A, T2Z, T1E, T1P, T2Y, T50, T51;
				   {
					E T1y, T81, T8q, T1D, T7Y, T8r;
					{
					     E T7X, Tt, T1C, T6l, T7W;
					     T4z = T1u + T1x;
					     T1y = T1u - T1x;
					     T7X = Tp - Ts;
					     Tt = Tp + Ts;
					     T1C = T1A + T1B;
					     T6l = T1A - T1B;
					     T81 = T7Z + T80;
					     T8q = T7Z - T80;
					     T6m = T6k + T6l;
					     T7W = T6k - T6l;
					     T4A = T1z + T1C;
					     T1D = T1z - T1C;
					     T6P = Tm - Tt;
					     Tu = Tm + Tt;
					     T7Y = T7W - T7X;
					     T8r = T7X + T7W;
					}
					T2Z = FMA(KP414213562, T1y, T1D);
					T1E = FNMS(KP414213562, T1D, T1y);
					T9a = T81 + T7Y;
					T82 = T7Y - T81;
					T8X = T8q + T8r;
					T8s = T8q - T8r;
					T1P = FMA(KP414213562, T1O, T1J);
					T2Y = FNMS(KP414213562, T1J, T1O);
				   }
				   T4y = FNMS(KP414213562, T4x, T4w);
				   T50 = FMA(KP414213562, T4w, T4x);
				   T40 = T1P + T1E;
				   T1Q = T1E - T1P;
				   T3M = T2Y + T2Z;
				   T30 = T2Y - T2Z;
				   T51 = FMA(KP414213562, T4z, T4A);
				   T4B = FNMS(KP414213562, T4A, T4z);
				   T5w = T50 + T51;
				   T52 = T50 - T51;
				   T7u = T6p + T6m;
				   T6q = T6m - T6p;
			      }
			 }
		    }
		    {
			 E T7D, T7K, T7J, T5K, T4C, T7E, T83, T8w, T8t, T8i, T6r, T70, T6X, T6K;
			 {
			      E T8Y, T9e, T9b, T95, T8F, T8G, T8L, T8M;
			      {
				   E T7v, T7p, T7y, Tv, T10;
				   T7D = Tf - Tu;
				   Tv = Tf + Tu;
				   T10 = TK + TZ;
				   T7K = TK - TZ;
				   T7J = T7t - T7u;
				   T7v = T7t + T7u;
				   T5K = T4B - T4y;
				   T4C = T4y + T4B;
				   T7p = Tv - T10;
				   T7E = T7x - T7w;
				   T7y = T7w + T7x;
				   Rp[0] = Tv + T10;
				   {
					E T9p, T9x, T9z, T9v;
					{
					     E T9n, T7A, T7q, T7z, T9o, T9t, T9u;
					     T8Y = FNMS(KP707106781, T8X, T8W);
					     T9n = FMA(KP707106781, T8X, T8W);
					     T7A = T7s * T7p;
					     T7q = T7o * T7p;
					     Rm[0] = T7v + T7y;
					     T7z = T7v - T7y;
					     T9o = T9c + T9d;
					     T9e = T9c - T9d;
					     T9b = FNMS(KP707106781, T9a, T99);
					     T9t = FMA(KP707106781, T9a, T99);
					     T9u = T94 + T91;
					     T95 = T91 - T94;
					     Rm[WS(rs, 8)] = FMA(T7o, T7z, T7A);
					     Rp[WS(rs, 8)] = FNMS(T7s, T7z, T7q);
					     T9p = FNMS(KP923879532, T9o, T9n);
					     T9x = FMA(KP923879532, T9o, T9n);
					     T9z = FMA(KP923879532, T9u, T9t);
					     T9v = FNMS(KP923879532, T9u, T9t);
					}
					{
					     E T9y, T9q, T9w, T9A;
					     T9y = T3v * T9x;
					     T9q = T9m * T9p;
					     T9w = T9m * T9v;
					     T9A = T3v * T9z;
					     Rp[WS(rs, 1)] = FNMS(T3x, T9z, T9y);
					     Rp[WS(rs, 9)] = FNMS(T9s, T9v, T9q);
					     Rm[WS(rs, 9)] = FMA(T9s, T9p, T9w);
					     Rm[WS(rs, 1)] = FMA(T3x, T9x, T9A);
					}
				   }
				   T83 = FMA(KP707106781, T82, T7V);
				   T8F = FNMS(KP707106781, T82, T7V);
				   T8G = T8u + T8v;
				   T8w = T8u - T8v;
				   T8t = FMA(KP707106781, T8s, T8p);
				   T8L = FNMS(KP707106781, T8s, T8p);
				   T8M = T8h + T8a;
				   T8i = T8a - T8h;
			      }
			      {
				   E T79, T7a, T7f, T7g;
				   T6r = T6j + T6q;
				   T79 = T6j - T6q;
				   {
					E T8Q, T8H, T8T, T8N;
					T8Q = FMA(KP923879532, T8G, T8F);
					T8H = FNMS(KP923879532, T8G, T8F);
					T8T = FMA(KP923879532, T8M, T8L);
					T8N = FNMS(KP923879532, T8M, T8L);
					{
					     E T8R, T8I, T8U, T8O;
					     T8R = T8P * T8Q;
					     T8I = T8E * T8H;
					     T8U = T8P * T8T;
					     T8O = T8E * T8N;
					     Rp[WS(rs, 15)] = FNMS(T8S, T8T, T8R);
					     Rp[WS(rs, 7)] = FNMS(T8K, T8N, T8I);
					     Rm[WS(rs, 15)] = FMA(T8S, T8Q, T8U);
					     Rm[WS(rs, 7)] = FMA(T8K, T8H, T8O);
					     T7a = T6Z - T6Y;
					     T70 = T6Y + T6Z;
					}
				   }
				   T6X = T6P + T6W;
				   T7f = T6W - T6P;
				   T7g = T6A - T6J;
				   T6K = T6A + T6J;
				   {
					E T7j, T7b, T7l, T7h;
					T7j = FMA(KP707106781, T7a, T79);
					T7b = FNMS(KP707106781, T7a, T79);
					T7l = FMA(KP707106781, T7g, T7f);
					T7h = FNMS(KP707106781, T7g, T7f);
					{
					     E T7k, T7c, T7m, T7i;
					     T7k = T5X * T7j;
					     T7c = T78 * T7b;
					     T7m = T5X * T7l;
					     T7i = T78 * T7h;
					     Rp[WS(rs, 6)] = FNMS(T5Z, T7l, T7k);
					     Rp[WS(rs, 14)] = FNMS(T7e, T7h, T7c);
					     Rm[WS(rs, 6)] = FMA(T5Z, T7j, T7m);
					     Rm[WS(rs, 14)] = FMA(T7e, T7b, T7i);
					}
				   }
				   {
					E T9h, T96, T9j, T9f;
					T9h = FMA(KP923879532, T95, T8Y);
					T96 = FNMS(KP923879532, T95, T8Y);
					T9j = FMA(KP923879532, T9e, T9b);
					T9f = FNMS(KP923879532, T9e, T9b);
					{
					     E T9k, T9i, T9g, T97;
					     T9k = T3J * T9h;
					     T9i = T3G * T9h;
					     T9g = T98 * T96;
					     T97 = T8V * T96;
					     Rm[WS(rs, 5)] = FMA(T3G, T9j, T9k);
					     Rp[WS(rs, 5)] = FNMS(T3J, T9j, T9i);
					     Rm[WS(rs, 13)] = FMA(T8V, T9f, T9g);
					     Rp[WS(rs, 13)] = FNMS(T98, T9f, T97);
					}
				   }
			      }
			 }
			 {
			      E T31, T3r, T1R, T3m, T33, T32, T3s, T2K, T8z, T8j;
			      {
				   E T73, T6L, T75, T71;
				   T73 = FMA(KP707106781, T6K, T6r);
				   T6L = FNMS(KP707106781, T6K, T6r);
				   T75 = FMA(KP707106781, T70, T6X);
				   T71 = FNMS(KP707106781, T70, T6X);
				   {
					E T76, T74, T72, T6M;
					T76 = T3d * T73;
					T74 = T3a * T73;
					T72 = T6O * T6L;
					T6M = T6i * T6L;
					Rm[WS(rs, 2)] = FMA(T3a, T75, T76);
					Rp[WS(rs, 2)] = FNMS(T3d, T75, T74);
					Rm[WS(rs, 10)] = FMA(T6i, T71, T72);
					Rp[WS(rs, 10)] = FNMS(T6O, T71, T6M);
				   }
			      }
			      {
				   E T7N, T7F, T7P, T7L;
				   T7N = T7D + T7E;
				   T7F = T7D - T7E;
				   T7P = T7K + T7J;
				   T7L = T7J - T7K;
				   {
					E T7O, T7G, T7Q, T7M;
					T7O = T4p * T7N;
					T7G = T7C * T7F;
					T7Q = T4p * T7P;
					T7M = T7C * T7L;
					Rp[WS(rs, 4)] = FNMS(T4r, T7P, T7O);
					Rp[WS(rs, 12)] = FNMS(T7I, T7L, T7G);
					Rm[WS(rs, 4)] = FMA(T4r, T7N, T7Q);
					Rm[WS(rs, 12)] = FMA(T7I, T7F, T7M);
				   }
			      }
			      T31 = FMA(KP923879532, T30, T2X);
			      T3r = FNMS(KP923879532, T30, T2X);
			      T8z = FMA(KP923879532, T8i, T83);
			      T8j = FNMS(KP923879532, T8i, T83);
			      {
				   E T8B, T8x, T8C, T8A;
				   T8B = FMA(KP923879532, T8w, T8t);
				   T8x = FNMS(KP923879532, T8w, T8t);
				   T8C = T1a * T8z;
				   T8A = T16 * T8z;
				   {
					E T8y, T8k, T2i, T2J;
					T8y = T8m * T8j;
					T8k = T7S * T8j;
					Rm[WS(rs, 3)] = FMA(T16, T8B, T8C);
					Rp[WS(rs, 3)] = FNMS(T1a, T8B, T8A);
					Rm[WS(rs, 11)] = FMA(T7S, T8x, T8y);
					Rp[WS(rs, 11)] = FNMS(T8m, T8x, T8k);
					T1R = FMA(KP923879532, T1Q, T1t);
					T3m = FNMS(KP923879532, T1Q, T1t);
					T33 = FNMS(KP668178637, T28, T2h);
					T2i = FMA(KP668178637, T2h, T28);
					T2J = FNMS(KP668178637, T2I, T2z);
					T32 = FMA(KP668178637, T2z, T2I);
					T3s = T2J + T2i;
					T2K = T2i - T2J;
				   }
			      }
			      {
				   E T5l, T53, T5f, T4D, T4K, T4R, T56, T5g;
				   T5l = FNMS(KP923879532, T52, T4Z);
				   T53 = FMA(KP923879532, T52, T4Z);
				   {
					E T3t, T3D, T3f, T2L;
					T3t = FNMS(KP831469612, T3s, T3r);
					T3D = FMA(KP831469612, T3s, T3r);
					T3f = FMA(KP831469612, T2K, T1R);
					T2L = FNMS(KP831469612, T2K, T1R);
					{
					     E T3n, T34, T3g, T2M;
					     T3n = T32 + T33;
					     T34 = T32 - T33;
					     T3g = T3e * T3f;
					     T2M = T1c * T2L;
					     {
						  E T3o, T3z, T3j, T35;
						  T3o = FNMS(KP831469612, T3n, T3m);
						  T3z = FMA(KP831469612, T3n, T3m);
						  T3j = FMA(KP831469612, T34, T31);
						  T35 = FNMS(KP831469612, T34, T31);
						  {
						       E T3u, T3p, T3E, T3A;
						       T3u = T3q * T3o;
						       T3p = T3l * T3o;
						       T3E = T3C * T3z;
						       T3A = T3y * T3z;
						       {
							    E T3k, T36, T54, T55;
							    T3k = T3e * T3j;
							    Ip[WS(rs, 2)] = FNMS(T3i, T3j, T3g);
							    T36 = T1c * T35;
							    Ip[WS(rs, 10)] = FNMS(T2O, T35, T2M);
							    Im[WS(rs, 6)] = FMA(T3l, T3t, T3u);
							    Ip[WS(rs, 6)] = FNMS(T3q, T3t, T3p);
							    Im[WS(rs, 14)] = FMA(T3y, T3D, T3E);
							    Ip[WS(rs, 14)] = FNMS(T3C, T3D, T3A);
							    Im[WS(rs, 2)] = FMA(T3i, T3f, T3k);
							    Im[WS(rs, 10)] = FMA(T2O, T2L, T36);
							    T5f = FMA(KP923879532, T4C, T4v);
							    T4D = FNMS(KP923879532, T4C, T4v);
							    T4K = FNMS(KP668178637, T4J, T4G);
							    T54 = FMA(KP668178637, T4G, T4J);
							    T55 = FMA(KP668178637, T4N, T4Q);
							    T4R = FNMS(KP668178637, T4Q, T4N);
							    T56 = T54 - T55;
							    T5g = T54 + T55;
						       }
						  }
					     }
					}
				   }
				   {
					E T4h, T41, T4c, T3N, T3Q, T3T, T44, T4d;
					T4h = FNMS(KP923879532, T40, T3Z);
					T41 = FMA(KP923879532, T40, T3Z);
					{
					     E T57, T5b, T5h, T5p;
					     T57 = FNMS(KP831469612, T56, T53);
					     T5b = FMA(KP831469612, T56, T53);
					     T5h = FNMS(KP831469612, T5g, T5f);
					     T5p = FMA(KP831469612, T5g, T5f);
					     {
						  E T5m, T4S, T5i, T5q;
						  T5m = T4K - T4R;
						  T4S = T4K + T4R;
						  T5i = T5e * T5h;
						  T5q = T17 * T5p;
						  {
						       E T5n, T5r, T59, T4T;
						       T5n = FMA(KP831469612, T5m, T5l);
						       T5r = FNMS(KP831469612, T5m, T5l);
						       T59 = FMA(KP831469612, T4S, T4D);
						       T4T = FNMS(KP831469612, T4S, T4D);
						       {
							    E T5o, T5s, T5c, T5a;
							    T5o = T5e * T5n;
							    Ip[WS(rs, 5)] = FNMS(T5k, T5n, T5i);
							    T5s = T17 * T5r;
							    Ip[WS(rs, 13)] = FNMS(T1b, T5r, T5q);
							    T5c = T14 * T59;
							    T5a = T11 * T59;
							    {
								 E T58, T4U, T42, T43;
								 T58 = T4W * T4T;
								 T4U = T4s * T4T;
								 Im[WS(rs, 5)] = FMA(T5k, T5h, T5o);
								 Im[WS(rs, 13)] = FMA(T1b, T5p, T5s);
								 Im[WS(rs, 1)] = FMA(T11, T5b, T5c);
								 Ip[WS(rs, 1)] = FNMS(T14, T5b, T5a);
								 Im[WS(rs, 9)] = FMA(T4s, T57, T58);
								 Ip[WS(rs, 9)] = FNMS(T4W, T57, T4U);
								 T4c = FNMS(KP923879532, T3M, T3L);
								 T3N = FMA(KP923879532, T3M, T3L);
								 T3Q = FNMS(KP198912367, T3P, T3O);
								 T42 = FMA(KP198912367, T3O, T3P);
								 T43 = FNMS(KP198912367, T3R, T3S);
								 T3T = FMA(KP198912367, T3S, T3R);
								 T44 = T42 + T43;
								 T4d = T43 - T42;
							    }
						       }
						  }
					     }
					}
					T67 = FNMS(KP923879532, T5K, T5J);
					T5L = FMA(KP923879532, T5K, T5J);
					{
					     E T45, T49, T4e, T4l;
					     T45 = FNMS(KP980785280, T44, T41);
					     T49 = FMA(KP980785280, T44, T41);
					     T4e = FNMS(KP980785280, T4d, T4c);
					     T4l = FMA(KP980785280, T4d, T4c);
					     {
						  E T4i, T3U, T4f, T4m;
						  T4i = T3Q - T3T;
						  T3U = T3Q + T3T;
						  T4f = T4b * T4e;
						  T4m = T12 * T4l;
						  {
						       E T4j, T4n, T47, T3V;
						       T4j = FNMS(KP980785280, T4i, T4h);
						       T4n = FMA(KP980785280, T4i, T4h);
						       T47 = FMA(KP980785280, T3U, T3N);
						       T3V = FNMS(KP980785280, T3U, T3N);
						       {
							    E T4k, T4o, T4a, T48;
							    T4k = T4b * T4j;
							    Ip[WS(rs, 12)] = FNMS(T4g, T4j, T4f);
							    T4o = T12 * T4n;
							    Ip[WS(rs, 4)] = FNMS(T15, T4n, T4m);
							    T4a = T39 * T47;
							    T48 = T37 * T47;
							    {
								 E T46, T3W, T5M, T5N;
								 T46 = T3Y * T3V;
								 T3W = T3K * T3V;
								 Im[WS(rs, 12)] = FMA(T4g, T4e, T4k);
								 Im[WS(rs, 4)] = FMA(T15, T4l, T4o);
								 Im[0] = FMA(T37, T49, T4a);
								 Ip[0] = FNMS(T39, T49, T48);
								 Im[WS(rs, 8)] = FMA(T3K, T45, T46);
								 Ip[WS(rs, 8)] = FNMS(T3Y, T45, T3W);
								 T61 = FMA(KP923879532, T5w, T5v);
								 T5x = FNMS(KP923879532, T5w, T5v);
								 T5A = FNMS(KP198912367, T5z, T5y);
								 T5M = FMA(KP198912367, T5y, T5z);
								 T5N = FMA(KP198912367, T5B, T5C);
								 T5D = FNMS(KP198912367, T5C, T5B);
								 T5O = T5M - T5N;
								 T62 = T5M + T5N;
							    }
						       }
						  }
					     }
					}
				   }
			      }
			 }
		    }
	       }
	  }
	  T5V = FMA(KP980785280, T5O, T5L);
	  T5P = FNMS(KP980785280, T5O, T5L);
	  {
	       E T6c, T63, T5E, T68;
	       T6c = FMA(KP980785280, T62, T61);
	       T63 = FNMS(KP980785280, T62, T61);
	       T5E = T5A + T5D;
	       T68 = T5D - T5A;
	       {