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/SNESVISetVariableBounds.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/externalpackages > >> > --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 > >> >
