I'll try tomorrow - thanks Jed!<div><br></div><div>Derek<br><br>On Saturday, November 29, 2014, Jed Brown <<a href="mailto:jed@jedbrown.org">jed@jedbrown.org</a>> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Derek Gaston <<a href="javascript:;" onclick="_e(event, 'cvml', 'friedmud@gmail.com')">friedmud@gmail.com</a>> writes:<br>
<br>
> I am changing the stencil of my preconditioning matrix in the middle of a<br>
> solve a solve... and I'm running into this:<br>
<br>
Can you PCReset() after changing the stencil?<br>
<br>
The implementation of MatGetSubMatrices_MPIAIJ caches some intermediate<br>
structures to speed up subsequent calls.  I'll grant that Mat should do<br>
a better job of remembering what has changed so that it can<br>
automatically invalidate caches in this case.  There are only two things<br>
hard in computer science...<br>
</blockquote></div><br><br>-- <br>Sent from my iPhone<br>