[MOAB-dev] DMMoab error w/ SNES

[Gerd Heber]

Vijay, thanks a lot for those comments. I did not call SNESSetDM, which seems
to have caused the problem. It's working fine now. I was playing with a linear
problem earlier and didn't call KSPSetDM either, and there it didn't seem
to mind... Thanks again for your help.

Best, G.

-----Original Message-----
From: Vijay S. Mahadevan [mailto:vijay.m at gmail.com] 
Sent: Monday, February 9, 2015 6:21 PM
To: Gerd Heber
Cc: MOAB dev
Subject: Re: [MOAB-dev] DMMoab error w/ SNES

Gerd,

The only reason I can think of (without seeing your code) is that DMSetUp is not being called before the call to SNESSolve although I think if you have called SNESSetDM with an appropriate DMMoab object, that should get called internally anyway.

After looking closely at the line it's failing, the issue seems to be due to a reduce operation. Is this being run in parallel ? I should be able to change how this is done if this fails for a non-mpi builds.
Without any other useful code or test case, that's all I have for now.
Let me know if any of that is helpful.

Vijay

On Mon, Feb 9, 2015 at 3:07 PM, Gerd Heber <gheber at hdfgroup.org> wrote:
> I'm probably doing something stupid, but here's what I get when I try 
> to get an SNES solve going. I created the vectors involved with DMCreateGlobalVector.
> Maybe that's the problem.
>
> [0]PETSC ERROR: #1 DMCreateVector_Moab_Private() line 460 in 
> /tmp/petsc-3.5.3/src/dm/impls/moab/dmmbvec.cxx
> [0]PETSC ERROR: #2 DMVecDuplicate_Moab() line 658 in 
> /tmp/petsc-3.5.3/src/dm/impls/moab/dmmbvec.cxx
> [0]PETSC ERROR: #3 VecDuplicate() line 399 in 
> /tmp/petsc-3.5.3/src/vec/vec/interface/vector.c
> [0]PETSC ERROR: #4 VecDuplicateVecs_Default() line 827 in 
> /tmp/petsc-3.5.3/src/vec/vec/interface/vector.c
> [0]PETSC ERROR: #5 VecDuplicateVecs() line 473 in 
> /tmp/petsc-3.5.3/src/vec/vec/interface/vector.c
> [0]PETSC ERROR: #6 SNESSetWorkVecs() line 561 in 
> /tmp/petsc-3.5.3/src/snes/interface/snesut.c
> [0]PETSC ERROR: #7 SNESSetUp_NEWTONLS() line 317 in 
> /tmp/petsc-3.5.3/src/snes/impls/ls/ls.c
> [0]PETSC ERROR: #8 SNESSetUp() line 2616 in 
> /tmp/petsc-3.5.3/src/snes/interface/snes.c
> [0]PETSC ERROR: #9 SNESSolve() line 3731 in 
> /tmp/petsc-3.5.3/src/snes/interface/snes.c
> [0]PETSC ERROR: #10 main() line 144 in 
> /home/gerd/FEAWD/Drivers/src/ifeawd.cc
>
> Any thoughts/comments?
>
> Thanks, G.