hc2cf2_32.c

快速fft变换

					     {
						  E T2z, T4Z, T5p, T2D, T51, T7v, T5q;
						  T3H = Ip[WS(rs, 1)];
						  {
						       E T2y, T5o, T2B, T50;
						       T2y = FMA(T2w, T2x, T2u);
						       T5o = T2s * T2x;
						       T2B = T8 * T2A;
						       T50 = T8 * T2C;
						       T2z = T2q + T2y;
						       T4Z = T2q - T2y;
						       T5p = FNMS(T2w, T2t, T5o);
						       T2D = FMA(Tc, T2C, T2B);
						       T51 = FNMS(Tc, T2A, T50);
						       T3I = T3 * T3H;
						  }
						  T7v = T5n + T5p;
						  T5q = T5n - T5p;
						  {
						       E T2M, T5r, T7w, T54;
						       T2M = T2D + T2L;
						       T5r = T2D - T2L;
						       T7w = T51 + T53;
						       T54 = T51 - T53;
						       T5s = T5q - T5r;
						       T6I = T5q + T5r;
						       T2N = T2z + T2M;
						       T7A = T2z - T2M;
						       T55 = T4Z + T54;
						       T6F = T4Z - T54;
						       T7x = T7v - T7w;
						       T8i = T7v + T7w;
						       T3J = Im[WS(rs, 1)];
						  }
					     }
					     T3Y = Ip[WS(rs, 5)];
					     T40 = Im[WS(rs, 5)];
					     T3N = Ip[WS(rs, 9)];
					     T3K = FMA(T6, T3J, T3I);
					     T5G = T3 * T3J;
					     T3Z = Td * T3Y;
					     T5P = Td * T40;
					     T3O = T3M * T3N;
					     T5H = FNMS(T6, T3H, T5G);
					     T41 = FMA(Th, T40, T3Z);
					     T5Q = FNMS(Th, T3Y, T5P);
					     T3R = Im[WS(rs, 9)];
					     T3U = Ip[WS(rs, 13)];
					     T3W = Im[WS(rs, 13)];
					}
					{
					     E T2O, T2P, T2Q, T37, T39, T57, T38, T5g;
					     {
						  E T3T, T5F, T5J, T3X, T5O, T7M, T5K;
						  T2O = Ip[WS(rs, 2)];
						  {
						       E T3S, T5I, T3V, T5N;
						       T3S = FMA(T3Q, T3R, T3O);
						       T5I = T3M * T3R;
						       T3V = Te * T3U;
						       T5N = Te * T3W;
						       T3T = T3K + T3S;
						       T5F = T3K - T3S;
						       T5J = FNMS(T3Q, T3N, T5I);
						       T3X = FMA(Ti, T3W, T3V);
						       T5O = FNMS(Ti, T3U, T5N);
						       T2P = T29 * T2O;
						  }
						  T7M = T5H + T5J;
						  T5K = T5H - T5J;
						  {
						       E T42, T5M, T7N, T5R;
						       T42 = T3X + T41;
						       T5M = T3X - T41;
						       T7N = T5O + T5Q;
						       T5R = T5O - T5Q;
						       T5L = T5F + T5K;
						       T62 = T5K - T5F;
						       T43 = T3T + T42;
						       T7G = T42 - T3T;
						       T5S = T5M - T5R;
						       T63 = T5M + T5R;
						       T7O = T7M - T7N;
						       T8o = T7M + T7N;
						       T2Q = Im[WS(rs, 2)];
						  }
					     }
					     T37 = Ip[WS(rs, 6)];
					     T39 = Im[WS(rs, 6)];
					     T2U = Ip[WS(rs, 10)];
					     T2R = FMA(T2b, T2Q, T2P);
					     T57 = T29 * T2Q;
					     T38 = T1R * T37;
					     T5g = T1R * T39;
					     T2V = T2T * T2U;
					     T58 = FNMS(T2b, T2O, T57);
					     T3a = FMA(T1U, T39, T38);
					     T5h = FNMS(T1U, T37, T5g);
					     T2Y = Im[WS(rs, 10)];
					     T32 = Ip[WS(rs, 14)];
					     T35 = Im[WS(rs, 14)];
					}
				   }
			      }
			 }
		    }
		    {
			 E T5c, T5t, T5j, T5u, T88, T90, T8Z, T8b;
			 {
			      E T7e, T8T, T7y, T7D, T7h, T8U, T8S, T8R;
			      {
				   E T8c, T1i, T8A, T8z, T8O, T8J, T8N, T2l, T8L, T45, T8t, T8l, T8u, T8q, T3c;
				   E T8k, T8p, T8w, T2m;
				   {
					E T8x, T8y, T8j, T8C, T8I;
					{
					     E TI, T30, T56, T5a, T36, T5f, T1h, T7B, T5b;
					     TI = Tq + TH;
					     T7e = Tq - TH;
					     {
						  E T2Z, T59, T33, T5e;
						  T2Z = FMA(T2X, T2Y, T2V);
						  T59 = T2T * T2Y;
						  T33 = T31 * T32;
						  T5e = T31 * T35;
						  T30 = T2R + T2Z;
						  T56 = T2R - T2Z;
						  T5a = FNMS(T2X, T2U, T59);
						  T36 = FMA(T34, T35, T33);
						  T5f = FNMS(T34, T32, T5e);
						  T1h = TZ + T1g;
						  T8T = T1g - TZ;
					     }
					     T7B = T58 + T5a;
					     T5b = T58 - T5a;
					     {
						  E T3b, T5d, T7C, T5i;
						  T3b = T36 + T3a;
						  T5d = T36 - T3a;
						  T7C = T5f + T5h;
						  T5i = T5f - T5h;
						  T5c = T56 + T5b;
						  T5t = T5b - T56;
						  T3c = T30 + T3b;
						  T7y = T3b - T30;
						  T5j = T5d - T5i;
						  T5u = T5d + T5i;
						  T7D = T7B - T7C;
						  T8j = T7B + T7C;
						  T8c = TI - T1h;
						  T1i = TI + T1h;
					     }
					}
					T8k = T8i - T8j;
					T8x = T8i + T8j;
					T8y = T8n + T8o;
					T8p = T8n - T8o;
					T7h = T7f - T7g;
					T8C = T7f + T7g;
					T8I = T8D + T8H;
					T8U = T8H - T8D;
					T8A = T8x + T8y;
					T8z = T8x - T8y;
					T8O = T8I - T8C;
					T8J = T8C + T8I;
				   }
				   {
					E T8h, T8m, T3d, T44;
					T8h = T2N - T3c;
					T3d = T2N + T3c;
					T44 = T3G + T43;
					T8m = T3G - T43;
					T8N = T2k - T1J;
					T2l = T1J + T2k;
					T8L = T44 - T3d;
					T45 = T3d + T44;
					T8t = T8k - T8h;
					T8l = T8h + T8k;
					T8u = T8m + T8p;
					T8q = T8m - T8p;
				   }
				   T8w = T1i - T2l;
				   T2m = T1i + T2l;
				   {
					E T8s, T8P, T8Q, T8v;
					{
					     E T8r, T8M, T8K, T8g, T8B, T8f;
					     T8S = T8q - T8l;
					     T8r = T8l + T8q;
					     T8B = T8d + T8e;
					     T8f = T8d - T8e;
					     Rp[0] = T2m + T45;
					     Rm[WS(rs, 15)] = T2m - T45;
					     Rp[WS(rs, 8)] = T8w + T8z;
					     Rm[WS(rs, 7)] = T8w - T8z;
					     T8M = T8J - T8B;
					     T8K = T8B + T8J;
					     T8g = T8c + T8f;
					     T8s = T8c - T8f;
					     T8R = T8O - T8N;
					     T8P = T8N + T8O;
					     Ip[WS(rs, 8)] = T8L + T8M;
					     Im[WS(rs, 7)] = T8L - T8M;
					     Ip[0] = T8A + T8K;
					     Im[WS(rs, 15)] = T8A - T8K;
					     Rp[WS(rs, 4)] = FMA(KP707106781, T8r, T8g);
					     Rm[WS(rs, 11)] = FNMS(KP707106781, T8r, T8g);
					     T8Q = T8t + T8u;
					     T8v = T8t - T8u;
					}
					Ip[WS(rs, 4)] = FMA(KP707106781, T8Q, T8P);
					Im[WS(rs, 11)] = FMS(KP707106781, T8Q, T8P);
					Rp[WS(rs, 12)] = FMA(KP707106781, T8v, T8s);
					Rm[WS(rs, 3)] = FNMS(KP707106781, T8v, T8s);
				   }
			      }
			      {
				   E T7P, T7W, T7i, T7K, T8a, T86, T91, T8V, T8W, T7t, T7T, T7F, T92, T7Z, T89;
				   E T83;
				   {
					E T7X, T7n, T7s, T7Y, T84, T85;
					T7P = T7L - T7O;
					T84 = T7L + T7O;
					Ip[WS(rs, 12)] = FMA(KP707106781, T8S, T8R);
					Im[WS(rs, 3)] = FMS(KP707106781, T8S, T8R);
					T7W = T7e + T7h;
					T7i = T7e - T7h;
					T85 = T7J + T7G;
					T7K = T7G - T7J;
					T7X = T7m + T7l;
					T7n = T7l - T7m;
					T8a = FMA(KP414213562, T84, T85);
					T86 = FNMS(KP414213562, T85, T84);
					T91 = T8U - T8T;
					T8V = T8T + T8U;
					T7s = T7o + T7r;
					T7Y = T7o - T7r;
					{
					     E T82, T81, T7z, T7E;
					     T82 = T7x + T7y;
					     T7z = T7x - T7y;
					     T7E = T7A - T7D;
					     T81 = T7A + T7D;
					     T8W = T7n + T7s;
					     T7t = T7n - T7s;
					     T7T = FNMS(KP414213562, T7z, T7E);
					     T7F = FMA(KP414213562, T7E, T7z);
					     T92 = T7Y - T7X;
					     T7Z = T7X + T7Y;
					     T89 = FNMS(KP414213562, T81, T82);
					     T83 = FMA(KP414213562, T82, T81);
					}
				   }
				   {
					E T7S, T7u, T93, T95, T7U, T7Q;
					T7S = FNMS(KP707106781, T7t, T7i);
					T7u = FMA(KP707106781, T7t, T7i);
					T93 = FMA(KP707106781, T92, T91);
					T95 = FNMS(KP707106781, T92, T91);
					T7U = FNMS(KP414213562, T7K, T7P);
					T7Q = FMA(KP414213562, T7P, T7K);
					{
					     E T80, T87, T8X, T8Y;
					     T88 = FNMS(KP707106781, T7Z, T7W);
					     T80 = FMA(KP707106781, T7Z, T7W);
					     {
						  E T7V, T94, T96, T7R;
						  T7V = T7T + T7U;
						  T94 = T7U - T7T;
						  T96 = T7Q - T7F;
						  T7R = T7F + T7Q;
						  Rm[WS(rs, 1)] = FMA(KP923879532, T7V, T7S);
						  Rp[WS(rs, 14)] = FNMS(KP923879532, T7V, T7S);
						  Ip[WS(rs, 6)] = FMA(KP923879532, T94, T93);
						  Im[WS(rs, 9)] = FMS(KP923879532, T94, T93);
						  Ip[WS(rs, 14)] = FMA(KP923879532, T96, T95);
						  Im[WS(rs, 1)] = FMS(KP923879532, T96, T95);
						  Rp[WS(rs, 6)] = FMA(KP923879532, T7R, T7u);
						  Rm[WS(rs, 9)] = FNMS(KP923879532, T7R, T7u);
						  T87 = T83 + T86;
						  T90 = T86 - T83;
					     }
					     T8Z = FNMS(KP707106781, T8W, T8V);
					     T8X = FMA(KP707106781, T8W, T8V);
					     T8Y = T89 + T8a;
					     T8b = T89 - T8a;
					     Rp[WS(rs, 2)] = FMA(KP923879532, T87, T80);
					     Rm[WS(rs, 13)] = FNMS(KP923879532, T87, T80);
					     Ip[WS(rs, 2)] = FMA(KP923879532, T8Y, T8X);
					     Im[WS(rs, 13)] = FMS(KP923879532, T8Y, T8X);
					}
				   }
			      }
			 }
			 {
			      E T6s, T9o, T9n, T6v, T6Q, T6N, T6J, T6G, T9k, T9j;
			      {
				   E T6c, T4s, T9i, T4X, T9h, T9b, T9c, T6f, T5U, T6k, T64, T5k, T5v;
				   {
					E T6d, T6e, T99, T9a, T5T;
					{
					     E T4c, T4r, T4H, T4W;
					     T6s = T46 - T4b;
					     T4c = T46 + T4b;
					     Rp[WS(rs, 10)] = FMA(KP923879532, T8b, T88);
					     Rm[WS(rs, 5)] = FNMS(KP923879532, T8b, T88);
					     Ip[WS(rs, 10)] = FMA(KP923879532, T90, T8Z);
					     Im[WS(rs, 5)] = FMS(KP923879532, T90, T8Z);
					     T4r = T4j + T4q;
					     T9o = T4q - T4j;
					     T6d = FNMS(KP414213562, T4z, T4G);
					     T4H = FMA(KP414213562, T4G, T4z);
					     T4W = FNMS(KP414213562, T4V, T4O);
					     T6e = FMA(KP414213562, T4O, T4V);
					     T9n = T98 + T97;
					     T99 = T97 - T98;
					     T6c = FNMS(KP707106781, T4r, T4c);
					     T4s = FMA(KP707106781, T4r, T4c);
					     T9i = T4W - T4H;
					     T4X = T4H + T4W;
					     T9a = T6t + T6u;
					     T6v = T6t - T6u;
					}
					T6Q = T5S - T5L;
					T5T = T5L + T5S;
					T9h = FNMS(KP707106781, T9a, T99);
					T9b = FMA(KP707106781, T9a, T99);
					T9c = T6d + T6e;
					T6f = T6d - T6e;
					T5U = FMA(KP707106781, T5T, T5E);
					T6k = FNMS(KP707106781, T5T, T5E);
					T64 = T62 + T63;
					T6N = T63 - T62;
					T6J = T5c - T5j;
					T5k = T5c + T5j;
					T5v = T5t + T5u;
					T6G = T5u - T5t;
				   }
				   {
					E T6m, T6q, T6j, T6p, T9f, T9g;
					{
					     E T68, T4Y, T6a, T66, T69, T5x, T9d, T6l, T65, T9e, T6b, T67;
					     T68 = FNMS(KP923879532, T4X, T4s);
					     T4Y = FMA(KP923879532, T4X, T4s);
					     T6l = FNMS(KP707106781, T64, T61);
					     T65 = FMA(KP707106781, T64, T61);
					     {
						  E T6h, T5l, T6i, T5w;
						  T6h = FNMS(KP707106781, T5k, T55);
						  T5l = FMA(KP707106781, T5k, T55);
						  T6i = FNMS(KP707106781, T5v, T5s);
						  T5w = FMA(KP707106781, T5v, T5s);
						  T6m = FMA(KP668178637, T6l, T6k);
						  T6q = FNMS(KP668178637, T6k, T6l);
						  T6a = FMA(KP198912367, T5U, T65);
						  T66 = FNMS(KP198912367, T65, T5U);
						  T6j = FNMS(KP668178637, T6i, T6h);
						  T6p = FMA(KP668178637, T6h, T6i);
						  T69 = FNMS(KP198912367, T5l, T5w);
						  T5x = FMA(KP198912367, T5w, T5l);
					     }
					     T9d = FMA(KP923879532, T9c, T9b);
					     T9f = FNMS(KP923879532, T9c, T9b);
					     T9e = T69 + T6a;
					     T6b = T69 - T6a;
					     T9g = T66 - T5x;
					     T67 = T5x + T66;
					     Ip[WS(rs, 1)] = FMA(KP980785280, T9e, T9d);
					     Im[WS(rs, 14)] = FMS(KP980785280, T9e, T9d);
					     Rp[WS(rs, 1)] = FMA(KP980785280, T67, T4Y);
					     Rm[WS(rs, 14)] = FNMS(KP980785280, T67, T4Y);
					     Rp[WS(rs, 9)] = FMA(KP980785280, T6b, T68);
					     Rm[WS(rs, 6)] = FNMS(KP980785280, T6b, T68);
					}
					{
					     E T6o, T9l, T9m, T6r, T6g, T6n;
					     T6o = FMA(KP923879532, T6f, T6c);
					     T6g = FNMS(KP923879532, T6f, T6c);
					     T6n = T6j + T6m;
					     T9k = T6m - T6j;
					     T9j = FMA(KP923879532, T9i, T9h);
					     T9l = FNMS(KP923879532, T9i, T9h);
					     Ip[WS(rs, 9)] = FMA(KP980785280, T9g, T9f);
					     Im[WS(rs, 6)] = FMS(KP980785280, T9g, T9f);
					     Rm[WS(rs, 2)] = FMA(KP831469612, T6n, T6g);
					     Rp[WS(rs, 13)] = FNMS(KP831469612, T6n, T6g);
					     T9m = T6p + T6q;
					     T6r = T6p - T6q;
					     Ip[WS(rs, 13)] = FNMS(KP831469612, T9m, T9l);
					     Im[WS(rs, 2)] = -(FMA(KP831469612, T9m, T9l));
					     Rp[WS(rs, 5)] = FMA(KP831469612, T6r, T6o);
					     Rm[WS(rs, 10)] = FNMS(KP831469612, T6r, T6o);
					}
				   }
			      }
			      {
				   E T6Y, T6w, T9w, T6D, T9v, T9p, T9q, T71, T6H, T74, T78, T7c, T6W, T6S;
				   {
					E T6Z, T6z, T6C, T70;
					T6Z = FNMS(KP414213562, T6x, T6y);
					T6z = FMA(KP414213562, T6y, T6x);
					Ip[WS(rs, 5)] = FMA(KP831469612, T9k, T9j);
					Im[WS(rs, 10)] = FMS(KP831469612, T9k, T9j);
					T6Y = FNMS(KP707106781, T6v, T6s);
					T6w = FMA(KP707106781, T6v, T6s);
					T6C = FNMS(KP414213562, T6B, T6A);
					T70 = FMA(KP414213562, T6A, T6B);
					T9w = T6z + T6C;
					T6D = T6z - T6C;
					T9v = FNMS(KP707106781, T9o, T9n);
					T9p = FMA(KP707106781, T9o, T9n);
					{
					     E T77, T6O, T76, T6R;
					     T9q = T70 - T6Z;
					     T71 = T6Z + T70;
					     T77 = FMA(KP707106781, T6N, T6M);
					     T6O = FNMS(KP707106781, T6N, T6M);
					     T76 = FMA(KP707106781, T6Q, T6P);
					     T6R = FNMS(KP707106781, T6Q, T6P);
					     T6H = FNMS(KP707106781, T6G, T6F);
					     T74 = FMA(KP707106781, T6G, T6F);
					     T78 = FMA(KP198912367, T77, T76);
					     T7c = FNMS(KP198912367, T76, T77);
					     T6W = FNMS(KP668178637, T6O, T6R);
					     T6S = FMA(KP668178637, T6R, T6O);
					}
				   }
				   {
					E T6U, T6E, T9r, T9t, T73, T6K;
					T6U = FNMS(KP923879532, T6D, T6w);
					T6E = FMA(KP923879532, T6D, T6w);
					T9r = FMA(KP923879532, T9q, T9p);
					T9t = FNMS(KP923879532, T9q, T9p);
					T73 = FMA(KP707106781, T6J, T6I);
					T6K = FNMS(KP707106781, T6J, T6I);
					{
					     E T7a, T9x, T9y, T7d;
					     {
						  E T72, T7b, T6V, T6L, T79, T75;
						  T7a = FMA(KP923879532, T71, T6Y);
						  T72 = FNMS(KP923879532, T71, T6Y);
						  T75 = FMA(KP198912367, T74, T73);
						  T7b = FNMS(KP198912367, T73, T74);
						  T6V = FNMS(KP668178637, T6H, T6K);
						  T6L = FMA(KP668178637, T6K, T6H);
						  T79 = T75 + T78;
						  T9A = T78 - T75;
						  T9z = FMA(KP923879532, T9w, T9v);
						  T9x = FNMS(KP923879532, T9w, T9v);
						  {
						       E T6X, T9s, T9u, T6T;
						       T6X = T6V + T6W;
						       T9s = T6V - T6W;
						       T9u = T6S - T6L;
						       T6T = T6L + T6S;
						       Rp[WS(rs, 7)] = FMA(KP980785280, T79, T72);
						       Rm[WS(rs, 8)] = FNMS(KP980785280, T79, T72);
						       Rp[WS(rs, 11)] = FMA(KP831469612, T6X, T6U);
						       Rm[WS(rs, 4)] = FNMS(KP831469612, T6X, T6U);
						       Ip[WS(rs, 3)] = FMA(KP831469612, T9s, T9r);