hf2_32.c

快速fft变换

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