Very nice answer, that I am still digesting. However I think there's a minor mistake concerning a minus sign, in that MINUS the divergence is the adjoint of the gradient, and that this leads to the nice formula for the discrete divergence that you give... right?