diff options
author | Jitse Niesen <jitse@maths.leeds.ac.uk> | 2010-04-12 18:14:32 +0100 |
---|---|---|
committer | Jitse Niesen <jitse@maths.leeds.ac.uk> | 2010-04-12 18:14:32 +0100 |
commit | 73d3a27667aff4506a20693140dc603110e48cbc (patch) | |
tree | 070dbcddaa7c9f5a6dee2a20400c3d361314573b /doc/snippets/RealSchur_RealSchur_MatrixType.cpp | |
parent | 872df22ca4cb5a44ad2e0be5c16e50f8877ebb14 (diff) |
RealSchur: Make sure zeros are really zero (cont'd); add default ctor, docs.
Diffstat (limited to 'doc/snippets/RealSchur_RealSchur_MatrixType.cpp')
-rw-r--r-- | doc/snippets/RealSchur_RealSchur_MatrixType.cpp | 10 |
1 files changed, 10 insertions, 0 deletions
diff --git a/doc/snippets/RealSchur_RealSchur_MatrixType.cpp b/doc/snippets/RealSchur_RealSchur_MatrixType.cpp new file mode 100644 index 000000000..a5530dcc8 --- /dev/null +++ b/doc/snippets/RealSchur_RealSchur_MatrixType.cpp @@ -0,0 +1,10 @@ +MatrixXd A = MatrixXd::Random(6,6); +cout << "Here is a random 6x6 matrix, A:" << endl << A << endl << endl; + +RealSchur<MatrixXd> schur(A); +cout << "The orthogonal matrix U is:" << endl << schur.matrixU() << endl; +cout << "The quasi-triangular matrix T is:" << endl << schur.matrixT() << endl << endl; + +MatrixXd U = schur.matrixU(); +MatrixXd T = schur.matrixT(); +cout << "U * T * U^T = " << endl << U * T * U.transpose() << endl; |