On Tue, 1 Jun 2010, Boyce Griffith wrote: > Would it work simply to append the send_list associated with the DofMap to my > list of ghosted dofs, e.g., > > const vector<unsigned int>& send_list = dof_map.get_send_list(); > vector<unsigned int> ghost_dofs = ...; > ghost_dofs.insert(ghost_dofs.end(), send_list.begin(), > send_list.end()); > AutoPtr<NumericVector<double> > ghosted_solution = > system.solution->clone(); > ghosted_solution->init(solution->size(), solution->local_size(), > ghost_dofs, GHOSTED);
That should be sufficient, although possibly inefficient if your non-standard vectors don't really need all the data on the traditional ghost dofs, which will get pulled in along with the constraint dependencies. You also might need to sort and remove duplicates from ghost_dofs before using it in init(). --- Roy ------------------------------------------------------------------------------ _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
