Regarding applying inhomogenous boundary conditions on RHS

[Balu sharma]

Dear everyone,

I am doing a load-controlled non-linear simulation in which I am solving a linear system for the newton update. While assembling the system before solve, I am trying have separate functions to assemble matrix and RHS. I am using AffineConstraints class to take care of the Dirichlet boundary conditions. From the documentation(distribute_local_to_global()) and my experience, I am not able to apply inhomogenous boundary condition on the RHS without already having an assemble system matrix. If I use the function with only RHS vectors, homogenous boundary condition is being applied by default. Is there any solution to this?

With regards 

Subramania Sharma