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 <[email protected]> wrote: > Zin Lin <[email protected]> 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
