
<html><head><meta http-equiv="Content-Type" content="text/html; charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">Hello Yann,<div class="">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. <a href="https://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/KSP/KSPMatSolve.html" class="">https://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/KSP/KSPMatSolve.html</a>.</div><div class="">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.</div><div class="">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: <a href="http://jolivet.perso.enseeiht.fr/article.pdf" class="">http://jolivet.perso.enseeiht.fr/article.pdf</a></div><div class="">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.</div><div class="">Thanks,</div><div class="">Pierre</div></body></html>