Nonlinear PDE + scalar constraint with PETSc?

Hi there! Check out this thread.