Dear all, I am trying to extend my current code to add Lagrange multipliers, leading to extra terms that couplet certain DOFs in a fixed range of the grid. To do so, I added another variable to my system, leading to the system-matrix of block-structure:
/ H L^T \ \ L 0 / where H is the matrix corresponding to the actual problem and L contains the Lagrange constrains. Thus, my system.solution contains in the first part the coefficients of the particular solution and in the lower part the Lagrange multipliers themselves. Now my question: is there some function to extract the Lagrange multipliers and the actual solution vector (dof-based) separately? Or how is it supposed to be done? Many thanks in advance, Tobias ------------------------------------------------------------------------------ Check out the vibrant tech community on one of the world's most engaging tech sites, Slashdot.org! http://sdm.link/slashdot _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
