>> 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.
Is there a mechanism to access the data stored in _dof_constraints, i.e., to loop over the list of constrained rows? > You also might need to sort and remove duplicates from ghost_dofs > before using it in init(). Yeah, I am already doing this out of paranoia. -- Boyce ------------------------------------------------------------------------------ _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
