[petsc-users] Issues with subsequent GMRES Solves

[Jed Brown]

On Thu, 13 May 2010 10:20:49 -0400, "C.A.(Sandy) Mader" <mader at utias.utoronto.ca> wrote:
> I was able to get a version of the code working last night by destroying
> all of the PETSc objects between subsequent adjoint setup calls. Based
> on some of my previous results I suspect that I am having an issue with
> resetting the matrices used in the linear solve. I am running some more
> tests now to confirm this...

Try reusing all the data structures, but place a MatZeroEntries() before
your reassembly.  Perhaps the reassembly isn't setting all the same
entries, or is using ADD_VALUES?

Jed