On Thu, 5 Dec 2019, Prashant Kumar Jha wrote: > I am not setting boundary id of the element (can we set them?).
Yes, with BoundaryInfo::add_side(). > For Dirichlet boundary condition on 1 dof, I am not sure how to put a > constraint row. With DofMap::add_constraint_row() (using the overload with the constraint_rhs argument, unless you only ever have homogeneous Dirichlet boundaries). --- Roy _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users