diff options
| author | Stanislaw Halik <sthalik@misaki.pl> | 2017-03-25 14:17:07 +0100 |
|---|---|---|
| committer | Stanislaw Halik <sthalik@misaki.pl> | 2017-03-25 14:17:07 +0100 |
| commit | 35f7829af10c61e33dd2e2a7a015058e11a11ea0 (patch) | |
| tree | 7135010dcf8fd0a49f3020d52112709bcb883bd6 /eigen/doc/snippets/BiCGSTAB_step_by_step.cpp | |
| parent | 6e8724193e40a932faf9064b664b529e7301c578 (diff) | |
update
Diffstat (limited to 'eigen/doc/snippets/BiCGSTAB_step_by_step.cpp')
| -rw-r--r-- | eigen/doc/snippets/BiCGSTAB_step_by_step.cpp | 14 |
1 files changed, 14 insertions, 0 deletions
diff --git a/eigen/doc/snippets/BiCGSTAB_step_by_step.cpp b/eigen/doc/snippets/BiCGSTAB_step_by_step.cpp new file mode 100644 index 0000000..06147bb --- /dev/null +++ b/eigen/doc/snippets/BiCGSTAB_step_by_step.cpp @@ -0,0 +1,14 @@ + int n = 10000; + VectorXd x(n), b(n); + SparseMatrix<double> A(n,n); + /* ... fill A and b ... */ + BiCGSTAB<SparseMatrix<double> > solver(A); + // start from a random solution + x = VectorXd::Random(n); + solver.setMaxIterations(1); + int i = 0; + do { + x = solver.solveWithGuess(b,x); + std::cout << i << " : " << solver.error() << std::endl; + ++i; + } while (solver.info()!=Success && i<100);
\ No newline at end of file |