RHMC with Block Solvers and Multiple Pseudofermions
arXiv:1808.01829 · doi:10.1103/PhysRevE.98.043306
Abstract
The dominant cost of most lattice QCD simulations is the inversion of the Dirac operator required to calculate the force term in the RHMC update. One way to improve this situation is to use multiple pseudofermions, which reduces the size and variance of this force and hence allows a larger integration step size to be used. This means fewer force term calculations are required, but at the cost of having to invert the Dirac operator for each pseudofermion field. This bottleneck can be addressed: recently there has been renewed interest in the use of block Krylov solvers, which can solve multiple right hand side vectors with significantly fewer iterations than are required if each vector is solved using a separate Krylov solver. We combine these two ideas, achieving a significant speed-up of RHMC lattice QCD simulations.
arXiv comment12 pages, 12 figures, published version
References in corpus (21)
- Speeding up the Hybrid-Monte-Carlo algorithm for dynamical fermions
- Accelerating Dynamical Fermion Computations using the Rational Hybrid Monte Carlo (RHMC) Algorithm with Multiple Pseudofermion Fields
- Testing and tuning symplectic integrators for Hybrid Monte Carlo algorithm in lattice QCD
- Speeding up Lattice QCD simulations with clover-improved Wilson Fermions
- Many Masses on One Stroke: Economic Computation of Quark Propagators
- Local coherence and deflation of the low quark modes in lattice QCD
- A New Exact Method for Dynamical Fermion Computations with Non-Local Actions
- The RHMC Algorithm for 2 Flavours of Dynamical Staggered Fermions
- Krylov space solvers for shifted linear systems
- Solution of the Dirac equation in lattice QCD using a domain decomposition method
- Fast Fermion Monte Carlo
- Exploring the HMC trajectory-length dependence of autocorrelation times in lattice QCD
- Application of block Krylov subspace algorithms to the Wilson-Dirac equation with multiple right-hand sides in lattice QCD
- UV-filtered fermionic Monte Carlo
- Shadow Hamiltonians, Poisson Brackets, and Gauge Theories
- Multigrid for Staggered Lattice Fermions
- Modified Block BiCGSTAB for Lattice QCD
- Pushing Memory Bandwidth Limitations Through Efficient Implementations of Block-Krylov Space Solvers on GPUs
- Application of preconditioned block BiCGGR to the Wilson-Dirac equation with multiple right-hand sides in lattice QCD
- Deflation and Flexible SAP-Preconditioning of GMRES in Lattice QCD Simulation
- Tuning the Hybrid Monte Carlo algorithm using molecular dynamics forces' variances