bylore

Hi @bzjan ,have you fixed this issue?. I got similar results in a very simple linear elastic problem. The dirichlet boundary values are output as zero even when they are...

> @bylore: Yes, I fixed it by modifying the vtk output routine and using a custom build_parallel_solution_vector(). I can look for the code in the next few days. Many thanks...

> Are you calling equation_systems.init() after adding the bc object to the dof map ? Yes，equation_systems.init() is called after adding the bc object to the dof map.

> Are you calling equation_systems.init() after adding the bc object to the dof map ? It seems that the solutions on the constrained boundaries are not saved.

> Are you calling equation_systems.init() after adding the bc object to the dof map ? Issue #2452 may share some relations with mine

> Are you solving a linear system that you build with `constrain_element_matrix_and_vector` (or one of the other homogeneous-only APIs) rather than `heterogeneously_constrain_element_matrix_and_vector`? I solved a linear system，but using nonlinearimplicitsystem with...

> Are you solving a linear system that you build with `constrain_element_matrix_and_vector` (or one of the other homogeneous-only APIs) rather than `heterogeneously_constrain_element_matrix_and_vector`? Examples supplied by libmesh, such as /systems_of_equations/ex5 ex6...

> Except that would leave you without the BC contribution to the problem at all ... I assume that interior solution isn't purely coming from a convective term plus forcing...

> @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)....