diff options
Diffstat (limited to 'doc/snippets/BiCGSTAB_simple.cpp')
| -rw-r--r-- | doc/snippets/BiCGSTAB_simple.cpp | 11 |
1 files changed, 11 insertions, 0 deletions
diff --git a/doc/snippets/BiCGSTAB_simple.cpp b/doc/snippets/BiCGSTAB_simple.cpp new file mode 100644 index 000000000..5520f4f1f --- /dev/null +++ b/doc/snippets/BiCGSTAB_simple.cpp @@ -0,0 +1,11 @@ + int n = 10000; + VectorXd x(n), b(n); + SparseMatrix<double> A(n,n); + /* ... fill A and b ... */ + BiCGSTAB<SparseMatrix<double> > solver; + solver.compute(A); + x = solver.solve(b); + std::cout << "#iterations: " << solver.iterations() << std::endl; + std::cout << "estimated error: " << solver.error() << std::endl; + /* ... update b ... */ + x = solver.solve(b); // solve again
\ No newline at end of file |