I'll try tomorrow - thanks Jed! Derek
On Saturday, November 29, 2014, Jed Brown <[email protected]> wrote: > Derek Gaston <[email protected] <javascript:;>> writes: > > > I am changing the stencil of my preconditioning matrix in the middle of a > > solve a solve... and I'm running into this: > > Can you PCReset() after changing the stencil? > > The implementation of MatGetSubMatrices_MPIAIJ caches some intermediate > structures to speed up subsequent calls. I'll grant that Mat should do > a better job of remembering what has changed so that it can > automatically invalidate caches in this case. There are only two things > hard in computer science... > -- Sent from my iPhone
