I am doing a lagrange-multipliers problem similar to what is in the page
but I am not sure how to set up normalization constrant and the single value variable l. Can someone show me the code piece with decomposed mesh.
it gives error when I throw in the single column into the Mat block matrix.
Please see the example from https://github.com/FreeFem/FreeFem-sources/blob/master/examples/hpddm/laplace-lagrange-PETSc.edp.
thanks. Very helpful