aboutsummaryrefslogtreecommitdiffhomepage
path: root/doc/snippets/Tridiagonalization_packedMatrix.cpp
diff options
context:
space:
mode:
authorGravatar Jitse Niesen <jitse@maths.leeds.ac.uk>2010-05-01 17:52:16 +0100
committerGravatar Jitse Niesen <jitse@maths.leeds.ac.uk>2010-05-01 17:52:16 +0100
commitafed0ef90d8f6c6c62dc5b1173bbbc63bee3c5c1 (patch)
tree78d3b648ce1ac4c65c8df49226688cda6ddbbc2c /doc/snippets/Tridiagonalization_packedMatrix.cpp
parentd9c1224133153045fb92df351bdd3660a801f8d7 (diff)
Document Tridiagonalization class, remove unused types.
Diffstat (limited to 'doc/snippets/Tridiagonalization_packedMatrix.cpp')
-rw-r--r--doc/snippets/Tridiagonalization_packedMatrix.cpp8
1 files changed, 8 insertions, 0 deletions
diff --git a/doc/snippets/Tridiagonalization_packedMatrix.cpp b/doc/snippets/Tridiagonalization_packedMatrix.cpp
new file mode 100644
index 000000000..0f55d0c28
--- /dev/null
+++ b/doc/snippets/Tridiagonalization_packedMatrix.cpp
@@ -0,0 +1,8 @@
+Matrix4d X = Matrix4d::Random(4,4);
+Matrix4d A = X + X.transpose();
+cout << "Here is a random symmetric 4x4 matrix:" << endl << A << endl;
+Tridiagonalization<Matrix4d> triOfA(A);
+Matrix4d pm = triOfA.packedMatrix();
+cout << "The packed matrix M is:" << endl << pm << endl;
+cout << "The diagonal and subdiagonal corresponds to the matrix T, which is:"
+ << endl << triOfA.matrixT() << endl;