hf2_32.c

快速傅立叶变换库函数

					E T3N, T3K, T3O, T5H, T41, T5Q, T3R, T3U, T3W;
					{
					     E T3H, T3I, T3J, T3Y, T40, T5G, T3Z, T5P;
					     {
						  E T30, T56, T5a, T36, T5f, T7B, T5b;
						  T3H = rio[WS(ios, 3)];
						  {
						       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);
						       T3I = T3 * T3H;
						  }
						  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;
						       T3J = iio[-WS(ios, 28)];
						  }
					     }
					     T3Y = rio[WS(ios, 11)];
					     T40 = iio[-WS(ios, 20)];
					     T3N = rio[WS(ios, 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 = iio[-WS(ios, 12)];
					     T3U = rio[WS(ios, 27)];
					     T3W = iio[-WS(ios, 4)];
					}
					{
					     E T3g, T3h, T3k, T3A, T3D, T5V, T3B, T5B;
					     {
						  E T3T, T5F, T5J, T3X, T5O, T7M, T5K;
						  T3g = rio[WS(ios, 31)];
						  {
						       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);
						       T3h = T3f * T3g;
						  }
						  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;
						       T3k = iio[0];
						  }
					     }
					     T3A = rio[WS(ios, 23)];
					     T3D = iio[-WS(ios, 8)];
					     T3o = rio[WS(ios, 15)];
					     T3l = FMA(T3j, T3k, T3h);
					     T5V = T3f * T3k;
					     T3B = T3z * T3A;
					     T5B = T3z * T3D;
					     T3p = T3n * T3o;
					     T5W = FNMS(T3j, T3g, T5V);
					     T3E = FMA(T3C, T3D, T3B);
					     T5C = FNMS(T3C, T3A, T5B);
					     T3s = iio[-WS(ios, 16)];
					     T3v = rio[WS(ios, 7)];
					     T3x = iio[-WS(ios, 24)];
					}
				   }
			      }
			 }
		    }
		    {
			 E T61, T6P, T5E, T6M, T88, T90, T8Z, T8b;
			 {
			      E T7e, T8T, T7L, T7J, T7h, T8U, T8S, T8R;
			      {
				   E T8c, T1i, T8A, T8z, T8O, T8J, T8N, T2l, T8L, T45, T8t, T8l, T8u, T8q, T3G;
				   E T8k, T8p, T8w, T2m;
				   {
					E T8x, T8y, T8n, T8C, T8I;
					{
					     E TI, T3u, T5y, T5Y, T3y, T5A, T1h, T7H, T5Z;
					     TI = Tq + TH;
					     T7e = Tq - TH;
					     {
						  E T3t, T5X, T3w, T5z;
						  T3t = FMA(T3r, T3s, T3p);
						  T5X = T3n * T3s;
						  T3w = TP * T3v;
						  T5z = TP * T3x;
						  T3u = T3l + T3t;
						  T5y = T3l - T3t;
						  T5Y = FNMS(T3r, T3o, T5X);
						  T3y = FMA(TR, T3x, T3w);
						  T5A = FNMS(TR, T3v, T5z);
						  T1h = TZ + T1g;
						  T8T = T1g - TZ;
					     }
					     T7H = T5W + T5Y;
					     T5Z = T5W - T5Y;
					     {
						  E T3F, T60, T7I, T5D;
						  T3F = T3y + T3E;
						  T60 = T3E - T3y;
						  T7I = T5A + T5C;
						  T5D = T5A - T5C;
						  T61 = T5Z + T60;
						  T6P = T60 - T5Z;
						  T3G = T3u + T3F;
						  T7L = T3u - T3F;
						  T5E = T5y + T5D;
						  T6M = T5y - T5D;
						  T7J = T7H - T7I;
						  T8n = T7H + T7I;
						  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;
					     rio[0] = T2m + T45;
					     iio[-WS(ios, 16)] = T2m - T45;
					     rio[WS(ios, 8)] = T8w + T8z;
					     iio[-WS(ios, 24)] = T8w - T8z;
					     T8M = T8J - T8B;
					     T8K = T8B + T8J;
					     T8g = T8c + T8f;
					     T8s = T8c - T8f;
					     T8R = T8O - T8N;
					     T8P = T8N + T8O;
					     iio[-WS(ios, 8)] = T8L + T8M;
					     rio[WS(ios, 24)] = T8L - T8M;
					     iio[0] = T8A + T8K;
					     rio[WS(ios, 16)] = T8A - T8K;
					     rio[WS(ios, 4)] = FMA(KP707106781, T8r, T8g);
					     iio[-WS(ios, 20)] = FNMS(KP707106781, T8r, T8g);
					     T8Q = T8t + T8u;
					     T8v = T8t - T8u;
					}
					iio[-WS(ios, 4)] = FMA(KP707106781, T8Q, T8P);
					rio[WS(ios, 20)] = FMS(KP707106781, T8Q, T8P);
					rio[WS(ios, 12)] = FMA(KP707106781, T8v, T8s);
					iio[-WS(ios, 28)] = 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;
					iio[-WS(ios, 12)] = FMA(KP707106781, T8S, T8R);
					rio[WS(ios, 28)] = 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;
						  iio[-WS(ios, 30)] = FMA(KP923879532, T7V, T7S);
						  rio[WS(ios, 14)] = FNMS(KP923879532, T7V, T7S);
						  iio[-WS(ios, 6)] = FMA(KP923879532, T94, T93);
						  rio[WS(ios, 22)] = FMS(KP923879532, T94, T93);
						  iio[-WS(ios, 14)] = FMA(KP923879532, T96, T95);
						  rio[WS(ios, 30)] = FMS(KP923879532, T96, T95);
						  rio[WS(ios, 6)] = FMA(KP923879532, T7R, T7u);
						  iio[-WS(ios, 22)] = 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;
					     rio[WS(ios, 2)] = FMA(KP923879532, T87, T80);
					     iio[-WS(ios, 18)] = FNMS(KP923879532, T87, T80);
					     iio[-WS(ios, 2)] = FMA(KP923879532, T8Y, T8X);
					     rio[WS(ios, 18)] = 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;
					     rio[WS(ios, 10)] = FMA(KP923879532, T8b, T88);
					     iio[-WS(ios, 26)] = FNMS(KP923879532, T8b, T88);
					     iio[-WS(ios, 10)] = FMA(KP923879532, T90, T8Z);
					     rio[WS(ios, 26)] = 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;
					     iio[-WS(ios, 1)] = FMA(KP980785280, T9e, T9d);
					     rio[WS(ios, 17)] = FMS(KP980785280, T9e, T9d);
					     rio[WS(ios, 1)] = FMA(KP980785280, T67, T4Y);
					     iio[-WS(ios, 17)] = FNMS(KP980785280, T67, T4Y);
					     rio[WS(ios, 9)] = FMA(KP980785280, T6b, T68);
					     iio[-WS(ios, 25)] = 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);
					     iio[-WS(ios, 9)] = FMA(KP980785280, T9g, T9f);
					     rio[WS(ios, 25)] = FMS(KP980785280, T9g, T9f);
					     iio[-WS(ios, 29)] = FMA(KP831469612, T6n, T6g);
					     rio[WS(ios, 13)] = FNMS(KP831469612, T6n, T6g);
					     T9m = T6p + T6q;
					     T6r = T6p - T6q;
					     iio[-WS(ios, 13)] = FNMS(KP831469612, T9m, T9l);
					     rio[WS(ios, 29)] = -(FMA(KP831469612, T9m, T9l));
					     rio[WS(ios, 5)] = FMA(KP831469612, T6r, T6o);
					     iio[-WS(ios, 21)] = 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);
					iio[-WS(ios, 5)] = FMA(KP831469612, T9k, T9j);
					rio[WS(ios, 21)] = 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;
						       rio[WS(ios, 7)] = FMA(KP980785280, T79, T72);
						       iio[-WS(ios, 23)] = FNMS(KP980785280, T79, T72);
						       rio[WS(ios, 11)] = FMA(KP831469612, T6X, T6U);
						       iio[-WS(ios, 27)] = FNMS(KP831469612, T6X, T6U);
						       iio[-WS(ios, 3)] = FMA(KP831469612, T9s, T9r);
						       rio[WS(ios, 19)] = FMS(KP831469612, T9s, T9r);