[petsc-users] Error with SNESVISetVariableBounds

[Que Cat]

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 <jed at jedbrown.org> wrote:

> Que Cat <quecat001 at gmail.com> 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 <jed at jedbrown.org> wrote:
> >
> >> Que Cat <quecat001 at gmail.com> 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
> >>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20140729/6e8b1c8f/attachment-0001.html>