[petsc-users] Error with SNESVISetVariableBounds

Tue Jul 29 14:06:49 CDT 2014

This error is also seen in src/snes/examples/tutorials/ex53.c.

[Shri at 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 petsc-maint at mcs.anl.gov----------
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 <knepley at gmail.com>
Date:  Tue, 29 Jul 2014 13:21:45 -0500
To:  Que Cat <quecat001 at gmail.com>
Cc:  petsc-users <petsc-users at mcs.anl.gov>
Subject:  Re: [petsc-users] Error with SNESVISetVariableBounds

>On Tue, Jul 29, 2014 at 1:19 PM, Que Cat <quecat001 at gmail.com> 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 <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/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