[petsc-users] Problem with MATSOLVERPASTIX using KSPSolveTranspose

[Zin Lin]

Is there any chance that I could implement the transpose solve with
PCApplyTranspose etc on my own (rather than using a new KSP and PC context
which re-factorizes the transpose matrix)?

Thanks
Zin



On Fri, Jul 25, 2014 at 4:36 PM, Jed Brown <jed at jedbrown.org> wrote:

> 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.
>



-- 
Zin Lin
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20140807/8ec0e6f1/attachment.html>