diff options
author | Benoit Jacob <jacob.benoit.1@gmail.com> | 2009-06-25 00:57:51 +0200 |
---|---|---|
committer | Benoit Jacob <jacob.benoit.1@gmail.com> | 2009-06-25 00:57:51 +0200 |
commit | 903acf0d5c2855fbcd0cfd6913e2cd8312f06eb1 (patch) | |
tree | 5edbfc60adce6c2bdf17dbe32b321d014fb9b8f0 /doc | |
parent | 03ad303d1417e53b8bcec8dd8d8e404cd4e13af4 (diff) |
add missing code snippets for newer Matrix methods and PartialLU::solve()
Diffstat (limited to 'doc')
-rw-r--r-- | doc/snippets/Matrix_setConstant_int.cpp | 3 | ||||
-rw-r--r-- | doc/snippets/Matrix_setConstant_int_int.cpp | 3 | ||||
-rw-r--r-- | doc/snippets/Matrix_setIdentity_int_int.cpp | 3 | ||||
-rw-r--r-- | doc/snippets/Matrix_setOnes_int.cpp | 3 | ||||
-rw-r--r-- | doc/snippets/Matrix_setOnes_int_int.cpp | 3 | ||||
-rw-r--r-- | doc/snippets/Matrix_setRandom_int.cpp | 3 | ||||
-rw-r--r-- | doc/snippets/Matrix_setRandom_int_int.cpp | 3 | ||||
-rw-r--r-- | doc/snippets/Matrix_setZero_int.cpp | 3 | ||||
-rw-r--r-- | doc/snippets/Matrix_setZero_int_int.cpp | 3 | ||||
-rw-r--r-- | doc/snippets/PartialLU_solve.cpp | 8 |
10 files changed, 35 insertions, 0 deletions
diff --git a/doc/snippets/Matrix_setConstant_int.cpp b/doc/snippets/Matrix_setConstant_int.cpp new file mode 100644 index 000000000..ff5a86c98 --- /dev/null +++ b/doc/snippets/Matrix_setConstant_int.cpp @@ -0,0 +1,3 @@ +VectorXf v; +v.setConstant(3, 5); +cout << v << endl; diff --git a/doc/snippets/Matrix_setConstant_int_int.cpp b/doc/snippets/Matrix_setConstant_int_int.cpp new file mode 100644 index 000000000..32b950cfd --- /dev/null +++ b/doc/snippets/Matrix_setConstant_int_int.cpp @@ -0,0 +1,3 @@ +MatrixXf m; +m.setConstant(3, 3, 5); +cout << m << endl; diff --git a/doc/snippets/Matrix_setIdentity_int_int.cpp b/doc/snippets/Matrix_setIdentity_int_int.cpp new file mode 100644 index 000000000..a65967199 --- /dev/null +++ b/doc/snippets/Matrix_setIdentity_int_int.cpp @@ -0,0 +1,3 @@ +MatrixXf m; +m.setIdentity(3, 3); +cout << m << endl; diff --git a/doc/snippets/Matrix_setOnes_int.cpp b/doc/snippets/Matrix_setOnes_int.cpp new file mode 100644 index 000000000..752cb35b2 --- /dev/null +++ b/doc/snippets/Matrix_setOnes_int.cpp @@ -0,0 +1,3 @@ +VectorXf v; +v.setOnes(3); +cout << v << endl; diff --git a/doc/snippets/Matrix_setOnes_int_int.cpp b/doc/snippets/Matrix_setOnes_int_int.cpp new file mode 100644 index 000000000..1ffb66bbd --- /dev/null +++ b/doc/snippets/Matrix_setOnes_int_int.cpp @@ -0,0 +1,3 @@ +MatrixXf m; +m.setOnes(3, 3); +cout << m << endl; diff --git a/doc/snippets/Matrix_setRandom_int.cpp b/doc/snippets/Matrix_setRandom_int.cpp new file mode 100644 index 000000000..e160dd7df --- /dev/null +++ b/doc/snippets/Matrix_setRandom_int.cpp @@ -0,0 +1,3 @@ +VectorXf v; +v.setRandom(3); +cout << v << endl; diff --git a/doc/snippets/Matrix_setRandom_int_int.cpp b/doc/snippets/Matrix_setRandom_int_int.cpp new file mode 100644 index 000000000..80cda11d7 --- /dev/null +++ b/doc/snippets/Matrix_setRandom_int_int.cpp @@ -0,0 +1,3 @@ +MatrixXf m; +m.setRandom(3, 3); +cout << m << endl; diff --git a/doc/snippets/Matrix_setZero_int.cpp b/doc/snippets/Matrix_setZero_int.cpp new file mode 100644 index 000000000..0fb16c1f3 --- /dev/null +++ b/doc/snippets/Matrix_setZero_int.cpp @@ -0,0 +1,3 @@ +VectorXf v; +v.setZero(3); +cout << v << endl; diff --git a/doc/snippets/Matrix_setZero_int_int.cpp b/doc/snippets/Matrix_setZero_int_int.cpp new file mode 100644 index 000000000..ad883b916 --- /dev/null +++ b/doc/snippets/Matrix_setZero_int_int.cpp @@ -0,0 +1,3 @@ +MatrixXf m; +m.setZero(3, 3); +cout << m << endl; 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; |