Re: [Libmesh-users] Dirichlet constraints in NonlinearImplicitSystem

[John P. <jwp...@gm...>]

On Thu, Aug 12, 2021 at 3:59 PM Renato Poli <re...@gm...> 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