Related papers: Adaptive multigrid algorithm for the lattice Wilso…
The microscopic spectral density of the Hermitian Wilson-Dirac operator is computed numerically in quenched lattice QCD. We demonstrate that the results given for fixed index of the Wilson-Dirac operator can be matched by the predictions…
We present a family of spacetree-based multigrid realizations using the tree's multiscale nature to derive coarse grids. They align with matrix-free geometric multigrid solvers as they never assemble the system matrices which is cumbersome…
We present a multi-level algorithm for the solution of five dimensional chiral fermion formulations, including domain wall and Mobius Fermions. The algorithm operates on the red-black preconditioned Hermitian operator, and directly…
The simulation of lattice QCD on massively parallel computers stimulated the development of scalable algorithms for the solution of sparse linear systems. We tackle the problem of the Wilson-Dirac operator inversion by combining a Schwarz…
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 Dirac operator, which satisfies the Ginsparg-Wilson relation, realizes exact chiral symmetry on the lattice without any unphysical doubler modes. To perform the path integrals, one should, however, note that the overlap fermion…
We propose a method to improve lattice operators composed of Wilson fermions which allows the removal of all corrections of $O(a)$, including those proportional to the quark mass, leaving only errors of $O(a^2)$. The method exploits the…
Starting from the chiral Lagrangian for Wilson fermions at nonzero lattice spacing we have obtained compact expressions for all spectral correlation functions of the Hermitian Wilson Dirac operator in the $\epsilon$-domain of QCD with…
Critical slowing down in Krylov methods for the Dirac operator presents a major obstacle to further advances in lattice field theory as it approaches the continuum solution. Here we formulate a multi-grid algorithm for the Kogut-Susskind…
We introduce Random Matrix Models for the Hermitian Wilson-Dirac operator of QCD-like theories. We show that they are equivalent to the $\epsilon$-limit of the chiral Lagrangian for Wilson chiral perturbation theory. Results are obtained…
We consider two-flavor QCD in the lattice regularization with improved Wilson fermions. In this formulation chiral symmetry is explicitly broken at order a and hence the isovector axial currents require improvement as well as a finite…
It has been a big challenge for lattice QCD to simulate dynamical quarks near the chiral limit. Theoretically, it is well-known that the naive chiral symmetry cannot be realized on the lattice (the Nielsen-Ninomiya theorem). Also…
Lattice QCD solvers encounter critical slowing down for fine lattice spacings and small quark mass. Traditional matrix eigenvalue deflation is one approach to mitigating this problem. However, to improve scaling we study the effects of…
We present a progress report on a new class of multigrid solver algorithm suitable for the solution of 5d chiral fermions such as Domain Wall fermions and the Continued Fraction overlap. Unlike HDCG \cite{Boyle:2014rwa}, the algorithm works…
We consider the effect of discretization errors on the microscopic spectrum of the Wilson Dirac operator using both chiral Perturbation Theory and chiral Random Matrix Theory. A graded chiral Lagrangian is used to evaluate the microscopic…
Computing the trace of the inverse of large matrices is typically addressed through statistical methods. Deflating out the lowest eigenvectors or singular vectors of the matrix reduces the variance of the trace estimator. This work…
We report on an ongoing project to parametrize the Fixed-Point Dirac operator for massless quarks, using a very general construction which has arbitrarily many fermion offsets and gauge paths, the complete Clifford algebra and satisfies all…
We show that it is possible to improve the chiral behaviour and the approach to the continuum limit of correlation functions in lattice QCD with Wilson fermions by taking arithmetic averages of correlators computed in theories regularized…
We study the global symmetries of naive lattices Dirac operators in QCD-like theories in any dimension larger than two. In particular we investigate how the chosen number of lattice sites in each direction affects the global symmetries of…
With the ever-growing number of computing architectures, performance portability is an important aspect of (Lattice QCD) software. The Grid library provides a good framework for writing such code, as it thoroughly separates…