hc2r_128.c

最新的FFT程序

     DK(KP390180644, +0.390180644032256535696569736954044481855383236);
     DK(KP1_961570560, +1.961570560806460898252364472268478073947867462);
     DK(KP382683432, +0.382683432365089771728459984030398866761344562);
     DK(KP923879532, +0.923879532511286756128183189396788286822416626);
     DK(KP707106781, +0.707106781186547524400844362104849039284835938);
     DK(KP765366864, +0.765366864730179543456919968060797733522689125);
     DK(KP1_847759065, +1.847759065022573512256366378793576573644833252);
     DK(KP1_414213562, +1.414213562373095048801688724209698078569671875);
     DK(KP2_000000000, +2.000000000000000000000000000000000000000000000);
     INT i;
     for (i = v; i > 0; i = i - 1, ri = ri + ivs, ii = ii + ivs, O = O + ovs, MAKE_VOLATILE_STRIDE(ris), MAKE_VOLATILE_STRIDE(iis), MAKE_VOLATILE_STRIDE(os)) {
	  E Ta, T6q, T2a, T5k, T8x, Tbx, TcF, Ten, Th, T6r, T2j, T5l, T8E, Tby, TcI;
	  E Teo, Tx, T6t, TcM, Teq, TcP, Ter, T2t, T5n, T2C, T5o, T8Q, TbA, T8X, TbB;
	  E T6w, T7L, T1j, T6L, Tde, TeC, TdL, TeR, T3v, T5z, T4I, T5O, T9O, TbM, TaV;
	  E Tc1, T78, T7Z, TN, T6z, TcU, Teu, Td8, Tey, T2N, T5r, T3j, T5v, T9a, TbE;
	  E T9A, TbI, T6H, T7O, T1O, T7V, T48, T4u, Tds, TeG, T5E, T5K, Taf, TbP, Tdp;
	  E TeF, T6U, T72, Tam, TbQ, T23, T7U, T4r, T4v, Tdz, TeJ, T5H, T5L, Tay, TbS;
	  E Tdw, TeI, T6Z, T73, TaF, TbT, T1y, T75, Tdl, TeQ, TdI, TeD, T3O, T5N, T4z;
	  E T5A, Ta3, Tc0, TaO, TbN, T6O, T80, T12, T6E, Td1, Tex, Td5, Tev, T36, T5u;
	  E T3a, T5s, T9p, TbH, T9t, TbF, T6C, T7P;
	  {
	       E T5, T8s, T3, T8q, T9, T8u, T29, T8v, T6, T26;
	       {
		    E T4, T8r, T1, T2;
		    T4 = ri[WS(ris, 32)];
		    T5 = KP2_000000000 * T4;
		    T8r = ii[WS(iis, 32)];
		    T8s = KP2_000000000 * T8r;
		    T1 = ri[0];
		    T2 = ri[WS(ris, 64)];
		    T3 = T1 + T2;
		    T8q = T1 - T2;
		    {
			 E T7, T8, T27, T28;
			 T7 = ri[WS(ris, 16)];
			 T8 = ri[WS(ris, 48)];
			 T9 = KP2_000000000 * (T7 + T8);
			 T8u = T7 - T8;
			 T27 = ii[WS(iis, 16)];
			 T28 = ii[WS(iis, 48)];
			 T29 = KP2_000000000 * (T27 - T28);
			 T8v = T27 + T28;
		    }
	       }
	       T6 = T3 + T5;
	       Ta = T6 + T9;
	       T6q = T6 - T9;
	       T26 = T3 - T5;
	       T2a = T26 - T29;
	       T5k = T26 + T29;
	       {
		    E T8t, T8w, TcD, TcE;
		    T8t = T8q - T8s;
		    T8w = KP1_414213562 * (T8u - T8v);
		    T8x = T8t + T8w;
		    Tbx = T8t - T8w;
		    TcD = T8q + T8s;
		    TcE = KP1_414213562 * (T8u + T8v);
		    TcF = TcD - TcE;
		    Ten = TcD + TcE;
	       }
	  }
	  {
	       E Td, T8y, T2e, T8C, Tg, T8B, T2h, T8z, T2b, T2i;
	       {
		    E Tb, Tc, T2c, T2d;
		    Tb = ri[WS(ris, 8)];
		    Tc = ri[WS(ris, 56)];
		    Td = Tb + Tc;
		    T8y = Tb - Tc;
		    T2c = ii[WS(iis, 8)];
		    T2d = ii[WS(iis, 56)];
		    T2e = T2c - T2d;
		    T8C = T2c + T2d;
	       }
	       {
		    E Te, Tf, T2f, T2g;
		    Te = ri[WS(ris, 40)];
		    Tf = ri[WS(ris, 24)];
		    Tg = Te + Tf;
		    T8B = Te - Tf;
		    T2f = ii[WS(iis, 40)];
		    T2g = ii[WS(iis, 24)];
		    T2h = T2f - T2g;
		    T8z = T2f + T2g;
	       }
	       Th = KP2_000000000 * (Td + Tg);
	       T6r = KP2_000000000 * (T2h + T2e);
	       T2b = Td - Tg;
	       T2i = T2e - T2h;
	       T2j = KP1_414213562 * (T2b - T2i);
	       T5l = KP1_414213562 * (T2b + T2i);
	       {
		    E T8A, T8D, TcG, TcH;
		    T8A = T8y - T8z;
		    T8D = T8B + T8C;
		    T8E = FNMS(KP765366864, T8D, KP1_847759065 * T8A);
		    Tby = FMA(KP765366864, T8A, KP1_847759065 * T8D);
		    TcG = T8y + T8z;
		    TcH = T8C - T8B;
		    TcI = FNMS(KP1_847759065, TcH, KP765366864 * TcG);
		    Teo = FMA(KP1_847759065, TcG, KP765366864 * TcH);
	       }
	  }
	  {
	       E Tl, T8G, T2x, T8V, To, T8U, T2A, T8H, Tv, T8S, T2o, T8O, Ts, T8R, T2r;
	       E T8L;
	       {
		    E Tj, Tk, T2y, T2z;
		    Tj = ri[WS(ris, 4)];
		    Tk = ri[WS(ris, 60)];
		    Tl = Tj + Tk;
		    T8G = Tj - Tk;
		    {
			 E T2v, T2w, Tm, Tn;
			 T2v = ii[WS(iis, 4)];
			 T2w = ii[WS(iis, 60)];
			 T2x = T2v - T2w;
			 T8V = T2v + T2w;
			 Tm = ri[WS(ris, 36)];
			 Tn = ri[WS(ris, 28)];
			 To = Tm + Tn;
			 T8U = Tm - Tn;
		    }
		    T2y = ii[WS(iis, 36)];
		    T2z = ii[WS(iis, 28)];
		    T2A = T2y - T2z;
		    T8H = T2y + T2z;
		    {
			 E Tt, Tu, T8M, T2m, T2n, T8N;
			 Tt = ri[WS(ris, 12)];
			 Tu = ri[WS(ris, 52)];
			 T8M = Tt - Tu;
			 T2m = ii[WS(iis, 52)];
			 T2n = ii[WS(iis, 12)];
			 T8N = T2n + T2m;
			 Tv = Tt + Tu;
			 T8S = T8M + T8N;
			 T2o = T2m - T2n;
			 T8O = T8M - T8N;
		    }
		    {
			 E Tq, Tr, T8J, T2p, T2q, T8K;
			 Tq = ri[WS(ris, 20)];
			 Tr = ri[WS(ris, 44)];
			 T8J = Tq - Tr;
			 T2p = ii[WS(iis, 20)];
			 T2q = ii[WS(iis, 44)];
			 T8K = T2p + T2q;
			 Ts = Tq + Tr;
			 T8R = T8J + T8K;
			 T2r = T2p - T2q;
			 T8L = T8J - T8K;
		    }
	       }
	       {
		    E Tp, Tw, TcK, TcL;
		    Tp = Tl + To;
		    Tw = Ts + Tv;
		    Tx = KP2_000000000 * (Tp + Tw);
		    T6t = Tp - Tw;
		    TcK = T8G + T8H;
		    TcL = KP707106781 * (T8R + T8S);
		    TcM = TcK - TcL;
		    Teq = TcK + TcL;
	       }
	       {
		    E TcN, TcO, T2l, T2s;
		    TcN = KP707106781 * (T8L - T8O);
		    TcO = T8V - T8U;
		    TcP = TcN + TcO;
		    Ter = TcO - TcN;
		    T2l = Tl - To;
		    T2s = T2o - T2r;
		    T2t = T2l + T2s;
		    T5n = T2l - T2s;
	       }
	       {
		    E T2u, T2B, T8I, T8P;
		    T2u = Ts - Tv;
		    T2B = T2x - T2A;
		    T2C = T2u + T2B;
		    T5o = T2B - T2u;
		    T8I = T8G - T8H;
		    T8P = KP707106781 * (T8L + T8O);
		    T8Q = T8I + T8P;
		    TbA = T8I - T8P;
	       }
	       {
		    E T8T, T8W, T6u, T6v;
		    T8T = KP707106781 * (T8R - T8S);
		    T8W = T8U + T8V;
		    T8X = T8T + T8W;
		    TbB = T8W - T8T;
		    T6u = T2A + T2x;
		    T6v = T2r + T2o;
		    T6w = T6u - T6v;
		    T7L = KP2_000000000 * (T6v + T6u);
	       }
	  }
	  {
	       E T17, T9E, T4D, TaT, T1a, TaS, T4G, T9F, T1h, TaQ, T3q, T9M, T1e, TaP, T3t;
	       E T9J;
	       {
		    E T15, T16, T4E, T4F;
		    T15 = ri[WS(ris, 1)];
		    T16 = ri[WS(ris, 63)];
		    T17 = T15 + T16;
		    T9E = T15 - T16;
		    {
			 E T4B, T4C, T18, T19;
			 T4B = ii[WS(iis, 1)];
			 T4C = ii[WS(iis, 63)];
			 T4D = T4B - T4C;
			 TaT = T4B + T4C;
			 T18 = ri[WS(ris, 33)];
			 T19 = ri[WS(ris, 31)];
			 T1a = T18 + T19;
			 TaS = T18 - T19;
		    }
		    T4E = ii[WS(iis, 33)];
		    T4F = ii[WS(iis, 31)];
		    T4G = T4E - T4F;
		    T9F = T4E + T4F;
		    {
			 E T1f, T1g, T9K, T3o, T3p, T9L;
			 T1f = ri[WS(ris, 15)];
			 T1g = ri[WS(ris, 49)];
			 T9K = T1f - T1g;
			 T3o = ii[WS(iis, 49)];
			 T3p = ii[WS(iis, 15)];
			 T9L = T3p + T3o;
			 T1h = T1f + T1g;
			 TaQ = T9K + T9L;
			 T3q = T3o - T3p;
			 T9M = T9K - T9L;
		    }
		    {
			 E T1c, T1d, T9H, T3r, T3s, T9I;
			 T1c = ri[WS(ris, 17)];
			 T1d = ri[WS(ris, 47)];
			 T9H = T1c - T1d;
			 T3r = ii[WS(iis, 17)];
			 T3s = ii[WS(iis, 47)];
			 T9I = T3r + T3s;
			 T1e = T1c + T1d;
			 TaP = T9H + T9I;
			 T3t = T3r - T3s;
			 T9J = T9H - T9I;
		    }
	       }
	       {
		    E T1b, T1i, Tdc, Tdd;
		    T1b = T17 + T1a;
		    T1i = T1e + T1h;
		    T1j = T1b + T1i;
		    T6L = T1b - T1i;
		    Tdc = T9E + T9F;
		    Tdd = KP707106781 * (TaP + TaQ);
		    Tde = Tdc - Tdd;
		    TeC = Tdc + Tdd;
	       }
	       {
		    E TdJ, TdK, T3n, T3u;
		    TdJ = KP707106781 * (T9J - T9M);
		    TdK = TaT - TaS;
		    TdL = TdJ + TdK;
		    TeR = TdK - TdJ;
		    T3n = T17 - T1a;
		    T3u = T3q - T3t;
		    T3v = T3n + T3u;
		    T5z = T3n - T3u;
	       }
	       {
		    E T4A, T4H, T9G, T9N;
		    T4A = T1e - T1h;
		    T4H = T4D - T4G;
		    T4I = T4A + T4H;
		    T5O = T4H - T4A;
		    T9G = T9E - T9F;
		    T9N = KP707106781 * (T9J + T9M);
		    T9O = T9G + T9N;
		    TbM = T9G - T9N;
	       }
	       {
		    E TaR, TaU, T76, T77;
		    TaR = KP707106781 * (TaP - TaQ);
		    TaU = TaS + TaT;
		    TaV = TaR + TaU;
		    Tc1 = TaU - TaR;
		    T76 = T4G + T4D;
		    T77 = T3t + T3q;
		    T78 = T76 - T77;
		    T7Z = T77 + T76;
	       }
	  }
	  {
	       E TB, T90, T3e, T9y, TE, T9x, T3h, T91, TL, T9v, T2I, T98, TI, T9u, T2L;
	       E T95;
	       {
		    E Tz, TA, T3f, T3g;
		    Tz = ri[WS(ris, 2)];
		    TA = ri[WS(ris, 62)];
		    TB = Tz + TA;
		    T90 = Tz - TA;
		    {
			 E T3c, T3d, TC, TD;
			 T3c = ii[WS(iis, 2)];
			 T3d = ii[WS(iis, 62)];
			 T3e = T3c - T3d;
			 T9y = T3c + T3d;
			 TC = ri[WS(ris, 34)];
			 TD = ri[WS(ris, 30)];
			 TE = TC + TD;
			 T9x = TC - TD;
		    }
		    T3f = ii[WS(iis, 34)];
		    T3g = ii[WS(iis, 30)];
		    T3h = T3f - T3g;
		    T91 = T3f + T3g;
		    {
			 E TJ, TK, T96, T2G, T2H, T97;
			 TJ = ri[WS(ris, 14)];
			 TK = ri[WS(ris, 50)];
			 T96 = TJ - TK;
			 T2G = ii[WS(iis, 50)];
			 T2H = ii[WS(iis, 14)];
			 T97 = T2H + T2G;
			 TL = TJ + TK;
			 T9v = T96 + T97;
			 T2I = T2G - T2H;
			 T98 = T96 - T97;
		    }
		    {
			 E TG, TH, T93, T2J, T2K, T94;
			 TG = ri[WS(ris, 18)];
			 TH = ri[WS(ris, 46)];
			 T93 = TG - TH;
			 T2J = ii[WS(iis, 18)];
			 T2K = ii[WS(iis, 46)];
			 T94 = T2J + T2K;
			 TI = TG + TH;
			 T9u = T93 + T94;
			 T2L = T2J - T2K;
			 T95 = T93 - T94;
		    }
	       }
	       {
		    E TF, TM, TcS, TcT;
		    TF = TB + TE;
		    TM = TI + TL;
		    TN = TF + TM;
		    T6z = TF - TM;
		    TcS = T90 + T91;
		    TcT = KP707106781 * (T9u + T9v);
		    TcU = TcS - TcT;
		    Teu = TcS + TcT;
	       }
	       {
		    E Td6, Td7, T2F, T2M;
		    Td6 = KP707106781 * (T95 - T98);
		    Td7 = T9y - T9x;
		    Td8 = Td6 + Td7;
		    Tey = Td7 - Td6;
		    T2F = TB - TE;
		    T2M = T2I - T2L;
		    T2N = T2F + T2M;
		    T5r = T2F - T2M;
	       }
	       {
		    E T3b, T3i, T92, T99;
		    T3b = TI - TL;
		    T3i = T3e - T3h;
		    T3j = T3b + T3i;
		    T5v = T3i - T3b;
		    T92 = T90 - T91;
		    T99 = KP707106781 * (T95 + T98);
		    T9a = T92 + T99;
		    TbE = T92 - T99;
	       }
	       {
		    E T9w, T9z, T6F, T6G;
		    T9w = KP707106781 * (T9u - T9v);
		    T9z = T9x + T9y;
		    T9A = T9w + T9z;
		    TbI = T9z - T9w;
		    T6F = T3h + T3e;
		    T6G = T2L + T2I;
		    T6H = T6F - T6G;
		    T7O = T6G + T6F;
	       }
	  }
	  {
	       E T1G, Taj, T3Q, Ta5, T46, Tak, T6R, Ta6, T1N, Tag, Tah, T3X, T3Z, Taa, Tad;
	       E T6S, Tdn, Tdo;
	       {
		    E T1A, T1B, T1C, T1D, T1E, T1F;
		    T1A = ri[WS(ris, 5)];
		    T1B = ri[WS(ris, 59)];
		    T1C = T1A + T1B;
		    T1D = ri[WS(ris, 37)];
		    T1E = ri[WS(ris, 27)];
		    T1F = T1D + T1E;
		    T1G = T1C + T1F;