Related papers: Linear systems solvers - recent developments and i…
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 present a parallelizable SSOR preconditioning scheme for Krylov subspace iterative solvers which proves to be efficient in lattice QCD applications involving Wilson fermions. Our preconditioner is based on a locally lexicographic…
The stabilized biconjugate gradient algorithm BiCGStab recently presented by van der Vorst is applied to the inversion of the lattice fermion operator in the Wilson formulation of lattice Quantum Chromodynamics. Its computational efficiency…
We discuss different cases of dissipative Hamiltonian differential-algebraic equations and the linear algebraic systems that arise in their linearization or discretization. For each case we give examples from practical applications. An…
Linear systems arise in generating samples and in calculating observables in lattice quantum chromodynamics~(QCD). Solving the Hermitian positive definite systems, which are sparse but ill-conditioned, involves using iterative methods, such…
We investigate the application of Krylov space methods to the solution of shifted linear systems of the form (A+\sigma) x - b = 0 for several values of \sigma simultaneously, using only as many matrix-vector operations as the solution of a…
We present numerical methods to solve the Generalized Hartree-Fock theory for fermionic systems in lattices, both in thermal equilibrium and out of equilibrium. Specifically, we show how to determine the covariance matrix corresponding to…
We give a full account of the Numerical Stochastic Perturbation Theory method for Lattice Gauge Theories. Particular relevance is given to the inclusion of dynamical fermions, which turns out to be surprisingly cheap in this context. We…
This paper studies the solution of nonsymmetric linear systems by preconditioned Krylov methods based on the normal equations, LSQR in particular. On some examples, preconditioned LSQR is seen to produce errors many orders of magnitude…
Some results of test runs on a $6^3\times 12$ lattice with Wilson quarks and gauge group SU(2) for a previously proposed fermion algorithm by A. Slavnov are presented.
The solution of systems of linear(ized) equations lies at the heart of many problems in Scientific Computing. In particular for systems of large dimension, iterative methods are a primary approach. Stationary iterative methods are generally…
A parallelizable SSOR preconditioning scheme for Krylov subspace iterative solvers in lattice QCD applications involving Wilson fermions is presented. In actual Hybrid Monte Carlo and quark propagator calculations it helps to reduce the…
We study and simulate N=2 supersymmetric Wess-Zumino models in one and two dimensions. For any choice of the lattice derivative, the theories can be made manifestly supersymmetric by adding appropriate improvement terms corresponding to…
In this work I apply a recently proposed improvement procedure, originally conceived to reduce finite lattice spacing effects in transfer matrices for dilute Fermi systems, to tuning operators for the calculation of observables. I…
We present analytical results to guide numerical simulations with Wilson fermions in higher representations of the colour group. The ratio of $\Lambda$ parameters, the additive renormalization of the fermion mass, and the renormalization of…
In this paper we use some basic facts from the theory of (matrix) Lie groups and algebras to show that many of the classical matrix splittings used to construct stationary iterative methods and preconditioniers for Krylov subspace methods…
The availability of efficient Krylov subspace solvers play a vital role for the solution of a variety of numerical problems in computational science. Here we consider lattice field theory. We present a new general numerical method to…
We extend the epsilon-expansion of continuum chiral perturbation theory to nonzero lattice spacing in the framework of Wilson Chiral Perturbation Theory. We distinguish various regimes by defining the relative power counting of the quark…
New convergence bounds are presented for weighted, preconditioned, and deflated GMRES for the solution of large, sparse, non-Hermitian linear systems. These bounds are given for the case when the Hermitian part of the coefficient matrix is…
We present results for the renormalization of gauge invariant nonlocal fermion operators which contain a Wilson line, to one loop level in lattice perturbation theory. Our calculations have been performed for Wilson/clover fermions and a…