[petsc-users] 1D parallel split of dense matrix

[Pierre Jolivet]

Hello Yann,
This is probably not fully answering your question, but the proper way to solve a system with N RHS is _not_ to use KSPSolve(), but instead KSPMatSolve(), cf. https://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/KSP/KSPMatSolve.html <https://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/KSP/KSPMatSolve.html>.
If you are tracking master (from the GitLab repository), it’s available out of the box. If you are using the release tarballs, it will be available in 3.14.0 scheduled to be released in a couple of days.
If you want to know more about the current status of block solvers in PETSc, please feel free to have a look at this preprint: http://jolivet.perso.enseeiht.fr/article.pdf <http://jolivet.perso.enseeiht.fr/article.pdf>
If you are using a specific PC which is not “multi-RHS ready”, see the list at the top of page 5, please let me know and I’ll tell you how easy to support it.
Thanks,
Pierre
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20200917/713924a3/attachment.html>