search

JuliaLinearAlgebra / IterativeSolvers.jl

Iterative algorithms for solving linear systems, eigensystems, and singular value problems
MIT License
399 stars 105 forks source link

No x0 argument for the out-of-place solvers #235

Open antoine-levitt opened 5 years ago

antoine-levitt commented 5 years ago

It'd be nice to be able to do cg(A,b,x0=x0).

haampie commented 5 years ago

Alternatively cg!(copy(x0), A, b) is quite short as well. I'm not sure I like x0 as a kwarg because it's not an obvious name for an initial guess.