Related papers: Adaptive multigrid algorithm for the lattice Wilso…
New exact upper and lower bounds are derived on the spectrum of the square of the hermitian Wilson Dirac operator. It is hoped that the derivations and the results will be of help in the search for ways to reduce the cost of simulations…
In this work, an efficient blackbox-type multigrid method is proposed for solving multipoint flux approximations of the Darcy problem on logically rectangular grids. The approach is based on a cell-centered multigrid algorithm, which…
A new class of lattice Dirac operators which satisfy the index theorem have been recently proposed on the basis of the algebraic relation $\gamma_{5}(\gamma_{5}D) + (\gamma_{5}D)\gamma_{5} = 2a^{2k+1}(\gamma_{5}D)^{2k+2}$. Here $k$ stands…
We review the numerical analysis' understanding of Krylov subspace methods for solving (non-hermitian) systems of equations and discuss its implications for lattice gauge theory computations using the example of the Wilson fermion matrix.…
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 dominant cutoff effects in lattice QCD with Wilson quarks are proportional to the lattice spacing a. In particular, the isovector axial current satisfies the PCAC relation only up to such effects. Following a suggestion of Symanzik,…
We discuss the implementation of a Sheikholeslami-Wohlert term for simulations of lattice QCD with dynamical Wilson fermions as required by Symanzik's improvement program. We show that for the Hybrid Monte Carlo or Kramers equation…
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…
It is shown that the 21264 Alpha processor can reach about 20% sustained efficiency for the inversion of the Wilson-Dirac operator. Since fast ethernet is not sufficient to get balancing between computation and communication on reasonable…
The fermion determinant is a highly non-local object and its logarithm is an extensive quantity. For these reasons it is widely believed that the determinant cannot be treated in acceptance steps of gauge link configurations that differ in…
Lattice simulations of Quantum Chromodynamics (QCD) enable one to calculate the low-energy properties of the strong interaction among quarks and gluons based on the first principle. The most time-consuming part of the numerical simulations…
The Adaptive Aggregation-based Domain Decomposition Multigrid method (arXiv:1303.1377) is extended for two degenerate flavors of twisted mass fermions. By fine-tuning the parameters we achieve a speed-up of the order of hundred times…
We apply the topology conserving gauge action proposed by Luescher to the four-dimensional lattice QCD simulation in the quenched approximation. With this gauge action the topological charge is stabilized along the hybrid Monte Carlo…
In lattice QCD it is possible, in principle, to determine the parameters in the effective chiral lagrangian (including weak interaction couplings) by performing numerical simulations in the $\epsilon$--regime, i.e. at quark masses where the…
We study the spectrum of the hermitian Wilson Dirac operator in the epsilon-regime of QCD in the quenched approximation and compare it to predictions from Wilson Random Matrix Theory. Using the distributions of single eigenvalues in the…
The non-normality of Wilson-type lattice Dirac operators has important consequences - the application of the usual concepts from the textbook (hermitian) quantum mechanics should be reconsidered. This includes an appropriate definition of…
It is well known that the block Krylov subspace solvers work efficiently for some cases of the solution of differential equations with multiple right-hand sides. In lattice QCD calculation of physical quantities on a given configuration…
Lattice formulations of QCD with Wilson fermions and a chirally twisted quark mass matrix provide an attractive framework for non-perturbative numerical studies. Owing to reparameterization invariance, the limiting continuum theory is just…
We compute fermionic observables relevant to the study of chiral symmetry in quenched QCD using the Overlap-Dirac operator for a wide range of the fermion mass. We use analytical results to disentangle the contribution from exact zero modes…
Lattice QCD using fermions whose Dirac operator obeys the Ginsparg-Wilson relation, is perhaps the best known formulation of QCD with a finite cutoff. It reproduces all the low energy QCD phenomenology associated with chiral symmetry at…