[petsc-users] A number of questions about DMDA with SNES and Quasi-Newton methods

[Barry Smith]

> On Aug 28, 2017, at 12:13 AM, zakaryah . <zakaryah at gmail.com> wrote:
> 
> Any clues on what I should test with the debugger?  Just running with "-snes_type test -snes_test_display -start_in_debugger noxterm" the debugger reports a PETSc error:
> 
> MatSetValues_SeqAIJ(), Inserting a new nonzero at (7,0) in the matrix
> 
> (7,0) is certainly within the allocated space, considering the call to DMDACreate3d() that I sent before.  I'm suspicious that the -snes_type test is somehow incompatible with the DMDA, at least without somehow initializing the Jacobian.

   No definitely not. We use it all the time in many tests and never have problems with new nonzeros. Can you send your code to petsc-maint at mcs.anl.gov and we can run it?

  Barry