本文整理汇总了C++中GenericMatrix::t方法的典型用法代码示例。如果您正苦于以下问题:C++ GenericMatrix::t方法的具体用法?C++ GenericMatrix::t怎么用?C++ GenericMatrix::t使用的例子?那么, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类GenericMatrix
的用法示例。
在下文中一共展示了GenericMatrix::t方法的2个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的C++代码示例。
示例1: trymat1
void trymat1()
{
// cout << "\nFirst test of Matrix package\n\n";
Tracer et("First test of Matrix package");
Tracer::PrintTrace();
{
Tracer et1("Stage 1");
int i,j;
LowerTriangularMatrix L(10);
for (i=1;i<=10;i++) for (j=1;j<=i;j++) L(i,j)=2.0+i*i+j;
SymmetricMatrix S(10);
for (i=1;i<=10;i++) for (j=1;j<=i;j++) S(i,j)=i*j+1.0;
SymmetricMatrix S1 = S / 2.0;
S = S1 * 2.0;
UpperTriangularMatrix U=L.t()*2.0;
Print(LowerTriangularMatrix(L-U.t()*0.5));
DiagonalMatrix D(10);
for (i=1;i<=10;i++) D(i,i)=(i-4)*(i-5)*(i-6);
Matrix M=(S+U-D+L)*(L+U-D+S);
DiagonalMatrix DD=D*D;
LowerTriangularMatrix LD=L*D;
// expressions split for Turbo C
Matrix M1 = S*L + U*L - D*L + L*L + 10.0;
{ M1 = M1 + S*U + U*U - D*U + L*U - S*D; }
{ M1 = M1 - U*D + DD - LD + S*S; }
{ M1 = M1 + U*S - D*S + L*S - 10.0; }
M=M1-M;
Print(M);
}
{
Tracer et1("Stage 2");
int i,j;
LowerTriangularMatrix L(9);
for (i=1;i<=9;i++) for (j=1;j<=i;j++) L(i,j)=1.0+j;
UpperTriangularMatrix U1(9);
for (j=1;j<=9;j++) for (i=1;i<=j;i++) U1(i,j)=1.0+i;
LowerTriangularMatrix LX(9);
for (i=1;i<=9;i++) for (j=1;j<=i;j++) LX(i,j)=1.0+i*i;
UpperTriangularMatrix UX(9);
for (j=1;j<=9;j++) for (i=1;i<=j;i++) UX(i,j)=1.0+j*j;
{
L=L+LX/0.5; L=L-LX*3.0; L=LX*2.0+L;
U1=U1+UX*2.0; U1=U1-UX*3.0; U1=UX*2.0+U1;
}
SymmetricMatrix S(9);
for (i=1;i<=9;i++) for (j=1;j<=i;j++) S(i,j)=i*i+j;
{
SymmetricMatrix S1 = S;
S=S1+5.0;
S=S-3.0;
}
DiagonalMatrix D(9);
for (i=1;i<=9;i++) D(i,i)=S(i,i);
UpperTriangularMatrix U=L.t()*2.0;
{
U1=U1*2.0 - U; Print(U1);
L=L*2.0-D; U=U-D;
}
Matrix M=U+L; S=S*2.0; M=S-M; Print(M);
}
{
Tracer et1("Stage 3");
int i,j;
Matrix M(10,3), N(10,3);
for (i = 1; i<=10; i++) for (j = 1; j<=3; j++)
{ M(i,j) = 2*i-j; N(i,j) = i*j + 20; }
Matrix MN = M + N, M1;
M1 = M; M1 += N; M1 -= MN; Print(M1);
M1 = M; M1 += M1; M1 = M1 - M * 2; Print(M1);
M1 = M; M1 += N * 2; M1 -= (MN + N); Print(M1);
M1 = M; M1 -= M1; Print(M1);
M1 = M; M1 -= MN + M1; M1 += N + M; Print(M1);
M1 = M; M1 -= 5; M1 -= M; M1 *= 0.2; M1 = M1 + 1; Print(M1);
Matrix NT = N.t();
M1 = M; M1 *= NT; M1 -= M * N.t(); Print(M1);
M = M * M.t();
DiagonalMatrix D(10); D = 2;
M1 = M; M1 += D; M1 -= M; M1 = M1 - D; Print(M1);
M1 = M; M1 -= D; M1 -= M; M1 = M1 + D; Print(M1);
M1 = M; M1 *= D; M1 /= 2; M1 -= M; Print(M1);
SymmetricMatrix SM; SM << M;
// UpperTriangularMatrix SM; SM << M;
SM += 10; M1 = SM - M; M1 /=10; M1 = M1 - 1; Print(M1);
}
{
Tracer et1("Stage 4");
int i,j;
Matrix M(10,3), N(10,5);
for (i = 1; i<=10; i++) for (j = 1; j<=3; j++) M(i,j) = 2*i-j;
for (i = 1; i<=10; i++) for (j = 1; j<=5; j++) N(i,j) = i*j + 20;
Matrix M1;
M1 = M; M1 |= N; M1 &= N | M;
M1 -= (M | N) & (N | M); Print(M1);
M1 = M; M1 |= M1; M1 &= M1;
//.........这里部分代码省略.........
示例2: Transposer
void Transposer(const GenericMatrix& GM1, GenericMatrix&GM2)
{ GM2 = GM1.t(); }