Re: [Libmesh-users] Dirichlet constraints in NonlinearImplicitSystem

[Renato P. <re...@gm...>]

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 <vik...@gm...>
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 <re...@gm...> 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 <re...@gm...> 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
>> Lib...@li...
>> https://lists.sourceforge.net/lists/listinfo/libmesh-users
>>
>