# Parallel computations for Metropolis Markov chains with Picard maps

arXiv: https://arxiv.org/abs/2506.09762 · HTML: https://mlsanity.com/paper/2506.09762/

## Key facts

- Cited by 1 papers in the map (data through 14 Oct 2026).
- First posted on arXiv 11 Jun 2025.

## Built on by (most cited first, 1 of 1)

- [A Unifying Framework for Parallelizing Sequential Models with Linear Dynamical Systems](https://mlsanity.com/paper/2509.21716/) (2025)

## Abstract (arXiv)

We develop parallel algorithms for simulating zeroth-order (aka gradient-free) Metropolis Markov chains based on the Picard map. For Random Walk Metropolis Markov chains targeting log-concave distributions π on ℝd, our algorithm generates samples close to π in O(√d) parallel iterations with O(√d) processors, therefore speeding up the convergence of the corresponding sequential implementation by a factor √d. Furthermore, a modification of our algorithm generates samples from an approximate measure  πr in O(1) parallel iterations and O(d) processors. We empirically assess the performance of the proposed algorithms in high-dimensional regression problems, an epidemic model where the gradient is unavailable and a real-word application in precision medicine. Our algorithms are straightforward to implement and may constitute a useful tool for practitioners seeking to sample from a prescribed distribution π using only point-wise evaluations of logπ and parallel computing.