<div dir="ltr">Hi,<br><br>Does anybody have an idea how to implement a Lagrange multiplier in GetDP instead of a constraint? In this case for a mechanical problem. I've tried adding the following:<br><br>To the FunctionSpace add the variables for the Lagrange multiplier [ lxn, lyn, lzn ] that support the nodes defined in region "Fixation":<br><br>  { Name H_l_Mec3D ; Type Vector ;<br>    BasisFunction {<br>      { Name txn ; NameOfCoef lxn ; Function BF_NodeX ;<br>        Support Fixation ; Entity NodesOf[ All ] ; }<br>      { Name tyn ; NameOfCoef lyn ; Function BF_NodeY ;<br>        Support Fixation ; Entity NodesOf[ All ] ; }<br>      { Name tzn ; NameOfCoef lzn ; Function BF_NodeZ ;<br>        Support Fixation ; Entity NodesOf[ All ] ; }<br>    }<br>    Constraint {<br>    }<br>  }  <br><br>Then in the Formulation I add the equations (identity matrix for now) on the off-diagonals, i.e. [ 0 I; I 0 ] * [ u; l ] = [ 0; 0 ]:<br><br>       // Constraints as Lagrange Multipliers<br>       Galerkin { [ TensorDiag[ 1, 1, 1 ] * Dof{D1 l} , {D1 u} ] ;<br>         In Fixation ;<br>         Jacobian Vol ; Integration GradGrad ; }<br>       Galerkin { [ TensorDiag[ 1, 1, 1 ] * Dof{D1 u} , {D1 l} ] ;<br>         In Fixation ;<br>         Jacobian Vol ; Integration GradGrad ; }<br><br>But somehow this does not work. I have checked the matrices and the rows are added, but somehow they are all zero, so the identity matrix does not show up.<div><br></div><div>Any help would be greatly appreciated.</div><div><br></div><div>Best regards,</div><div><br></div><div>Jasper</div></div>