[petsc-users] Question about speed of MatMatSolveTranspose with MUMPS

Junchao Zhang junchao.zhang at gmail.com
Thu Aug 4 17:21:02 CDT 2022

It seems yes, from mump manual

ICNTL(9) computes the solution using A or AT

Phase: accessed by the host during the solve phase.

Possible values :

1    : AX = B is solved.

!= 1 : ATX = B is solved.

On Thu, Aug 4, 2022 at 4:20 PM Barry Smith <bsmith at petsc.dev> wrote:

>   Yes, it appears no MUMPS matmatsolvetranspose is wired in. Does MUMPS
> provide such a function? If so it can be added with a few lines of code.
>   Barry
> On Aug 4, 2022, at 3:55 PM, Lucas Banting <bantingl at myumanitoba.ca> wrote:
> Hello,
> I am using MUMPS to solve a system with multiple right hand sides.
> I have noticed that the MatMatSolve function appears to be quite a bit
> faster than the MatMatSolveTranspose function.
> I was trying to look through the petsc source to see where the
> MatMatSolveTranspose_MUMPS function was but couldn't find it.
> Is PETSc just doing MatSolveTranspose_MUMPS for each right hand side when
> doing MatMatSolveTranspose with MUMPS?
> If an accelerated version of MatMatSolveTranspose existed, it would be
> very useful for my application.
> Regards,
> Lucas
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20220804/dab9c6a5/attachment.html>

More information about the petsc-users mailing list