[MOAB-dev] DMMoab error w/ SNES

[Vijay S. Mahadevan]

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.