Add conjugate gradient functionality to the library. Things to think about: - From #23, should we define natural gradients through passing through a `solve` argument? - Or should we have its own method on each object `linop.solve_cg(...)` - Or ??? Either way, we should ensure that the rest of the upstream code can be minimal/static as possible to facilitate cg solves in GP models, if desired.
Add conjugate gradient functionality to the library.
Things to think about:
solveargument?linop.solve_cg(...)Either way, we should ensure that the rest of the upstream code can be minimal/static as possible to facilitate cg solves in GP models, if desired.