Related papers: Minimizing storage in implementations of the overl…
We investigate a number of algorithms that calculate the quark propagators for the overlap-Dirac fermion operator. The QCD simulations were performed at beta = 5.9 with a lattice volume of 16**3*32.
We investigate the algorithms for dynamical overlap fermions aiming at improving the performance for large-scale simulations. We look for the best combination of Hybrid Monte Carlo options and iterative quark solvers with respect to the…
We present a novel method to compute the overlap Dirac operator at zero and nonzero quark chemical potential. To approximate the sign function of large, sparse matrices, standard methods project the operator on a much smaller Krylov…
Rational approximations of the matrix sign function lead to multishift methods. For non-Hermitian matrices long recurrences can cause storage problems, which can be circumvented with restarts. Together with deflation we obtain efficient…
Results of porting parts of the Lattice Quantum Chromodynamics code to modern FPGA devices are presented. A single-node, double precision implementation of the Conjugate Gradient algorithm is used to invert numerically the Dirac-Wilson…
Current dynamical overlap fermion hybrid Monte Carlo simulations encounter large fermionic forces when there is mixing between near zero-eigenvectors of the kernel operator. This leads to low acceptance rates when there is a large density…
A modern Fortran implementation of three Dirac operators (Wilson, Brillouin, Susskind) in lattice QCD is presented, based on OpenMP shared-memory parallelization and SIMD pragmas. The main idea is to apply a Dirac operator to $N_v$ vectors…
We report on our progress in using the overlap-Dirac fermion operator in simulations of lattice QCD. We have investigated the Lanczos based method of Borici, as well as various rational approximations, to calculate the step function in the…
The overlap operator in lattice QCD requires the computation of the sign function of a matrix, which is non-Hermitian in the presence of a quark chemical potential. In previous work we introduced an Arnoldi-based Krylov subspace…
In QCD chiral symmetry is explicitly broken by quark masses, the effect of which can be described reliably by chiral perturbation theory. Effects of explicit chiral symmetry breaking by the lattice regularisation of the Dirac operator,…
We derive the vector-like four dimensional overlap Dirac operator starting from a five dimensional Dirac action in the presence of a delta-function space-time defect. The effective operator is obtained by first integrating out all the…
We propose a novel algorithm for calculating multi-baryon correlation functions on the lattice. By considering the permutation of quarks (Wick contractions) and color/spinor contractions simultaneously, we construct a unified index list for…
The Overlap-Dirac operator provides a lattice regularization of massless vector gauge theories with an exact chiral symmetry. Practical implementations of this operator and recent results in quenched QCD using this Overlap-Dirac operator…
We present relaxation and preconditioning techniques which accelerate the inversion of the overlap operator by a factor of four on small lattices, with larger gains as the lattice size increases. These improvements can be used in both…
I apply a recently developed algorithm for reweighting simulations of lattice QCD from one quark mass to another to simulations performed with overlap fermions in the epsilon regime. I test it by computing the condensate from distributions…
The lattice QCD simulation with the lattice chiral symmetry is very attractive, however, it is difficult to maintain the symmetry at a modest numerical computation cost. A candidate to reduce the computational cost during the configuration…
We investigate the distribution of low-lying eigenmodes of hermitian Wilson-Dirac operator, $H_W$, with the gauge action whose form is designed to avoid topology change, as well as with the standard plaquette action. On the quenched…
We report on the two-flavor QCD simulation in the epsilon-regime using the overlap fermion formulation. Sea quark mass is reduced to ~ 2 MeV on a 16^3x32 lattice with the lattice spacing ~ 0.11fm. Topological charge is fixed at Q=0. We…
A recently developed numerical method for the calculation of derivatives of functions of general complex matrices, which can also be combined with implicit matrix function approximations such as Krylov-Ritz type algorithms, is presented. An…
Generating 2-by-2 unitary matrices in floating-precision arithmetic is a delicate task. One way to reduce the accumulation error is to use less floating-point operations to compute each of the entries in the 2-by-2 unitary matrix. This…