This error is also seen in src/snes/examples/tutorials/ex53.c. [Shri@Shrirangs-MacBook-Pro tutorials (master)]$ ./ex53 Vec Object: 1 MPI processes type: seq 1 1 1 1 1 -1 -1 0 -1 -1 [0]PETSC ERROR: --------------------- Error Message -------------------------------------------------------------- [0]PETSC ERROR: [0]PETSC ERROR: Must call DMShellSetGlobalVector() or DMShellSetCreateGlobalVector() [0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html for trouble shooting. [0]PETSC ERROR: Petsc Development GIT revision: v3.5-25-g0ace994 GIT Date: 2014-07-07 14:50:15 +0200 [0]PETSC ERROR: ./ex53 on a debug-master named Shrirangs-MacBook-Pro.local by Shri Tue Jul 29 14:00:33 2014 [0]PETSC ERROR: Configure options PETSC_ARCH=debug-master [0]PETSC ERROR: #1 DMCreateGlobalVector_Shell() line 245 in /Users/Shri/Documents/petsc/src/dm/impls/shell/dmshell.c [0]PETSC ERROR: #2 DMCreateGlobalVector() line 681 in /Users/Shri/Documents/petsc/src/dm/interface/dm.c [0]PETSC ERROR: #3 DMGetGlobalVector() line 154 in /Users/Shri/Documents/petsc/src/dm/interface/dmget.c [0]PETSC ERROR: --------------------- Error Message -------------------------------------------------------------- [0]PETSC ERROR: Object is in wrong state [0]PETSC ERROR: Matrix is missing diagonal entry 0 [0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html for trouble shooting. [0]PETSC ERROR: Petsc Development GIT revision: v3.5-25-g0ace994 GIT Date: 2014-07-07 14:50:15 +0200 [0]PETSC ERROR: ./ex53 on a debug-master named Shrirangs-MacBook-Pro.local by Shri Tue Jul 29 14:00:33 2014 [0]PETSC ERROR: Configure options PETSC_ARCH=debug-master [0]PETSC ERROR: #4 MatILUFactorSymbolic_SeqAIJ() line 1738 in /Users/Shri/Documents/petsc/src/mat/impls/aij/seq/aijfact.c [0]PETSC ERROR: #5 MatILUFactorSymbolic() line 6203 in /Users/Shri/Documents/petsc/src/mat/interface/matrix.c [0]PETSC ERROR: #6 PCSetUp_ILU() line 206 in /Users/Shri/Documents/petsc/src/ksp/pc/impls/factor/ilu/ilu.c [0]PETSC ERROR: #7 PCSetUp() line 902 in /Users/Shri/Documents/petsc/src/ksp/pc/interface/precon.c [0]PETSC ERROR: #8 KSPSetUp() line 305 in /Users/Shri/Documents/petsc/src/ksp/ksp/interface/itfunc.c [0]PETSC ERROR: #9 SNESSolve_VINEWTONRSLS() line 506 in /Users/Shri/Documents/petsc/src/snes/impls/vi/rs/virs.c [0]PETSC ERROR: #10 SNESSolve() line 3743 in /Users/Shri/Documents/petsc/src/snes/interface/snes.c [0]PETSC ERROR: #11 main() line 49 in /Users/Shri/Documents/petsc/src/snes/examples/tutorials/ex53.c [0]PETSC ERROR: ----------------End of Error Message -------send entire error message to [email protected] application called MPI_Abort(MPI_COMM_WORLD, 73) - process 0 [unset]: aborting job: application called MPI_Abort(MPI_COMM_WORLD, 73) - process 0
From: Matthew Knepley <[email protected]> Date: Tue, 29 Jul 2014 13:21:45 -0500 To: Que Cat <[email protected]> Cc: petsc-users <[email protected]> Subject: Re: [petsc-users] Error with SNESVISetVariableBounds >On Tue, Jul 29, 2014 at 1:19 PM, Que Cat <[email protected]> wrote: > >I also used the ghosted vector (VecCreateGhost). Does it have any >relation to the DMGetGlobalVector() ? > > > >We clearly have the problem that a full stack in not being printed. I am >open to the idea that it is a problem >in our SNESVI, but we need a way to track it down. Could you produce a >test example small enough to send >us? > > Thanks, > > Matt > > >Que > > > >On Tue, Jul 29, 2014 at 12:42 PM, Que Cat <[email protected]> wrote: > >Yes, I don't use DM at all and no information about who is calling >DMGetGlobalVector(). When I use snesnewtonls (without constraint on the >solution), the code worked fine. I only changed one of code : > > >ierr = SNESSetType(snes,SNESNEWTONTR );CHKERRQ(ierr); >--- > ierr = SNESSetType(snes,SNESVINEWTONSSLS );CHKERRQ(ierr); >and added : >ierr = SNESVISetVariableBounds(snes, low_b, up_b);CHKERRQ(ierr); > > >and the errors appeared. > >Que > > >On Tue, Jul 29, 2014 at 12:29 PM, Jed Brown <[email protected]> wrote: > >Que Cat <[email protected]> writes: > >> Hi Jed, >> >> I removed only one line from the error message that indicated the >> unnecessary information like directory path, program name, processor >>type. >> >> In this problem, I created the vector using VecCreateMPI, not with DM. >>Does >> it cause such problem? I have the unstructured grid and can not use DM > > >There should be more lines, like who is calling DMGetGlobalVector(). >Maybe that is your code and you don't check errors, but why are you >calling that function if you don't use DM (so the DMShell is created >automatically)? > >> >> Que >> >> >> On Tue, Jul 29, 2014 at 12:09 PM, Jed Brown <[email protected]> wrote: >> >>> Que Cat <[email protected]> writes: >>> >>> > Dear Petsc-Users, >>> > >>> > I called the SNESVISetVariableBounds >>> > < >>> >>>http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/SNES/SNESVIS >>>etVariableBounds.html#SNESVISetVariableBounds > > > >>> > >>> > and >>> > received the following error: >>> > >>> > [0]PETSC ERROR: Must call DMShellSetGlobalVector() or >>> > DMShellSetCreateGlobalVector() >>> > [0]PETSC ERROR: See >>>http://www.mcs.anl.gov/petsc/documentation/faq.html >>> for >>> > trouble shooting. >>> > [0]PETSC ERROR: Petsc Release Version 3.5.0, Jun, 30, 2014 >>> > .... >>> >>> It looks like you are truncating the error message. Please don't ever >>> do that. >>> >>> > [0]PETSC ERROR: Configure options --with-x=0 -with-pic >>> > >>> >>>--with-external-packages-dir=/opt/apps/intel13/mvapich2_1_9/petsc/3.5/ex >>>ternalpackages >>> > --with-mpi-compilers=1 --with-mpi-dir=/opt/apps/intel13/mvapich2/1.9 >>> > --with-scalar-type=real --with-shared-libraries=1 --with-matlab >>> > --with-matlab-dir=/work/apps/matlab/2013a --with-spai=1 >>>--download-spai >>> > --with-hypre=1 --download-hypre --with-mumps=1 --download-mumps >>> > --with-scalapack=1 --download-scalapack --with-blacs=1 >>>--download-blacs >>> > --with-spooles=1 --download-spooles --with-superlu=1 >>>--download-superlu >>> > --with-superlu_dist=1 --download-superlu_dist --with-parmetis=1 >>> > --download-parmetis --with-metis=1 --download-metis --with-hdf5=1 >>> > --with-hdf5-dir=/opt/apps/intel13/mvapich2_1_9/phdf5/1.8.9 >>> > --with-debugging=no >>> > >>> >>>--with-blas-lapack-dir=/opt/apps/intel/13/composer_xe_2013.3.163/mkl/lib >>>/intel64 >>> > --with-mpiexec=mpirun_rsh --COPTFLAGS= --CXXOPTFLAGS= --FOPTFLAGS= >>> > [ >>> > 0]PETSC ERROR: #1 DMCreateGlobalVector_Shell() line 245 in >>> > /opt/apps/intel13/mvapich2_1_9/petsc/3.5/src/dm/impls/shell/dmshell.c >>> > [0]PETSC ERROR: #2 DMCreateGlobalVector() line 681 in >>> > /opt/apps/intel13/mvapich2_1_9/petsc/3.5/src/dm/interface/dm.c >>> > [0]PETSC ERROR: #3 DMGetGlobalVector() line 154 in >>> > /opt/apps/intel13/mvapich2_1_9/petsc/3.5/src/dm/interface/dmget.c >>> > >>> > Could anyone give me a hint to the source of the error? Thank you. >>> > >>> > Que >>> > > > > > > > > > > > > > > > > > > > > > >-- >What most experimenters take for granted before they begin their >experiments is infinitely more interesting than any results to which >their experiments lead. >-- Norbert Wiener
