Processing math: 100%
5
views
0
recommends
+1 Recommend
0 collections
    0
    shares
      • Record: found
      • Abstract: found
      • Article: found
      Is Open Access

      Faster Sampling from Log-Concave Densities over Polytopes via Efficient Linear Solvers

      Preprint
      ,

      Read this article at

      Bookmark
          There is no author summary for this article yet. Authors can add summaries to their articles on ScienceOpen to make them more accessible to a non-specialist audience.

          Abstract

          We consider the problem of sampling from a log-concave distribution π(θ)ef(θ) constrained to a polytope K:={θRd:Aθb}, where ARm×d and bRm.The fastest-known algorithm \cite{mangoubi2022faster} for the setting when f is O(1)-Lipschitz or O(1)-smooth runs in roughly O(md×mdω1) arithmetic operations, where the mdω1 term arises because each Markov chain step requires computing a matrix inversion and determinant (here ω2.37 is the matrix multiplication constant). We present a nearly-optimal implementation of this Markov chain with per-step complexity which is roughly the number of non-zero entries of A while the number of Markov chain steps remains the same. The key technical ingredients are 1) to show that the matrices that arise in this Dikin walk change slowly, 2) to deploy efficient linear solvers that can leverage this slow change to speed up matrix inversion by using information computed in previous steps, and 3) to speed up the computation of the determinantal term in the Metropolis filter step via a randomized Taylor series-based estimator.

          Related collections

          Author and article information

          Journal
          06 September 2024
          Article
          2409.04320
          f0cd649c-cc69-417d-8579-29c98d349103

          http://creativecommons.org/licenses/by/4.0/

          History
          Custom metadata
          The conference version of this paper appears in ICLR 2024
          cs.DS cs.LG stat.ML

          Data structures & Algorithms,Machine learning,Artificial intelligence
          Data structures & Algorithms, Machine learning, Artificial intelligence

          Comments

          Comment on this article