Hi all --
I second the thanks to AlephZero for useful replies! :)
I am also in a situation where it's desirable to solve a large Ax = b problem, where A is large (N x N with N ~ 10^6) and sparse (nonzero entries are typically < 10 N in number). The thing is, my A has some extra known...