r2cb_64.c

快速fft变换

		    T4b = TG - TH;
		    T1H = Ci[WS(csi, 9)];
		    T1I = Ci[WS(csi, 23)];
		    T1J = T1H - T1I;
		    T4c = T1H + T1I;
	       }
	       {
		    E TJ, TK, T1E, T1F;
		    TJ = Cr[WS(csr, 7)];
		    TK = Cr[WS(csr, 25)];
		    TL = TJ + TK;
		    T4e = TJ - TK;
		    T1E = Ci[WS(csi, 25)];
		    T1F = Ci[WS(csi, 7)];
		    T1G = T1E - T1F;
		    T4f = T1F + T1E;
	       }
	       {
		    E TF, TM, T1D, T1K;
		    T4a = T48 - T49;
		    T5y = T48 + T49;
		    T5N = T4G - T4F;
		    T4H = T4F + T4G;
		    TF = TB + TE;
		    TM = TI + TL;
		    TN = TF + TM;
		    T31 = TF - TM;
		    {
			 E T4C, T4D, T37, T38;
			 T4C = T4b + T4c;
			 T4D = T4e + T4f;
			 T4E = KP707106781 * (T4C - T4D);
			 T5z = KP707106781 * (T4C + T4D);
			 T37 = T2f + T2c;
			 T38 = T1J + T1G;
			 T39 = T37 - T38;
			 T3q = T38 + T37;
		    }
		    T1D = TB - TE;
		    T1K = T1G - T1J;
		    T1L = T1D + T1K;
		    T2B = T1D - T1K;
		    {
			 E T4d, T4g, T29, T2g;
			 T4d = T4b - T4c;
			 T4g = T4e - T4f;
			 T4h = KP707106781 * (T4d + T4g);
			 T5M = KP707106781 * (T4d - T4g);
			 T29 = TI - TL;
			 T2g = T2c - T2f;
			 T2h = T29 + T2g;
			 T2F = T2g - T29;
		    }
	       }
	  }
	  {
	       E TQ, T4j, T1P, T4n, TT, T4m, T1S, T4k, TX, T4q, T1Y, T4u, T10, T4t, T21;
	       E T4r;
	       {
		    E TO, TP, T1N, T1O;
		    TO = Cr[WS(csr, 5)];
		    TP = Cr[WS(csr, 27)];
		    TQ = TO + TP;
		    T4j = TO - TP;
		    T1N = Ci[WS(csi, 5)];
		    T1O = Ci[WS(csi, 27)];
		    T1P = T1N - T1O;
		    T4n = T1N + T1O;
	       }
	       {
		    E TR, TS, T1Q, T1R;
		    TR = Cr[WS(csr, 21)];
		    TS = Cr[WS(csr, 11)];
		    TT = TR + TS;
		    T4m = TR - TS;
		    T1Q = Ci[WS(csi, 21)];
		    T1R = Ci[WS(csi, 11)];
		    T1S = T1Q - T1R;
		    T4k = T1Q + T1R;
	       }
	       {
		    E TV, TW, T1W, T1X;
		    TV = Cr[WS(csr, 3)];
		    TW = Cr[WS(csr, 29)];
		    TX = TV + TW;
		    T4q = TV - TW;
		    T1W = Ci[WS(csi, 29)];
		    T1X = Ci[WS(csi, 3)];
		    T1Y = T1W - T1X;
		    T4u = T1X + T1W;
	       }
	       {
		    E TY, TZ, T1Z, T20;
		    TY = Cr[WS(csr, 13)];
		    TZ = Cr[WS(csr, 19)];
		    T10 = TY + TZ;
		    T4t = TY - TZ;
		    T1Z = Ci[WS(csi, 13)];
		    T20 = Ci[WS(csi, 19)];
		    T21 = T1Z - T20;
		    T4r = T1Z + T20;
	       }
	       {
		    E TU, T11, T5B, T5C;
		    TU = TQ + TT;
		    T11 = TX + T10;
		    T12 = TU + T11;
		    T36 = TU - T11;
		    T5B = T4j + T4k;
		    T5C = T4n - T4m;
		    T5D = FNMS(KP923879532, T5C, KP382683432 * T5B);
		    T5J = FMA(KP923879532, T5B, KP382683432 * T5C);
	       }
	       {
		    E T5E, T5F, T1M, T1T;
		    T5E = T4q + T4r;
		    T5F = T4t + T4u;
		    T5G = FNMS(KP923879532, T5F, KP382683432 * T5E);
		    T5K = FMA(KP923879532, T5E, KP382683432 * T5F);
		    T1M = TQ - TT;
		    T1T = T1P - T1S;
		    T1U = T1M - T1T;
		    T26 = T1M + T1T;
	       }
	       {
		    E T1V, T22, T4l, T4o;
		    T1V = TX - T10;
		    T22 = T1Y - T21;
		    T23 = T1V + T22;
		    T27 = T22 - T1V;
		    T4l = T4j - T4k;
		    T4o = T4m + T4n;
		    T4p = FNMS(KP382683432, T4o, KP923879532 * T4l);
		    T4z = FMA(KP382683432, T4l, KP923879532 * T4o);
	       }
	       {
		    E T4s, T4v, T32, T33;
		    T4s = T4q - T4r;
		    T4v = T4t - T4u;
		    T4w = FMA(KP923879532, T4s, KP382683432 * T4v);
		    T4A = FNMS(KP382683432, T4s, KP923879532 * T4v);
		    T32 = T21 + T1Y;
		    T33 = T1S + T1P;
		    T34 = T32 - T33;
		    T3r = T33 + T32;
	       }
	  }
	  {
	       E T13, T3x, Ty, T3w, Ti, Tx;
	       T13 = KP2_000000000 * (TN + T12);
	       T3x = KP2_000000000 * (T3r + T3q);
	       Ti = Ta + Th;
	       Tx = KP2_000000000 * (Tp + Tw);
	       Ty = Ti + Tx;
	       T3w = Ti - Tx;
	       R0[WS(rs, 16)] = Ty - T13;
	       R0[WS(rs, 24)] = T3w + T3x;
	       R0[0] = Ty + T13;
	       R0[WS(rs, 8)] = T3w - T3x;
	  }
	  {
	       E T3g, T3k, T3j, T3l;
	       {
		    E T3e, T3f, T3h, T3i;
		    T3e = T2S + T2T;
		    T3f = KP1_414213562 * (T2V + T2Y);
		    T3g = T3e - T3f;
		    T3k = T3e + T3f;
		    T3h = T31 - T34;
		    T3i = T39 - T36;
		    T3j = FNMS(KP1_847759065, T3i, KP765366864 * T3h);
		    T3l = FMA(KP1_847759065, T3h, KP765366864 * T3i);
	       }
	       R0[WS(rs, 22)] = T3g - T3j;
	       R0[WS(rs, 30)] = T3k + T3l;
	       R0[WS(rs, 6)] = T3g + T3j;
	       R0[WS(rs, 14)] = T3k - T3l;
	  }
	  {
	       E T3o, T3u, T3t, T3v;
	       {
		    E T3m, T3n, T3p, T3s;
		    T3m = Ta - Th;
		    T3n = KP2_000000000 * (T2X + T2W);
		    T3o = T3m - T3n;
		    T3u = T3m + T3n;
		    T3p = TN - T12;
		    T3s = T3q - T3r;
		    T3t = KP1_414213562 * (T3p - T3s);
		    T3v = KP1_414213562 * (T3p + T3s);
	       }
	       R0[WS(rs, 20)] = T3o - T3t;
	       R0[WS(rs, 28)] = T3u + T3v;
	       R0[WS(rs, 4)] = T3o + T3t;
	       R0[WS(rs, 12)] = T3u - T3v;
	  }
	  {
	       E T30, T3c, T3b, T3d;
	       {
		    E T2U, T2Z, T35, T3a;
		    T2U = T2S - T2T;
		    T2Z = KP1_414213562 * (T2V - T2Y);
		    T30 = T2U + T2Z;
		    T3c = T2U - T2Z;
		    T35 = T31 + T34;
		    T3a = T36 + T39;
		    T3b = FNMS(KP765366864, T3a, KP1_847759065 * T35);
		    T3d = FMA(KP765366864, T35, KP1_847759065 * T3a);
	       }
	       R0[WS(rs, 18)] = T30 - T3b;
	       R0[WS(rs, 26)] = T3c + T3d;
	       R0[WS(rs, 2)] = T30 + T3b;
	       R0[WS(rs, 10)] = T3c - T3d;
	  }
	  {
	       E T25, T2p, T2i, T2q, T1C, T2k, T2o, T2s, T24, T28;
	       T24 = KP707106781 * (T1U + T23);
	       T25 = T1L + T24;
	       T2p = T1L - T24;
	       T28 = KP707106781 * (T26 + T27);
	       T2i = T28 + T2h;
	       T2q = T2h - T28;
	       {
		    E T1i, T1B, T2m, T2n;
		    T1i = T18 + T1h;
		    T1B = FNMS(KP765366864, T1A, KP1_847759065 * T1r);
		    T1C = T1i + T1B;
		    T2k = T1i - T1B;
		    T2m = T18 - T1h;
		    T2n = FMA(KP765366864, T1r, KP1_847759065 * T1A);
		    T2o = T2m - T2n;
		    T2s = T2m + T2n;
	       }
	       {
		    E T2j, T2t, T2l, T2r;
		    T2j = FNMS(KP390180644, T2i, KP1_961570560 * T25);
		    R0[WS(rs, 17)] = T1C - T2j;
		    R0[WS(rs, 1)] = T1C + T2j;
		    T2t = FMA(KP1_662939224, T2p, KP1_111140466 * T2q);
		    R0[WS(rs, 13)] = T2s - T2t;
		    R0[WS(rs, 29)] = T2s + T2t;
		    T2l = FMA(KP390180644, T25, KP1_961570560 * T2i);
		    R0[WS(rs, 9)] = T2k - T2l;
		    R0[WS(rs, 25)] = T2k + T2l;
		    T2r = FNMS(KP1_662939224, T2q, KP1_111140466 * T2p);
		    R0[WS(rs, 21)] = T2o - T2r;
		    R0[WS(rs, 5)] = T2o + T2r;
	       }
	  }
	  {
	       E T2D, T2N, T2G, T2O, T2A, T2I, T2M, T2Q, T2C, T2E;
	       T2C = KP707106781 * (T27 - T26);
	       T2D = T2B + T2C;
	       T2N = T2B - T2C;
	       T2E = KP707106781 * (T1U - T23);
	       T2G = T2E + T2F;
	       T2O = T2F - T2E;
	       {
		    E T2w, T2z, T2K, T2L;
		    T2w = T2u - T2v;
		    T2z = FNMS(KP1_847759065, T2y, KP765366864 * T2x);
		    T2A = T2w + T2z;
		    T2I = T2w - T2z;
		    T2K = T2u + T2v;
		    T2L = FMA(KP1_847759065, T2x, KP765366864 * T2y);
		    T2M = T2K - T2L;
		    T2Q = T2K + T2L;
	       }
	       {
		    E T2H, T2R, T2J, T2P;
		    T2H = FNMS(KP1_111140466, T2G, KP1_662939224 * T2D);
		    R0[WS(rs, 19)] = T2A - T2H;
		    R0[WS(rs, 3)] = T2A + T2H;
		    T2R = FMA(KP1_961570560, T2N, KP390180644 * T2O);
		    R0[WS(rs, 15)] = T2Q - T2R;
		    R0[WS(rs, 31)] = T2Q + T2R;
		    T2J = FMA(KP1_111140466, T2D, KP1_662939224 * T2G);
		    R0[WS(rs, 11)] = T2I - T2J;
		    R0[WS(rs, 27)] = T2I + T2J;
		    T2P = FNMS(KP1_961570560, T2O, KP390180644 * T2N);
		    R0[WS(rs, 23)] = T2M - T2P;
		    R0[WS(rs, 7)] = T2M + T2P;
	       }
	  }
	  {
	       E T5p, T5T, T5w, T5U, T5I, T5W, T5P, T5X, T5s, T5v;
	       T5p = T5l + T5o;
	       T5T = T5l - T5o;
	       T5s = T5q - T5r;
	       T5v = T5t + T5u;
	       T5w = FNMS(KP1_111140466, T5v, KP1_662939224 * T5s);
	       T5U = FMA(KP1_111140466, T5s, KP1_662939224 * T5v);
	       {
		    E T5A, T5H, T5L, T5O;
		    T5A = T5y - T5z;
		    T5H = T5D + T5G;
		    T5I = T5A + T5H;
		    T5W = T5A - T5H;
		    T5L = T5J - T5K;
		    T5O = T5M + T5N;
		    T5P = T5L + T5O;
		    T5X = T5O - T5L;
	       }
	       {
		    E T5x, T5Q, T5Z, T60;
		    T5x = T5p + T5w;
		    T5Q = FNMS(KP580569354, T5P, KP1_913880671 * T5I);
		    R1[WS(rs, 17)] = T5x - T5Q;
		    R1[WS(rs, 1)] = T5x + T5Q;
		    T5Z = T5T + T5U;
		    T60 = FMA(KP1_763842528, T5W, KP942793473 * T5X);
		    R1[WS(rs, 13)] = T5Z - T60;
		    R1[WS(rs, 29)] = T5Z + T60;
	       }
	       {
		    E T5R, T5S, T5V, T5Y;
		    T5R = T5p - T5w;
		    T5S = FMA(KP580569354, T5I, KP1_913880671 * T5P);
		    R1[WS(rs, 9)] = T5R - T5S;
		    R1[WS(rs, 25)] = T5R + T5S;
		    T5V = T5T - T5U;
		    T5Y = FNMS(KP1_763842528, T5X, KP942793473 * T5W);
		    R1[WS(rs, 21)] = T5V - T5Y;
		    R1[WS(rs, 5)] = T5V + T5Y;
	       }
	  }
	  {
	       E T3N, T4N, T46, T4O, T4y, T4Q, T4J, T4R, T3Y, T45;
	       T3N = T3F + T3M;
	       T4N = T3F - T3M;
	       T3Y = T3Q + T3X;
	       T45 = T41 + T44;
	       T46 = FNMS(KP390180644, T45, KP1_961570560 * T3Y);
	       T4O = FMA(KP390180644, T3Y, KP1_961570560 * T45);
	       {
		    E T4i, T4x, T4B, T4I;
		    T4i = T4a + T4h;
		    T4x = T4p + T4w;
		    T4y = T4i + T4x;
		    T4Q = T4i - T4x;
		    T4B = T4z + T4A;
		    T4I = T4E + T4H;
		    T4J = T4B + T4I;
		    T4R = T4I - T4B;
	       }
	       {
		    E T47, T4K, T4T, T4U;
		    T47 = T3N + T46;
		    T4K = FNMS(KP196034280, T4J, KP1_990369453 * T4y);
		    R1[WS(rs, 16)] = T47 - T4K;
		    R1[0] = T47 + T4K;
		    T4T = T4N + T4O;
		    T4U = FMA(KP1_546020906, T4Q, KP1_268786568 * T4R);
		    R1[WS(rs, 12)] = T4T - T4U;
		    R1[WS(rs, 28)] = T4T + T4U;
	       }
	       {
		    E T4L, T4M, T4P, T4S;
		    T4L = T3N - T46;
		    T4M = FMA(KP196034280, T4y, KP1_990369453 * T4J);
		    R1[WS(rs, 8)] = T4L - T4M;
		    R1[WS(rs, 24)] = T4L + T4M;
		    T4P = T4N - T4O;
		    T4S = FNMS(KP1_546020906, T4R, KP1_268786568 * T4Q);
		    R1[WS(rs, 20)] = T4P - T4S;
		    R1[WS(rs, 4)] = T4P + T4S;
	       }
	  }
	  {
	       E T63, T6h, T66, T6i, T6a, T6k, T6d, T6l, T64, T65;
	       T63 = T61 - T62;
	       T6h = T61 + T62;
	       T64 = T5q + T5r;
	       T65 = T5u - T5t;
	       T66 = FNMS(KP1_961570560, T65, KP390180644 * T64);
	       T6i = FMA(KP1_961570560, T64, KP390180644 * T65);
	       {
		    E T68, T69, T6b, T6c;
		    T68 = T5y + T5z;
		    T69 = T5J + T5K;
		    T6a = T68 - T69;
		    T6k = T68 + T69;
		    T6b = T5D - T5G;
		    T6c = T5N - T5M;
		    T6d = T6b + T6c;
		    T6l = T6c - T6b;
	       }
	       {
		    E T67, T6e, T6n, T6o;
		    T67 = T63 + T66;
		    T6e = FNMS(KP1_268786568, T6d, KP1_546020906 * T6a);
		    R1[WS(rs, 19)] = T67 - T6e;
		    R1[WS(rs, 3)] = T67 + T6e;
		    T6n = T6h + T6i;
		    T6o = FMA(KP1_990369453, T6k, KP196034280 * T6l);
		    R1[WS(rs, 15)] = T6n - T6o;
		    R1[WS(rs, 31)] = T6n + T6o;
	       }
	       {
		    E T6f, T6g, T6j, T6m;
		    T6f = T63 - T66;
		    T6g = FMA(KP1_268786568, T6a, KP1_546020906 * T6d);
		    R1[WS(rs, 11)] = T6f - T6g;
		    R1[WS(rs, 27)] = T6f + T6g;
		    T6j = T6h - T6i;
		    T6m = FNMS(KP1_990369453, T6l, KP196034280 * T6k);
		    R1[WS(rs, 23)] = T6j - T6m;
		    R1[WS(rs, 7)] = T6j + T6m;
	       }
	  }
	  {
	       E T4X, T5b, T50, T5c, T54, T5e, T57, T5f, T4Y, T4Z;
	       T4X = T4V - T4W;
	       T5b = T4V + T4W;
	       T4Y = T3Q - T3X;
	       T4Z = T44 - T41;
	       T50 = FNMS(KP1_662939224, T4Z, KP1_111140466 * T4Y);
	       T5c = FMA(KP1_662939224, T4Y, KP1_111140466 * T4Z);
	       {
		    E T52, T53, T55, T56;
		    T52 = T4a - T4h;
		    T53 = T4A - T4z;
		    T54 = T52 + T53;
		    T5e = T52 - T53;
		    T55 = T4p - T4w;
		    T56 = T4H - T4E;
		    T57 = T55 + T56;
		    T5f = T56 - T55;
	       }
	       {
		    E T51, T58, T5h, T5i;
		    T51 = T4X + T50;
		    T58 = FNMS(KP942793473, T57, KP1_763842528 * T54);
		    R1[WS(rs, 18)] = T51 - T58;
		    R1[WS(rs, 2)] = T51 + T58;
		    T5h = T5b + T5c;
		    T5i = FMA(KP1_913880671, T5e, KP580569354 * T5f);
		    R1[WS(rs, 14)] = T5h - T5i;
		    R1[WS(rs, 30)] = T5h + T5i;
	       }
	       {
		    E T59, T5a, T5d, T5g;
		    T59 = T4X - T50;
		    T5a = FMA(KP942793473, T54, KP1_763842528 * T57);
		    R1[WS(rs, 10)] = T59 - T5a;
		    R1[WS(rs, 26)] = T59 + T5a;
		    T5d = T5b - T5c;
		    T5g = FNMS(KP1_913880671, T5f, KP580569354 * T5e);
		    R1[WS(rs, 22)] = T5d - T5g;
		    R1[WS(rs, 6)] = T5d + T5g;
	       }
	  }
     }
}

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

void X(codelet_r2cb_64) (planner *p) {
     X(kr2c_register) (p, r2cb_64, &desc);
}

#endif				/* HAVE_FMA */