Thank you! It helped me a lot. (By submesh, I mean the part of mesh on each processor (I use it in parallel). I need them to assemble constrain matrix for FETI type method.)
Dne 2016-09-01 16:50, John Peterson napsal: > On Thu, Sep 1, 2016 at 6:10 AM, Alena Vasatova <alena.vasat...@vsb.cz> > 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 Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
