hb2_32.c

快速fft变换

					     Tg = cr[WS(rs, 2)];
					     T5z = FMA(KP707106781, T4F, T4E);
					     T4G = FNMS(KP707106781, T4F, T4E);
					     T6Y = T6s + T6z;
					     T6A = T6s - T6z;
					     Th = ci[WS(rs, 13)];
					     T1L = ci[WS(rs, 21)];
					     T1M = cr[WS(rs, 26)];
					     Tj = cr[WS(rs, 10)];
					     T1K = Tg - Th;
					     Ti = Tg + Th;
					     T6o = T1L - T1M;
					     T1N = T1L + T1M;
					     Tk = ci[WS(rs, 5)];
					     T1G = ci[WS(rs, 29)];
					     T1H = cr[WS(rs, 18)];
					}
					{
					     E T1F, Tl, T6n, T1I;
					     T4x = T1K + T1N;
					     T1O = T1K - T1N;
					     T1F = Tj - Tk;
					     Tl = Tj + Tk;
					     T6n = T1G - T1H;
					     T1I = T1G + T1H;
					     Tm = Ti + Tl;
					     T7Z = Ti - Tl;
					     T80 = T6n - T6o;
					     T6p = T6n + T6o;
					     T4w = T1I - T1F;
					     T1J = T1F + T1I;
					}
				   }
				   {
					E Tq, T1z, Tp, T6l, T1C, Tr, T1v, T1w;
					{
					     E Tn, To, T1A, T1B;
					     Tn = ci[WS(rs, 1)];
					     To = cr[WS(rs, 14)];
					     T1A = ci[WS(rs, 25)];
					     T1B = cr[WS(rs, 22)];
					     Tq = cr[WS(rs, 6)];
					     T1z = Tn - To;
					     Tp = Tn + To;
					     T6l = T1A - T1B;
					     T1C = T1A + T1B;
					     Tr = ci[WS(rs, 9)];
					     T1v = ci[WS(rs, 17)];
					     T1w = cr[WS(rs, 30)];
					}
					{
					     E T1u, Ts, T6k, T1x;
					     T4A = T1z + T1C;
					     T1D = T1z - T1C;
					     T1u = Tq - Tr;
					     Ts = Tq + Tr;
					     T6k = T1v - T1w;
					     T1x = T1v + T1w;
					     Tt = Tp + Ts;
					     T7X = Tp - Ts;
					     T7W = T6k - T6l;
					     T6m = T6k + T6l;
					     T4z = T1u + T1x;
					     T1y = T1u - T1x;
					}
				   }
				   {
					E T8r, T8q, T2Z, T1E, T1P, T2Y, T7Y, T81, T50, T51;
					T8r = T7X + T7W;
					T7Y = T7W - T7X;
					T81 = T7Z + T80;
					T8q = T7Z - T80;
					T6P = Tm - Tt;
					Tu = Tm + Tt;
					T9a = T81 + T7Y;
					T82 = T7Y - T81;
					T2Z = FMA(KP414213562, T1y, T1D);
					T1E = FNMS(KP414213562, T1D, T1y);
					T1P = FMA(KP414213562, T1O, T1J);
					T2Y = FNMS(KP414213562, T1J, T1O);
					T8X = T8q + T8r;
					T8s = T8q - T8r;
					T40 = T1P + T1E;
					T1Q = T1E - T1P;
					T4y = FNMS(KP414213562, T4x, T4w);
					T50 = FMA(KP414213562, T4w, T4x);
					T51 = FMA(KP414213562, T4z, T4A);
					T4B = FNMS(KP414213562, T4A, T4z);
					T3M = T2Y + T2Z;
					T30 = T2Y - T2Z;
					T5w = T50 + T51;
					T52 = T50 - T51;
				   }
			      }
			 }
		    }
		    {
			 E T7D, T5K, T4C, T7K, T7J, T7E, T83, T8w, T8t, T8i, T6r, T70, T6X, T6K;
			 {
			      E T6q, T8Y, T9e, T9b, T95, T8L, T8Q, T8H, T8M, T8I, T8R;
			      {
				   E Tv, T10, T7v, T7y, T7u;
				   T7D = Tf - Tu;
				   Tv = Tf + Tu;
				   T7u = T6p + T6m;
				   T6q = T6m - T6p;
				   T5K = T4B - T4y;
				   T4C = T4y + T4B;
				   T10 = TK + TZ;
				   T7K = TK - TZ;
				   T7J = T7t - T7u;
				   T7v = T7t + T7u;
				   T7y = T7w + T7x;
				   T7E = T7x - T7w;
				   {
					E T9t, T9x, T9p, T9u, T9q, T9y;
					{
					     E T9n, T7z, T9o, T7A, T7q, T7p;
					     T8Y = FNMS(KP707106781, T8X, T8W);
					     T9n = FMA(KP707106781, T8X, T8W);
					     cr[0] = Tv + T10;
					     T7p = Tv - T10;
					     ci[0] = T7v + T7y;
					     T7z = T7v - T7y;
					     T9o = T9c + T9d;
					     T9e = T9c - T9d;
					     T7A = T7s * T7p;
					     T7q = T7o * T7p;
					     T9b = FNMS(KP707106781, T9a, T99);
					     T9t = FMA(KP707106781, T9a, T99);
					     T9x = FMA(KP923879532, T9o, T9n);
					     T9p = FNMS(KP923879532, T9o, T9n);
					     ci[WS(rs, 16)] = FMA(T7o, T7z, T7A);
					     cr[WS(rs, 16)] = FNMS(T7s, T7z, T7q);
					     T9u = T94 + T91;
					     T95 = T91 - T94;
					}
					T9q = T9m * T9p;
					T9y = T3v * T9x;
					{
					     E T8F, T9z, T9v, T8G, T9A, T9w;
					     T83 = FMA(KP707106781, T82, T7V);
					     T8F = FNMS(KP707106781, T82, T7V);
					     T9z = FMA(KP923879532, T9u, T9t);
					     T9v = FNMS(KP923879532, T9u, T9t);
					     T8G = T8u + T8v;
					     T8w = T8u - T8v;
					     T8t = FMA(KP707106781, T8s, T8p);
					     T8L = FNMS(KP707106781, T8s, T8p);
					     T9A = T3v * T9z;
					     cr[WS(rs, 2)] = FNMS(T3x, T9z, T9y);
					     T9w = T9m * T9v;
					     cr[WS(rs, 18)] = FNMS(T9s, T9v, T9q);
					     T8Q = FMA(KP923879532, T8G, T8F);
					     T8H = FNMS(KP923879532, T8G, T8F);
					     ci[WS(rs, 2)] = FMA(T3x, T9x, T9A);
					     ci[WS(rs, 18)] = FMA(T9s, T9p, T9w);
					     T8M = T8h + T8a;
					     T8i = T8a - T8h;
					}
					T8I = T8E * T8H;
					T8R = T8P * T8Q;
				   }
			      }
			      {
				   E T7f, T7j, T7b, T7g, T7c, T7k;
				   {
					E T79, T8T, T8N, T7a, T8U, T8O;
					T6r = T6j + T6q;
					T79 = T6j - T6q;
					T8T = FMA(KP923879532, T8M, T8L);
					T8N = FNMS(KP923879532, T8M, T8L);
					T7a = T6Z - T6Y;
					T70 = T6Y + T6Z;
					T6X = T6P + T6W;
					T7f = T6W - T6P;
					T8U = T8P * T8T;
					cr[WS(rs, 30)] = FNMS(T8S, T8T, T8R);
					T8O = T8E * T8N;
					cr[WS(rs, 14)] = FNMS(T8K, T8N, T8I);
					T7j = FMA(KP707106781, T7a, T79);
					T7b = FNMS(KP707106781, T7a, T79);
					ci[WS(rs, 30)] = FMA(T8S, T8Q, T8U);
					ci[WS(rs, 14)] = FMA(T8K, T8H, T8O);
					T7g = T6A - T6J;
					T6K = T6A + T6J;
				   }
				   T7c = T78 * T7b;
				   T7k = T5X * T7j;
				   {
					E T97, T9g, T9i, T9j, T9f, T9k, T9h, T96;
					{
					     E T7l, T7h, T7m, T7i;
					     T7l = FMA(KP707106781, T7g, T7f);
					     T7h = FNMS(KP707106781, T7g, T7f);
					     T7m = T5X * T7l;
					     cr[WS(rs, 12)] = FNMS(T5Z, T7l, T7k);
					     T7i = T78 * T7h;
					     cr[WS(rs, 28)] = FNMS(T7e, T7h, T7c);
					     T9h = FMA(KP923879532, T95, T8Y);
					     T96 = FNMS(KP923879532, T95, T8Y);
					     ci[WS(rs, 12)] = FMA(T5Z, T7j, T7m);
					     ci[WS(rs, 28)] = FMA(T7e, T7b, T7i);
					}
					T97 = T8V * T96;
					T9g = T98 * T96;
					T9i = T3G * T9h;
					T9j = FMA(KP923879532, T9e, T9b);
					T9f = FNMS(KP923879532, T9e, T9b);
					T9k = T3J * T9h;
					cr[WS(rs, 10)] = FNMS(T3J, T9j, T9i);
					ci[WS(rs, 26)] = FMA(T8V, T9f, T9g);
					cr[WS(rs, 26)] = FNMS(T98, T9f, T97);
					ci[WS(rs, 10)] = FMA(T3G, T9j, T9k);
				   }
			      }
			 }
			 {
			      E T31, T3r, T1R, T3m, T33, T32, T3s, T2K, T8z, T8j;
			      {
				   E T73, T6L, T75, T71;
				   T73 = FMA(KP707106781, T6K, T6r);
				   T6L = FNMS(KP707106781, T6K, T6r);
				   T75 = FMA(KP707106781, T70, T6X);
				   T71 = FNMS(KP707106781, T70, T6X);
				   {
					E T76, T74, T72, T6M;
					T76 = T3d * T73;
					T74 = T3a * T73;
					T72 = T6O * T6L;
					T6M = T6i * T6L;
					ci[WS(rs, 4)] = FMA(T3a, T75, T76);
					cr[WS(rs, 4)] = FNMS(T3d, T75, T74);
					ci[WS(rs, 20)] = FMA(T6i, T71, T72);
					cr[WS(rs, 20)] = FNMS(T6O, T71, T6M);
				   }
			      }
			      {
				   E T7N, T7F, T7P, T7L;
				   T7N = T7D + T7E;
				   T7F = T7D - T7E;
				   T7P = T7K + T7J;
				   T7L = T7J - T7K;
				   {
					E T7O, T7G, T7Q, T7M;
					T7O = T4p * T7N;
					T7G = T7C * T7F;
					T7Q = T4p * T7P;
					T7M = T7C * T7L;
					cr[WS(rs, 8)] = FNMS(T4r, T7P, T7O);
					cr[WS(rs, 24)] = FNMS(T7I, T7L, T7G);
					ci[WS(rs, 8)] = FMA(T4r, T7N, T7Q);
					ci[WS(rs, 24)] = FMA(T7I, T7F, T7M);
				   }
			      }
			      T31 = FMA(KP923879532, T30, T2X);
			      T3r = FNMS(KP923879532, T30, T2X);
			      T8z = FMA(KP923879532, T8i, T83);
			      T8j = FNMS(KP923879532, T8i, T83);
			      {
				   E T8B, T8x, T8C, T8A;
				   T8B = FMA(KP923879532, T8w, T8t);
				   T8x = FNMS(KP923879532, T8w, T8t);
				   T8C = T1a * T8z;
				   T8A = T16 * T8z;
				   {
					E T8y, T8k, T2i, T2J;
					T8y = T8m * T8j;
					T8k = T7S * T8j;
					ci[WS(rs, 6)] = FMA(T16, T8B, T8C);
					cr[WS(rs, 6)] = FNMS(T1a, T8B, T8A);
					ci[WS(rs, 22)] = FMA(T7S, T8x, T8y);
					cr[WS(rs, 22)] = FNMS(T8m, T8x, T8k);
					T1R = FMA(KP923879532, T1Q, T1t);
					T3m = FNMS(KP923879532, T1Q, T1t);
					T33 = FNMS(KP668178637, T28, T2h);
					T2i = FMA(KP668178637, T2h, T28);
					T2J = FNMS(KP668178637, T2I, T2z);
					T32 = FMA(KP668178637, T2z, T2I);
					T3s = T2J + T2i;
					T2K = T2i - T2J;
				   }
			      }
			      {
				   E T5l, T53, T5f, T4D, T4K, T4R, T56, T5g;
				   T5l = FNMS(KP923879532, T52, T4Z);
				   T53 = FMA(KP923879532, T52, T4Z);
				   {
					E T3t, T3D, T3f, T2L;
					T3t = FNMS(KP831469612, T3s, T3r);
					T3D = FMA(KP831469612, T3s, T3r);
					T3f = FMA(KP831469612, T2K, T1R);
					T2L = FNMS(KP831469612, T2K, T1R);
					{
					     E T3n, T34, T3g, T2M;
					     T3n = T32 + T33;
					     T34 = T32 - T33;
					     T3g = T3e * T3f;
					     T2M = T1c * T2L;
					     {
						  E T3o, T3z, T3j, T35;
						  T3o = FNMS(KP831469612, T3n, T3m);
						  T3z = FMA(KP831469612, T3n, T3m);
						  T3j = FMA(KP831469612, T34, T31);
						  T35 = FNMS(KP831469612, T34, T31);
						  {
						       E T3u, T3p, T3E, T3A;
						       T3u = T3q * T3o;
						       T3p = T3l * T3o;
						       T3E = T3C * T3z;
						       T3A = T3y * T3z;
						       {
							    E T3k, T36, T54, T55;
							    T3k = T3e * T3j;
							    cr[WS(rs, 5)] = FNMS(T3i, T3j, T3g);
							    T36 = T1c * T35;
							    cr[WS(rs, 21)] = FNMS(T2O, T35, T2M);
							    ci[WS(rs, 13)] = FMA(T3l, T3t, T3u);
							    cr[WS(rs, 13)] = FNMS(T3q, T3t, T3p);
							    ci[WS(rs, 29)] = FMA(T3y, T3D, T3E);
							    cr[WS(rs, 29)] = FNMS(T3C, T3D, T3A);
							    ci[WS(rs, 5)] = FMA(T3i, T3f, T3k);
							    ci[WS(rs, 21)] = FMA(T2O, T2L, T36);
							    T5f = FMA(KP923879532, T4C, T4v);
							    T4D = FNMS(KP923879532, T4C, T4v);
							    T4K = FNMS(KP668178637, T4J, T4G);
							    T54 = FMA(KP668178637, T4G, T4J);
							    T55 = FMA(KP668178637, T4N, T4Q);
							    T4R = FNMS(KP668178637, T4Q, T4N);
							    T56 = T54 - T55;
							    T5g = T54 + T55;
						       }
						  }
					     }
					}
				   }
				   {
					E T4h, T41, T4c, T3N, T3Q, T3T, T44, T4d;
					T4h = FNMS(KP923879532, T40, T3Z);
					T41 = FMA(KP923879532, T40, T3Z);
					{
					     E T57, T5b, T5h, T5p;
					     T57 = FNMS(KP831469612, T56, T53);
					     T5b = FMA(KP831469612, T56, T53);
					     T5h = FNMS(KP831469612, T5g, T5f);
					     T5p = FMA(KP831469612, T5g, T5f);
					     {
						  E T5m, T4S, T5i, T5q;
						  T5m = T4K - T4R;
						  T4S = T4K + T4R;
						  T5i = T5e * T5h;
						  T5q = T17 * T5p;
						  {
						       E T5n, T5r, T59, T4T;
						       T5n = FMA(KP831469612, T5m, T5l);
						       T5r = FNMS(KP831469612, T5m, T5l);
						       T59 = FMA(KP831469612, T4S, T4D);
						       T4T = FNMS(KP831469612, T4S, T4D);
						       {
							    E T5o, T5s, T5c, T5a;
							    T5o = T5e * T5n;
							    cr[WS(rs, 11)] = FNMS(T5k, T5n, T5i);
							    T5s = T17 * T5r;
							    cr[WS(rs, 27)] = FNMS(T1b, T5r, T5q);
							    T5c = T14 * T59;
							    T5a = T11 * T59;
							    {
								 E T58, T4U, T42, T43;
								 T58 = T4W * T4T;
								 T4U = T4s * T4T;
								 ci[WS(rs, 11)] = FMA(T5k, T5h, T5o);
								 ci[WS(rs, 27)] = FMA(T1b, T5p, T5s);
								 ci[WS(rs, 3)] = FMA(T11, T5b, T5c);
								 cr[WS(rs, 3)] = FNMS(T14, T5b, T5a);
								 ci[WS(rs, 19)] = FMA(T4s, T57, T58);
								 cr[WS(rs, 19)] = FNMS(T4W, T57, T4U);
								 T4c = FNMS(KP923879532, T3M, T3L);
								 T3N = FMA(KP923879532, T3M, T3L);
								 T3Q = FNMS(KP198912367, T3P, T3O);
								 T42 = FMA(KP198912367, T3O, T3P);
								 T43 = FNMS(KP198912367, T3R, T3S);
								 T3T = FMA(KP198912367, T3S, T3R);
								 T44 = T42 + T43;
								 T4d = T43 - T42;
							    }
						       }
						  }
					     }
					}
					T67 = FNMS(KP923879532, T5K, T5J);
					T5L = FMA(KP923879532, T5K, T5J);
					{
					     E T45, T49, T4e, T4l;
					     T45 = FNMS(KP980785280, T44, T41);
					     T49 = FMA(KP980785280, T44, T41);
					     T4e = FNMS(KP980785280, T4d, T4c);
					     T4l = FMA(KP980785280, T4d, T4c);
					     {
						  E T4i, T3U, T4f, T4m;
						  T4i = T3Q - T3T;
						  T3U = T3Q + T3T;
						  T4f = T4b * T4e;
						  T4m = T12 * T4l;
						  {
						       E T4j, T4n, T47, T3V;
						       T4j = FNMS(KP980785280, T4i, T4h);
						       T4n = FMA(KP980785280, T4i, T4h);
						       T47 = FMA(KP980785280, T3U, T3N);
						       T3V = FNMS(KP980785280, T3U, T3N);
						       {
							    E T4k, T4o, T4a, T48;
							    T4k = T4b * T4j;
							    cr[WS(rs, 25)] = FNMS(T4g, T4j, T4f);
							    T4o = T12 * T4n;
							    cr[WS(rs, 9)] = FNMS(T15, T4n, T4m);
							    T4a = T39 * T47;
							    T48 = T37 * T47;
							    {
								 E T46, T3W, T5M, T5N;
								 T46 = T3Y * T3V;
								 T3W = T3K * T3V;
								 ci[WS(rs, 25)] = FMA(T4g, T4e, T4k);
								 ci[WS(rs, 9)] = FMA(T15, T4l, T4o);
								 ci[WS(rs, 1)] = FMA(T37, T49, T4a);
								 cr[WS(rs, 1)] = FNMS(T39, T49, T48);
								 ci[WS(rs, 17)] = FMA(T3K, T45, T46);
								 cr[WS(rs, 17)] = FNMS(T3Y, T45, T3W);
								 T61 = FMA(KP923879532, T5w, T5v);
								 T5x = FNMS(KP923879532, T5w, T5v);
								 T5A = FNMS(KP198912367, T5z, T5y);
								 T5M = FMA(KP198912367, T5y, T5z);
								 T5N = FMA(KP198912367, T5B, T5C);
								 T5D = FNMS(KP198912367, T5C, T5B);
								 T5O = T5M - T5N;
								 T62 = T5M + T5N;
							    }
						       }
						  }
					     }
					}
				   }
			      }
			 }
		    }
	       }
	  }
	  T5V = FMA(KP980785280, T5O, T5L);
	  T5P = FNMS(KP980785280, T5O, T5L);
	  {
	       E T6c, T63, T5E, T68;
	       T6c = FMA(KP980785280, T62, T61);
	       T63 = FNMS(KP980785280, T62, T61);
	       T5E = T5A + T5D;
	       T68 = T5D - T5A;
	       {
		    E T64, T6d, T6f, T69;
		    T64 = T60 * T63;
		    T6d = T6b * T6c;
		    T6f = FNMS(KP980785280, T68, T67);
		    T69 = FMA(KP980785280, T68, T67);
		    {