diff options
author | Jim Guggemos <jimg@google.com> | 2012-12-04 16:46:41 -0700 |
---|---|---|
committer | Jim Guggemos <jimg@google.com> | 2012-12-04 17:08:38 -0700 |
commit | b015e75e8c7ba1ab4ddb91e9372a57e76f3fd159 (patch) | |
tree | 54d1c7d66098154c1d7c5bd414394ef4cf255810 /doc/special_examples/Tutorial_sparse_example.cpp | |
parent | 63f67d748682b46d58be31235a0a2d64d81b998c (diff) | |
parent | c981c48f5bc9aefeffc0bcb0cc3934c2fae179dd (diff) | |
download | eigen-b015e75e8c7ba1ab4ddb91e9372a57e76f3fd159.tar.gz |
Merge 'goog/jb-mr1.1-dev' into platform/external/eigen.android-sdk-4.4.2_r1.0.1android-sdk-4.4.2_r1android-cts-4.4_r1android-4.4w_r1android-4.4_r0.7android-4.4.4_r2.0.1android-4.4.4_r2android-4.4.4_r1.0.1android-4.4.4_r1android-4.4.3_r1.1.0.1android-4.4.3_r1.1android-4.4.3_r1.0.1android-4.4.3_r1android-4.4.2_r2.0.1android-4.4.2_r2android-4.4.2_r1.0.1android-4.4.2_r1android-4.4.1_r1.0.1android-4.4.1_r1android-4.3_r3.1android-4.3_r3android-4.3.1_r1tools_r22.2tools_r22kitkat-wearkitkat-mr2.2-releasekitkat-mr2.1-releasekitkat-mr2-releasekitkat-mr1.1-releasekitkat-mr1-releasekitkat-devkitkat-cts-devjb-mr2.0.0-releasejb-mr2-devjb-mr1.1-dev-plus-aospjb-mr1-dev-plus-aospidea133-weekly-releaseidea133
Change-Id: Ic9004531328145ea36ba513bb96a23595427f6a4
Diffstat (limited to 'doc/special_examples/Tutorial_sparse_example.cpp')
-rw-r--r-- | doc/special_examples/Tutorial_sparse_example.cpp | 32 |
1 files changed, 32 insertions, 0 deletions
diff --git a/doc/special_examples/Tutorial_sparse_example.cpp b/doc/special_examples/Tutorial_sparse_example.cpp new file mode 100644 index 000000000..002f19f01 --- /dev/null +++ b/doc/special_examples/Tutorial_sparse_example.cpp @@ -0,0 +1,32 @@ +#include <Eigen/Sparse> +#include <vector> + +typedef Eigen::SparseMatrix<double> SpMat; // declares a column-major sparse matrix type of double +typedef Eigen::Triplet<double> T; + +void buildProblem(std::vector<T>& coefficients, Eigen::VectorXd& b, int n); +void saveAsBitmap(const Eigen::VectorXd& x, int n, const char* filename); + +int main(int argc, char** argv) +{ + int n = 300; // size of the image + int m = n*n; // number of unknows (=number of pixels) + + // Assembly: + std::vector<T> coefficients; // list of non-zeros coefficients + Eigen::VectorXd b(m); // the right hand side-vector resulting from the constraints + buildProblem(coefficients, b, n); + + SpMat A(m,m); + A.setFromTriplets(coefficients.begin(), coefficients.end()); + + // Solving: + Eigen::SimplicialCholesky<SpMat> chol(A); // performs a Cholesky factorization of A + Eigen::VectorXd x = chol.solve(b); // use the factorization to solve for the given right hand side + + // Export the result to a file: + saveAsBitmap(x, n, argv[1]); + + return 0; +} + |