<div dir="ltr">True, also I should have read Chapter 4 of the manual before asking the question. I will try the various options there and see if I get any improvement.<div><br></div><div>One thing I'm very confused about though: as I mentioned in the first post, I am solving the matrix multiple times and changing the elements, but not the nonzero structure in between solves. We've established that the numerical LU is done on the first solve only. Also, I forgot to mention that I am *not* calling KSPSetOperator before the subsequent solves, which is why the subsequent solves are so fast. </div>
<div><br></div><div>According to section 4.2 of the manual, this is actually not allowed if I change anything about the matrix (which I am doing since I didn't know about this). Then why are my solves still getting the correct results that I am expecting? </div>
</div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Jan 29, 2014 at 4:54 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="im">David Liu <<a href="mailto:daveliu@mit.edu">daveliu@mit.edu</a>> writes:<br>
> Not too sure if there's anything I can do about this. Also, I am definately<br>
> willing to use other preconditioners. If direct solvers are expensive, what<br>
> are the cheaper alternatives?<br>
<br>
</div>This is a vast question, a bit like asking what kind of "math" is<br>
important for "physics".  If you want help, you'll have to explain more<br>
about your problem, discretization, and regime.<br>
</blockquote></div><br></div>