Anderson Mixing in Bures Wasserstein Space of Gaussian Measures
Fuente:
arXiv
Saved in:
| Main Authors: | , , |
|---|---|
| Format: | Preprint |
| Published: |
2026
|
| Subjects: | |
| Online Access: | |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866914291843596288 |
|---|---|
| author | Aksenov, Vitalii Eigel, Martin Oster, Mathias |
| author_facet | Aksenov, Vitalii Eigel, Martin Oster, Mathias |
| contents | Various statistical tasks, including sampling or computing Wasserstein barycenters, can be reformulated as fixed-point problems for operators on probability distributions.
Accelerating standard fixed-point iteration schemes provides a promising novel approach to the design of efficient numerical methods for these problems.
The Wasserstein geometry on the space of probability measures, although not precisely Riemannian, allows us to define various useful Riemannian notions, such as tangent spaces, exponential maps and parallel transport, motivating the adaptation of Riemannian numerical methods.
We demonstrate this by developing and implementing the Riemannian Anderson Mixing (RAM) method for Gaussian distributions.
The method reuses the history of the residuals and improves the iteration complexity, and we argue that the additional costs, compared to Picard method, are negligible.
We show that certain open balls in the Bures-Wasserstein manifold satisfy the requirements for convergence of RAM.
The numerical experiments show a significant acceleration compared to a Picard iteration, and performance on par with Riemannian Gradient Descent and Conjugate Gradient methods. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2601_22038 |
| institution | arXiv |
| publishDate | 2026 |
| record_format | arxiv |
| spellingShingle | Anderson Mixing in Bures Wasserstein Space of Gaussian Measures Aksenov, Vitalii Eigel, Martin Oster, Mathias Optimization and Control Various statistical tasks, including sampling or computing Wasserstein barycenters, can be reformulated as fixed-point problems for operators on probability distributions. Accelerating standard fixed-point iteration schemes provides a promising novel approach to the design of efficient numerical methods for these problems. The Wasserstein geometry on the space of probability measures, although not precisely Riemannian, allows us to define various useful Riemannian notions, such as tangent spaces, exponential maps and parallel transport, motivating the adaptation of Riemannian numerical methods. We demonstrate this by developing and implementing the Riemannian Anderson Mixing (RAM) method for Gaussian distributions. The method reuses the history of the residuals and improves the iteration complexity, and we argue that the additional costs, compared to Picard method, are negligible. We show that certain open balls in the Bures-Wasserstein manifold satisfy the requirements for convergence of RAM. The numerical experiments show a significant acceleration compared to a Picard iteration, and performance on par with Riemannian Gradient Descent and Conjugate Gradient methods. |
| title | Anderson Mixing in Bures Wasserstein Space of Gaussian Measures |
| topic | Optimization and Control |
| url | https://arxiv.org/abs/2601.22038 |