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

Barry Smith bsmith at mcs.anl.gov
Mon Aug 28 06:46:54 CDT 2017

> 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?


More information about the petsc-users mailing list