On Mon, Jul 23, 2012 at 6:56 PM, Mark F. Adams <mark.adams at columbia.edu>wrote:
> On Mon, Jul 23, 2012 at 2:58 PM, Mark F. Adams <mark.adams at > columbia.edu>wrote: > >> 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 >> >> pc->setupcalled=0; >> > > This is at the end of PCReset(). > > Where are you calling this from? > > > PCSetUp_GAMG > > Why use the private interface directly? > > > PCReset() resets the matrix apparently. So the base class resets too much. > Then call PCSetUp_MG and then pc->setupcalled = 0. -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120723/c484b441/attachment-0001.html>
