Document Tridiagonalization class, remove unused types.

1 files changed, 11 insertions, 0 deletions

diff --git a/doc/snippets/Tridiagonalization_Tridiagonalization_MatrixType.cpp b/doc/snippets/Tridiagonalization_Tridiagonalization_MatrixType.cpp
new file mode 100644
index 000000000..109650041
--- /dev/null
+++ b/doc/snippets/Tridiagonalization_Tridiagonalization_MatrixType.cpp
@@ -0,0 +1,11 @@
+MatrixXd X = MatrixXd::Random(5,5);
+MatrixXd A = X + X.transpose();
+cout << "Here is a random symmetric 5x5 matrix:" << endl << A << endl << endl;
+
+Tridiagonalization<MatrixXd> triOfA(A);
+cout << "The orthogonal matrix Q is:" << endl << triOfA.matrixQ() << endl;
+cout << "The tridiagonal matrix T is:" << endl << triOfA.matrixT() << endl << endl;
+
+MatrixXd Q = triOfA.matrixQ();
+MatrixXd T = triOfA.matrixT();
+cout << "Q * T * Q^T = " << endl << Q * T * Q.transpose() << endl;