Polynomial Filtered HMC -- an algorithm for lattice QCD with dynamical quarks
arXiv:1106.5625 · doi:10.1016/j.cpc.2012.05.002
Abstract
Polynomial approximations to the inverse of the fermion matrix are used to filter the dynamics of the upper energy scales in HMC simulations. The use of a multiple time-scale integration scheme allows the filtered pseudofermions to be evolved using a coarse step size. We introduce a novel generalisation of the nested leapfrog which allows for far greater flexibility in the choice of time scales. We observe a reduction in the computational expense of the molecular dynamics integration of between 3--5 which improves as the quark mass decreases.
10 pages, 7 figures
References in corpus (5)
- Trivializing maps, the Wilson flow and the HMC algorithm
- Accelerating Dynamical Fermion Computations using the Rational Hybrid Monte Carlo (RHMC) Algorithm with Multiple Pseudofermion Fields
- Physical Point Simulation in 2+1 Flavor Lattice QCD
- Non-renormalizability of the HMC algorithm
- A Novel Multiple-Time Scale Integrator for the Hybrid Monte Carlo Algorithm
Cited by in corpus (7)
- A new class of variance reduction techniques using lattice symmetries
- Scale setting for large- SUSY Yang-Mills on the lattice
- Applying polynomial filtering to mass preconditioned Hybrid Monte Carlo
- Single flavour optimisations to Hybrid Monte Carlo
- Single flavour filtering for RHMC in BQCD
- Exploiting the hopping parameter expansion in the hybrid Monte Carlo (HMC) simulation of lattice QCD with two degenerate flavours of Wilson fermions
- Accuracy of Symmetric Partitioned Runge-Kutta Methods for Differential Equations on Lie-Groups