[petsc-users] Error with SNESVISetVariableBounds

[Que Cat]

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 <quecat001 at gmail.com> 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 <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/f802c1f7/attachment.html>