Related papers: An iterative method to compute the overlap Dirac o…
We review a procedure of factorizing the Minkowski space Dirac operator over a~suitable superspace, discuss its Euclidean space version and apply the worked out formalism in the case od an almost-commutative Dirac operator. The presented…
In this study, we consider the numerical solution of large systems of linear equations obtained from the stochastic Galerkin formulation of stochastic partial differential equations. We propose an iterative algorithm that exploits the…
In the continuum, a topological obstruction to the vanishing of the non-abelian anomaly in 2n dimensions is given by the index of a certain Dirac operator in 2n+2 dimensions, or equivalently, the index of a 2-parameter family of Dirac…
We launched a project to perform dymanical fermion simulations using the overlap fermion formulation for sea quarks. In order to avoid the appearace of near-zero modes of the hermitian Wilson-Dirac operator $H_W$, we introduce a pair of…
We consider the low-rank alternating directions implicit (ADI) iteration for approximately solving large-scale algebraic Sylvester equations. Inside every iteration step of this iterative process a pair of linear systems of equations has to…
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 investigate the algorithms for dynamical overlap fermions aiming at improving the performance for large-scale simulations. We look for the best combination of Hybrid Monte Carlo options and iterative quark solvers with respect to the…
A common way to approximate $F(A)b$ -- the action of a matrix function on a vector -- is to use the Arnoldi approximation. Since a new vector needs to be generated and stored in every iteration, one is often forced to rely on restart…
Krylov subspace methods are a powerful tool for efficiently solving high-dimensional linear algebra problems. In this work, we study the approximation quality that a Krylov subspace provides for estimating the numerical range of a matrix.…
We prove that the Neumann, Dirichlet and regularity problems for divergence form elliptic equations in the half space are well posed in $L_2$ for small complex $L_\infty$ perturbations of a coefficient matrix which is either real symmetric,…
In this paper, a reduced-rank scheme with joint iterative optimization is presented for direction of arrival estimation. A rank-reduction matrix and an auxiliary reduced-rank parameter vector are jointly optimized to calculate the output…
Vlastimil Pt\'ak's method of nondiscrete induction is based on the idea that in the analysis of iterative processes one should aim at rates of convergence as functions rather than just numbers, because functions may give convergence…
A specific algebraic realization of the Ginsparg-Wilson relation in the form $\gamma_{5}(\gamma_{5}D)+(\gamma_{5}D)\gamma_{5} = 2a^{2k+1}(\gamma_{5}D)^{2k+2}$ is discussed, where $k$ stands for a non-negative integer and $k=0$ corresponds…
The rational Krylov subspace method (RKSM) and the low-rank alternating directions implicit (LR-ADI) iteration are established numerical tools for computing low-rank solution factors of large-scale Lyapunov equations. In order to generate…
Algorithms and implementations for computing the sign function of a triangular matrix are fundamental building blocks in algorithms for computing the sign of arbitrary square real or complex matrices. We present novel recursive and cache…
Applying domain decomposition to the lattice Dirac operator and the associated quark propagator, we arrive at expressions which, with the proper insertion of random sources therein, can provide improvement to the estimation of the…
We study the spectrum and dynamics of a one-dimensional discrete Dirac operator in a random potential obtained by damping an i.i.d. environment with an envelope of type $n^{-\alpha}$ for $\alpha>0$. We recover all the spectral regimes…
Starting from the chiral Lagrangian for Wilson fermions at nonzero lattice spacing we have obtained compact expressions for all spectral correlation functions of the Hermitian Wilson Dirac operator in the $\epsilon$-domain of QCD with…
The lower moments of the unpolarized and polarized deep-inelastic structure functions of the nucleon are calculated on the lattice. The calculation is done with Wilson fermions and for three values of the hopping parameter $\kappa$, so that…
An efficient and robust restart strategy is important for any Krylov-based method for eigenvalue problems. The tensor infinite Arnoldi method (TIAR) is a Krylov-based method for solving nonlinear eigenvalue problems (NEPs). This method can…