On Thu, Aug 12, 2021 at 3:59 PM Renato Poli <[email protected]> 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 [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
