From 73d3a27667aff4506a20693140dc603110e48cbc Mon Sep 17 00:00:00 2001 From: Jitse Niesen Date: Mon, 12 Apr 2010 18:14:32 +0100 Subject: RealSchur: Make sure zeros are really zero (cont'd); add default ctor, docs. --- doc/snippets/RealSchur_RealSchur_MatrixType.cpp | 10 ++++++++++ doc/snippets/RealSchur_compute.cpp | 6 ++++++ 2 files changed, 16 insertions(+) create mode 100644 doc/snippets/RealSchur_RealSchur_MatrixType.cpp create mode 100644 doc/snippets/RealSchur_compute.cpp (limited to 'doc') 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 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; diff --git a/doc/snippets/RealSchur_compute.cpp b/doc/snippets/RealSchur_compute.cpp new file mode 100644 index 000000000..4dcfaf0f2 --- /dev/null +++ b/doc/snippets/RealSchur_compute.cpp @@ -0,0 +1,6 @@ +MatrixXf A = MatrixXf::Random(4,4); +RealSchur schur(4); +schur.compute(A); +cout << "The matrix T in the decomposition of A is:" << endl << schur.matrixT() << endl; +schur.compute(A.inverse()); +cout << "The matrix T in the decomposition of A^(-1) is:" << endl << schur.matrixT() << endl; -- cgit v1.2.3