diff options
author | 2009-06-25 00:57:51 +0200 | |
---|---|---|
committer | 2009-06-25 00:57:51 +0200 | |
commit | 903acf0d5c2855fbcd0cfd6913e2cd8312f06eb1 (patch) | |
tree | 5edbfc60adce6c2bdf17dbe32b321d014fb9b8f0 /doc/snippets/PartialLU_solve.cpp | |
parent | 03ad303d1417e53b8bcec8dd8d8e404cd4e13af4 (diff) |
add missing code snippets for newer Matrix methods and PartialLU::solve()
Diffstat (limited to 'doc/snippets/PartialLU_solve.cpp')
-rw-r--r-- | doc/snippets/PartialLU_solve.cpp | 8 |
1 files changed, 8 insertions, 0 deletions
diff --git a/doc/snippets/PartialLU_solve.cpp b/doc/snippets/PartialLU_solve.cpp new file mode 100644 index 000000000..12441437f --- /dev/null +++ b/doc/snippets/PartialLU_solve.cpp @@ -0,0 +1,8 @@ +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; +if(A.lu().solve(B, &X)) +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; |