diff options
author | Jitse Niesen <jitse@maths.leeds.ac.uk> | 2009-08-22 20:12:47 +0100 |
---|---|---|
committer | Jitse Niesen <jitse@maths.leeds.ac.uk> | 2009-08-22 20:12:47 +0100 |
commit | 90735b6a9cc4c91eb41ef1d2b6828f15a7f4fba3 (patch) | |
tree | f466895074e2e67fa734fe975da6fe6a984b37db /doc/snippets/Tutorial_solve_reuse_decomposition.cpp | |
parent | 37dede6077250156be084e55e629c68535456d2a (diff) |
Rewrite tutorial section on solving linear systems
Diffstat (limited to 'doc/snippets/Tutorial_solve_reuse_decomposition.cpp')
-rw-r--r-- | doc/snippets/Tutorial_solve_reuse_decomposition.cpp | 13 |
1 files changed, 13 insertions, 0 deletions
diff --git a/doc/snippets/Tutorial_solve_reuse_decomposition.cpp b/doc/snippets/Tutorial_solve_reuse_decomposition.cpp new file mode 100644 index 000000000..b4112adc4 --- /dev/null +++ b/doc/snippets/Tutorial_solve_reuse_decomposition.cpp @@ -0,0 +1,13 @@ +Matrix3f A(3,3); +A << 1,2,3, 4,5,6, 7,8,10; +PartialLU<Matrix3f> luOfA(A); // compute LU decomposition of A +Vector3f b; +b << 3,3,4; +Vector3f x; +luOfA.solve(b, &x); +cout << "The solution with right-hand side (3,3,4) is:" << endl; +cout << x << endl; +b << 1,1,1; +luOfA.solve(b, &x); +cout << "The solution with right-hand side (1,1,1) is:" << endl; +cout << x << endl; |