Does anyone has experience with inequality contraints in dolfin? I know this has to do with Lagrange multipliers and rewriting the weak form somehow. Let's say you have an energy functional
E[u]*dx, such that E >=0 and u >=0 everywhere on the mesh
that you are trying to minimize. How do you add these constraints to the weak form?
