diff options
author | Jitse Niesen <jitse@maths.leeds.ac.uk> | 2010-03-31 11:59:11 +0100 |
---|---|---|
committer | Jitse Niesen <jitse@maths.leeds.ac.uk> | 2010-03-31 11:59:11 +0100 |
commit | 1b3f7f2feef34d0506e66e7514130a3994589be8 (patch) | |
tree | edcd3881af84b59f23eb608f4aeb1a5ee97255ab /doc/snippets/EigenSolver_pseudoEigenvectors.cpp | |
parent | 338ec0390fc89d9f202d2eb2a41725b7fc738a0a (diff) |
Extend documentation and add examples for EigenSolver class.
Diffstat (limited to 'doc/snippets/EigenSolver_pseudoEigenvectors.cpp')
-rw-r--r-- | doc/snippets/EigenSolver_pseudoEigenvectors.cpp | 9 |
1 files changed, 9 insertions, 0 deletions
diff --git a/doc/snippets/EigenSolver_pseudoEigenvectors.cpp b/doc/snippets/EigenSolver_pseudoEigenvectors.cpp new file mode 100644 index 000000000..85e2569df --- /dev/null +++ b/doc/snippets/EigenSolver_pseudoEigenvectors.cpp @@ -0,0 +1,9 @@ +MatrixXd A = MatrixXd::Random(6,6); +cout << "Here is a random 6x6 matrix, A:" << endl << A << endl << endl; + +EigenSolver<MatrixXd> es(A); +MatrixXd D = es.pseudoEigenvalueMatrix(); +MatrixXd V = es.pseudoEigenvectors(); +cout << "The pseudo-eigenvalue matrix D is:" << endl << D << endl; +cout << "The pseudo-eigenvector matrix V is:" << endl << V << endl; +cout << "Finally, V * D * V^(-1) = " << endl << V * D * V.inverse() << endl; |