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

[Stefano Zampini]

Yann

 you want to have the number of columns equal to the number of rhs. Not the
number of rows.
This is also consistent in terms of logical layouts. Simply create and dump
the matrix this way and you can read in parallel splitting yhe rows ( ie
the logical size per process of the vectors)

Il Gio 17 Set 2020, 18:09 Pierre Jolivet <pierre.jolivet at enseeiht.fr> ha
scritto:

> 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
> .
> 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
> 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/a5eff054/attachment-0001.html>