added lagrange multiplier to impose pressure mean value