[petsc-users] wrong no of nonzeros

[Derek Gaston]

I am changing the stencil of my preconditioning matrix in the middle of a
solve a solve... and I'm running into this:

[0]PETSC ERROR: Nonconforming object sizes
[0]PETSC ERROR: Cannot reuse matrix. wrong no of nonzeros
[0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html for
trouble shooting.
[0]PETSC ERROR: Petsc Release Version 3.5.1, Jul, 24, 2014
[0]PETSC ERROR: ../bison-opt on a arch-darwin-c-opt named dereksmacpro.home
by gastdr Sat Nov 29 21:15:28 2014
[0]PETSC ERROR: Configure options
--prefix=/opt/moose/petsc/openmpi_petsc-3.5.1/clang-opt-superlu
--with-hypre-dir=/opt/moose/hypre/openmpi_hypre-2.8.0b/clang-opt
--with-debugging=no --with-pic=1 --with-shared-libraries=1 --with-cc=mpicc
--with-cxx=mpicxx --with-fc=mpif90 --download-fblaslapack=1
--download-metis=1 --download-parmetis=1 --download-superlu_dist=1 CC=mpicc
CXX=mpicxx FC=mpif90 F77=mpif77 F90=mpif90 CFLAGS="-fPIC -fopenmp"
CXXFLAGS="-fPIC -fopenmp" FFLAGS="-fPIC -fopenmp" FCFLAGS="-fPIC -fopenmp"
F90FLAGS="-fPIC -fopenmp" F77FLAGS="-fPIC -fopenmp"
PETSC_DIR=/opt/moose/stack_src/petsc-3.5.1
[0]PETSC ERROR: #1 MatGetSubMatrices_MPIAIJ_Local() line 1312 in
/opt/moose/stack_src/petsc-3.5.1/src/mat/impls/aij/mpi/mpiov.c
[0]PETSC ERROR: #2 MatGetSubMatrices_MPIAIJ() line 790 in
/opt/moose/stack_src/petsc-3.5.1/src/mat/impls/aij/mpi/mpiov.c
[0]PETSC ERROR: #3 MatGetSubMatrices() line 6352 in
/opt/moose/stack_src/petsc-3.5.1/src/mat/interface/matrix.c
[0]PETSC ERROR: #4 PCSetUp_ASM() line 369 in
/opt/moose/stack_src/petsc-3.5.1/src/ksp/pc/impls/asm/asm.c
[0]PETSC ERROR: #5 PCSetUp() line 902 in
/opt/moose/stack_src/petsc-3.5.1/src/ksp/pc/interface/precon.c
[0]PETSC ERROR: #6 KSPSetUp() line 305 in
/opt/moose/stack_src/petsc-3.5.1/src/ksp/ksp/interface/itfunc.c
[0]PETSC ERROR: #7 KSPSolve() line 417 in
/opt/moose/stack_src/petsc-3.5.1/src/ksp/ksp/interface/itfunc.c
[0]PETSC ERROR: #8 SNESSolve_NEWTONLS() line 232 in
/opt/moose/stack_src/petsc-3.5.1/src/snes/impls/ls/ls.c
[0]PETSC ERROR: #9 SNESSolve() line 3743 in
/opt/moose/stack_src/petsc-3.5.1/src/snes/interface/snes.c
[0]PETSC ERROR: #10 solve() line 559 in src/solvers/petsc_nonlinear_solver.C


I know that in the past we used to have to specify SAME_NONZERO_PATTERN or
DIFFERENT_NONZERO_PATTERN... but I thought those times were all past us ;-)

I'm using PETSc 3.5.1 on OSX ... underneath libMesh and MOOSE of course...

Any ideas where to start digging?

Thanks!
Derek
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20141129/358f1148/attachment.html>