This paper proposes and studies a numerical method for approximation of posterior expectations based on interpolation with Stein reproducing kernel. Finite-sample-size bounds the error are established distributions supported compact Riemannian manifold, we relate these to kernel discrepancy (KSD). Moreover, prove in our setting that KSD is equivalent Sobolev and, doing so, completely characteri...