q1_8.c

快速fft变换

	       T8O = T7N - T7Q;
	       T90 = T8I + T8J;
	       T84 = T80 + T83;
	       T8m = T8i - T8l;
	       T8A = T83 - T80;
	       T8K = T8I - T8J;
	       T8w = T8i + T8l;
	  }
	  {
	       E T9c, T9s, T9v, Ta7, T9f, T9n, T9q, Ta8;
	       {
		    E T9a, T9b, T9t, T9u;
		    T9a = rio[WS(vs, 7) + WS(rs, 1)];
		    T9b = rio[WS(vs, 7) + WS(rs, 5)];
		    T9c = T9a + T9b;
		    T9s = T9a - T9b;
		    T9t = iio[WS(vs, 7) + WS(rs, 1)];
		    T9u = iio[WS(vs, 7) + WS(rs, 5)];
		    T9v = T9t - T9u;
		    Ta7 = T9t + T9u;
	       }
	       {
		    E T9d, T9e, T9o, T9p;
		    T9d = rio[WS(vs, 7) + WS(rs, 7)];
		    T9e = rio[WS(vs, 7) + WS(rs, 3)];
		    T9f = T9d + T9e;
		    T9n = T9d - T9e;
		    T9o = iio[WS(vs, 7) + WS(rs, 7)];
		    T9p = iio[WS(vs, 7) + WS(rs, 3)];
		    T9q = T9o - T9p;
		    Ta8 = T9o + T9p;
	       }
	       T9g = T9c + T9f;
	       Ta9 = Ta7 - Ta8;
	       Taj = Ta7 + Ta8;
	       T9r = T9n - T9q;
	       T9w = T9s + T9v;
	       T9G = T9n + T9q;
	       Ta3 = T9f - T9c;
	       T9F = T9v - T9s;
	  }
	  {
	       E Ta, Tq, Tt, T15, Td, Tl, To, T16;
	       {
		    E T8, T9, Tr, Ts;
		    T8 = rio[WS(rs, 1)];
		    T9 = rio[WS(rs, 5)];
		    Ta = T8 + T9;
		    Tq = T8 - T9;
		    Tr = iio[WS(rs, 1)];
		    Ts = iio[WS(rs, 5)];
		    Tt = Tr - Ts;
		    T15 = Tr + Ts;
	       }
	       {
		    E Tb, Tc, Tm, Tn;
		    Tb = rio[WS(rs, 7)];
		    Tc = rio[WS(rs, 3)];
		    Td = Tb + Tc;
		    Tl = Tb - Tc;
		    Tm = iio[WS(rs, 7)];
		    Tn = iio[WS(rs, 3)];
		    To = Tm - Tn;
		    T16 = Tm + Tn;
	       }
	       Te = Ta + Td;
	       T17 = T15 - T16;
	       T1h = T15 + T16;
	       Tp = Tl - To;
	       Tu = Tq + Tt;
	       TE = Tl + To;
	       T11 = Td - Ta;
	       TD = Tt - Tq;
	  }
	  {
	       E T1l, T1Q, T1B, T2g, T1o, T1y, T1T, T2h;
	       {
		    E T1j, T1k, T1z, T1A;
		    T1j = rio[WS(vs, 1)];
		    T1k = rio[WS(vs, 1) + WS(rs, 4)];
		    T1l = T1j + T1k;
		    T1Q = T1j - T1k;
		    T1z = iio[WS(vs, 1)];
		    T1A = iio[WS(vs, 1) + WS(rs, 4)];
		    T1B = T1z - T1A;
		    T2g = T1z + T1A;
	       }
	       {
		    E T1m, T1n, T1R, T1S;
		    T1m = rio[WS(vs, 1) + WS(rs, 2)];
		    T1n = rio[WS(vs, 1) + WS(rs, 6)];
		    T1o = T1m + T1n;
		    T1y = T1m - T1n;
		    T1R = iio[WS(vs, 1) + WS(rs, 2)];
		    T1S = iio[WS(vs, 1) + WS(rs, 6)];
		    T1T = T1R - T1S;
		    T2h = T1R + T1S;
	       }
	       T1p = T1l + T1o;
	       T2m = T1l - T1o;
	       T2y = T2g + T2h;
	       T1C = T1y + T1B;
	       T1U = T1Q - T1T;
	       T28 = T1B - T1y;
	       T2i = T2g - T2h;
	       T24 = T1Q + T1T;
	  }
	  {
	       E T7U, T8a, T8d, T8P, T7X, T85, T88, T8Q;
	       {
		    E T7S, T7T, T8b, T8c;
		    T7S = rio[WS(vs, 6) + WS(rs, 1)];
		    T7T = rio[WS(vs, 6) + WS(rs, 5)];
		    T7U = T7S + T7T;
		    T8a = T7S - T7T;
		    T8b = iio[WS(vs, 6) + WS(rs, 1)];
		    T8c = iio[WS(vs, 6) + WS(rs, 5)];
		    T8d = T8b - T8c;
		    T8P = T8b + T8c;
	       }
	       {
		    E T7V, T7W, T86, T87;
		    T7V = rio[WS(vs, 6) + WS(rs, 7)];
		    T7W = rio[WS(vs, 6) + WS(rs, 3)];
		    T7X = T7V + T7W;
		    T85 = T7V - T7W;
		    T86 = iio[WS(vs, 6) + WS(rs, 7)];
		    T87 = iio[WS(vs, 6) + WS(rs, 3)];
		    T88 = T86 - T87;
		    T8Q = T86 + T87;
	       }
	       T7Y = T7U + T7X;
	       T8R = T8P - T8Q;
	       T91 = T8P + T8Q;
	       T89 = T85 - T88;
	       T8e = T8a + T8d;
	       T8o = T85 + T88;
	       T8L = T7X - T7U;
	       T8n = T8d - T8a;
	  }
	  {
	       E T95, T9A, T9l, Ta0, T98, T9i, T9D, Ta1;
	       {
		    E T93, T94, T9j, T9k;
		    T93 = rio[WS(vs, 7)];
		    T94 = rio[WS(vs, 7) + WS(rs, 4)];
		    T95 = T93 + T94;
		    T9A = T93 - T94;
		    T9j = iio[WS(vs, 7)];
		    T9k = iio[WS(vs, 7) + WS(rs, 4)];
		    T9l = T9j - T9k;
		    Ta0 = T9j + T9k;
	       }
	       {
		    E T96, T97, T9B, T9C;
		    T96 = rio[WS(vs, 7) + WS(rs, 2)];
		    T97 = rio[WS(vs, 7) + WS(rs, 6)];
		    T98 = T96 + T97;
		    T9i = T96 - T97;
		    T9B = iio[WS(vs, 7) + WS(rs, 2)];
		    T9C = iio[WS(vs, 7) + WS(rs, 6)];
		    T9D = T9B - T9C;
		    Ta1 = T9B + T9C;
	       }
	       T99 = T95 + T98;
	       Ta6 = T95 - T98;
	       Tai = Ta0 + Ta1;
	       T9m = T9i + T9l;
	       T9E = T9A - T9D;
	       T9S = T9l - T9i;
	       Ta2 = Ta0 - Ta1;
	       T9O = T9A + T9D;
	  }
	  {
	       E T2D, T38, T2T, T3y, T2G, T2Q, T3b, T3z;
	       {
		    E T2B, T2C, T2R, T2S;
		    T2B = rio[WS(vs, 2)];
		    T2C = rio[WS(vs, 2) + WS(rs, 4)];
		    T2D = T2B + T2C;
		    T38 = T2B - T2C;
		    T2R = iio[WS(vs, 2)];
		    T2S = iio[WS(vs, 2) + WS(rs, 4)];
		    T2T = T2R - T2S;
		    T3y = T2R + T2S;
	       }
	       {
		    E T2E, T2F, T39, T3a;
		    T2E = rio[WS(vs, 2) + WS(rs, 2)];
		    T2F = rio[WS(vs, 2) + WS(rs, 6)];
		    T2G = T2E + T2F;
		    T2Q = T2E - T2F;
		    T39 = iio[WS(vs, 2) + WS(rs, 2)];
		    T3a = iio[WS(vs, 2) + WS(rs, 6)];
		    T3b = T39 - T3a;
		    T3z = T39 + T3a;
	       }
	       T2H = T2D + T2G;
	       T3E = T2D - T2G;
	       T3Q = T3y + T3z;
	       T2U = T2Q + T2T;
	       T3c = T38 - T3b;
	       T3q = T2T - T2Q;
	       T3A = T3y - T3z;
	       T3m = T38 + T3b;
	  }
	  {
	       E T42, T4i, T4l, T4X, T45, T4d, T4g, T4Y;
	       {
		    E T40, T41, T4j, T4k;
		    T40 = rio[WS(vs, 3) + WS(rs, 1)];
		    T41 = rio[WS(vs, 3) + WS(rs, 5)];
		    T42 = T40 + T41;
		    T4i = T40 - T41;
		    T4j = iio[WS(vs, 3) + WS(rs, 1)];
		    T4k = iio[WS(vs, 3) + WS(rs, 5)];
		    T4l = T4j - T4k;
		    T4X = T4j + T4k;
	       }
	       {
		    E T43, T44, T4e, T4f;
		    T43 = rio[WS(vs, 3) + WS(rs, 7)];
		    T44 = rio[WS(vs, 3) + WS(rs, 3)];
		    T45 = T43 + T44;
		    T4d = T43 - T44;
		    T4e = iio[WS(vs, 3) + WS(rs, 7)];
		    T4f = iio[WS(vs, 3) + WS(rs, 3)];
		    T4g = T4e - T4f;
		    T4Y = T4e + T4f;
	       }
	       T46 = T42 + T45;
	       T4Z = T4X - T4Y;
	       T59 = T4X + T4Y;
	       T4h = T4d - T4g;
	       T4m = T4i + T4l;
	       T4w = T4d + T4g;
	       T4T = T45 - T42;
	       T4v = T4l - T4i;
	  }
	  {
	       E T5d, T5I, T5t, T68, T5g, T5q, T5L, T69;
	       {
		    E T5b, T5c, T5r, T5s;
		    T5b = rio[WS(vs, 4)];
		    T5c = rio[WS(vs, 4) + WS(rs, 4)];
		    T5d = T5b + T5c;
		    T5I = T5b - T5c;
		    T5r = iio[WS(vs, 4)];
		    T5s = iio[WS(vs, 4) + WS(rs, 4)];
		    T5t = T5r - T5s;
		    T68 = T5r + T5s;
	       }
	       {
		    E T5e, T5f, T5J, T5K;
		    T5e = rio[WS(vs, 4) + WS(rs, 2)];
		    T5f = rio[WS(vs, 4) + WS(rs, 6)];
		    T5g = T5e + T5f;
		    T5q = T5e - T5f;
		    T5J = iio[WS(vs, 4) + WS(rs, 2)];
		    T5K = iio[WS(vs, 4) + WS(rs, 6)];
		    T5L = T5J - T5K;
		    T69 = T5J + T5K;
	       }
	       T5h = T5d + T5g;
	       T6e = T5d - T5g;
	       T6q = T68 + T69;
	       T5u = T5q + T5t;
	       T5M = T5I - T5L;
	       T60 = T5t - T5q;
	       T6a = T68 - T69;
	       T5W = T5I + T5L;
	  }
	  {
	       E T6C, T6S, T6V, T7x, T6F, T6N, T6Q, T7y;
	       {
		    E T6A, T6B, T6T, T6U;
		    T6A = rio[WS(vs, 5) + WS(rs, 1)];
		    T6B = rio[WS(vs, 5) + WS(rs, 5)];
		    T6C = T6A + T6B;
		    T6S = T6A - T6B;
		    T6T = iio[WS(vs, 5) + WS(rs, 1)];
		    T6U = iio[WS(vs, 5) + WS(rs, 5)];
		    T6V = T6T - T6U;
		    T7x = T6T + T6U;
	       }
	       {
		    E T6D, T6E, T6O, T6P;
		    T6D = rio[WS(vs, 5) + WS(rs, 7)];
		    T6E = rio[WS(vs, 5) + WS(rs, 3)];
		    T6F = T6D + T6E;
		    T6N = T6D - T6E;
		    T6O = iio[WS(vs, 5) + WS(rs, 7)];
		    T6P = iio[WS(vs, 5) + WS(rs, 3)];
		    T6Q = T6O - T6P;
		    T7y = T6O + T6P;
	       }
	       T6G = T6C + T6F;
	       T7z = T7x - T7y;
	       T7J = T7x + T7y;
	       T6R = T6N - T6Q;
	       T6W = T6S + T6V;
	       T76 = T6N + T6Q;
	       T7t = T6F - T6C;
	       T75 = T6V - T6S;
	  }
	  {
	       E T2K, T30, T33, T3F, T2N, T2V, T2Y, T3G;
	       {
		    E T2I, T2J, T31, T32;
		    T2I = rio[WS(vs, 2) + WS(rs, 1)];
		    T2J = rio[WS(vs, 2) + WS(rs, 5)];
		    T2K = T2I + T2J;
		    T30 = T2I - T2J;
		    T31 = iio[WS(vs, 2) + WS(rs, 1)];
		    T32 = iio[WS(vs, 2) + WS(rs, 5)];
		    T33 = T31 - T32;
		    T3F = T31 + T32;
	       }
	       {
		    E T2L, T2M, T2W, T2X;
		    T2L = rio[WS(vs, 2) + WS(rs, 7)];
		    T2M = rio[WS(vs, 2) + WS(rs, 3)];
		    T2N = T2L + T2M;
		    T2V = T2L - T2M;
		    T2W = iio[WS(vs, 2) + WS(rs, 7)];
		    T2X = iio[WS(vs, 2) + WS(rs, 3)];
		    T2Y = T2W - T2X;
		    T3G = T2W + T2X;
	       }
	       T2O = T2K + T2N;
	       T3H = T3F - T3G;
	       T3R = T3F + T3G;
	       T2Z = T2V - T2Y;
	       T34 = T30 + T33;
	       T3e = T2V + T2Y;
	       T3B = T2N - T2K;
	       T3d = T33 - T30;
	  }
	  {
	       E T3V, T4q, T4b, T4Q, T3Y, T48, T4t, T4R;
	       {
		    E T3T, T3U, T49, T4a;
		    T3T = rio[WS(vs, 3)];
		    T3U = rio[WS(vs, 3) + WS(rs, 4)];
		    T3V = T3T + T3U;
		    T4q = T3T - T3U;
		    T49 = iio[WS(vs, 3)];
		    T4a = iio[WS(vs, 3) + WS(rs, 4)];
		    T4b = T49 - T4a;
		    T4Q = T49 + T4a;
	       }
	       {
		    E T3W, T3X, T4r, T4s;
		    T3W = rio[WS(vs, 3) + WS(rs, 2)];
		    T3X = rio[WS(vs, 3) + WS(rs, 6)];
		    T3Y = T3W + T3X;
		    T48 = T3W - T3X;
		    T4r = iio[WS(vs, 3) + WS(rs, 2)];
		    T4s = iio[WS(vs, 3) + WS(rs, 6)];
		    T4t = T4r - T4s;
		    T4R = T4r + T4s;
	       }
	       T3Z = T3V + T3Y;
	       T4W = T3V - T3Y;
	       T58 = T4Q + T4R;
	       T4c = T48 + T4b;
	       T4u = T4q - T4t;
	       T4I = T4b - T48;
	       T4S = T4Q - T4R;
	       T4E = T4q + T4t;
	  }
	  {
	       E T5k, T5A, T5D, T6f, T5n, T5v, T5y, T6g;
	       {
		    E T5i, T5j, T5B, T5C;
		    T5i = rio[WS(vs, 4) + WS(rs, 1)];
		    T5j = rio[WS(vs, 4) + WS(rs, 5)];
		    T5k = T5i + T5j;
		    T5A = T5i - T5j;
		    T5B = iio[WS(vs, 4) + WS(rs, 1)];
		    T5C = iio[WS(vs, 4) + WS(rs, 5)];
		    T5D = T5B - T5C;
		    T6f = T5B + T5C;
	       }
	       {
		    E T5l, T5m, T5w, T5x;
		    T5l = rio[WS(vs, 4) + WS(rs, 7)];
		    T5m = rio[WS(vs, 4) + WS(rs, 3)];
		    T5n = T5l + T5m;
		    T5v = T5l - T5m;
		    T5w = iio[WS(vs, 4) + WS(rs, 7)];
		    T5x = iio[WS(vs, 4) + WS(rs, 3)];
		    T5y = T5w - T5x;
		    T6g = T5w + T5x;
	       }
	       T5o = T5k + T5n;
	       T6h = T6f - T6g;
	       T6r = T6f + T6g;
	       T5z = T5v - T5y;
	       T5E = T5A + T5D;
	       T5O = T5v + T5y;
	       T6b = T5n - T5k;
	       T5N = T5D - T5A;
	  }
	  {
	       E T6v, T70, T6L, T7q, T6y, T6I, T73, T7r;
	       {
		    E T6t, T6u, T6J, T6K;
		    T6t = rio[WS(vs, 5)];
		    T6u = rio[WS(vs, 5) + WS(rs, 4)];
		    T6v = T6t + T6u;
		    T70 = T6t - T6u;
		    T6J = iio[WS(vs, 5)];
		    T6K = iio[WS(vs, 5) + WS(rs, 4)];
		    T6L = T6J - T6K;
		    T7q = T6J + T6K;
	       }
	       {
		    E T6w, T6x, T71, T72;
		    T6w = rio[WS(vs, 5) + WS(rs, 2)];
		    T6x = rio[WS(vs, 5) + WS(rs, 6)];
		    T6y = T6w + T6x;
		    T6I = T6w - T6x;
		    T71 = iio[WS(vs, 5) + WS(rs, 2)];
		    T72 = iio[WS(vs, 5) + WS(rs, 6)];
		    T73 = T71 - T72;
		    T7r = T71 + T72;
	       }
	       T6z = T6v + T6y;
	       T7w = T6v - T6y;
	       T7I = T7q + T7r;
	       T6M = T6I + T6L;
	       T74 = T70 - T73;
	       T7i = T6L - T6I;
	       T7s = T7q - T7r;
	       T7e = T70 + T73;
	  }
	  rio[0] = T7 + Te;
	  iio[0] = T1g + T1h;
	  rio[WS(rs, 1)] = T1p + T1w;
	  iio[WS(rs, 1)] = T2y + T2z;
	  rio[WS(rs, 3)] = T3Z + T46;
	  rio[WS(rs, 2)] = T2H + T2O;
	  iio[WS(rs, 2)] = T3Q + T3R;
	  iio[WS(rs, 3)] = T58 + T59;
	  rio[WS(rs, 6)] = T7R + T7Y;
	  iio[WS(rs, 6)] = T90 + T91;
	  iio[WS(rs, 5)] = T7I + T7J;
	  rio[WS(rs, 5)] = T6z + T6G;
	  iio[WS(rs, 4)] = T6q + T6r;
	  rio[WS(rs, 4)] = T5h + T5o;
	  rio[WS(rs, 7)] = T99 + T9g;
	  iio[WS(rs, 7)] = Tai + Taj;
	  {
	       E T12, T18, TX, T13;
	       T12 = T10 - T11;
	       T18 = T14 - T17;
	       TX = W[10];
	       T13 = W[11];
	       iio[WS(vs, 6)] = FNMS(T13, T18, TX * T12);
	       rio[WS(vs, 6)] = FMA(T13, T12, TX * T18);
	  }
	  {