Saved in:
| Main Authors: | , , , , |
|---|---|
| Format: | Preprint |
| Published: |
2025
|
| Subjects: | |
| Online Access: | https://arxiv.org/abs/2509.11424 |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866911153958944768 |
|---|---|
| author | Frommer, Andreas Jimenez-Merchan, Jose Knechtli, Francesco Korzec, Tomasz Ramirez-Hidalgo, Gustavo |
| author_facet | Frommer, Andreas Jimenez-Merchan, Jose Knechtli, Francesco Korzec, Tomasz Ramirez-Hidalgo, Gustavo |
| contents | We introduce a multigrid multilevel Monte Carlo method for stochastic trace estimation in lattice QCD based on orthogonal projections. This formulation extends the previously proposed oblique decomposition and it is assessed on three representative problems: the connected pseudoscalar correlator, the trace of the full Dirac operator's inverse $\mathrm{tr}(D^{-1})$, and disconnected fermion loops. For the connected correlator, variance reductions grow systematically with the time separation and lead to cost savings of up to a factor of 30 at large separations, outperforming both the plain Hutchinson's estimator and the oblique formulation. For $\mathrm{tr}(D^{-1})$, reductions are more modest but remain systematic, with stronger effects on more ill-conditioned systems. Disconnected loops show no improvement, since their variance is dominated by local same-slice contributions not targeted by the decomposition. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2509_11424 |
| institution | arXiv |
| publishDate | 2025 |
| record_format | arxiv |
| spellingShingle | Using orthogonal projectors in multigrid multilevel Monte Carlo for trace estimation in lattice QCD Frommer, Andreas Jimenez-Merchan, Jose Knechtli, Francesco Korzec, Tomasz Ramirez-Hidalgo, Gustavo High Energy Physics - Lattice We introduce a multigrid multilevel Monte Carlo method for stochastic trace estimation in lattice QCD based on orthogonal projections. This formulation extends the previously proposed oblique decomposition and it is assessed on three representative problems: the connected pseudoscalar correlator, the trace of the full Dirac operator's inverse $\mathrm{tr}(D^{-1})$, and disconnected fermion loops. For the connected correlator, variance reductions grow systematically with the time separation and lead to cost savings of up to a factor of 30 at large separations, outperforming both the plain Hutchinson's estimator and the oblique formulation. For $\mathrm{tr}(D^{-1})$, reductions are more modest but remain systematic, with stronger effects on more ill-conditioned systems. Disconnected loops show no improvement, since their variance is dominated by local same-slice contributions not targeted by the decomposition. |
| title | Using orthogonal projectors in multigrid multilevel Monte Carlo for trace estimation in lattice QCD |
| topic | High Energy Physics - Lattice |
| url | https://arxiv.org/abs/2509.11424 |