Related papers: An iterative method to compute the overlap Dirac o…
In the presence of a non-vanishing chemical potential the eigenvalues of the Dirac operator become complex. We use a random matrix model approach to calculate analytically all correlation functions at weak and strong non-Hermiticity for…
Optimization of quadratic functions and the quotient of those are relevant in subspace and iterative optimization methods. In this paper, the calculation of the generalized operator norm and extremal generalized Rayleigh quotient is…
Solving symmetric positive definite linear problems is a fundamental computational task in machine learning. The exact solution, famously, is cubicly expensive in the size of the matrix. To alleviate this problem, several linear-time…
In this paper, we present a method for the computation of the low-lying real eigenvalues of the Wilson-Dirac operator based on the Arnoldi algorithm. These eigenvalues contain information about several observables. We used them to calculate…
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 recent years, a great deal of attention has been focused on numerically solving exponential integrators. The important ingredient to the implementation of exponential integrators is the efficient and accurate evaluation of the so called…
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 =…
We develop an accelerated gradient descent algorithm on the Grassmann manifold to compute the subspace spanned by a number of leading eigenvectors of a symmetric positive semi-definite matrix. This has a constant cost per iteration and a…
A class of second-order algorithms is proposed for minimizing smooth nonconvex functions that alternates between regularized Newton and negative curvature steps in an iteration-dependent subspace. In most cases, the Hessian matrix is…
In this paper, we get estimates on the higher eigenvalues of the Dirac operator on locally reducible Riemannian manifolds, in terms of the eigenvalues of the Laplace-Beltrami operator and the scalar curvature. These estimates are sharp, in…
Given a set of matrices, modeled as samples of a matrix-valued function, we suggest a method to approximate the underline function using a product approximation operator. This operator extends known approximation methods by exploiting the…
In this preliminary study, we examine the chiral properties of the parametrized Fixed-Point Dirac operator D^FP, see how to improve its chirality via the Overlap construction, measure the renormalized quark condensate Sigma and the…
The efficient and accurate QR decomposition for matrices with hierarchical low-rank structures, such as HODLR and hierarchical matrices, has been challenging. Existing structure-exploiting algorithms are prone to numerical instability as…
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…
We apply a complex chiral random matrix model as an effective model to QCD with a small chemical potential at zero temperature. In our model the correlation functions of complex eigenvalues can be determined analytically in two different…
We calculate non-singlet quark operator matrix elements of deep-inelastic scattering in the chiral limit including operators with total derivatives. This extends previous calculations with zero-momentum transfer through the operator vertex…
We consider an elliptic self-adjoint first order differential operator acting on pairs (2-columns) of complex-valued half-densities over a connected compact 3-dimensional manifold without boundary. The principal symbol of our operator is…
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 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.
Recent studies of the topological properties of a general class of lattice Dirac operators are reported. This is based on a specific algebraic realization of the Ginsparg-Wilson relation in the form…