hc2cfdft_32.c

快速fft变换

				   T1U = FNMS(T1Q, T1T, T1P);
				   T4O = FMA(T1Q, T1O, T4N);
				   T20 = W[23];
				   T31 = FNMS(T2Z, T30, T2Y);
				   T37 = FNMS(T35, T36, T34);
				   {
					E T5N, T4Q, T24, T5P;
					T5N = FMA(T2W, T30, T5M);
					T4Q = FMA(T20, T1Y, T4P);
					T24 = FNMS(T20, T23, T1Z);
					T8f = T31 - T37;
					T38 = T31 + T37;
					T5P = FMA(T32, T36, T5O);
					T7v = T4O - T4Q;
					T4R = T4O + T4Q;
					T7w = T1U - T24;
					T25 = T1U + T24;
					T8g = T5N - T5P;
					T5Q = T5N + T5P;
				   }
			      }
			 }
			 {
			      E Tp, Tu, Tq, T4Z, Tw, T3X, Tm, T3R, Tr;
			      {
				   E Tn, To, Ts, Tt;
				   Tn = Ip[WS(rs, 5)];
				   T8h = T8f + T8g;
				   T8o = T8g - T8f;
				   To = Im[WS(rs, 5)];
				   Ts = Rp[WS(rs, 5)];
				   Tt = Rm[WS(rs, 5)];
				   Tm = W[18];
				   T3V = Tn + To;
				   Tp = Tn - To;
				   T3S = Ts - Tt;
				   Tu = Ts + Tt;
				   T3R = W[20];
				   Tq = Tm * Tp;
			      }
			      {
				   E Tx, Ty, TC, TD;
				   Tx = Ip[WS(rs, 13)];
				   T4Z = Tm * Tu;
				   T5p = T3R * T3V;
				   T3T = T3R * T3S;
				   Ty = Im[WS(rs, 13)];
				   TC = Rp[WS(rs, 13)];
				   TD = Rm[WS(rs, 13)];
				   Tw = W[50];
				   T41 = Tx + Ty;
				   Tz = Tx - Ty;
				   T3Y = TC - TD;
				   TE = TC + TD;
				   T3X = W[52];
				   TA = Tw * Tz;
			      }
			      Tr = W[19];
			      T51 = Tw * TE;
			      T5r = T3X * T41;
			      T3Z = T3X * T3Y;
			      Tv = FNMS(Tr, Tu, Tq);
			      T50 = FMA(Tr, Tp, T4Z);
			      TB = W[51];
			      T3U = W[21];
			      T40 = W[53];
			 }
		    }
	       }
	       {
		    E T6y, T7B, T7E, T6u, T6S, T85, T7Y, T6s, T6v, T6x, T6R, T6r, T6F, T6D, T6C;
		    E T6G, T6M, T6K, T6J, T6N, T6l, T6o, T7j, T7m;
		    {
			 E T6i, T1p, T68, T2H, T67, T5g, T6h, T4T, T4w, T5G, T6d, T3D, T6c, T6m, T63;
			 E T6e;
			 {
			      E T5t, T43, T26, T2G, T54, T5f, T4H, T4S;
			      {
				   E T1o, T53, T7W, T7X, TH, T52, TF, T5q;
				   T6y = T12 - T1n;
				   T1o = T12 + T1n;
				   T52 = FMA(TB, Tz, T51);
				   TF = FNMS(TB, TE, TA);
				   T5q = FNMS(T3U, T3S, T5p);
				   {
					E T3W, T5s, T42, TG;
					T3W = FMA(T3U, T3V, T3T);
					T5s = FNMS(T40, T3Y, T5r);
					T42 = FMA(T40, T41, T3Z);
					T7B = T50 - T52;
					T53 = T50 + T52;
					T7E = Tv - TF;
					TG = Tv + TF;
					T7W = T5s - T5q;
					T5t = T5q + T5s;
					T7X = T3W - T42;
					T43 = T3W + T42;
					TH = Tl + TG;
					T6u = Tl - TG;
				   }
				   T6S = T1K - T25;
				   T26 = T1K + T25;
				   T85 = T7W - T7X;
				   T7Y = T7W + T7X;
				   T6i = TH - T1o;
				   T1p = TH + T1o;
				   T2G = T2r + T2F;
				   T6s = T2F - T2r;
				   T6v = T4Y - T53;
				   T54 = T4Y + T53;
				   T5f = T59 + T5e;
				   T6x = T59 - T5e;
			      }
			      T6R = T4B - T4G;
			      T4H = T4B + T4G;
			      T68 = T2G - T26;
			      T2H = T26 + T2G;
			      T67 = T5f - T54;
			      T5g = T54 + T5f;
			      T4S = T4M + T4R;
			      T6r = T4R - T4M;
			      {
				   E T5u, T6b, T5F, T44, T4v;
				   T6F = T43 - T3Q;
				   T44 = T3Q + T43;
				   T4v = T4h + T4u;
				   T6D = T4u - T4h;
				   T6C = T5t - T5o;
				   T5u = T5o + T5t;
				   T6h = T4H - T4S;
				   T4T = T4H + T4S;
				   T6b = T44 - T4v;
				   T4w = T44 + T4v;
				   T6G = T5z - T5E;
				   T5F = T5z + T5E;
				   {
					E T5R, T62, T39, T3C, T6a;
					T6M = T2V - T38;
					T39 = T2V + T38;
					T3C = T3m + T3B;
					T6K = T3B - T3m;
					T6a = T5F - T5u;
					T5G = T5u + T5F;
					T6J = T5Q - T5L;
					T5R = T5L + T5Q;
					T6d = T3C - T39;
					T3D = T39 + T3C;
					T6N = T61 - T5W;
					T62 = T5W + T61;
					T6c = T6a + T6b;
					T6m = T6a - T6b;
					T63 = T5R + T62;
					T6e = T62 - T5R;
				   }
			      }
			 }
			 {
			      E T5j, T6n, T6f, T64;
			      {
				   E T5i, T5h, T65, T66, T2I, T4x;
				   T5j = T2H - T1p;
				   T2I = T1p + T2H;
				   T4x = T3D - T4w;
				   T5i = T4w + T3D;
				   T6n = T6d + T6e;
				   T6f = T6d - T6e;
				   T5h = T4T - T5g;
				   T65 = T4T + T5g;
				   Im[WS(rs, 15)] = KP500000000 * (T4x - T2I);
				   Ip[0] = KP500000000 * (T2I + T4x);
				   T66 = T5G + T63;
				   T64 = T5G - T63;
				   Rp[0] = KP500000000 * (T65 + T66);
				   Rm[WS(rs, 15)] = KP500000000 * (T65 - T66);
				   Rp[WS(rs, 8)] = KP500000000 * (T5h + T5i);
				   Rm[WS(rs, 7)] = KP500000000 * (T5h - T5i);
			      }
			      {
				   E T6k, T6j, T6p, T6q, T69, T6g;
				   T6l = T68 - T67;
				   T69 = T67 + T68;
				   T6g = T6c + T6f;
				   T6k = T6f - T6c;
				   T6j = T6h - T6i;
				   T6p = T6h + T6i;
				   Im[WS(rs, 7)] = KP500000000 * (T64 - T5j);
				   Ip[WS(rs, 8)] = KP500000000 * (T5j + T64);
				   Im[WS(rs, 11)] = -(KP500000000 * (FNMS(KP707106781, T6g, T69)));
				   Ip[WS(rs, 4)] = KP500000000 * (FMA(KP707106781, T6g, T69));
				   T6q = T6m + T6n;
				   T6o = T6m - T6n;
				   Rp[WS(rs, 4)] = KP500000000 * (FMA(KP707106781, T6q, T6p));
				   Rm[WS(rs, 11)] = KP500000000 * (FNMS(KP707106781, T6q, T6p));
				   Rp[WS(rs, 12)] = KP500000000 * (FMA(KP707106781, T6k, T6j));
				   Rm[WS(rs, 3)] = KP500000000 * (FNMS(KP707106781, T6k, T6j));
			      }
			 }
		    }
		    {
			 E T75, T6t, T7f, T6T, T76, T6W, T7g, T6A, T7b, T6L, T7a, T7k, T70, T6I, T6U;
			 E T6w;
			 Im[WS(rs, 3)] = -(KP500000000 * (FNMS(KP707106781, T6o, T6l)));
			 Ip[WS(rs, 12)] = KP500000000 * (FMA(KP707106781, T6o, T6l));
			 T75 = T6s - T6r;
			 T6t = T6r + T6s;
			 T7f = T6R - T6S;
			 T6T = T6R + T6S;
			 T6U = T6v + T6u;
			 T6w = T6u - T6v;
			 {
			      E T78, T6E, T6V, T6z, T79, T6H;
			      T6V = T6x - T6y;
			      T6z = T6x + T6y;
			      T78 = T6C - T6D;
			      T6E = T6C + T6D;
			      T76 = T6V - T6U;
			      T6W = T6U + T6V;
			      T7g = T6w - T6z;
			      T6A = T6w + T6z;
			      T79 = T6G - T6F;
			      T6H = T6F + T6G;
			      T7b = T6K - T6J;
			      T6L = T6J + T6K;
			      T7a = FMA(KP414213562, T79, T78);
			      T7k = FNMS(KP414213562, T78, T79);
			      T70 = FNMS(KP414213562, T6E, T6H);
			      T6I = FMA(KP414213562, T6H, T6E);
			 }
			 {
			      E T6Z, T6B, T73, T6X, T7c, T6O;
			      T6Z = FNMS(KP707106781, T6A, T6t);
			      T6B = FMA(KP707106781, T6A, T6t);
			      T73 = FMA(KP707106781, T6W, T6T);
			      T6X = FNMS(KP707106781, T6W, T6T);
			      T7c = T6N - T6M;
			      T6O = T6M + T6N;
			      {
				   E T7i, T7h, T7n, T7o;
				   {
					E T77, T7l, T71, T6P, T7e, T7d;
					T7j = FMA(KP707106781, T76, T75);
					T77 = FNMS(KP707106781, T76, T75);
					T7d = FMA(KP414213562, T7c, T7b);
					T7l = FNMS(KP414213562, T7b, T7c);
					T71 = FMA(KP414213562, T6L, T6O);
					T6P = FNMS(KP414213562, T6O, T6L);
					T7e = T7a - T7d;
					T7i = T7a + T7d;
					T7h = FMA(KP707106781, T7g, T7f);
					T7n = FNMS(KP707106781, T7g, T7f);
					{
					     E T72, T74, T6Y, T6Q;
					     T72 = T70 - T71;
					     T74 = T70 + T71;
					     T6Y = T6P - T6I;
					     T6Q = T6I + T6P;
					     Im[WS(rs, 1)] = -(KP500000000 * (FNMS(KP923879532, T7e, T77)));
					     Ip[WS(rs, 14)] = KP500000000 * (FMA(KP923879532, T7e, T77));
					     Im[WS(rs, 5)] = -(KP500000000 * (FNMS(KP923879532, T72, T6Z)));
					     Ip[WS(rs, 10)] = KP500000000 * (FMA(KP923879532, T72, T6Z));
					     Rp[WS(rs, 2)] = KP500000000 * (FMA(KP923879532, T74, T73));
					     Rm[WS(rs, 13)] = KP500000000 * (FNMS(KP923879532, T74, T73));
					     Rp[WS(rs, 10)] = KP500000000 * (FMA(KP923879532, T6Y, T6X));
					     Rm[WS(rs, 5)] = KP500000000 * (FNMS(KP923879532, T6Y, T6X));
					     Im[WS(rs, 13)] = -(KP500000000 * (FNMS(KP923879532, T6Q, T6B)));
					     Ip[WS(rs, 2)] = KP500000000 * (FMA(KP923879532, T6Q, T6B));
					     T7o = T7k + T7l;
					     T7m = T7k - T7l;
					}
				   }
				   Rm[WS(rs, 1)] = KP500000000 * (FMA(KP923879532, T7o, T7n));
				   Rp[WS(rs, 14)] = KP500000000 * (FNMS(KP923879532, T7o, T7n));
				   Rp[WS(rs, 6)] = KP500000000 * (FMA(KP923879532, T7i, T7h));
				   Rm[WS(rs, 9)] = KP500000000 * (FNMS(KP923879532, T7i, T7h));
			      }
			 }
		    }
		    {
			 E T9x, T9T, T8L, T7z, T97, T9J, T8V, T8z, T8M, T8C, T8W, T7O, T9O, T9Y, T9E;
			 E T9t, T8Q, T90, T8G, T88, T8p, T8m, T9K, T9A, T9U, T9e, T8R, T8j, T9R, T9Z;
			 E T9F, T9m;
			 {
			      E T9c, T9b, T99, T98, T7S, T86, T83, T9q, T9M, T9p, T9r, T7Z, T9z, T9a;
			      {
				   E T95, T7r, T9v, T8v, T8w, T8x, T9w, T7y, T7u, T7x;
				   T95 = T7q + T7p;
				   T7r = T7p - T7q;
				   T9v = T8t - T8u;
				   T8v = T8t + T8u;
				   T8w = T7t + T7s;
				   T7u = T7s - T7t;
				   Im[WS(rs, 9)] = -(KP500000000 * (FNMS(KP923879532, T7m, T7j)));
				   Ip[WS(rs, 6)] = KP500000000 * (FMA(KP923879532, T7m, T7j));
				   T7x = T7v + T7w;
				   T8x = T7v - T7w;
				   T9w = T7u - T7x;
				   T7y = T7u + T7x;
				   {
					E T7J, T8A, T7G, T7M;
					{
					     E T7C, T96, T8y, T7F;
					     T9c = T7A + T7B;
					     T7C = T7A - T7B;
					     T9x = FMA(KP707106781, T9w, T9v);
					     T9T = FNMS(KP707106781, T9w, T9v);
					     T8L = FNMS(KP707106781, T7y, T7r);
					     T7z = FMA(KP707106781, T7y, T7r);
					     T96 = T8x - T8w;
					     T8y = T8w + T8x;
					     T7F = T7D + T7E;
					     T9b = T7D - T7E;
					     T99 = T7H + T7I;
					     T7J = T7H - T7I;
					     T97 = FMA(KP707106781, T96, T95);
					     T9J = FNMS(KP707106781, T96, T95);
					     T8V = FNMS(KP707106781, T8y, T8v);
					     T8z = FMA(KP707106781, T8y, T8v);
					     T8A = FMA(KP414213562, T7C, T7F);
					     T7G = FNMS(KP414213562, T7F, T7C);
					     T7M = T7K + T7L;
					     T98 = T7K - T7L;
					}
					{
					     E T9n, T9o, T8B, T7N;
					     T7S = T7Q + T7R;
					     T9n = T7R - T7Q;
					     T9o = T85 - T84;
					     T86 = T84 + T85;
					     T83 = T81 + T82;
					     T9q = T81 - T82;
					     T8B = FNMS(KP414213562, T7J, T7M);
					     T7N = FMA(KP414213562, T7M, T7J);
					     T9M = FMA(KP707106781, T9o, T9n);
					     T9p = FNMS(KP707106781, T9o, T9n);
					     T8M = T8B - T8A;
					     T8C = T8A + T8B;
					     T8W = T7G - T7N;
					     T7O = T7G + T7N;
					     T9r = T7Y - T7V;
					     T7Z = T7V + T7Y;
					}
				   }
			      }
			      {
				   E T8O, T80, T9N, T9s, T8P, T87;
				   T9N = FMA(KP707106781, T9r, T9q);
				   T9s = FNMS(KP707106781, T9r, T9q);
				   T8O = FNMS(KP707106781, T7Z, T7S);
				   T80 = FMA(KP707106781, T7Z, T7S);
				   T9O = FMA(KP198912367, T9N, T9M);
				   T9Y = FNMS(KP198912367, T9M, T9N);
				   T9E = FMA(KP668178637, T9p, T9s);
				   T9t = FNMS(KP668178637, T9s, T9p);
				   T8P = FNMS(KP707106781, T86, T83);
				   T87 = FMA(KP707106781, T86, T83);
				   T9z = FNMS(KP414213562, T98, T99);
				   T9a = FMA(KP414213562, T99, T98);
				   T8Q = FNMS(KP668178637, T8P, T8O);
				   T90 = FMA(KP668178637, T8O, T8P);
				   T8G = FNMS(KP198912367, T80, T87);
				   T88 = FMA(KP198912367, T87, T80);
			      }
			      {
				   E T8b, T9j, T9P, T9i, T9k, T8i, T9Q, T9l;
				   {
					E T9g, T9h, T9y, T9d;
					T8b = T89 - T8a;
					T9g = T8a + T89;
					T9h = T8n - T8o;
					T8p = T8n + T8o;
					T8m = T8k + T8l;
					T9j = T8l - T8k;
					T9y = FMA(KP414213562, T9b, T9c);
					T9d = FNMS(KP414213562, T9c, T9b);
					T9P = FMA(KP707106781, T9h, T9g);
					T9i = FNMS(KP707106781, T9h, T9g);
					T9K = T9y + T9z;
					T9A = T9y - T9z;
					T9U = T9d + T9a;
					T9e = T9a - T9d;
					T9k = T8h - T8e;
					T8i = T8e + T8h;
				   }
				   T9Q = FMA(KP707106781, T9k, T9j);
				   T9l = FNMS(KP707106781, T9k, T9j);
				   T8R = FNMS(KP707106781, T8i, T8b);
				   T8j = FMA(KP707106781, T8i, T8b);
				   T9R = FMA(KP198912367, T9Q, T9P);
				   T9Z = FNMS(KP198912367, T9P, T9Q);
				   T9F = FMA(KP668178637, T9i, T9l);
				   T9m = FNMS(KP668178637, T9l, T9i);
			      }
			 }
			 {
			      E T8Z, T92, T9D, T9G;
			      {
				   E T8F, T7P, T8J, T8D, T8S, T8q;
				   T8F = FNMS(KP923879532, T7O, T7z);
				   T7P = FMA(KP923879532, T7O, T7z);
				   T8J = FMA(KP923879532, T8C, T8z);
				   T8D = FNMS(KP923879532, T8C, T8z);
				   T8S = FNMS(KP707106781, T8p, T8m);
				   T8q = FMA(KP707106781, T8p, T8m);
				   {
					E T8Y, T8X, T93, T94;
					{
					     E T8N, T91, T8H, T8r, T8U, T8T;
					     T8Z = FMA(KP923879532, T8M, T8L);
					     T8N = FNMS(KP923879532, T8M, T8L);
					     T8T = FMA(KP668178637, T8S, T8R);
					     T91 = FNMS(KP668178637, T8R, T8S);
					     T8H = FMA(KP198912367, T8j, T8q);
					     T8r = FNMS(KP198912367, T8q, T8j);
					     T8U = T8Q + T8T;
					     T8Y = T8T - T8Q;
					     T8X = FMA(KP923879532, T8W, T8V);
					     T93 = FNMS(KP923879532, T8W, T8V);
					     {
						  E T8I, T8K, T8E, T8s;
						  T8I = T8G - T8H;
						  T8K = T8G + T8H;
						  T8E = T8r - T88;
						  T8s = T88 + T8r;
						  Im[WS(rs, 2)] = -(KP500000000 * (FMA(KP831469612, T8U, T8N)));
						  Ip[WS(rs, 13)] = KP500000000 * (FNMS(KP831469612, T8U, T8N));
						  Im[WS(rs, 6)] = -(KP500000000 * (FNMS(KP980785280, T8I, T8F)));
						  Ip[WS(rs, 9)] = KP500000000 * (FMA(KP980785280, T8I, T8F));
						  Rp[WS(rs, 1)] = KP500000000 * (FMA(KP980785280, T8K, T8J));
						  Rm[WS(rs, 14)] = KP500000000 * (FNMS(KP980785280, T8K, T8J));
						  Rp[WS(rs, 9)] = KP500000000 * (FMA(KP980785280, T8E, T8D));
						  Rm[WS(rs, 6)] = KP500000000 * (FNMS(KP980785280, T8E, T8D));
						  Im[WS(rs, 14)] = -(KP500000000 * (FNMS(KP980785280, T8s, T7P)));
						  Ip[WS(rs, 1)] = KP500000000 * (FMA(KP980785280, T8s, T7P));
						  T94 = T90 + T91;
						  T92 = T90 - T91;
					     }
					}
					Rm[WS(rs, 2)] = KP500000000 * (FMA(KP831469612, T94, T93));
					Rp[WS(rs, 13)] = KP500000000 * (FNMS(KP831469612, T94, T93));
					Rp[WS(rs, 5)] = KP500000000 * (FMA(KP831469612, T8Y, T8X));
					Rm[WS(rs, 10)] = KP500000000 * (FNMS(KP831469612, T8Y, T8X));
				   }
			      }
			      {
				   E T9C, T9B, T9H, T9I, T9f, T9u;
				   T9D = FNMS(KP923879532, T9e, T97);
				   T9f = FMA(KP923879532, T9e, T97);
				   T9u = T9m - T9t;
				   T9C = T9t + T9m;
				   T9B = FNMS(KP923879532, T9A, T9x);
				   T9H = FMA(KP923879532, T9A, T9x);
				   Im[WS(rs, 10)] = -(KP500000000 * (FNMS(KP831469612, T92, T8Z)));
				   Ip[WS(rs, 5)] = KP500000000 * (FMA(KP831469612, T92, T8Z));
				   Im[WS(rs, 12)] = -(KP500000000 * (FNMS(KP831469612, T9u, T9f)));
				   Ip[WS(rs, 3)] = KP500000000 * (FMA(KP831469612, T9u, T9f));
				   T9I = T9E + T9F;
				   T9G = T9E - T9F;
				   Rp[WS(rs, 3)] = KP500000000 * (FMA(KP831469612, T9I, T9H));
				   Rm[WS(rs, 12)] = KP500000000 * (FNMS(KP831469612, T9I, T9H));
				   Rp[WS(rs, 11)] = KP500000000 * (FMA(KP831469612, T9C, T9B));
				   Rm[WS(rs, 4)] = KP500000000 * (FNMS(KP831469612, T9C, T9B));
			      }
			      {
				   E T9W, T9V, Ta1, Ta2, T9L, T9S;
				   T9X = FNMS(KP923879532, T9K, T9J);
				   T9L = FMA(KP923879532, T9K, T9J);
				   T9S = T9O - T9R;
				   T9W = T9O + T9R;
				   T9V = FNMS(KP923879532, T9U, T9T);
				   Ta1 = FMA(KP923879532, T9U, T9T);
				   Im[WS(rs, 4)] = -(KP500000000 * (FNMS(KP831469612, T9G, T9D)));
				   Ip[WS(rs, 11)] = KP500000000 * (FMA(KP831469612, T9G, T9D));
				   Im[0] = -(KP500000000 * (FNMS(KP980785280, T9S, T9L)));
				   Ip[WS(rs, 15)] = KP500000000 * (FMA(KP980785280, T9S, T9L));
				   Ta2 = T9Y + T9Z;
				   Ta0 = T9Y - T9Z;
				   Rm[0] = KP500000000 * (FMA(KP980785280, Ta2, Ta1));
				   Rp[WS(rs, 15)] = KP500000000 * (FNMS(KP980785280, Ta2, Ta1));
				   Rp[WS(rs, 7)] = KP500000000 * (FMA(KP980785280, T9W, T9V));
				   Rm[WS(rs, 8)] = KP500000000 * (FNMS(KP980785280, T9W, T9V));
			      }
			 }
		    }
	       }
	  }