<div dir="ltr">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)? <div><br></div><div>Thanks</div>
<div>Zin<br><div><br></div></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Fri, Jul 25, 2014 at 4:36 PM, Jed Brown <span dir="ltr"><<a href="mailto:jed@jedbrown.org" target="_blank">jed@jedbrown.org</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="">Zin Lin <<a href="mailto:zinlin.zinlin@gmail.com">zinlin.zinlin@gmail.com</a>> writes:<br>
<br>
> Hi,<br>
><br>
> I am having trouble calling KSPSolveTranspose under MATSOLVERPASTIX.<br>
> The following is a very simple C code solving A^T x = b to illustrate the<br>
> issue. The matrix A and vector b are built following the PETSC tutorial.<br>
> KSPSolve works fine but when it comes to KSPSolveTranspose, it spits out<br>
> the error:<br>
><br>
> [0]PETSC ERROR: --------------------- Error Message<br>
> ------------------------------------<br>
> [0]PETSC ERROR: No support for this operation for this object type!<br>
> [0]PETSC ERROR: Matrix type mpiaij!<br>
><br>
> The code runs fine if I change the solver package MATSOLVERPASTIX to<br>
> MATSOLVERMUMPS.<br>
> However, I need the PASTIX solver and the transpose solution for solving a<br>
> large problem.<br>
> Could you please let me know what is wrong with the PASTIX solver or with<br>
> my usage?<br>
<br>
</div>The contributed PaStiX interface does not include<br>
MatSolveTranspose_PaStiX.  If there is a way to do this in PaStiX, it<br>
should be added to the interface.<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div dir="ltr">Zin Lin<br><br></div>
</div>