<div dir="auto">MatMatSolve  with sparse rhs works with MUMPS </div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, May 21, 2022, 19:41 Mateo José Moinelo <<a href="mailto:mateojmoinelo@gmail.com">mateojmoinelo@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">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.<div><br></div><div>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?.</div><div><br></div><div>Thanks in advance.</div></div>
</blockquote></div>