Related papers: The eigSUMR inverter for overlap fermion
Low-rank matrix estimation is a canonical problem that finds numerous applications in signal processing, machine learning and imaging science. A popular approach in practice is to factorize the matrix into two compact low-rank factors, and…
We explore a scaled spectral preconditioner for the efficient solution of sequences of symmetric and positive-definite linear systems. We design the scaled preconditioner not only as an approximation of the inverse of the linear system but…
This thesis is dedicated to explore the feasibility of numerical calculations in the $\epsilon$--regime of QCD for the extraction of physical information. We apply two formulations of the Ginsparg-Wilson fermions the Neuberger operator and…
We consider a quadrature-based eigensolver to find eigenpairs of Hermitian matrices arising in lattice quantum chromodynamics. To reduce the computational cost for finding eigenpairs of such Hermitian matrices, we propose a new technique…
Because of the significant increase in size and complexity of the networks, the distributed computation of eigenvalues and eigenvectors of graph matrices has become very challenging and yet it remains as important as before. In this paper…
In recent years, a new kind of accelerated hardware has gained popularity in the Artificial Intelligence (AI) and Machine Learning (ML) communities which enables extremely high-performance tensor contractions in reduced precision for deep…
We report a lattice calculation of nucleon forward matrix elements on a $48^3 \times 96$ lattice at the physical pion mass and a spatial size of 5.5 fm. The $2+1$ flavor dynamical fermion configurations are generated with domain-wall…
An algorithm for the numerical inversion of large matrices, the biconjugate gradient algorithm (BGA), is investigated in view of its use for Monte Carlo simulations of fermionic field theories. It is compared with the usual conjugate…
We construct fast, structure-preserving iterations for computing the sign decomposition of a unitary matrix $A$ with no eigenvalues equal to $\pm i$. This decomposition factorizes $A$ as the product of an involutory matrix $S =…
This paper proposes Inverse Gram Matrix (IGM) methods to prioritize the Pairwise Reciprocal Matrix (PRM) in the Analytic Hierarchy Process. The IGM methods include Pseudo-IGM, Normalized-IGM, and Lagrange-IGM. Interestingly, the proposed…
A new O(N) algorithm based on a recursion method, in which the computational effort is proportional to the number of atoms N, is presented for calculating the inverse of an overlap matrix which is needed in electronic structure calculations…
We study SU(2) gluodynamics at finite temperature near the deconfining phase transition. We create the lattice ensembles using the tadpole improved Luscher-Weisz action. The overlap Dirac operator is used to determine the following three…
Versions of GMRES with deflation of eigenvalues are applied to lattice QCD problems. Approximate eigenvectors corresponding to the smallest eigenvalues are generated at the same time that linear equations are solved. The eigenvectors…
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 present a systematic comparison of various iterative methods to obtain the fermion propagator with both overlap and twisted mass fermions at fixed pion mass in the quenched approximation. Taking the best available algorithm in each case…
We define a sparse hermitian lattice Dirac matrix, $H$, coupling $2n+1$ Dirac fermions. When $2n$ fermions are integrated out the induced action for the last fermion is a rational approximation to the hermitian overlap Dirac operator. We…
We present an acceleration of the well-established Krylov-Ritz methods to compute the sign function of large complex matrices, as needed in lattice QCD simulations involving the overlap Dirac operator at both zero and nonzero baryon…
Measurements of the lowest-lying eigenvalues of the staggered fermion Dirac operator are made on ensembles of equilibrium gauge field configurations in quenched SU(3) lattice gauge theory. The results are compared with exact analytical…
The solution of large scale Sylvester matrix equation plays an important role in control and large scientific computations. A popular approach is to use the global GMRES algorithm. In this work, we first consider the global GMRES algorithm…
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…