<div class="gmail_quote">On Mon, Aug 29, 2011 at 14:30, John Fettig <span dir="ltr"><<a href="mailto:john.fettig@gmail.com">john.fettig@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div id=":1t6">Here's the patch.  I left in the serial code that creates the matrices<br>
in PETSc rather than letting ML take care of it, just in case you<br>
wanted that for any reason.</div></blockquote><div><br></div><div>Thanks. I pushed a modified version of this patch. Run with -pc_ml_reuse_interpolation and SAME_NONZERO_PATTERN to activate. The savings in setup time are sometimes significant, so perhaps we should make it the default. We usually choose safe defaults, but wanting to reuse interpolants might be sufficiently common to warrant. Thoughts?</div>
<div><br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div id=":1t6">  One thing that I notice is that when I<br>
use UMFPACK for the coarse grid solve, the preconditioner is different<br>
(e.g. minres takes more iterations) when I do:<br>
<div class="im"><br>
KSPSetOperators(gridctx[fine_level].ksp,gridctx[level].A,gridctx[fine_level].A,SAME_NONZERO_PATTERN);<br>
<br>
</div>versus<br>
<br>
KSPSetOperators(gridctx[fine_level].ksp,gridctx[level].A,gridctx[fine_level].A,DIFFERENT_NONZERO_PATTERN);<br>
<br>
This shouldn't be the case, should it?  It doesn't do this with the<br>
PETSc builtin lu, so I suspect a bug in the UMFPACK interface.</div></blockquote></div><br><div>Not if the symbolic factorization is still valid (e.g. nothing caused the nonzero pattern to change). Do you have a test case for this?</div>