Doesn't appear to be... -- Boyce
On 4/14/10 10:02 AM, Derek Gaston wrote: > 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
