Wait... all you want is a ghosted solution vector? I thought that should be the default when using Petsc now... can anyone correct me?
Derek On Apr 13, 2010, at 7:05 PM, Boyce Griffith wrote: > 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 ------------------------------------------------------------------------------ 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
