Unified Interface for Linear Solving

I opened an issue to discuss this with IterativeSolvers:

https://github.com/JuliaMath/IterativeSolvers.jl/issues/99

Also, I put out a blog post explaining the benefits this kind of interface would yield: