本文整理汇总了C++中FMS函数的典型用法代码示例。如果您正苦于以下问题：C++ FMS函数的具体用法？C++ FMS怎么用？C++ FMS使用的例子？那么恭喜您, 这里精选的函数代码示例或许可以为您提供帮助。
在下文中一共展示了FMS函数的15个代码示例，这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞，您的评价将有助于系统推荐出更棒的C++代码示例。
示例1: r2cbIII_5
static void r2cbIII_5(R *R0, R *R1, R *Cr, R *Ci, stride rs, stride csr, stride csi, INT v, INT ivs, INT ovs)
{
DK(KP1_902113032, +1.902113032590307144232878666758764286811397268);
DK(KP1_118033988, +1.118033988749894848204586834365638117720309180);
DK(KP500000000, +0.500000000000000000000000000000000000000000000);
DK(KP2_000000000, +2.000000000000000000000000000000000000000000000);
DK(KP618033988, +0.618033988749894848204586834365638117720309180);
{
INT i;
for (i = v; i > 0; i = i - 1, R0 = R0 + ovs, R1 = R1 + ovs, Cr = Cr + ivs, Ci = Ci + ivs, MAKE_VOLATILE_STRIDE(rs), MAKE_VOLATILE_STRIDE(csr), MAKE_VOLATILE_STRIDE(csi)) {
E T1, T2, T3, Tc, Ta, T8, T9;
T8 = Ci[WS(csi, 1)];
T9 = Ci[0];
T1 = Cr[WS(csr, 2)];
T2 = Cr[WS(csr, 1)];
T3 = Cr[0];
Tc = FMS(KP618033988, T8, T9);
Ta = FMA(KP618033988, T9, T8);
{
E T6, T4, T5, T7, Tb;
T6 = T3 - T2;
T4 = T2 + T3;
R0[0] = FMA(KP2_000000000, T4, T1);
T5 = FNMS(KP500000000, T4, T1);
T7 = FNMS(KP1_118033988, T6, T5);
Tb = FMA(KP1_118033988, T6, T5);
R0[WS(rs, 2)] = FNMS(KP1_902113032, Ta, T7);
R1[0] = -(FMA(KP1_902113032, Ta, T7));
R1[WS(rs, 1)] = FMS(KP1_902113032, Tc, Tb);
R0[WS(rs, 1)] = FMA(KP1_902113032, Tc, Tb);
}
}
}
}示例2: hc2rIII_5
static void hc2rIII_5(const R *ri, const R *ii, R *O, stride ris, stride iis, stride os, INT v, INT ivs, INT ovs)
{
DK(KP1_902113032, +1.902113032590307144232878666758764286811397268);
DK(KP1_118033988, +1.118033988749894848204586834365638117720309180);
DK(KP500000000, +0.500000000000000000000000000000000000000000000);
DK(KP2_000000000, +2.000000000000000000000000000000000000000000000);
DK(KP618033988, +0.618033988749894848204586834365638117720309180);
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 T1, T2, T3, Tc, Ta, T8, T9;
T8 = ii[WS(iis, 1)];
T9 = ii[0];
T1 = ri[WS(ris, 2)];
T2 = ri[WS(ris, 1)];
T3 = ri[0];
Tc = FMS(KP618033988, T8, T9);
Ta = FMA(KP618033988, T9, T8);
{
E T6, T4, T5, T7, Tb;
T6 = T3 - T2;
T4 = T2 + T3;
O[0] = FMA(KP2_000000000, T4, T1);
T5 = FNMS(KP500000000, T4, T1);
T7 = FNMS(KP1_118033988, T6, T5);
Tb = FMA(KP1_118033988, T6, T5);
O[WS(os, 4)] = FNMS(KP1_902113032, Ta, T7);
O[WS(os, 1)] = -(FMA(KP1_902113032, Ta, T7));
O[WS(os, 3)] = FMS(KP1_902113032, Tc, Tb);
O[WS(os, 2)] = FMA(KP1_902113032, Tc, Tb);
}
}
}示例3: findMedianSortedArrays
double findMedianSortedArrays(int A[], int m, int B[], int n) {
int total = m + n ;
if(total%2 == 1)
return FMS(A,m,B,n,total/2+1);
else
return (FMS(A,m,B,n,total/2+1) + FMS(A,m,B,n,total/2))/2;
}示例4: r2cfII_8
static void r2cfII_8(R *R0, R *R1, R *Cr, R *Ci, stride rs, stride csr, stride csi, INT v, INT ivs, INT ovs)
{
DK(KP923879532, +0.923879532511286756128183189396788286822416626);
DK(KP707106781, +0.707106781186547524400844362104849039284835938);
DK(KP414213562, +0.414213562373095048801688724209698078569671875);
{
INT i;
for (i = v; i > 0; i = i - 1, R0 = R0 + ivs, R1 = R1 + ivs, Cr = Cr + ovs, Ci = Ci + ovs, MAKE_VOLATILE_STRIDE(rs), MAKE_VOLATILE_STRIDE(csr), MAKE_VOLATILE_STRIDE(csi)) {
E Te, T8, Td, T5, Tj, Tl, Tf, Tb;
{
E T1, Th, T9, Ti, T4, Ta;
T1 = R0[0];
Th = R0[WS(rs, 2)];
{
E T2, T3, T6, T7;
T2 = R0[WS(rs, 1)];
T3 = R0[WS(rs, 3)];
T6 = R1[0];
T7 = R1[WS(rs, 2)];
T9 = R1[WS(rs, 3)];
Ti = T2 + T3;
T4 = T2 - T3;
Te = FMA(KP414213562, T6, T7);
T8 = FNMS(KP414213562, T7, T6);
Ta = R1[WS(rs, 1)];
}
Td = FNMS(KP707106781, T4, T1);
T5 = FMA(KP707106781, T4, T1);
Tj = FMA(KP707106781, Ti, Th);
Tl = FNMS(KP707106781, Ti, Th);
Tf = FMA(KP414213562, T9, Ta);
Tb = FMS(KP414213562, Ta, T9);
}
{
E Tk, Tg, Tc, Tm;
Tk = Te + Tf;
Tg = Te - Tf;
Tc = T8 + Tb;
Tm = Tb - T8;
Cr[WS(csr, 1)] = FMA(KP923879532, Tg, Td);
Cr[WS(csr, 2)] = FNMS(KP923879532, Tg, Td);
Ci[WS(csi, 3)] = FNMS(KP923879532, Tk, Tj);
Ci[0] = -(FMA(KP923879532, Tk, Tj));
Ci[WS(csi, 1)] = FMA(KP923879532, Tm, Tl);
Ci[WS(csi, 2)] = FMS(KP923879532, Tm, Tl);
Cr[0] = FMA(KP923879532, Tc, T5);
Cr[WS(csr, 3)] = FNMS(KP923879532, Tc, T5);
}
}
}
}示例5: FMS
double FMS(int a[], int m , int b[], int n, const int k)
{
if(m > n )
return FMS(b,n,a,m,k);
if(m == 0)
return b[k-1];
if(k == 1)
return min(a[0],b[0]);
int pa = min(m,k/2);
int pb = k - pa;
if(a[pa-1] < b[pb-1])
return FMS(a+pa,m-pa,b,n,k-pa);
else
return FMS(a,m,b+pb,n-pb,k-pb);
}示例6: hc2rIII_5
static void hc2rIII_5(const R *ri, const R *ii, R *O, stride ris, stride iis, stride os, int v, int ivs, int ovs)
{
DK(KP2_000000000, +2.000000000000000000000000000000000000000000000);
DK(KP1_118033988, +1.118033988749894848204586834365638117720309180);
DK(KP500000000, +0.500000000000000000000000000000000000000000000);
DK(KP1_175570504, +1.175570504584946258337411909278145537195304875);
DK(KP1_902113032, +1.902113032590307144232878666758764286811397268);
int i;
for (i = v; i > 0; i = i - 1, ri = ri + ivs, ii = ii + ivs, O = O + ovs) {
E Ta, Tc, T1, T4, T5, T6, Tb, T7;
{
E T8, T9, T2, T3;
T8 = ii[WS(iis, 1)];
T9 = ii[0];
Ta = FMA(KP1_902113032, T8, KP1_175570504 * T9);
Tc = FNMS(KP1_902113032, T9, KP1_175570504 * T8);
T1 = ri[WS(ris, 2)];
T2 = ri[WS(ris, 1)];
T3 = ri[0];
T4 = T2 + T3;
T5 = FMS(KP500000000, T4, T1);
T6 = KP1_118033988 * (T3 - T2);
}
O[0] = FMA(KP2_000000000, T4, T1);
Tb = T6 - T5;
O[WS(os, 2)] = Tb + Tc;
O[WS(os, 3)] = Tc - Tb;
T7 = T5 + T6;
O[WS(os, 1)] = T7 - Ta;
O[WS(os, 4)] = -(T7 + Ta);
}
}示例7: r2cfII_6
static void r2cfII_6(R *R0, R *R1, R *Cr, R *Ci, stride rs, stride csr, stride csi, INT v, INT ivs, INT ovs)
{
DK(KP866025403, +0.866025403784438646763723170752936183471402627);
DK(KP500000000, +0.500000000000000000000000000000000000000000000);
{
INT i;
for (i = v; i > 0; i = i - 1, R0 = R0 + ivs, R1 = R1 + ivs, Cr = Cr + ovs, Ci = Ci + ovs, MAKE_VOLATILE_STRIDE(24, rs), MAKE_VOLATILE_STRIDE(24, csr), MAKE_VOLATILE_STRIDE(24, csi)) {
E T1, T9, T2, T3, T6, T7;
T1 = R0[0];
T9 = R1[WS(rs, 1)];
T2 = R0[WS(rs, 2)];
T3 = R0[WS(rs, 1)];
T6 = R1[WS(rs, 2)];
T7 = R1[0];
{
E Tc, T4, Ta, T8, T5, Tb;
Cr[WS(csr, 1)] = T1 + T2 - T3;
Tc = T2 + T3;
T4 = T3 - T2;
Ta = T6 + T7;
T8 = T6 - T7;
T5 = FMA(KP500000000, T4, T1);
Tb = FMA(KP500000000, Ta, T9);
Ci[WS(csi, 1)] = T9 - Ta;
Cr[WS(csr, 2)] = FMA(KP866025403, T8, T5);
Cr[0] = FNMS(KP866025403, T8, T5);
Ci[WS(csi, 2)] = FMS(KP866025403, Tc, Tb);
Ci[0] = -(FMA(KP866025403, Tc, Tb));
}
}
}
}示例8: r2hcII_6
static void r2hcII_6(const R *I, R *ro, R *io, stride is, stride ros, stride ios, INT v, INT ivs, INT ovs)
{
DK(KP866025403, +0.866025403784438646763723170752936183471402627);
DK(KP500000000, +0.500000000000000000000000000000000000000000000);
INT i;
for (i = v; i > 0; i = i - 1, I = I + ivs, ro = ro + ovs, io = io + ovs, MAKE_VOLATILE_STRIDE(is), MAKE_VOLATILE_STRIDE(ros), MAKE_VOLATILE_STRIDE(ios)) {
E T1, T9, T2, T3, T6, T7;
T1 = I[0];
T9 = I[WS(is, 3)];
T2 = I[WS(is, 4)];
T3 = I[WS(is, 2)];
T6 = I[WS(is, 5)];
T7 = I[WS(is, 1)];
{
E Tc, T4, Ta, T8, T5, Tb;
ro[WS(ros, 1)] = T1 + T2 - T3;
Tc = T2 + T3;
T4 = T3 - T2;
Ta = T6 + T7;
T8 = T6 - T7;
T5 = FMA(KP500000000, T4, T1);
Tb = FMA(KP500000000, Ta, T9);
io[WS(ios, 1)] = T9 - Ta;
ro[WS(ros, 2)] = FMA(KP866025403, T8, T5);
ro[0] = FNMS(KP866025403, T8, T5);
io[WS(ios, 2)] = FMS(KP866025403, Tc, Tb);
io[0] = -(FMA(KP866025403, Tc, Tb));
}
}
}示例9: r2cbIII_7
static void r2cbIII_7(R *R0, R *R1, R *Cr, R *Ci, stride rs, stride csr, stride csi, INT v, INT ivs, INT ovs)
{
DK(KP1_949855824, +1.949855824363647214036263365987862434465571601);
DK(KP1_801937735, +1.801937735804838252472204639014890102331838324);
DK(KP2_000000000, +2.000000000000000000000000000000000000000000000);
DK(KP692021471, +0.692021471630095869627814897002069140197260599);
DK(KP801937735, +0.801937735804838252472204639014890102331838324);
DK(KP356895867, +0.356895867892209443894399510021300583399127187);
DK(KP554958132, +0.554958132087371191422194871006410481067288862);
{
INT i;
for (i = v; i > 0; i = i - 1, R0 = R0 + ovs, R1 = R1 + ovs, Cr = Cr + ivs, Ci = Ci + ivs, MAKE_VOLATILE_STRIDE(rs), MAKE_VOLATILE_STRIDE(csr), MAKE_VOLATILE_STRIDE(csi)) {
E Tn, Td, Tg, Ti, Tl, T8;
{
E T1, T9, Tb, Ta, T2, T4, Th, Tm, Tc, T3, Te;
T1 = Cr[WS(csr, 3)];
T9 = Ci[WS(csi, 1)];
Tb = Ci[0];
Ta = Ci[WS(csi, 2)];
T2 = Cr[WS(csr, 2)];
T4 = Cr[0];
Th = FMA(KP554958132, T9, Tb);
Tm = FNMS(KP554958132, Ta, T9);
Tc = FMA(KP554958132, Tb, Ta);
T3 = Cr[WS(csr, 1)];
Te = FNMS(KP356895867, T2, T4);
Tn = FNMS(KP801937735, Tm, Tb);
{
E Tf, Tk, T7, T5, Tj, T6;
Td = FMA(KP801937735, Tc, T9);
T5 = T2 + T3 + T4;
Tj = FNMS(KP356895867, T4, T3);
T6 = FNMS(KP356895867, T3, T2);
Tf = FNMS(KP692021471, Te, T3);
R0[0] = FMA(KP2_000000000, T5, T1);
Tk = FNMS(KP692021471, Tj, T2);
T7 = FNMS(KP692021471, T6, T4);
Tg = FNMS(KP1_801937735, Tf, T1);
Ti = FNMS(KP801937735, Th, Ta);
Tl = FNMS(KP1_801937735, Tk, T1);
T8 = FNMS(KP1_801937735, T7, T1);
}
}
R1[WS(rs, 2)] = FMS(KP1_949855824, Ti, Tg);
R0[WS(rs, 1)] = FMA(KP1_949855824, Ti, Tg);
R0[WS(rs, 2)] = FNMS(KP1_949855824, Tn, Tl);
R1[WS(rs, 1)] = -(FMA(KP1_949855824, Tn, Tl));
R0[WS(rs, 3)] = FNMS(KP1_949855824, Td, T8);
R1[0] = -(FMA(KP1_949855824, Td, T8));
}
}
}示例10: r2cf_8
static void r2cf_8(R *R0, R *R1, R *Cr, R *Ci, stride rs, stride csr, stride csi, INT v, INT ivs, INT ovs)
{
DK(KP707106781, +0.707106781186547524400844362104849039284835938);
{
INT i;
for (i = v; i > 0; i = i - 1, R0 = R0 + ivs, R1 = R1 + ivs, Cr = Cr + ovs, Ci = Ci + ovs, MAKE_VOLATILE_STRIDE(rs), MAKE_VOLATILE_STRIDE(csr), MAKE_VOLATILE_STRIDE(csi)) {
E T4, T7, T3, Tj, Td, T5, T8, T9;
{
E T1, T2, Tb, Tc;
T1 = R0[0];
T2 = R0[WS(rs, 2)];
Tb = R1[WS(rs, 3)];
Tc = R1[WS(rs, 1)];
T4 = R0[WS(rs, 1)];
T7 = T1 - T2;
T3 = T1 + T2;
Tj = Tb + Tc;
Td = Tb - Tc;
T5 = R0[WS(rs, 3)];
T8 = R1[0];
T9 = R1[WS(rs, 2)];
}
{
E T6, Tf, Ta, Ti;
T6 = T4 + T5;
Tf = T4 - T5;
Ta = T8 - T9;
Ti = T8 + T9;
{
E Th, Tk, Te, Tg;
Th = T3 + T6;
Cr[WS(csr, 2)] = T3 - T6;
Tk = Ti + Tj;
Ci[WS(csi, 2)] = Tj - Ti;
Te = Ta + Td;
Tg = Td - Ta;
Cr[0] = Th + Tk;
Cr[WS(csr, 4)] = Th - Tk;
Ci[WS(csi, 3)] = FMA(KP707106781, Tg, Tf);
Ci[WS(csi, 1)] = FMS(KP707106781, Tg, Tf);
Cr[WS(csr, 1)] = FMA(KP707106781, Te, T7);
Cr[WS(csr, 3)] = FNMS(KP707106781, Te, T7);
}
}
}
}
}示例11: hc2cf2_20
//.........这里部分代码省略.........
T4u = FNMS(KP559016994, T4t, T4s);
T4C = FMA(KP559016994, T4t, T4s);
}
{
E T2x, T4Q, T4p, T4r, T4R, T2A, T51, T55;
{
E T4h, T50, T4X, T4o, T4Z;
T4E = T4d + T4g;
T4h = T4d - T4g;
Rm[WS(rs, 1)] = FMA(KP951056516, T4B, T4u);
Rp[WS(rs, 2)] = FNMS(KP951056516, T4B, T4u);
Rp[WS(rs, 6)] = FMA(KP951056516, T4D, T4C);
Rm[WS(rs, 5)] = FNMS(KP951056516, T4D, T4C);
T50 = T4W - T4V;
T4X = T4V + T4W;
T4o = T4k - T4n;
T4F = T4k + T4n;
T2x = T2v + T2w;
T4Q = T2v - T2w;
Im[WS(rs, 9)] = T4X - T4Y;
T4Z = FMA(KP250000000, T4X, T4Y);
T4p = FMA(KP618033988, T4o, T4h);
T4r = FNMS(KP618033988, T4h, T4o);
T4R = T2z - T2y;
T2A = T2y + T2z;
T51 = FNMS(KP559016994, T50, T4Z);
T55 = FMA(KP559016994, T50, T4Z);
}
{
E T49, T48, T2B, T4a, T4q;
T2B = T2x + T2A;
T49 = T2x - T2A;
Ip[WS(rs, 2)] = FMA(KP951056516, T54, T51);
Im[WS(rs, 1)] = FMS(KP951056516, T54, T51);
Ip[WS(rs, 6)] = FMA(KP951056516, T56, T55);
Im[WS(rs, 5)] = FMS(KP951056516, T56, T55);
Rp[0] = T2u + T2B;
T48 = FNMS(KP250000000, T2B, T2u);
T4a = FMA(KP559016994, T49, T48);
T4q = FNMS(KP559016994, T49, T48);
T4U = FMA(KP618033988, T4Q, T4R);
T4S = FNMS(KP618033988, T4R, T4Q);
Rm[WS(rs, 3)] = FMA(KP951056516, T4p, T4a);
Rp[WS(rs, 4)] = FNMS(KP951056516, T4p, T4a);
Rp[WS(rs, 8)] = FMA(KP951056516, T4r, T4q);
Rm[WS(rs, 7)] = FNMS(KP951056516, T4r, T4q);
}
}
}
{
E T3O, T5u, T5w, T5o, T5q, T5n;
{
E T5m, T5l, T2I, T4O, T3N, T3L, T2X, T5s, T4N, T5t, T3c, T3v, T3K, T4G;
T5m = T3u + T3n;
T3v = T3n - T3u;
T3K = T3C - T3J;
T5l = T3J + T3C;
T3O = T2C + T2H;
T2I = T2C - T2H;
T4O = T4E - T4F;
T4G = T4E + T4F;
T3N = FMA(KP618033988, T3v, T3K);
T3L = FNMS(KP618033988, T3K, T3v);
T2X = T2P + T2W;
T5s = T2P - T2W;
Ip[0] = T4G + T4M;示例12: DK
//.........这里部分代码省略.........
T1V = FNMS(KP866025403, T1U, T1P);
T3c = FNMS(KP618033988, T3b, T38);
T3e = FMA(KP618033988, T38, T3b);
T3q = FNMS(KP618033988, T3p, T3o);
T3s = FMA(KP618033988, T3o, T3p);
T1A = TU + T1z;
T34 = TU - T1z;
}
}
{
E T2W, T33, T3m, T3h, T2Z, T3d, T35, T3l;
T3J = T2U + T2V;
T2W = T2U - T2V;
rio[0] = Tf + T1A;
T33 = FNMS(KP250000000, T1A, Tf);
T3m = T3f - T3g;
T3h = T3f + T3g;
T2Z = T2X - T2Y;
T3I = T2X + T2Y;
T3d = FMA(KP559016994, T34, T33);
T35 = FNMS(KP559016994, T34, T33);
iio[0] = T3h + T3k;
T3l = FNMS(KP250000000, T3h, T3k);
rio[WS(ios, 3)] = FMA(KP951056516, T3c, T35);
iio[-WS(ios, 12)] = FNMS(KP951056516, T3c, T35);
rio[WS(ios, 6)] = FMA(KP951056516, T3e, T3d);
iio[-WS(ios, 9)] = FNMS(KP951056516, T3e, T3d);
T3r = FMA(KP559016994, T3m, T3l);
T3n = FNMS(KP559016994, T3m, T3l);
T32 = FMA(KP618033988, T2W, T2Z);
T30 = FNMS(KP618033988, T2Z, T2W);
}
iio[-WS(ios, 3)] = FMA(KP951056516, T3q, T3n);
rio[WS(ios, 12)] = FMS(KP951056516, T3q, T3n);
iio[-WS(ios, 6)] = FMA(KP951056516, T3s, T3r);
rio[WS(ios, 9)] = FMS(KP951056516, T3s, T3r);
T2d = FMA(KP866025403, T1G, T1B);
T1H = FNMS(KP866025403, T1G, T1B);
T1W = T1O + T1V;
T3O = T1O - T1V;
T3P = T23 - T2a;
T2b = T23 + T2a;
}
}
{
E T3H, T3y, T2S, T3Q, T3S, T2R, T2c;
T3H = FNMS(KP866025403, T3x, T3w);
T3y = FMA(KP866025403, T3x, T3w);
T2c = T1W + T2b;
T2S = T1W - T2b;
T3Q = FMA(KP618033988, T3P, T3O);
T3S = FNMS(KP618033988, T3O, T3P);
rio[WS(ios, 5)] = T1H + T2c;
T2R = FNMS(KP250000000, T2c, T1H);
{
E T2g, T2j, T3G, T3E, T2A, T2N, T3v, T3A, T3M, T3L, T3z, T3F, T3B;
{
E T3C, T3D, T31, T2T, T3K;
T2g = T2e + T2f;
T3C = T2e - T2f;
T3D = T2h - T2i;
T2j = T2h + T2i;
T31 = FMA(KP559016994, T2S, T2R);
T2T = FNMS(KP559016994, T2S, T2R);
T3K = T3I + T3J;
T3M = T3J - T3I;示例13: hc2cf2_16
//.........这里部分代码省略.........
E T3m, T3a, T3J, T3H;
Ip[0] = T3u + T3C;
Im[WS(rs, 7)] = T3u - T3C;
Rp[WS(rs, 4)] = T3q + T3t;
Rm[WS(rs, 3)] = T3q - T3t;
T3m = T36 - T39;
T3a = T36 + T39;
T3J = T3G - T3F;
T3H = T3F + T3G;
{
E T2Q, T20, T3N, T3T, T2J, T2C, T3O, T2f, T34, T30, T2W, T2V, T3U, T2T, T2N;
E T2v;
{
E T2R, T27, T2e, T2S;
{
E T3n, T3f, T3o, T3k;
T2Q = T1U + T1Z;
T20 = T1U - T1Z;
T3n = T3e - T3b;
T3f = T3b + T3e;
T3o = T3g + T3j;
T3k = T3g - T3j;
T3N = T3L - T3M;
T3T = T3M + T3L;
{
E T3p, T3I, T3K, T3l;
T3p = T3n - T3o;
T3I = T3n + T3o;
T3K = T3k - T3f;
T3l = T3f + T3k;
Rp[WS(rs, 6)] = FMA(KP707106781, T3p, T3m);
Rm[WS(rs, 1)] = FNMS(KP707106781, T3p, T3m);
Ip[WS(rs, 2)] = FMA(KP707106781, T3I, T3H);
Im[WS(rs, 5)] = FMS(KP707106781, T3I, T3H);
Ip[WS(rs, 6)] = FMA(KP707106781, T3K, T3J);
Im[WS(rs, 1)] = FMS(KP707106781, T3K, T3J);
Rp[WS(rs, 2)] = FMA(KP707106781, T3l, T3a);
Rm[WS(rs, 5)] = FNMS(KP707106781, T3l, T3a);
T2R = T26 + T25;
T27 = T25 - T26;
T2e = T28 + T2d;
T2S = T28 - T2d;
}
}
{
E T2Y, T2Z, T2n, T2u;
T2J = T2D - T2I;
T2Y = T2D + T2I;
T2Z = T2B + T2w;
T2C = T2w - T2B;
T3O = T27 + T2e;
T2f = T27 - T2e;
T34 = FMA(KP414213562, T2Y, T2Z);
T30 = FNMS(KP414213562, T2Z, T2Y);
T2W = T2l - T2m;
T2n = T2l + T2m;
T2u = T2o - T2t;
T2V = T2o + T2t;
T3U = T2S - T2R;
T2T = T2R + T2S;
T2N = FNMS(KP414213562, T2n, T2u);
T2v = FMA(KP414213562, T2u, T2n);
}
}
{
E T33, T2X, T3X, T3Y;示例14: DK
//.........这里部分代码省略.........
TQ = iio[-WS(ios, 4)];
Td = FMA(Tb, Tc, T9);
T1q = T7 * Tc;
TL = Tl * TK;
T1b = T3 * TQ;
TM = iio[0];
T1r = FNMS(Tb, T8, T1q);
TO = rio[WS(ios, 3)];
{
E TN, T1c, T1a, TR, T19, TP;
Tx = rio[WS(ios, 1)];
TN = FMA(Tn, TM, TL);
T19 = Tl * TM;
T1c = FNMS(T6, TO, T1b);
TP = T3 * TO;
Ty = T2 * Tx;
T1a = FNMS(Tn, TK, T19);
TH = iio[-WS(ios, 2)];
TR = FMA(T6, TQ, TP);
Tz = iio[-WS(ios, 6)];
T1m = T1a + T1c;
T1d = T1a - T1c;
T14 = TC * TH;
TS = TN + TR;
T18 = TN - TR;
T12 = T2 * Tz;
TD = rio[WS(ios, 5)];
}
TA = FMA(T5, Tz, Ty);
}
}
{
E Te, T1p, T1t, T10, T1g, Tv, TJ, T1i, T1e, T1z, T1B, T17, T1h;
{
E T1x, T16, T1y, T11;
{
E TU, T13, T15, TE, TZ, TI;
Te = T1 + Td;
TU = T1 - Td;
T13 = FNMS(T5, Tx, T12);
T15 = FNMS(TG, TD, T14);
TE = TC * TD;
T1p = TW + TY;
TZ = TW - TY;
T1x = T1s - T1r;
T1t = T1r + T1s;
T1l = T13 + T15;
T16 = T13 - T15;
TI = FMA(TG, TH, TE);
T10 = TU + TZ;
T1g = TU - TZ;
Tv = Tk + Tu;
T1y = Tk - Tu;
T11 = TA - TI;
TJ = TA + TI;
}
T1i = T18 + T1d;
T1e = T18 - T1d;
T1z = T1x - T1y;
T1B = T1y + T1x;
T17 = T11 + T16;
T1h = T16 - T11;
}
{
E T1j, T1A, T1f, T1C;
T1j = T1h - T1i;
T1A = T1h + T1i;
T1f = T17 + T1e;
T1C = T1e - T17;
iio[-WS(ios, 1)] = FMA(KP707106781, T1A, T1z);
rio[WS(ios, 5)] = FMS(KP707106781, T1A, T1z);
rio[WS(ios, 3)] = FMA(KP707106781, T1j, T1g);
iio[-WS(ios, 7)] = FNMS(KP707106781, T1j, T1g);
iio[-WS(ios, 3)] = FMA(KP707106781, T1C, T1B);
rio[WS(ios, 7)] = FMS(KP707106781, T1C, T1B);
rio[WS(ios, 1)] = FMA(KP707106781, T1f, T10);
iio[-WS(ios, 5)] = FNMS(KP707106781, T1f, T10);
T1k = Te - Tv;
Tw = Te + Tv;
}
T1w = T1t - T1p;
T1u = T1p + T1t;
T1v = TS - TJ;
TT = TJ + TS;
}
}
}
iio[-WS(ios, 2)] = T1v + T1w;
rio[WS(ios, 6)] = T1v - T1w;
rio[0] = Tw + TT;
iio[-WS(ios, 4)] = Tw - TT;
T1n = T1l - T1m;
T1o = T1l + T1m;
iio[0] = T1o + T1u;
rio[WS(ios, 4)] = T1o - T1u;
rio[WS(ios, 2)] = T1k + T1n;
iio[-WS(ios, 6)] = T1k - T1n;
}
return W;
}示例15: r2cfII_25
//.........这里部分代码省略.........
T1t = FMA(KP481753674, T1b, KP1_752613360 * T1e);
T1u = FMA(KP684547105, T1i, KP1_457937254 * T1l);
T1B = T1t + T1u;
T1r = T1p - T1q;
T1C = T1A + T1B;
T1v = T1t - T1u;
}
{
E T10, T17, T1f, T1m;
T10 = FNMS(KP497379774, TZ, KP968583161 * TW);
T17 = FNMS(KP1_688655851, T16, KP535826794 * T13);
T18 = T10 + T17;
T1f = FNMS(KP963507348, T1e, KP876306680 * T1b);
T1m = FNMS(KP1_369094211, T1l, KP728968627 * T1i);
T1n = T1f + T1m;
T1o = T18 + T1n;
T1G = T10 - T17;
T1D = T1f - T1m;
}
{
E T1R, T1Q, T20, T1Z;
Cr[0] = TT + T1o;
Ci[0] = -(T1z + T1C);
T1R = KP559016994 * (T1P + T1M);
T1Q = FMA(KP250000000, T1M - T1P, TT);
Cr[WS(csr, 4)] = FMA(KP951056516, T1J, T1Q) + FMA(KP587785252, T1U, T1R);
Cr[WS(csr, 9)] = FMA(KP951056516, T1U, T1Q) + FNMA(KP587785252, T1J, T1R);
T20 = KP559016994 * (T1Y + T1X);
T1Z = FMA(KP250000000, T1X - T1Y, T1z);
Ci[WS(csi, 9)] = FMA(KP587785252, T1V, KP951056516 * T1W) + T1Z - T20;
Ci[WS(csi, 4)] = FMA(KP587785252, T1W, T1Z) + FNMS(KP951056516, T1V, T20);
{
E T1E, T1F, T1s, T1w;
T1E = FMS(KP250000000, T1C, T1z);
T1F = KP559016994 * (T1B - T1A);
Ci[WS(csi, 5)] = FMA(KP951056516, T1D, T1E) + FNMA(KP587785252, T1G, T1F);
Ci[WS(csi, 10)] = FMA(KP951056516, T1G, KP587785252 * T1D) + T1E + T1F;
T1s = FNMS(KP250000000, T1o, TT);
T1w = KP559016994 * (T18 - T1n);
Cr[WS(csr, 5)] = FMA(KP587785252, T1r, T1s) + FMS(KP951056516, T1v, T1w);
Cr[WS(csr, 10)] = T1w + FMA(KP587785252, T1v, T1s) - (KP951056516 * T1r);
}
}
}
{
E T21, T2z, T2L, T2K, T2M, T2F, T2P, T2C, T2Q, T2l, T2o, T2p, T2w, T2u, T28;
E T2f, T2g, T2s, T2h;
T21 = TS - TR;
{
E T2x, T2y, T2I, T2J;
T2x = FNMS(KP844327925, T29, KP1_071653589 * T2a);
T2y = FNMS(KP125581039, T2d, KP998026728 * T2c);
T2z = T2x + T2y;
T2L = T2y - T2x;
T2I = FNMS(KP481753674, T22, KP1_752613360 * T23);
T2J = FMA(KP904827052, T25, KP851558583 * T26);
T2K = T2I + T2J;
T2M = T2I - T2J;
}
{
E T2D, T2E, T2A, T2B;
T2D = FMA(KP535826794, T29, KP1_688655851 * T2a);
T2E = FMA(KP062790519, T2c, KP1_996053456 * T2d);
T2F = T2D + T2E;
T2P = T2E - T2D;
T2A = FMA(KP876306680, T22, KP963507348 * T23);