Hi, Folks --
I am trying to setup a ghosted solution vector for a LinearImplicitSystem.
My incomplete understanding of how this works is that, on each
processor, I make a list of the off-processor global DOFs which are
needed, and then I do something like
system.solution->init(system.solution->size(),
system.solution->local_size(),
ghost_dofs, GHOSTED);
When I do this, an assertion failure results:
[0] /Users/griffith/sfw/libmesh/include/numerics/petsc_vector.h, line
787, compiled Apr 13 2010 at 16:41:48
which corresponds to
// If the mesh is disjoint, the following assertion will fail.
// If the mesh is not disjoint, every processor will either have
// all the dofs, none of the dofs, or some non-zero dofs at the
// boundary between processors.
libmesh_assert(n_local == 0 || n_local == n || !ghost.empty());
If I comment out this assertion, everything seems to work OK.
Am I doing something wrong here?
Thanks,
-- Boyce
------------------------------------------------------------------------------
Download Intel® Parallel Studio Eval
Try the new software tools for yourself. Speed compiling, find bugs
proactively, and fine-tune applications for parallel performance.
See why Intel Parallel Studio got high marks during beta.
http://p.sf.net/sfu/intel-sw-dev
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
