Hi James,
Due to how the finite element method works, implementing such algebraic constraints only naturally comes in in the form of integral constraints (
https://featool.com/doc/physics#eqn_bdr_int_constr), that is the constraint is only applied in a integral sense (for any specific point it might not apply, as such pointwise constraints would be difficult to implement).