Hi, Hong and Jed,<br><br>Thanks a lot as usual. <br><br>Best wishes<br>Zhengyong <br><br><div class="gmail_quote">On Tue, Jun 5, 2012 at 8:01 PM, Jed Brown <span dir="ltr"><<a href="mailto:jedbrown@mcs.anl.gov" target="_blank">jedbrown@mcs.anl.gov</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="gmail_quote"><div class="im">On Tue, Jun 5, 2012 at 12:48 PM, RenZhengYong <span dir="ltr"><<a href="mailto:renzhengyong@gmail.com" target="_blank">renzhengyong@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Dear Petsc users, <br><font><br style="font-family:georgia,serif"><span style="font-family:georgia,serif">My purpose is to solve AmatX=b, A is a symmetric indefinite matrix. </span><br style="font-family:georgia,serif"><span style="font-family:georgia,serif">Here, I read the manual about the usage of preconditioned GMRES method:</span><br style="font-family:georgia,serif">


<span style="font-family:georgia,serif">Given a matrix Pmat which is a good approximation to Amat, and call</span><br style="font-family:georgia,serif"></font><pre style="font-family:georgia,serif"><font>PetscErrorCode  KSPSetOperators(KSP ksp,Mat Amat,Mat Pmat,MatStructure flag).<br>


<br>My question is: If I have an explicit good approximation (T) to the inverse of the Amat (that is T=Pmat^{_1}), <br>which call will be used to tell the KSP of having this explicit preconditioner. </font></pre></blockquote>

</div><div><a href="http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/PC/PCMAT.html" target="_blank">http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/PC/PCMAT.html</a></div></div>
</blockquote></div><br><br clear="all"><br>-- <br>Zhengyong Ren<br>AUG Group, Institute of Geophysics<br>Department of Geosciences, ETH Zurich<br>NO H 47 Sonneggstrasse 5<br>CH-8092, Zürich, Switzerland<br>Tel: +41 44 633 37561<br>
e-mail: <a href="mailto:zhengyong.ren@aug.ig.erdw.ethz.ch">zhengyong.ren@aug.ig.erdw.ethz.ch</a><br>Gmail: <a href="mailto:renzhengyong@gmail.com">renzhengyong@gmail.com</a><br>