Related papers: Minimizing storage in implementations of the overl…
A practical implementation of the Overlap-Dirac operator ${{1+\gamma_5\epsilon(H)}\over 2}$ is presented. The implementation exploits the sparseness of $H$ and does not require full storage. A simple application to parity invariant three…
The extreme computational costs of calculating the sign of the Wilson matrix within the overlap operator have so far prevented four dimensional dynamical overlap simulations on realistic lattice sizes, because the computational power…
We study three practical implementations of the Overlap-Dirac operator $D_o= (1/2) [1 + \gamma_5\epsilon(H_w)]$ in four dimensions. Two implementations are based on different representations of $\epsilon(H_w)$ as a sum over poles. One of…
The overlap Dirac operator at nonzero quark chemical potential involves the computation of the sign function of a non-Hermitian matrix. In this talk we present an iterative method, first proposed by us in Ref. [1], which allows for an…
Numerical evaluation of the overlap Dirac operator is difficult since it contains the sign function $\epsilon(H_w)$ of the Hermitian Wilson-Dirac operator $H_w$ with a negative mass term. The problems are due to $H_w$ having very small…
We use a continued fraction expansion of the sign-function in order to obtain a five dimensional formulation of the overlap lattice Dirac operator. Within this formulation the inverse of the overlap operator can be calculated by a single…
Improvements of various methods to compute the sign function of the hermitian Wilson-Dirac matrix within the overlap operator are presented. An optimal partial fraction expansion (PFE) based on a theorem of Zolotarev is given. Benchmarks…
We define a sparse hermitian lattice Dirac matrix, $H$, coupling $2n+1$ Dirac fermions. When $2n$ fermions are integrated out the induced action for the last fermion is a rational approximation to the hermitian overlap Dirac operator. We…
We use a continued fraction expansion of the sign-function in order to obtain a five dimensional formulation of the overlap lattice Dirac operator. Within this formulation the inverse of the overlap operator can be calculated by a single…
The action of the overlap-Dirac operator on a vector is typically implemented in directly through a multi-shift conjugate gradient solver. The compute-time this takes to evaluate depends upon the condition number $\kappa$ of the matrix that…
The low-lying eigenvalues of a (sparse) hermitian matrix can be computed with controlled numerical errors by a conjugate gradient (CG) method. This CG algorithm is accelerated by alternating it with exact diagonalisations in the subspace…
The Overlap operator fulfills the Ginsparg-Wilson relation exactly and therefore represents an optimal discretization of the QCD Dirac operator with respect to chiral symmetry. When computing propagators or in HMC simulations, where one has…
We propose a practical formulation of the overlap Dirac operator in lattice QCD that employs the diagonal Kenney-Laub rational iterates - expressed via their partial fraction decomposition - to approximate the matrix sign function. We…
The overlap Dirac operator at nonzero quark chemical potential involves the computation of the sign function of a non-Hermitian matrix. In this talk we present iterative Krylov subspace approximations, with deflation of critical…
We compare the behavior of different lattice Dirac operators in gauge backgrounds which are lattice discretizations of a classical instanton. In particular we analyze the standard Wilson operator, a chirally improved Dirac operator and the…
The overlap Dirac operator obeys the Ginsparg-Wilson equation and offers a possibility to introduce chiral symmetry on the lattice. Evaluating the overlap operator is numerically very expensive and one has to rely on approximation methods.…
We compare the efficiency of four different algorithms to compute the overlap Dirac operator, both for the speed, i.e., time required to reach a desired numerical accuracy, and for the adaptability, i.e., the scaling of speed with the…
The overlap Dirac operator in lattice QCD requires the computation of the sign function of a matrix. While this matrix is usually Hermitian, it becomes non-Hermitian in the presence of a quark chemical potential. We show how the action of…
Lattice QCD calculations were one of the first applications to show the potential of GPUs in the area of high performance computing. Our interest is to find ways to effectively use GPUs for lattice calculations using the overlap operator.…
The action of the overlap-Dirac operator on a vector is typically implemented indirectly through a multi-shift conjugate gradient solver. The compute-time required depends upon the condition number, $\kappa$, of the matrix that is used as…