Numerical Methods for the QCD Overlap Operator:III. Nested Iterations
arXiv:hep-lat/0405003 · doi:10.1016/j.cpc.2004.10.005
Abstract
The numerical and computational aspects of chiral fermions in lattice quantum chromodynamics are extremely demanding. In the overlap framework, the computation of the fermion propagator leads to a nested iteration where the matrix vector multiplications in each step of an outer iteration have to be accomplished by an inner iteration; the latter approximates the product of the sign function of the hermitian Wilson fermion matrix with a vector. In this paper we investigate aspects of this nested paradigm. We examine several Krylov subspace methods to be used as an outer iteration for both propagator computations and the Hybrid Monte-Carlo scheme. We establish criteria on the accuracy of the inner iteration which allow to preserve an a priori given precision for the overall computation. It will turn out that the accuracy of the sign function can be relaxed as the outer iteration proceeds. Furthermore, we consider preconditioning strategies, where the preconditioner is built upon an inaccurate approximation to the sign function. Relaxation combined with preconditioning allows for considerable savings in computational efforts up to a factor of 4 as our numerical experiments illustrate. We also discuss the possibility of projecting the squared overlap operator into one chiral sector.
33 Pages; citations added
References in corpus (36)
- Exactly massless quarks on the lattice
- More about exactly massless quarks on the lattice
- Locality properties of Neuberger's lattice Dirac operator
- Vector like gauge theories with almost massless fermions on the lattice
- Abelian chiral gauge theories on the lattice with exact gauge invariance
- Prospects for perfect actions
- A practical implementation of the Overlap-Dirac operator
- Numerical Methods for the QCD Overlap Operator: I. Sign-Function and Error Bounds
- Numerical techniques for lattice QCD in the --regime
- Dynamical overlap fermions, results with hybrid Monte-Carlo algorithm
- A Study of Practical Implementations of the Overlap-Dirac Operator in Four Dimensions
- Krylov space solvers for shifted linear systems
- Finite-size scaling of the quark condensate in quenched lattice QCD
- Light Quark Masses with Overlap Fermions in Quenched QCD
- Chiral Condensate in the Deconfined Phase of Quenched Gauge Theories
- Light quark masses, chiral condensate and quark-gluon condensate in quenched lattice QCD with exact chiral symmetry
- Optimal lattice domain-wall fermions
- A note on Zolotarev optimal rational approximation for the overlap Dirac operator
- Chiral Logs in Quenched QCD
- On the Neuberger overlap operator
- Numerical Methods for the QCD Overlap Operator IV: Hybrid Monte Carlo
- Light hadron spectroscopy in two-flavor QCD with small sea quark masses
- An alternative to domain wall fermions
- Locality of optimal lattice domain-wall fermions
- Lanczos Approach to the Inverse Square Root of a Large and Sparse Matrix
- Fractional Inversion in Krylov Space
- Alternative to Domain Wall Fermions
- A note on the practical feasibility of domain-wall fermions
- Chiral Fermions and Multigrid
- Efficient computation of low-lying eigenmodes of non-Hermitian Wilson-Dirac type matrices
- Speed and Adaptability of Overlap Fermion Algorithms
- Overlap and domain wall fermions: what is the price of chirality?
- Chiral symmetry breaking from Ginsparg-Wilson fermions
- Optimal lattice domain-wall fermions with finite Ns
- Fast methods for computing the Neuberger Operator
- Numerical Methods for the QCD Overlap Operator: II. Optimal Krylov Subspace Methods
Cited by in corpus (20)
- Deconfinement transition in two-flavour lattice QCD with dynamical overlap fermions in an external magnetic field
- Two-flavor lattice QCD simulation in the epsilon-regime with exact chiral symmetry
- Two-flavor QCD simulation with exact chiral symmetry
- Pion form factors from two-flavor lattice QCD with exact chiral symmetry
- Two-flavor lattice QCD in the epsilon-regime and chiral Random Matrix Theory
- Topology with Dynamical Overlap Fermions
- Numerical Methods for the QCD Overlap Operator IV: Hybrid Monte Carlo
- First-order restoration of SU(Nf) x SU(Nf) chiral symmetry with large Nf and Electroweak phase transition
- Iterative methods for overlap and twisted mass fermions
- Lattice QCD with fixed topology
- Dynamical overlap fermion at fixed topology
- Improvement of algorithms for dynamical overlap fermions
- Exploring the chiral regime with dynamical overlap fermions
- Electric polarizability of hadrons with overlap fermions on multi-GPUs
- A lattice Dirac operator for QCD with light dynamical quarks
- The eigSUMR inverter for overlap fermion
- Fast algorithms for simulating chiral fermions in U(1)lattice gauge theory
- Error Bounds for the Sign Function
- Dynamical fermions in lattice quantum chromodynamics
- Deconfinement transition in two-flavour lattice QCD with dynamical overlap fermions