[petsc-users] Solving a linear system with sparse matrices

Stefano Zampini stefano.zampini at gmail.com
Sat May 21 13:42:48 CDT 2022

MatMatSolve  with sparse rhs works with MUMPS

On Sat, May 21, 2022, 19:41 Mateo José Moinelo <mateojmoinelo at gmail.com>

> Hi. I am working on a Petsc program, and I have a problem. I have two
> sparse matrices A and B, and I want to compute inv(A)*B. I understand that
> computing the inverse of a matrix directly is not effective, and that in
> this case, the best way to do the operation is instead solving the system
> A*X = B, being X the result. I reviewed the documentation, and found some
> interesting options, like the function MatMatSolve, or using a KSP.
> The problem I found is that these options seem to work only for dense
> matrices, and I do not know how to convert my sparce matrices to dense.
> What can I do in this case?.
> Thanks in advance.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20220521/e7ddd1f0/attachment.html>

More information about the petsc-users mailing list