I also used the ghosted vector (VecCreateGhost). Does it have any relation to the DMGetGlobalVector() ?
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/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 >> >> >> > >
