<div class="gmail_quote">On Mon, Jul 23, 2012 at 6:56 PM, Mark F. Adams <span dir="ltr"><<a href="mailto:mark.adams@columbia.edu" target="_blank">mark.adams@columbia.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div><div class="im"><blockquote type="cite"><div class="gmail_quote">On Mon, Jul 23, 2012 at 2:58 PM, Mark F. Adams <span dir="ltr"><<a href="mailto:mark.adams@columbia.edu" target="_blank">mark.adams@columbia.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
I'm trying to have PC_GAMG redo its setup (of MG) on every solve.  I call PCReset_MG, but looks like PCReset_MG does not reset<br>
<br>
pc->setupcalled=0;<br></blockquote><div><br></div><div>This is at the end of PCReset().</div><div><br></div><div>Where are you calling this from? </div></div></blockquote><div><br></div></div>PCSetUp_GAMG</div><div><div class="im">
<br><blockquote type="cite"><div class="gmail_quote"><div>Why use the private interface directly?</div></div></blockquote><div><br></div></div><div>PCReset() resets the matrix apparently.  So the base class resets too much.</div>
</div></blockquote><div><br></div><div> Then call PCSetUp_MG and then pc->setupcalled = 0.</div></div>