| 1234567891011 |
- 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
|
