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

Reply via email to