Related papers: Error Bounds for the Sign Function
In this paper we propose and analyze an algorithm for identifying spectral gaps of a real symmetric matrix $A$ by simultaneously approximating the traces of spectral projectors associated with multiple different spectral slices. Our method…
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…
In this paper we discuss how the peculiar properties of twisted lattice QCD at maximal twist can be employed to set up a consistent computational scheme in which, despite the explicit breaking of chiral symmetry induced by the presence of…
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…
Approaches to finite baryon density lattice QCD usually suffer from uncontrolled systematic uncertainties in addition to the well-known sign problem. We test a method - sign reweighting - that works directly at finite chemical potential and…
The Lanczos method with implicit restarting is one of the most popular methods for finding a few exterior eigenpairs of a large symmetric matrix $A$. Usually based on polynomial filtering, restarting is crucial to limit memory and the cost…
The low-lying eigenvalue spectrum of the QCD Dirac operator in the epsilon-regime is expected to match with that of chiral Random Matrix Theory (ChRMT). We study this correspondence for the case including sea quarks by performing two-flavor…
I construct a Lanczos process on a large and sparse matrix and use the results of this iteration to compute the inverse square root of the same matrix. The algorithm is a stable version of an earlier proposal by the author. It can be used…
We outline our implementation of the adaptive thick-restart Lanczos algorithm ($a$-TRLan) for the projection of the low-lying eigenmodes of the overlap Dirac operator in lattice QCD, and compare the performances of our code and the widely…
In analogy to Neuberger's double-pass algorithm for the Conjugate Gradient inversion with multi-shifts we introduce a double-pass variant for BiCGstab(ell). One possible application is the overlap operator of QCD at non-zero chemical…
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…
Lattice fermions obeying the Ginsparg-Wilson relation do correctly represent the physical properties related to chirality. This can be achieved by local fermions, which involve an infinite number of couplings, however. For practical…
This paper considers the linear inverse problem where we wish to estimate a structured signal $x$ from its corrupted observations. When the problem is ill-posed, it is natural to make use of a convex function $f(\cdot)$ that exploits the…
We numerically evaluate the one-loop counterterms for the four-dimensional Wess-Zumino model formulated on the lattice using Ginsparg-Wilson fermions of the overlap (Neuberger) variety, together with an auxiliary fermion (plus…
Using non-perturbative techniques we have found the renormalisation factor, Z, in the RI-MOM scheme for quark bilinear operators in quenched QCD. We worked with overlap fermions using the Luescher-Weisz gauge action. Our calculation was…
The Golub-Welsch algorithm [ Math. Comp., 23: 221-230 (1969)] has long been assumed symmetric for estimating quadratic forms. Recent research indicates that asymmetric quadrature nodes may be more often and the existence of a practical…
A supersymmetric extension of the nonlinear O(3) sigma model in two spacetime dimensions is investigated by means of Monte Carlo simulations. We argue that it is impossible to construct a lattice action that implements both the O(3)…
At nonzero quark chemical potential dynamical lattice simulations of QCD are hindered by the sign problem caused by the complex fermion determinant. The severity of the sign problem can be assessed by the average phase of the fermion…
Lattice QCD should allow a derivation of the $\Delta I=1/2$ rule from first principles, but numerical calculations to date have been plagued by a variety of problems. After a brief review of these problems, we present several new methods…
We derive the effective action of the light fermion field of the domain-wall fermion, which is referred as $q(x)$ by Furman and Shamir. The inverse of the effective Dirac operator turns out to be identical to the inverse of the truncated…