On Tue, 13 Apr 2010, Boyce Griffith wrote: > system.solution->init(system.solution->size(), > system.solution->local_size(), > ghost_dofs, GHOSTED);
Why create your own ghost_dofs variable rather than just passing in system.get_dof_map().get_send_list()? How do you construct ghost_dofs? > [0] /Users/griffith/sfw/libmesh/include/numerics/petsc_vector.h, line > 787, compiled Apr 13 2010 at 16:41:48 > > 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? Possibly. If you're passing in an empty ghost_dofs in a non-serial computation, that's almost certainly a mistake - under what conditions, other than the degenerate cases in the assertion, does one processor really need *no* data from any other? It could also be we've got an overzealous assertion, if there's a possibility I missed. --- Roy ------------------------------------------------------------------------------ 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
