[petsc-users] Problem with SNES

[Milan Mitrovic]

sorry, it was the local to global mapping... forgot to turn it on
again... it seems that there is a problem with preallocation as I get
one of these:
[0]PETSC ERROR: --------------------- Error Message
------------------------------------
[0]PETSC ERROR: Argument out of range!
[0]PETSC ERROR: New nonzero at (93,351) caused a malloc!
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: Petsc Release Version 3.1.0, Patch 8, Thu Mar 17
13:37:48 CDT 2011
[0]PETSC ERROR: See docs/changes/index.html for recent updates.
[0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.
[0]PETSC ERROR: See docs/index.html for manual pages.
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: ./ppm_pf on a darwin10. named
milan-mitrovics-macbook-pro.local by milan Sat Aug 27 21:13:10 2011
[0]PETSC ERROR: Libraries linked from /Users/milan/Work/Library/petsc/lib
[0]PETSC ERROR: Configure run at Fri May  6 16:32:02 2011
[0]PETSC ERROR: Configure options --with-mpi=true
--with-mpi-dir=/Users/milan/Work/Library/openmpi.gcc/
--with-fortran=true --with-sundials=true
--with-sundials-dir=/Users/milan/Work/Library/sundials/
--prefix=/Users/milan/Work/Library/petsc
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: MatSetValues_MPIAIJ() line 358 in src/mat/impls/aij/mpi/mpiaij.c

even though in:

       call MatSetValuesLocal(Jac,1,(/i/),nzr(i),ind(i,:),val(i,:), &
            ADD_VALUES,info)

maxval(nzr) < nmax (the one used for preallocate)...

after that I get:

[1]PETSC ERROR: --------------------- Error Message
------------------------------------
[1]PETSC ERROR: Object is in wrong state!
[1]PETSC ERROR: Not for unassembled matrix!
[1]PETSC ERROR:
------------------------------------------------------------------------
[1]PETSC ERROR: Petsc Release Version 3.1.0, Patch 8, Thu Mar 17
13:37:48 CDT 2011
[1]PETSC ERROR: See docs/changes/index.html for recent updates.
[1]PETSC ERROR: See docs/faq.html for hints about trouble shooting.
[1]PETSC ERROR: See docs/index.html for manual pages.
[1]PETSC ERROR:
------------------------------------------------------------------------
[1]PETSC ERROR: ./ppm_pf on a darwin10. named
milan-mitrovics-macbook-pro.local by milan Sat Aug 27 21:13:10 2011
[1]PETSC ERROR: Libraries linked from /Users/milan/Work/Library/petsc/lib
[1]PETSC ERROR: Configure run at Fri May  6 16:32:02 2011
[1]PETSC ERROR: Configure options --with-mpi=true
--with-mpi-dir=/Users/milan/Work/Library/openmpi.gcc/
--with-fortran=true --with-sundials=true
--with-sundials-dir=/Users/milan/Work/Library/sundials/
--prefix=/Users/milan/Work/Library/petsc
[1]PETSC ERROR:
------------------------------------------------------------------------
[1]PETSC ERROR: MatMult() line 1883 in src/mat/interface/matrix.c
[1]PETSC ERROR: SNESLSCheckResidual_Private() line 60 in src/snes/impls/ls/ls.c
[1]PETSC ERROR: SNESSolve_LS() line 210 in src/snes/impls/ls/ls.c
[1]PETSC ERROR: SNESSolve() line 2255 in src/snes/interface/snes.c
[0] KSPDefaultConverged(): Linear solver has converged. Residual norm
0.0236436 is less than relative tolerance 800695 times initial right
hand side norm 0.0236436 at iteration 0
[0] SNESSolve_LS(): iter=0, linear solve iterations=0