search

JuliaLinearAlgebra / IterativeSolvers.jl

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

compatibility with MPI #102

Open nbren12 opened 7 years ago

nbren12 commented 7 years ago

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?

andreasnoack commented 7 years ago

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.