r2cf_64.c

快速fft变换

	       E Ti, T16, Ts, T1a, Tl, T17, Tp, T19, T4m, T4n;
	       {
		    E Tg, Th, Tq, Tr;
		    Tg = R0[WS(rs, 2)];
		    Th = R0[WS(rs, 18)];
		    Ti = Tg + Th;
		    T16 = Tg - Th;
		    Tq = R0[WS(rs, 6)];
		    Tr = R0[WS(rs, 22)];
		    Ts = Tq + Tr;
		    T1a = Tq - Tr;
	       }
	       {
		    E Tj, Tk, Tn, To;
		    Tj = R0[WS(rs, 10)];
		    Tk = R0[WS(rs, 26)];
		    Tl = Tj + Tk;
		    T17 = Tj - Tk;
		    Tn = R0[WS(rs, 30)];
		    To = R0[WS(rs, 14)];
		    Tp = Tn + To;
		    T19 = Tn - To;
	       }
	       Tm = Ti + Tl;
	       Tt = Tp + Ts;
	       Tu = Tm + Tt;
	       T4m = Ti - Tl;
	       T4n = Tp - Ts;
	       T4o = KP707106781 * (T4m + T4n);
	       T5b = KP707106781 * (T4n - T4m);
	       {
		    E T18, T1b, T2O, T2P;
		    T18 = FNMS(KP382683432, T17, KP923879532 * T16);
		    T1b = FMA(KP923879532, T19, KP382683432 * T1a);
		    T1c = T18 + T1b;
		    T3R = T1b - T18;
		    T2O = FNMS(KP923879532, T1a, KP382683432 * T19);
		    T2P = FMA(KP382683432, T16, KP923879532 * T17);
		    T2Q = T2O - T2P;
		    T3o = T2P + T2O;
	       }
	  }
	  {
	       E T1A, T4E, T1K, T4C, T1D, T4F, T1H, T4B;
	       {
		    E T1y, T1z, T1I, T1J;
		    T1y = R1[WS(rs, 30)];
		    T1z = R1[WS(rs, 14)];
		    T1A = T1y - T1z;
		    T4E = T1y + T1z;
		    T1I = R1[WS(rs, 10)];
		    T1J = R1[WS(rs, 26)];
		    T1K = T1I - T1J;
		    T4C = T1I + T1J;
	       }
	       {
		    E T1B, T1C, T1F, T1G;
		    T1B = R1[WS(rs, 6)];
		    T1C = R1[WS(rs, 22)];
		    T1D = T1B - T1C;
		    T4F = T1B + T1C;
		    T1F = R1[WS(rs, 2)];
		    T1G = R1[WS(rs, 18)];
		    T1H = T1F - T1G;
		    T4B = T1F + T1G;
	       }
	       {
		    E T1E, T1L, T5J, T5K;
		    T1E = FNMS(KP923879532, T1D, KP382683432 * T1A);
		    T1L = FMA(KP382683432, T1H, KP923879532 * T1K);
		    T1M = T1E - T1L;
		    T3z = T1L + T1E;
		    T5J = T4B + T4C;
		    T5K = T4E + T4F;
		    T5L = T5J + T5K;
		    T67 = T5K - T5J;
	       }
	       {
		    E T24, T25, T4D, T4G;
		    T24 = FNMS(KP382683432, T1K, KP923879532 * T1H);
		    T25 = FMA(KP923879532, T1A, KP382683432 * T1D);
		    T26 = T24 + T25;
		    T3C = T25 - T24;
		    T4D = T4B - T4C;
		    T4G = T4E - T4F;
		    T4H = KP707106781 * (T4D + T4G);
		    T4M = KP707106781 * (T4G - T4D);
	       }
	  }
	  {
	       E T2m, T4S, T2w, T4W, T2p, T4T, T2t, T4V;
	       {
		    E T2k, T2l, T2u, T2v;
		    T2k = R1[WS(rs, 1)];
		    T2l = R1[WS(rs, 17)];
		    T2m = T2k - T2l;
		    T4S = T2k + T2l;
		    T2u = R1[WS(rs, 5)];
		    T2v = R1[WS(rs, 21)];
		    T2w = T2u - T2v;
		    T4W = T2u + T2v;
	       }
	       {
		    E T2n, T2o, T2r, T2s;
		    T2n = R1[WS(rs, 9)];
		    T2o = R1[WS(rs, 25)];
		    T2p = T2n - T2o;
		    T4T = T2n + T2o;
		    T2r = R1[WS(rs, 29)];
		    T2s = R1[WS(rs, 13)];
		    T2t = T2r - T2s;
		    T4V = T2r + T2s;
	       }
	       {
		    E T2q, T2x, T5Q, T5R;
		    T2q = FNMS(KP382683432, T2p, KP923879532 * T2m);
		    T2x = FMA(KP923879532, T2t, KP382683432 * T2w);
		    T2y = T2q + T2x;
		    T3J = T2x - T2q;
		    T5Q = T4S + T4T;
		    T5R = T4V + T4W;
		    T5S = T5Q + T5R;
		    T6a = T5R - T5Q;
	       }
	       {
		    E T2A, T2B, T4U, T4X;
		    T2A = FNMS(KP923879532, T2w, KP382683432 * T2t);
		    T2B = FMA(KP382683432, T2m, KP923879532 * T2p);
		    T2C = T2A - T2B;
		    T3G = T2B + T2A;
		    T4U = T4S - T4T;
		    T4X = T4V - T4W;
		    T4Y = KP707106781 * (T4U + T4X);
		    T53 = KP707106781 * (T4X - T4U);
	       }
	  }
	  {
	       E Tv, T10, T5X, T5Y, T5Z, T60;
	       Tv = Tf + Tu;
	       T10 = TK + TZ;
	       T5X = Tv + T10;
	       T5Y = T5I + T5L;
	       T5Z = T5P + T5S;
	       T60 = T5Y + T5Z;
	       Cr[WS(csr, 16)] = Tv - T10;
	       Ci[WS(csi, 16)] = T5Z - T5Y;
	       Cr[WS(csr, 32)] = T5X - T60;
	       Cr[0] = T5X + T60;
	  }
	  {
	       E T5F, T5V, T5U, T5W, T5M, T5T;
	       T5F = Tf - Tu;
	       T5V = TZ - TK;
	       T5M = T5I - T5L;
	       T5T = T5P - T5S;
	       T5U = KP707106781 * (T5M + T5T);
	       T5W = KP707106781 * (T5T - T5M);
	       Cr[WS(csr, 24)] = T5F - T5U;
	       Ci[WS(csi, 24)] = T5W - T5V;
	       Cr[WS(csr, 8)] = T5F + T5U;
	       Ci[WS(csi, 8)] = T5V + T5W;
	  }
	  {
	       E T65, T6l, T6k, T6m, T6c, T6g, T6f, T6h;
	       {
		    E T61, T64, T6i, T6j;
		    T61 = T7 - Te;
		    T64 = KP707106781 * (T62 + T63);
		    T65 = T61 + T64;
		    T6l = T61 - T64;
		    T6i = FNMS(KP382683432, T66, KP923879532 * T67);
		    T6j = FMA(KP382683432, T69, KP923879532 * T6a);
		    T6k = T6i + T6j;
		    T6m = T6j - T6i;
	       }
	       {
		    E T68, T6b, T6d, T6e;
		    T68 = FMA(KP923879532, T66, KP382683432 * T67);
		    T6b = FNMS(KP382683432, T6a, KP923879532 * T69);
		    T6c = T68 + T6b;
		    T6g = T6b - T68;
		    T6d = KP707106781 * (T63 - T62);
		    T6e = Tt - Tm;
		    T6f = T6d - T6e;
		    T6h = T6e + T6d;
	       }
	       Cr[WS(csr, 28)] = T65 - T6c;
	       Ci[WS(csi, 28)] = T6k - T6h;
	       Cr[WS(csr, 4)] = T65 + T6c;
	       Ci[WS(csi, 4)] = T6h + T6k;
	       Ci[WS(csi, 12)] = T6f + T6g;
	       Cr[WS(csr, 12)] = T6l + T6m;
	       Ci[WS(csi, 20)] = T6g - T6f;
	       Cr[WS(csr, 20)] = T6l - T6m;
	  }
	  {
	       E T5n, T5D, T5x, T5z, T5q, T5A, T5t, T5B;
	       {
		    E T5l, T5m, T5v, T5w;
		    T5l = T4l - T4o;
		    T5m = T58 - T57;
		    T5n = T5l + T5m;
		    T5D = T5l - T5m;
		    T5v = T4v - T4s;
		    T5w = T5b - T5a;
		    T5x = T5v - T5w;
		    T5z = T5w + T5v;
	       }
	       {
		    E T5o, T5p, T5r, T5s;
		    T5o = T4A - T4H;
		    T5p = T4M - T4L;
		    T5q = FMA(KP831469612, T5o, KP555570233 * T5p);
		    T5A = FNMS(KP555570233, T5o, KP831469612 * T5p);
		    T5r = T4R - T4Y;
		    T5s = T53 - T52;
		    T5t = FNMS(KP555570233, T5s, KP831469612 * T5r);
		    T5B = FMA(KP555570233, T5r, KP831469612 * T5s);
	       }
	       {
		    E T5u, T5C, T5y, T5E;
		    T5u = T5q + T5t;
		    Cr[WS(csr, 26)] = T5n - T5u;
		    Cr[WS(csr, 6)] = T5n + T5u;
		    T5C = T5A + T5B;
		    Ci[WS(csi, 6)] = T5z + T5C;
		    Ci[WS(csi, 26)] = T5C - T5z;
		    T5y = T5t - T5q;
		    Ci[WS(csi, 10)] = T5x + T5y;
		    Ci[WS(csi, 22)] = T5y - T5x;
		    T5E = T5B - T5A;
		    Cr[WS(csr, 22)] = T5D - T5E;
		    Cr[WS(csr, 10)] = T5D + T5E;
	       }
	  }
	  {
	       E T4x, T5j, T5d, T5f, T4O, T5g, T55, T5h;
	       {
		    E T4p, T4w, T59, T5c;
		    T4p = T4l + T4o;
		    T4w = T4s + T4v;
		    T4x = T4p + T4w;
		    T5j = T4p - T4w;
		    T59 = T57 + T58;
		    T5c = T5a + T5b;
		    T5d = T59 - T5c;
		    T5f = T5c + T59;
	       }
	       {
		    E T4I, T4N, T4Z, T54;
		    T4I = T4A + T4H;
		    T4N = T4L + T4M;
		    T4O = FMA(KP980785280, T4I, KP195090322 * T4N);
		    T5g = FNMS(KP195090322, T4I, KP980785280 * T4N);
		    T4Z = T4R + T4Y;
		    T54 = T52 + T53;
		    T55 = FNMS(KP195090322, T54, KP980785280 * T4Z);
		    T5h = FMA(KP195090322, T4Z, KP980785280 * T54);
	       }
	       {
		    E T56, T5i, T5e, T5k;
		    T56 = T4O + T55;
		    Cr[WS(csr, 30)] = T4x - T56;
		    Cr[WS(csr, 2)] = T4x + T56;
		    T5i = T5g + T5h;
		    Ci[WS(csi, 2)] = T5f + T5i;
		    Ci[WS(csi, 30)] = T5i - T5f;
		    T5e = T55 - T4O;
		    Ci[WS(csi, 14)] = T5d + T5e;
		    Ci[WS(csi, 18)] = T5e - T5d;
		    T5k = T5h - T5g;
		    Cr[WS(csr, 18)] = T5j - T5k;
		    Cr[WS(csr, 14)] = T5j + T5k;
	       }
	  }
	  {
	       E T3p, T41, T4c, T3S, T3w, T4b, T49, T4h, T3P, T42, T3E, T3W, T46, T4g, T3L;
	       E T3X;
	       {
		    E T3s, T3v, T3A, T3D;
		    T3p = T3n + T3o;
		    T41 = T3n - T3o;
		    T4c = T3R - T3Q;
		    T3S = T3Q + T3R;
		    T3s = FMA(KP831469612, T3q, KP555570233 * T3r);
		    T3v = FNMS(KP555570233, T3u, KP831469612 * T3t);
		    T3w = T3s + T3v;
		    T4b = T3v - T3s;
		    {
			 E T47, T48, T3N, T3O;
			 T47 = T3F - T3G;
			 T48 = T3J - T3I;
			 T49 = FNMS(KP471396736, T48, KP881921264 * T47);
			 T4h = FMA(KP471396736, T47, KP881921264 * T48);
			 T3N = FNMS(KP555570233, T3q, KP831469612 * T3r);
			 T3O = FMA(KP555570233, T3t, KP831469612 * T3u);
			 T3P = T3N + T3O;
			 T42 = T3O - T3N;
		    }
		    T3A = T3y + T3z;
		    T3D = T3B + T3C;
		    T3E = FMA(KP956940335, T3A, KP290284677 * T3D);
		    T3W = FNMS(KP290284677, T3A, KP956940335 * T3D);
		    {
			 E T44, T45, T3H, T3K;
			 T44 = T3y - T3z;
			 T45 = T3C - T3B;
			 T46 = FMA(KP881921264, T44, KP471396736 * T45);
			 T4g = FNMS(KP471396736, T44, KP881921264 * T45);
			 T3H = T3F + T3G;
			 T3K = T3I + T3J;
			 T3L = FNMS(KP290284677, T3K, KP956940335 * T3H);
			 T3X = FMA(KP290284677, T3H, KP956940335 * T3K);
		    }
	       }
	       {
		    E T3x, T3M, T3V, T3Y;
		    T3x = T3p + T3w;
		    T3M = T3E + T3L;
		    Cr[WS(csr, 29)] = T3x - T3M;
		    Cr[WS(csr, 3)] = T3x + T3M;
		    T3V = T3S + T3P;
		    T3Y = T3W + T3X;
		    Ci[WS(csi, 3)] = T3V + T3Y;
		    Ci[WS(csi, 29)] = T3Y - T3V;
	       }
	       {
		    E T3T, T3U, T3Z, T40;
		    T3T = T3P - T3S;
		    T3U = T3L - T3E;
		    Ci[WS(csi, 13)] = T3T + T3U;
		    Ci[WS(csi, 19)] = T3U - T3T;
		    T3Z = T3p - T3w;
		    T40 = T3X - T3W;
		    Cr[WS(csr, 19)] = T3Z - T40;
		    Cr[WS(csr, 13)] = T3Z + T40;
	       }
	       {
		    E T43, T4a, T4f, T4i;
		    T43 = T41 + T42;
		    T4a = T46 + T49;
		    Cr[WS(csr, 27)] = T43 - T4a;
		    Cr[WS(csr, 5)] = T43 + T4a;
		    T4f = T4c + T4b;
		    T4i = T4g + T4h;
		    Ci[WS(csi, 5)] = T4f + T4i;
		    Ci[WS(csi, 27)] = T4i - T4f;
	       }
	       {
		    E T4d, T4e, T4j, T4k;
		    T4d = T4b - T4c;
		    T4e = T49 - T46;
		    Ci[WS(csi, 11)] = T4d + T4e;
		    Ci[WS(csi, 21)] = T4e - T4d;
		    T4j = T41 - T42;
		    T4k = T4h - T4g;
		    Cr[WS(csr, 21)] = T4j - T4k;
		    Cr[WS(csr, 11)] = T4j + T4k;
	       }
	  }
	  {
	       E T1d, T33, T3e, T2U, T1w, T3d, T3b, T3j, T2N, T34, T28, T2Y, T38, T3i, T2J;
	       E T2Z;
	       {
		    E T1m, T1v, T1Y, T27;
		    T1d = T15 - T1c;
		    T33 = T15 + T1c;
		    T3e = T2T + T2Q;
		    T2U = T2Q - T2T;
		    T1m = FMA(KP195090322, T1i, KP980785280 * T1l);
		    T1v = FNMS(KP195090322, T1u, KP980785280 * T1r);
		    T1w = T1m - T1v;
		    T3d = T1v + T1m;
		    {
			 E T39, T3a, T2L, T2M;
			 T39 = T2j + T2y;
			 T3a = T2H + T2C;
			 T3b = FNMS(KP098017140, T3a, KP995184726 * T39);
			 T3j = FMA(KP995184726, T3a, KP098017140 * T39);
			 T2L = FNMS(KP195090322, T1l, KP980785280 * T1i);
			 T2M = FMA(KP980785280, T1u, KP195090322 * T1r);
			 T2N = T2L - T2M;
			 T34 = T2M + T2L;
		    }
		    T1Y = T1M - T1X;
		    T27 = T23 - T26;
		    T28 = FMA(KP634393284, T1Y, KP773010453 * T27);
		    T2Y = FNMS(KP634393284, T27, KP773010453 * T1Y);
		    {
			 E T36, T37, T2z, T2I;
			 T36 = T1X + T1M;
			 T37 = T23 + T26;
			 T38 = FMA(KP098017140, T36, KP995184726 * T37);
			 T3i = FNMS(KP098017140, T37, KP995184726 * T36);
			 T2z = T2j - T2y;
			 T2I = T2C - T2H;
			 T2J = FNMS(KP634393284, T2I, KP773010453 * T2z);
			 T2Z = FMA(KP773010453, T2I, KP634393284 * T2z);
		    }
	       }
	       {
		    E T1x, T2K, T2X, T30;
		    T1x = T1d + T1w;
		    T2K = T28 + T2J;
		    Cr[WS(csr, 25)] = T1x - T2K;
		    Cr[WS(csr, 7)] = T1x + T2K;
		    T2X = T2U + T2N;
		    T30 = T2Y + T2Z;
		    Ci[WS(csi, 7)] = T2X + T30;
		    Ci[WS(csi, 25)] = T30 - T2X;
	       }
	       {
		    E T2V, T2W, T31, T32;
		    T2V = T2N - T2U;
		    T2W = T2J - T28;
		    Ci[WS(csi, 9)] = T2V + T2W;
		    Ci[WS(csi, 23)] = T2W - T2V;
		    T31 = T1d - T1w;
		    T32 = T2Z - T2Y;
		    Cr[WS(csr, 23)] = T31 - T32;
		    Cr[WS(csr, 9)] = T31 + T32;
	       }
	       {
		    E T35, T3c, T3h, T3k;
		    T35 = T33 + T34;
		    T3c = T38 + T3b;
		    Cr[WS(csr, 31)] = T35 - T3c;
		    Cr[WS(csr, 1)] = T35 + T3c;
		    T3h = T3e + T3d;
		    T3k = T3i + T3j;
		    Ci[WS(csi, 1)] = T3h + T3k;
		    Ci[WS(csi, 31)] = T3k - T3h;
	       }
	       {
		    E T3f, T3g, T3l, T3m;
		    T3f = T3d - T3e;
		    T3g = T3b - T38;
		    Ci[WS(csi, 15)] = T3f + T3g;
		    Ci[WS(csi, 17)] = T3g - T3f;
		    T3l = T33 - T34;
		    T3m = T3j - T3i;
		    Cr[WS(csr, 17)] = T3l - T3m;
		    Cr[WS(csr, 15)] = T3l + T3m;
	       }
	  }
     }
}

static const kr2c_desc desc = { 64, "r2cf_64", {342, 72, 52, 0}, &GENUS };

void X(codelet_r2cf_64) (planner *p) {
     X(kr2c_register) (p, r2cf_64, &desc);
}

#endif				/* HAVE_FMA */