Related papers: A nested Krylov subspace method for the overlap op…
A Krylov subspace recycling method for the efficient evaluation of a sequence of matrix functions acting on a set of vectors is developed. The method improves over the recycling methods presented in [Burke et al., arXiv:2209.14163, 2022] in…
We have recently given a construction of the overlap Dirac operator at nonzero quark chemical potential. Here, we introduce a quark chemical potential in the domain-wall fermion formalism and show that our earlier result is reproduced if…
An efficient Krylov subspace algorithm for computing actions of the $\varphi$ matrix function for large matrices is proposed. This matrix function is widely used in exponential time integration, Markov chains and network analysis and many…
Bivariate matrix functions provide a unified framework for various tasks in numerical linear algebra, including the solution of linear matrix equations and the application of the Fr\'echet derivative. In this work, we propose a novel…
This paper reviews the most popular methods which are used in lattice QCD to compute the determinant of the lattice Dirac operator: Gaussian integral representation and noisy methods. Both of them lead naturally to matrix function problems.…
The overlap operator is just the simplest of a class of Dirac operators with an exact chiral symmetry. I demonstrate how a general class of chiral Dirac operators can be constructed, show that they have no fermion doublers and that they are…
Thanks to its great potential in reducing both computational cost and memory requirements, combining sketching and Krylov subspace techniques has attracted a lot of attention in the recent literature on projection methods for linear…
In recent years two Krylov subspace methods have been proposed for solving skew symmetric linear systems, one based on the minimum residual condition, the other on the Galerkin condition. We give new, algorithm-independent proofs that in…
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…
We propose new techniques to implement numerically the overlap-Dirac operator which exploit the physical properties of the underlying theory to avoid nested algorithms. We test these procedures in the two-dimensional Schwinger model and the…
We propose new techniques for the numerical implementation of the overlap-Dirac operator, which exploit the physical properties of the underlying theory to avoid nested algorithms. We test these procedures in the two-dimensional Schwinger…
Most current prevalent iterative methods can be classified into the so-called extended Krylov subspace methods, a class of iterative methods which do not fall into this category are also proposed in this paper. Comparing with traditional…
We introduce an algorithm for estimating the trace of a matrix function $f(\mathbf{A})$ using implicit products with a symmetric matrix $\mathbf{A}$. Existing methods for implicit trace estimation of a matrix function tend to treat…
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…
I present several tricks to help implement the overlap Dirac operator numerically.
The computation of approximating e^tA B, where A is a large sparse matrix and B is a rectangular matrix, serves as a crucial element in numerous scientific and engineering calculations. A powerful way to consider this problem is to use…
The properties of the spectrum of the overlap Dirac operator and their relation to random matrix theory are studied. In particular, the predictions from chiral random matrix theory in topologically non-trivial gauge field sectors are…
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…
In QCD chiral symmetry is explicitly broken by quark masses, the effect of which can be described reliably by chiral perturbation theory. Effects of explicit chiral symmetry breaking by the lattice regularisation of the Dirac operator,…
This work is concerned with approximating matrix functions for banded matrices, hierarchically semiseparable matrices, and related structures. We develop a new divide-and-conquer method based on (rational) Krylov subspace methods for…