hb_32.c

最新的FFT程序

			 }
		    }
	       }
	       {
		    E T6D, T5y, T2C, T2H, T2E, T2z;
		    {
			 E T4o, T4v, T4y, T4r, T3C, T3Y, T3V, T4i, T3N, T4h, T4l, T4j, T4k;
			 {
			      E T4H, T4K, T4N, T4J, T4O;
			      {
				   E Tv, T10, T4L, T4I, T4M;
				   T4o = Tf - Tu;
				   Tv = Tf + Tu;
				   T10 = TK + TZ;
				   T4v = TZ - TK;
				   T4y = T4w - T4x;
				   T4L = T4x + T4w;
				   T6D = FNMS(KP707106781, T5x, T5u);
				   T5y = FMA(KP707106781, T5x, T5u);
				   T4I = Tv - T10;
				   T4r = T4p - T4q;
				   T4M = T4p + T4q;
				   rio[0] = Tv + T10;
				   T4H = W[30];
				   T4K = W[31];
				   iio[-WS(ios, 31)] = T4M + T4L;
				   T4N = T4L - T4M;
				   T4J = T4H * T4I;
				   T4O = T4K * T4I;
			      }
			      {
				   E T4d, T4a, T4e, T48, T49;
				   T3C = T3y + T3B;
				   T48 = T3y - T3B;
				   T49 = T3X - T3W;
				   T3Y = T3W + T3X;
				   rio[WS(ios, 16)] = FNMS(T4K, T4N, T4J);
				   iio[-WS(ios, 15)] = FMA(T4H, T4N, T4O);
				   T3V = T3R + T3U;
				   T4d = T3U - T3R;
				   T4i = FNMS(KP707106781, T49, T48);
				   T4a = FMA(KP707106781, T49, T48);
				   T4e = T3H - T3M;
				   T3N = T3H + T3M;
				   {
					E T47, T4c, T4f, T4b, T4g;
					T47 = W[6];
					T4c = W[7];
					T4h = W[38];
					T4l = FNMS(KP707106781, T4e, T4d);
					T4f = FMA(KP707106781, T4e, T4d);
					T4b = T47 * T4a;
					T4g = T4c * T4a;
					T4j = T4h * T4i;
					T4k = W[39];
					rio[WS(ios, 4)] = FNMS(T4c, T4f, T4b);
					iio[-WS(ios, 27)] = FMA(T47, T4f, T4g);
				   }
			      }
			 }
			 {
			      E T4s, T4C, T4F, T4z, T4m, T4n, T4u;
			      rio[WS(ios, 20)] = FNMS(T4k, T4l, T4j);
			      T4m = T4k * T4i;
			      iio[-WS(ios, 11)] = FMA(T4h, T4l, T4m);
			      T4s = T4o + T4r;
			      T4C = T4o - T4r;
			      T4F = T4y - T4v;
			      T4z = T4v + T4y;
			      T4n = W[46];
			      T4u = W[47];
			      {
				   E T4B, T4E, T4t, T4A, T4D, T4G;
				   T4B = W[14];
				   T4E = W[15];
				   T4t = T4n * T4s;
				   T4A = T4u * T4s;
				   T4D = T4B * T4C;
				   T4G = T4E * T4C;
				   rio[WS(ios, 24)] = FNMS(T4u, T4z, T4t);
				   iio[-WS(ios, 7)] = FMA(T4n, T4z, T4A);
				   rio[WS(ios, 8)] = FNMS(T4E, T4F, T4D);
				   iio[-WS(ios, 23)] = FMA(T4B, T4F, T4G);
			      }
			 }
			 {
			      E T45, T41, T44, T43, T46;
			      {
				   E T3x, T42, T3O, T3Z, T3Q, T3P, T40;
				   T3x = W[54];
				   T42 = FNMS(KP707106781, T3N, T3C);
				   T3O = FMA(KP707106781, T3N, T3C);
				   T45 = FNMS(KP707106781, T3Y, T3V);
				   T3Z = FMA(KP707106781, T3Y, T3V);
				   T3Q = W[55];
				   T3P = T3x * T3O;
				   T41 = W[22];
				   T44 = W[23];
				   T40 = T3Q * T3O;
				   rio[WS(ios, 28)] = FNMS(T3Q, T3Z, T3P);
				   T43 = T41 * T42;
				   T46 = T44 * T42;
				   iio[-WS(ios, 3)] = FMA(T3x, T3Z, T40);
			      }
			      {
				   E T2B, T2q, T2G, T2A, T2n, T2F, T28, T2u, T27, T1u;
				   T2B = T2p - T2o;
				   T2q = T2o + T2p;
				   rio[WS(ios, 12)] = FNMS(T44, T45, T43);
				   iio[-WS(ios, 19)] = FMA(T41, T45, T46);
				   T2G = T1N - T26;
				   T27 = T1N + T26;
				   T1u = FMA(KP707106781, T1t, T1a);
				   T2A = FNMS(KP707106781, T1t, T1a);
				   T2n = FMA(KP707106781, T2m, T2j);
				   T2F = FNMS(KP707106781, T2m, T2j);
				   T28 = FMA(KP923879532, T27, T1u);
				   T2u = FNMS(KP923879532, T27, T1u);
				   {
					E T2x, T2K, T2N, T2M, T2J, T2O;
					{
					     E T2r, T2a, T11, T2s, T29;
					     T2r = FMA(KP923879532, T2q, T2n);
					     T2x = FNMS(KP923879532, T2q, T2n);
					     T2a = W[59];
					     T11 = W[58];
					     T2C = FMA(KP923879532, T2B, T2A);
					     T2K = FNMS(KP923879532, T2B, T2A);
					     T2s = T2a * T28;
					     T29 = T11 * T28;
					     T2H = FMA(KP923879532, T2G, T2F);
					     T2N = FNMS(KP923879532, T2G, T2F);
					     T2M = W[43];
					     iio[-WS(ios, 1)] = FMA(T11, T2r, T2s);
					     rio[WS(ios, 30)] = FNMS(T2a, T2r, T29);
					     T2J = W[42];
					     T2O = T2M * T2K;
					}
					{
					     E T2w, T2t, T2y, T2L, T2v;
					     T2w = W[27];
					     T2L = T2J * T2K;
					     iio[-WS(ios, 9)] = FMA(T2J, T2N, T2O);
					     T2t = W[26];
					     T2y = T2w * T2u;
					     rio[WS(ios, 22)] = FNMS(T2M, T2N, T2L);
					     T2v = T2t * T2u;
					     iio[-WS(ios, 17)] = FMA(T2t, T2x, T2y);
					     T2E = W[11];
					     T2z = W[10];
					     rio[WS(ios, 14)] = FNMS(T2w, T2x, T2v);
					}
				   }
			      }
			 }
		    }
		    {
			 E T74, T6Z, T79, T6W, T6X, T70, T6Y, T71, T76;
			 {
			      E T3p, T3m, T3h, T3q, T3l;
			      {
				   E T38, T3j, T2Z, T3o, T3i, T2S, T3n, T35, T2I, T2D;
				   T38 = T36 + T37;
				   T3j = T36 - T37;
				   T2I = T2E * T2C;
				   T2D = T2z * T2C;
				   T2Z = T2V + T2Y;
				   T3o = T2Y - T2V;
				   iio[-WS(ios, 25)] = FMA(T2z, T2H, T2I);
				   rio[WS(ios, 6)] = FNMS(T2E, T2H, T2D);
				   T3i = FNMS(KP707106781, T2R, T2Q);
				   T2S = FMA(KP707106781, T2R, T2Q);
				   T3n = FNMS(KP707106781, T34, T33);
				   T35 = FMA(KP707106781, T34, T33);
				   {
					E T3f, T3c, T3k, T3v, T3u, T3t;
					{
					     E T39, T30, T3s, T32, T2P;
					     T3f = FNMS(KP923879532, T38, T35);
					     T39 = FMA(KP923879532, T38, T35);
					     T3c = FNMS(KP923879532, T2Z, T2S);
					     T30 = FMA(KP923879532, T2Z, T2S);
					     T3s = FNMS(KP923879532, T3j, T3i);
					     T3k = FMA(KP923879532, T3j, T3i);
					     T32 = W[3];
					     T2P = W[2];
					     {
						  E T3r, T3w, T3a, T31;
						  T3v = FNMS(KP923879532, T3o, T3n);
						  T3p = FMA(KP923879532, T3o, T3n);
						  T3u = W[19];
						  T3a = T32 * T30;
						  T31 = T2P * T30;
						  T3r = W[18];
						  T3w = T3u * T3s;
						  iio[-WS(ios, 29)] = FMA(T2P, T39, T3a);
						  rio[WS(ios, 2)] = FNMS(T32, T39, T31);
						  T3t = T3r * T3s;
						  iio[-WS(ios, 21)] = FMA(T3r, T3v, T3w);
					     }
					}
					{
					     E T3e, T3b, T3g, T3d;
					     T3e = W[35];
					     rio[WS(ios, 10)] = FNMS(T3u, T3v, T3t);
					     T3b = W[34];
					     T3g = T3e * T3c;
					     T3m = W[51];
					     T3d = T3b * T3c;
					     T3h = W[50];
					     iio[-WS(ios, 13)] = FMA(T3b, T3f, T3g);
					     T3q = T3m * T3k;
					     rio[WS(ios, 18)] = FNMS(T3e, T3f, T3d);
					     T3l = T3h * T3k;
					}
				   }
			      }
			      {
				   E T77, T6C, T6P, T72, T6Q, T6R, T73, T6J, T6F, T6I;
				   T77 = FNMS(KP923879532, T6B, T6A);
				   T6C = FMA(KP923879532, T6B, T6A);
				   iio[-WS(ios, 5)] = FMA(T3h, T3p, T3q);
				   T6P = FMA(KP923879532, T6O, T6N);
				   T72 = FNMS(KP923879532, T6O, T6N);
				   rio[WS(ios, 26)] = FNMS(T3m, T3p, T3l);
				   T6Q = FMA(KP668178637, T6D, T6E);
				   T6F = FNMS(KP668178637, T6E, T6D);
				   T6I = FMA(KP668178637, T6H, T6G);
				   T6R = FNMS(KP668178637, T6G, T6H);
				   T73 = T6I - T6F;
				   T6J = T6F + T6I;
				   {
					E T6z, T7f, T6K, T6M, T6L, T7c, T6T, T7b, T7e, T78, T6S;
					T6z = W[56];
					T74 = FMA(KP831469612, T73, T72);
					T7f = FNMS(KP831469612, T73, T72);
					T6Z = FNMS(KP831469612, T6J, T6C);
					T6K = FMA(KP831469612, T6J, T6C);
					T78 = T6Q - T6R;
					T6S = T6Q + T6R;
					T6M = W[57];
					T6L = T6z * T6K;
					T79 = FMA(KP831469612, T78, T77);
					T7c = FNMS(KP831469612, T78, T77);
					T6W = FNMS(KP831469612, T6S, T6P);
					T6T = FMA(KP831469612, T6S, T6P);
					T7b = W[40];
					T7e = W[41];
					{
					     E T6V, T6U, T7g, T7d;
					     T6V = W[24];
					     iio[-WS(ios, 2)] = FMA(T6M, T6T, T6L);
					     T6U = T6z * T6T;
					     T7g = T7b * T7f;
					     T7d = T7b * T7c;
					     T6X = T6V * T6W;
					     rio[WS(ios, 29)] = FNMS(T6M, T6K, T6U);
					     rio[WS(ios, 21)] = FNMS(T7e, T7c, T7g);
					     iio[-WS(ios, 10)] = FMA(T7e, T7f, T7d);
					     T70 = T6V * T6Z;
					}
					T6Y = W[25];
					T71 = W[8];
					T76 = W[9];
				   }
			      }
			 }
			 {
			      E T6m, T6h, T6r, T6e, T6f, T6i, T6g, T6j, T6o;
			      {
				   E T6p, T5g, T6k, T67, T68, T69, T6l, T5T;
				   {
					E T5z, T5S, T7a, T75;
					T6p = FNMS(KP923879532, T5f, T50);
					T5g = FMA(KP923879532, T5f, T50);
					iio[-WS(ios, 18)] = FMA(T6Y, T6W, T70);
					rio[WS(ios, 13)] = FNMS(T6Y, T6Z, T6X);
					T7a = T71 * T79;
					T75 = T71 * T74;
					T6k = FNMS(KP923879532, T66, T63);
					T67 = FMA(KP923879532, T66, T63);
					T68 = FMA(KP198912367, T5r, T5y);
					T5z = FNMS(KP198912367, T5y, T5r);
					iio[-WS(ios, 26)] = FMA(T76, T74, T7a);
					rio[WS(ios, 5)] = FNMS(T76, T79, T75);
					T5S = FMA(KP198912367, T5R, T5K);
					T69 = FNMS(KP198912367, T5K, T5R);
					T6l = T5S - T5z;
					T5T = T5z + T5S;
				   }
				   {
					E T4P, T6x, T5U, T5W, T5V, T6u, T6b, T6t, T6w, T6q, T6a;
					T4P = W[0];
					T6m = FMA(KP980785280, T6l, T6k);
					T6x = FNMS(KP980785280, T6l, T6k);
					T6h = FNMS(KP980785280, T5T, T5g);
					T5U = FMA(KP980785280, T5T, T5g);
					T6q = T68 - T69;
					T6a = T68 + T69;
					T5W = W[1];
					T5V = T4P * T5U;
					T6r = FMA(KP980785280, T6q, T6p);
					T6u = FNMS(KP980785280, T6q, T6p);
					T6e = FNMS(KP980785280, T6a, T67);
					T6b = FMA(KP980785280, T6a, T67);
					T6t = W[16];
					T6w = W[17];
					{
					     E T6d, T6c, T6y, T6v;
					     T6d = W[32];
					     rio[WS(ios, 1)] = FNMS(T5W, T6b, T5V);
					     T6c = T4P * T6b;
					     T6y = T6t * T6x;
					     T6v = T6t * T6u;
					     T6f = T6d * T6e;
					     iio[-WS(ios, 30)] = FMA(T5W, T5U, T6c);
					     iio[-WS(ios, 22)] = FMA(T6w, T6u, T6y);
					     rio[WS(ios, 9)] = FNMS(T6w, T6x, T6v);
					     T6i = T6d * T6h;
					}
					T6g = W[33];
					T6j = W[48];
					T6o = W[49];
				   }
			      }
			      {
				   E T8O, T8J, T8T, T8G, T8H, T8K, T8I, T8L, T8Q;
				   {
					E T8R, T8m, T8M, T8z, T8A, T8B, T8N, T8t;
					{
					     E T8p, T8s, T6s, T6n;
					     T8R = FMA(KP923879532, T8l, T8k);
					     T8m = FNMS(KP923879532, T8l, T8k);
					     rio[WS(ios, 17)] = FNMS(T6g, T6e, T6i);
					     iio[-WS(ios, 14)] = FMA(T6g, T6h, T6f);
					     T6s = T6j * T6r;
					     T6n = T6j * T6m;
					     T8M = FNMS(KP923879532, T8y, T8x);
					     T8z = FMA(KP923879532, T8y, T8x);
					     T8A = FMA(KP668178637, T8n, T8o);
					     T8p = FNMS(KP668178637, T8o, T8n);
					     rio[WS(ios, 25)] = FNMS(T6o, T6m, T6s);
					     iio[-WS(ios, 6)] = FMA(T6o, T6r, T6n);
					     T8s = FNMS(KP668178637, T8r, T8q);
					     T8B = FMA(KP668178637, T8q, T8r);
					     T8N = T8s - T8p;
					     T8t = T8p + T8s;
					}
					{
					     E T8j, T8Z, T8u, T8w, T8v, T8W, T8D, T8V, T8Y, T8S, T8C;
					     T8j = W[4];
					     T8O = FMA(KP831469612, T8N, T8M);
					     T8Z = FNMS(KP831469612, T8N, T8M);
					     T8J = FNMS(KP831469612, T8t, T8m);
					     T8u = FMA(KP831469612, T8t, T8m);
					     T8S = T8A + T8B;
					     T8C = T8A - T8B;
					     T8w = W[5];
					     T8v = T8j * T8u;
					     T8T = FMA(KP831469612, T8S, T8R);
					     T8W = FNMS(KP831469612, T8S, T8R);
					     T8G = FNMS(KP831469612, T8C, T8z);
					     T8D = FMA(KP831469612, T8C, T8z);
					     T8V = W[20];
					     T8Y = W[21];
					     {
						  E T8F, T8E, T90, T8X;
						  T8F = W[36];
						  rio[WS(ios, 3)] = FNMS(T8w, T8D, T8v);
						  T8E = T8j * T8D;
						  T90 = T8V * T8Z;
						  T8X = T8V * T8W;
						  T8H = T8F * T8G;
						  iio[-WS(ios, 28)] = FMA(T8w, T8u, T8E);
						  iio[-WS(ios, 20)] = FMA(T8Y, T8W, T90);
						  rio[WS(ios, 11)] = FNMS(T8Y, T8Z, T8X);
						  T8K = T8F * T8J;
					     }
					     T8I = W[37];
					     T8L = W[52];
					     T8Q = W[53];
					}
				   }
				   {
					E T89, T7s, T84, T7R, T7S, T7T, T85, T7H;
					{
					     E T7z, T7G, T8U, T8P;
					     T89 = FNMS(KP923879532, T7r, T7k);
					     T7s = FMA(KP923879532, T7r, T7k);
					     rio[WS(ios, 19)] = FNMS(T8I, T8G, T8K);
					     iio[-WS(ios, 12)] = FMA(T8I, T8J, T8H);
					     T8U = T8L * T8T;
					     T8P = T8L * T8O;
					     T84 = FNMS(KP923879532, T7Q, T7N);
					     T7R = FMA(KP923879532, T7Q, T7N);
					     T7S = FMA(KP198912367, T7v, T7y);
					     T7z = FNMS(KP198912367, T7y, T7v);
					     rio[WS(ios, 27)] = FNMS(T8Q, T8O, T8U);
					     iio[-WS(ios, 4)] = FMA(T8Q, T8T, T8P);
					     T7G = FNMS(KP198912367, T7F, T7C);
					     T7T = FMA(KP198912367, T7C, T7F);
					     T85 = T7z + T7G;
					     T7H = T7z - T7G;
					}
					{
					     E T7h, T8h, T7I, T7K, T7J, T8e, T7V, T8d, T8g, T8a, T7U;
					     T7h = W[60];
					     T86 = FNMS(KP980785280, T85, T84);
					     T8h = FMA(KP980785280, T85, T84);
					     T81 = FNMS(KP980785280, T7H, T7s);
					     T7I = FMA(KP980785280, T7H, T7s);
					     T8a = T7S - T7T;
					     T7U = T7S + T7T;
					     T7K = W[61];
					     T7J = T7h * T7I;
					     T8b = FMA(KP980785280, T8a, T89);
					     T8e = FNMS(KP980785280, T8a, T89);
					     T7Y = FNMS(KP980785280, T7U, T7R);
					     T7V = FMA(KP980785280, T7U, T7R);
					     T8d = W[44];
					     T8g = W[45];
					     {
						  E T7X, T7W, T8i, T8f;
						  T7X = W[28];
						  iio[0] = FMA(T7K, T7V, T7J);
						  T7W = T7h * T7V;
						  T8i = T8d * T8h;
						  T8f = T8d * T8e;
						  T7Z = T7X * T7Y;
						  rio[WS(ios, 31)] = FNMS(T7K, T7I, T7W);
						  rio[WS(ios, 23)] = FNMS(T8g, T8e, T8i);
						  iio[-WS(ios, 8)] = FMA(T8g, T8h, T8f);
						  T82 = T7X * T81;
					     }
					     T80 = W[29];
					     T83 = W[12];
					     T88 = W[13];
					}
				   }
			      }
			 }
		    }
	       }
	  }
	  iio[-WS(ios, 16)] = FMA(T80, T7Y, T82);