Functional Gradient Flows for Constrained Sampling

Shiyue Zhang,Longlin Yu,Ziheng Cheng,Cheng Zhang
2024-10-31
Abstract:Recently, through a unified gradient flow perspective of Markov chain Monte Carlo (MCMC) and variational inference (VI), particle-based variational inference methods (ParVIs) have been proposed that tend to combine the best of both worlds. While typical ParVIs such as Stein Variational Gradient Descent (SVGD) approximate the gradient flow within a reproducing kernel Hilbert space (RKHS), many attempts have been made recently to replace RKHS with more expressive function spaces, such as neural networks. While successful, these methods are mainly designed for sampling from unconstrained domains. In this paper, we offer a general solution to constrained sampling by introducing a boundary condition for the gradient flow which would confine the particles within the specific domain. This allows us to propose a new functional gradient ParVI method for constrained sampling, called constrained functional gradient flow (CFG), with provable continuous-time convergence in total variation (TV). We also present novel numerical strategies to handle the boundary integral term arising from the domain constraints. Our theory and experiments demonstrate the effectiveness of the proposed framework.
Machine Learning
What problem does this paper attempt to address?
This paper aims to solve the problem of efficiently sampling probability distributions from constrained domains. Specifically, existing Markov Chain Monte Carlo (MCMC) methods and Variational Inference (VI) methods perform well when dealing with sampling tasks on unconstrained domains, but face challenges when confronted with constrained domains. These challenges include the need to solve nonlinear equations or rely on complex implicit and symplectic schemes designed for specific constrained domains. The paper proposes a new functional - gradient - based method - Constrained Functional Gradient (CFG). By introducing boundary conditions to restrict particles to move within a specific domain, this problem is solved. This method can not only effectively handle domains with inequality constraints, but also provides Total Variation (TV) convergence guarantees in continuous time. In addition, the paper also proposes novel numerical strategies to handle the boundary integral terms caused by domain constraints, further enhancing the effectiveness and efficiency of the method.