vikramvgarg

Do the AMR element level rules need to apply across subdomains ? I dont think they should. If not, I think you can additionally mark all the elements in the...

Hello bzjan, looking at [miscellanous/example11](https://github.com/libMesh/libmesh/blob/master/examples/miscellaneous/miscellaneous_ex11/miscellaneous_ex11.C), where the subdivision element capability is demonstrated, it seems that the authors used a penalty method to enforce boundary conditions. It looks like that they...

@bzjan, this should be easy to test with example 11. If you replace 0 with some small value in the lines, const dof_id_type u_dof = nodes[n]->dof_number (system.number(), u_var, 0); const...

Oops, my bad. Yes those lines are returning dof_id_type, not setting the value of the boundary condition. I somehow assumed that this was doing the same thing we do while...

Are you calling equation_systems.init() after adding the bc object to the dof map ?

@bylore, can you try running adjoints/ex3 ? In that example, we apply a zero, const(1) as well as user defined parabolic Dirichlet boundary condition. See line 147 onwards [here](https://github.com/libMesh/libmesh/blob/master/examples/adjoints/adjoints_ex3/coupled_system.C). You...