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?
Thanks, Renato On Thu, Aug 12, 2021 at 11:07 AM Vikram Garg <vikram.v.g...@gmail.com> wrote: > Hello Renato, > I think example 8 in systems_of_equations should be > useful in this context. > > Vikram Garg > > > On Thu, Aug 12, 2021 at 8:48 AM Renato Poli <rebp...@gmail.com> wrote: > >> Hi all, >> >> I had no reply on this one, I am struggling with this design. >> I could not make this work, should I go for the Penalty method? >> >> Renato >> >> >> On Tue, Aug 10, 2021 at 10:49 PM Renato Poli <rebp...@gmail.com> wrote: >> >> > Hi, >> > >> > In linear systems, I was used to adding dirichlet constraints with >> > dogmap.add_dirichlet_constraints(). I needed to use >> > heterogeneously_constrain_matrix_and_vector after assembly. >> > >> > What would be the way to do that in NonlinearImplicitSystem? I >> understand >> > that NR solvers solve for the increment, so that "heterogeneously" would >> > make little sense. >> > >> > Should I expect Dirichlet constraints to work anyway? >> > >> > Renato >> > >> >> _______________________________________________ >> Libmesh-users mailing list >> Libmesh-users@lists.sourceforge.net >> https://lists.sourceforge.net/lists/listinfo/libmesh-users >> > _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
