On Fri, Dec 16, 2011 at 9:37 AM, Dave Nystrom <span dir="ltr"><<a href="mailto:dnystrom1@comcast.net">dnystrom1@comcast.net</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
I'm trying to figure out whether I can do a couple of things with petsc.<br>
<br>
1.  It looks like the preconditioning matrix can actually be different from<br>
the full problem matrix.  So I'm wondering if I could provide a different<br>
preconditioning matrix for my problem and then do an LU solve of the<br>
preconditioning matrix using the -pc_type lu as my preconditioner.<br></blockquote><div><br></div><div>Yes, that is what it is for.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

2.  When I build petsc, I use the --download-f-blas-lapack=yes option.  I'm<br>
wondering if petsc uses lapack under the hood or has the capability to use<br>
lapack under the hood when one uses the -pc_type lu option.  In particular,<br>
since my matrices are band matrices from doing a discretization on a 2d<br>
regular mesh, I'm wondering if the petsc lu solve has the ability to use the<br>
lapack band solver dgbsv or dgbsvx.  Or is it possible to use the lapack band<br>
solver through one of the external packages that petsc can interface with.<br>
I'm interested in this capability for smaller problem sizes that fit on a<br>
single node and that make sense.<br></blockquote><div><br></div><div>We do not have any banded matrix stuff. Its either dense or sparse right now.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

3.  I'm also wondering how I might be able to learn more about the petsc ilu<br>
capability.  My impression is that it does ilu(k) and I have tried it with<br>
k>0 but am wondering if one of the options might allow it to do ilut and<br>
whether as k gets big whether ilu(k) approximates lu.  I currently do not<br>
understand the petsc ilu well enough to know how much extra fill I get as I<br>
increase k and where that extra fill might be located for the case of a band<br>
matrix that one gets from discretization on a regular 2d mesh.<br></blockquote><div><br></div><div>We do not do ilu(dt). Its complicated, and we determined that it was not worth</div><div>the effort. You can get that from Hypre is you want. Certainly, for big enough</div>
<div>k, ilu(k) is lu but its a slow way to do it.</div><div><br></div><div>   Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Thanks,<br>
<br>
Dave<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener<br>