We introduce numerical study on the discrete counterpart of Gauss' theorem. The purpose is to seek and establish a third approach, beside the analytical and the kernel-independent approaches, for efficient dimension reduction and preconditioning of equations initially in differential form. Integration is done locally, or globally, using analytical/symbolic rules as well as numerical rules and utilizing geometric information.