[petsc-users] Problem with MATSOLVERPASTIX using KSPSolveTranspose

[Jed Brown]

Zin Lin <zinlin.zinlin at gmail.com> writes:

> Hi,
>
> I am having trouble calling KSPSolveTranspose under MATSOLVERPASTIX.
> The following is a very simple C code solving A^T x = b to illustrate the
> issue. The matrix A and vector b are built following the PETSC tutorial.
> KSPSolve works fine but when it comes to KSPSolveTranspose, it spits out
> the error:
>
> [0]PETSC ERROR: --------------------- Error Message
> ------------------------------------
> [0]PETSC ERROR: No support for this operation for this object type!
> [0]PETSC ERROR: Matrix type mpiaij!
>
> The code runs fine if I change the solver package MATSOLVERPASTIX to
> MATSOLVERMUMPS.
> However, I need the PASTIX solver and the transpose solution for solving a
> large problem.
> Could you please let me know what is wrong with the PASTIX solver or with
> my usage?

The contributed PaStiX interface does not include
MatSolveTranspose_PaStiX.  If there is a way to do this in PaStiX, it
should be added to the interface.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 818 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20140725/d72363bb/attachment.pgp>