On Thu, Aug 12, 2021 at 3:59 PM Renato Poli <rebp...@gmail.com> wrote:
> Thanks, > > As I understood, it just works without any major concern in residual and > jacobian calculations, right? > I must add the DirichletBoundary to the system dof_map and LibMesh handles > the matrix and vectors so Petsc will receive the constrained system. > Is that correct? > Yes, as long as you have a line similar to dof_map.constrain_element_matrix_and_vector (Ke, Re, dof_indices); in your assembly routine like in systems_of_equations_ex8, then I think it should "just work". -- John _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
