Discretization of positive harmonic functions on Riemannian