Monte Carlo Anti-Differentiation for Approximate Weighted Model Integration

Probabilistic inference in the hybrid domain, i.e. inference over discrete-continuous domains, requires tackling two well known #P-hard problems 1)~weighted model counting (WMC) over discrete variables and 2)~integration over continuous variables. For both of these problems inference techniques have been developed separately in order to manage their #P-hardness, such as knowledge compilation for WMC and Monte Carlo (MC) methods for (approximate) integration in the continuous domain. Weighted model integration (WMI), the extension of WMC to the hybrid domain, has been proposed as a formalism to study probabilistic inference over discrete and continuous variables alike. Recently developed WMI solvers have focused on exploiting structure in WMI problems, for which they rely on symbolic integration to find the primitive of an integrand, i.e. to perform anti-differentiation. To combine these advances with state-of-the-art Monte Carlo integration techniques, we introduce \textit{Monte Carlo anti-differentiation} (MCAD), which computes MC approximations of anti-derivatives. In our empirical evaluation we substitute the exact symbolic integration backend in an existing WMI solver with an MCAD backend. Our experiments show that that equipping existing WMI solvers with MCAD yields a fast yet reliable approximate inference scheme.

[1]  Sean R Eddy,et al.  What is dynamic programming? , 2004, Nature Biotechnology.

[2]  Guy Van den Broeck,et al.  Hashing-Based Approximate Probabilistic Inference in Hybrid Domains , 2015, UAI.

[3]  A. Kennedy,et al.  Hybrid Monte Carlo , 1988 .

[4]  Marian Verhelst,et al.  PROBLP: A framework for Iow-precision probabilistic inference , 2019, 2019 56th ACM/IEEE Design Automation Conference (DAC).

[5]  Guy Van den Broeck,et al.  Component Caching in Hybrid Domains with Piecewise Polynomial Densities , 2016, AAAI.

[6]  Timon Gehr,et al.  PSI: Exact Symbolic Inference for Probabilistic Programs , 2016, CAV.

[7]  M. Simonovits,et al.  Random walks and an O * ( n 5 ) volume algorithm for convex bodies , 1997 .

[8]  Jesús A. De Loera,et al.  Software for exact integration of polynomials over polyhedra , 2012, ACCA.

[9]  Guy Van den Broeck,et al.  Probabilistic Inference in Hybrid Domains by Weighted Model Integration , 2015, IJCAI.

[10]  Michael A. Osborne,et al.  Frank-Wolfe Bayesian Quadrature: Probabilistic Integration with Theoretical Guarantees , 2015, NIPS.

[11]  Scott Sanner,et al.  Efficient Symbolic Integration for Probabilistic Inference , 2018, IJCAI.

[12]  David J. C. Mackay,et al.  Introduction to Monte Carlo Methods , 1998, Learning in Graphical Models.

[13]  Komei Fukuda,et al.  Exact volume computation for polytopes: a practical study , 1996 .

[14]  Jesús A. De Loera,et al.  Effective lattice point counting in rational convex polytopes , 2004, J. Symb. Comput..

[15]  Scott Sanner,et al.  Closed-Form Gibbs Sampling for Graphical Models with Algebraic Constraints , 2016, AAAI.

[16]  Adnan Darwiche,et al.  On probabilistic inference by weighted model counting , 2008, Artif. Intell..

[17]  Guy Van den Broeck,et al.  Efficient Search-Based Weighted Model Integration , 2019, UAI.

[18]  Jean-Louis Imbert About Redundant Inequalities Generated by Fourier's Algorithm , 1990, AIMSA.

[19]  Santosh S. Vempala,et al.  A practical volume algorithm , 2016, Math. Program. Comput..

[20]  Andrea Passerini,et al.  Advanced SMT techniques for weighted model integration , 2019, Artif. Intell..

[21]  Miklós Simonovits,et al.  How to compute the volume in high dimension? , 2003, Math. Program..

[22]  Andrea Passerini,et al.  Efficient Weighted Model Integration via SMT-Based Predicate Abstraction , 2017, IJCAI.

[23]  Scott Sanner,et al.  Symbolic Dynamic Programming for Discrete and Continuous State MDPs , 2011, UAI.

[24]  Parikshit Ram,et al.  Density estimation trees , 2011, KDD.

[25]  Adnan Darwiche,et al.  Modeling and Reasoning with Bayesian Networks , 2009 .

[26]  M. Rosenblatt Remarks on Some Nonparametric Estimates of a Density Function , 1956 .

[27]  Kristian Kersting,et al.  Mixed Sum-Product Networks: A Deep Architecture for Hybrid Domains , 2018, AAAI.

[28]  Ioannis Z. Emiris,et al.  Efficient Random-Walk Methods for Approximating Polytope Volume , 2014, Symposium on Computational Geometry.

[29]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[30]  Miklós Simonovits,et al.  Random walks and an O*(n5) volume algorithm for convex bodies , 1997, Random Struct. Algorithms.

[31]  Luc De Raedt,et al.  How to Exploit Structure while Solving Weighted Model Integration Problems , 2019, UAI.

[32]  Martin E. Dyer,et al.  A Random Polynomial Time Algorithm for Approximating the Volume of Convex Bodies , 1989, STOC.

[33]  H. Bungartz,et al.  Sparse grids , 2004, Acta Numerica.

[34]  Luc De Raedt,et al.  Exact and Approximate Weighted Model Integration with Probability Density Functions Using Knowledge Compilation , 2019, AAAI.

[35]  E. Parzen On Estimation of a Probability Density Function and Mode , 1962 .

[36]  M. Betancourt Nested Sampling with Constrained Hamiltonian Monte Carlo , 2010, 1005.0157.

[37]  Leslie G. Valiant,et al.  The Complexity of Computing the Permanent , 1979, Theor. Comput. Sci..

[38]  Gizem Karaali Book Review: Vector Calculus, Linear Algebra, and Differential Forms: A Unified Approach , 2007 .

[39]  S. Weinzierl Introduction to Monte Carlo methods , 2000, hep-ph/0006269.