aboutsummaryrefslogtreecommitdiffhomepage
path: root/doc/snippets/PartialPivLU_solve.cpp
diff options
context:
space:
mode:
authorGravatar Jitse Niesen <jitse@maths.leeds.ac.uk>2009-12-08 15:12:27 +0000
committerGravatar Jitse Niesen <jitse@maths.leeds.ac.uk>2009-12-08 15:12:27 +0000
commit8bfa354ee3267691a7aec6fe576847296e074c90 (patch)
treeaddb0b77e7546ca1da4984f15b3d6bf52d82f1ce /doc/snippets/PartialPivLU_solve.cpp
parenta682a0eeb1246725a0f6bbbe02eeedf81d26b22b (diff)
Documentation clean up.
* remove non-existant reference to CwiseAll * define \householder_module (used in HouseholderSequence.h) * update I01_TopicLazyEvaluation.dox - Product is now called GeneralProduct * remove reference to list of examples which was deleted ages ago * rename PartialLU_solve.cpp snippet to PartialPivLU_solve.cpp
Diffstat (limited to 'doc/snippets/PartialPivLU_solve.cpp')
-rw-r--r--doc/snippets/PartialPivLU_solve.cpp7
1 files changed, 7 insertions, 0 deletions
diff --git a/doc/snippets/PartialPivLU_solve.cpp b/doc/snippets/PartialPivLU_solve.cpp
new file mode 100644
index 000000000..fa3570ab8
--- /dev/null
+++ b/doc/snippets/PartialPivLU_solve.cpp
@@ -0,0 +1,7 @@
+MatrixXd A = MatrixXd::Random(3,3);
+MatrixXd B = MatrixXd::Random(3,2);
+cout << "Here is the invertible matrix A:" << endl << A << endl;
+cout << "Here is the matrix B:" << endl << B << endl;
+MatrixXd X = A.lu().solve(B);
+cout << "Here is the (unique) solution X to the equation AX=B:" << endl << X << endl;
+cout << "Relative error: " << (A*X-B).norm() / B.norm() << endl;