| Gael Guennebaud | c6da881 | 2013-10-16 15:25:39 +0200 | [diff] [blame] | 1 | int n = 10000; |
| 2 | VectorXd x(n), b(n); | ||||
| 3 | SparseMatrix<double> A(n,n); | ||||
| 4 | /* ... fill A and b ... */ | ||||
| 5 | BiCGSTAB<SparseMatrix<double> > solver; | ||||
| 6 | solver.compute(A); | ||||
| 7 | x = solver.solve(b); | ||||
| 8 | std::cout << "#iterations: " << solver.iterations() << std::endl; | ||||
| 9 | std::cout << "estimated error: " << solver.error() << std::endl; | ||||
| 10 | /* ... update b ... */ | ||||
| Aaron Franke | 5c22c7a | 2020-03-21 21:30:06 -0400 | [diff] [blame] | 11 | x = solver.solve(b); // solve again |