On Thu, Sep 1, 2016 at 6:10 AM, Alena Vasatova <[email protected]> wrote:
> Hi, > > I'm using libMesh by extending introduction example 4 (Solving a 2D or > 3D Poisson Problem in Parallel), is there any easy way how to get nodes > with Dirichlet boundary condition for each submesh? I thought I can use > bnd_nodes_begin/end or bid_nodes_begin/end function, but I found out I > can't because Dirichlet condition is stored as "Side Boundary > condition". > I'm not sure exactly what you mean by "submesh", but if you just want a list of nodes that are on the Dirichlet boundaries, you can use the BoundaryInfo object to get this. std::vector<dof_id_type> node_id_list; std::vector<boundary_id_type> bc_id_list; mesh.get_boundary_info().build_node_list_from_side_list(); mesh.get_boundary_info().build_node_list (node_id_list, bc_id_list); The best answer to your question really depends on what you are planning to do with the Nodes that are on the Dirichlet boundary though... -- John ------------------------------------------------------------------------------ _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
