compatibility with MPI
Hey, I just wanted to ask about the feasibility of using this library as a replacement for a large-scale parallel solver like PETSc or hypre. Looking at the code for cg.jl it seems like all that is required would be to 1) implement a KrylovSubspace type with a nextvec method, and 2) implement a parallel version of dot. Has anyone used this library with MPI.jl or DistributedArrays.jl?
This should probably go on https://discourse.julialang.org/ instead of the issue tracker but it should be possible to do something along the lines of what you describe. See this example https://github.com/JuliaParallel/Elemental.jl#truncated-svd for an approach to this for a procedure very similar to the procedures in this package.