PoS - Proceedings of Science
Volume 256 - 34th annual International Symposium on Lattice Field Theory (LATTICE2016) - Algorithms and Machines
Domain Wall Fermion Simulations with the Exact One-Flavor Algorithm
D. Murphy
Full text: pdf
Pre-published on: December 27, 2016
Published on: March 24, 2017
As algorithmic developments have driven down the cost of simulating degenerate light quark pairs the relative cost of simulating single quark flavors with the Rational Hybrid Monte Carlo (RHMC) algorithm has become more expensive. TWQCD has proposed an exact one-flavor algorithm (EOFA) that allows for HMC simulations of a single quark flavor without taking a square root of the fermion determinant. We have independently implemented EOFA in the Columbia Physics System (CPS) and BAGEL Fermion Sparse-Matrix Library (BFM) for Shamir and M\"{o}bius domain wall fermions, and begun to optimize and test our implementation against RHMC. In this talk we discuss the derivation of the EOFA action, our tests of its equivalence to RHMC, and the current state of our implementation and optimization. We find, after introducing a novel preconditioning technique for the EOFA Dirac operator, that EOFA is a factor of 2.4 times faster than RHMC per molecular dynamics trajectory for the strange quark determinant on an $N_{f} = 2+1$ M\"{o}bius DWF ensemble with physical quark masses and a $24^{3} \times 64 \times 24$ volume. We expect that further improvement is possible by retuning the integrator parameters for EOFA and by continuing to optimize our code.
DOI: https://doi.org/10.22323/1.256.0272
How to cite

Metadata are provided both in "article" format (very similar to INSPIRE) as this helps creating very compact bibliographies which can be beneficial to authors and readers, and in "proceeding" format which is more detailed and complete.

Open Access
Creative Commons LicenseCopyright owned by the author(s) under the term of the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License.