diff options
Diffstat (limited to 'doc/snippets/Tridiagonalization_compute.cpp')
-rw-r--r-- | doc/snippets/Tridiagonalization_compute.cpp | 9 |
1 files changed, 9 insertions, 0 deletions
diff --git a/doc/snippets/Tridiagonalization_compute.cpp b/doc/snippets/Tridiagonalization_compute.cpp new file mode 100644 index 000000000..0062a99e8 --- /dev/null +++ b/doc/snippets/Tridiagonalization_compute.cpp @@ -0,0 +1,9 @@ +Tridiagonalization<MatrixXf> tri; +MatrixXf X = MatrixXf::Random(4,4); +MatrixXf A = X + X.transpose(); +tri.compute(A); +cout << "The matrix T in the tridiagonal decomposition of A is: " << endl; +cout << tri.matrixT() << endl; +tri.compute(2*A); // re-use tri to compute eigenvalues of 2A +cout << "The matrix T in the tridiagonal decomposition of 2A is: " << endl; +cout << tri.matrixT() << endl; |