Related papers: A nested Krylov subspace method for the overlap op…
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…
We discuss the salient features of Zolotarev optimal rational approximation for the inverse square root function, in particular, for its applications in lattice QCD with overlap Dirac quark. The theoretical error bound for the matrix-vector…
Many scientific applications require the evaluation of the action of the matrix function over a vector and the most common methods for this task are those based on the Krylov subspace. Since the orthogonalization cost and memory requirement…
One of the most computationally expensive steps of the low-rank ADI method for large-scale Lyapunov equations is the solution of a shifted linear system at each iteration. We propose the use of the extended Krylov subspace method for this…
Randomized iterative methods, such as the randomized Kaczmarz method, have gained significant attention for solving large-scale linear systems due to their simplicity and efficiency. Meanwhile, Krylov subspace methods have emerged as a…
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…
We propose a novel general approach to locality of lattice composite fields, which in case of QCD involves locality in both quark and gauge degrees of freedom. The method is applied to gauge operators based on the overlap Dirac matrix…
The overlap fermion offers the considerable advantage of exact chiral symmetry on the lattice, but is numerically intensive. This can be made affordable while still providing large lattice volumes, by using coarse lattice spacing, given…
For several classes of mathematical models that yield linear systems, the splitting of the matrix into its Hermitian and skew Hermitian parts is naturally related to properties of the underlying model. This is particularly so for…
We briefly review the overlap formalism for chiral gauge theories, the overlap Dirac operator for massless fermions and its connection to domain wall fermions. We describe properties of the overlap Dirac operator, and methods to implement…
In this paper we present an efficient active-set method for the solution of convex quadratic programming problems with general piecewise-linear terms in the objective, with applications to sparse approximations and risk-minimization. The…
An accurate residual--time (AccuRT) restarting for computing matrix exponential actions of nonsymmetric matrices by the shift-and-invert (SAI) Krylov subspace method is proposed. The proposed restarting method is an extension of the…
We determine the free field hypercubic Dirac operator which is optimally close to satisfying the Ginsparg-Wilson relation. Inserting this operator into the overlap formula, we show that the analytic locality bound on the resulting overlap…
Recently, a non-Hermitian chiral random matrix model was proposed to describe the eigenvalues of the QCD Dirac operator at nonzero chemical potential. This matrix model can be constructed from QCD by mapping it to an equivalent matrix model…
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.
In the present paper, we present some numerical methods for computing approximate solutions to some large differential linear matrix equations. In the first part of this work, we deal with differential generalized Sylvester matrix equations…
The overlap fermion offers the tremendous advantage of exact chiral symmetry on the lattice, but is numerically intensive. This can be made affordable while still providing large lattice volumes, by using coarse lattice spacing, given that…
Sparse matrices have recently played a significant and impactful role in scientific computing, including artificial intelligence-related fields. According to historical studies on sparse matrix--vector multiplication (SpMV), Krylov subspace…
On configurations with 2+1-flavor dynamical domain-wall fermions, we calculate the RI/(S)MOM renormalization constants (RC) of overlap quark bilinears. Hypercubic (HYP) smearing is used to construct the overlap Dirac operator. We…
This paper presents a new algorithm KIOPS for computing linear combinations of $\varphi$-functions that appear in exponential integrators. This algorithm is suitable for large-scale problems in computational physics where little or no…