[petsc-users] wrong no of nonzeros

[Derek Gaston]

I'll try tomorrow - thanks Jed!

Derek

On Saturday, November 29, 2014, Jed Brown <jed at jedbrown.org> wrote:

> Derek Gaston <friedmud at gmail.com <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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20141129/74f76743/attachment.html>