Related papers: Linear systems solvers - recent developments and i…
A preconditioning for the Wilson fermion matrix on the lattice is defined which is particularly suited to the case when the temporal lattice spacing is much smaller than the spatial one. Details on the implementation of the scheme are…
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.…
In this work we investigate theoretical and computational aspects of novel lattice fermion formulations for the simulation of lattice gauge theories. The lattice approach to quantum gauge theories is an important tool for studying quantum…
Multistep matrix splitting iterations serve as preconditioning for Krylov subspace methods for solving singular linear systems. The preconditioner is applied to the generalized minimal residual (GMRES) method and the flexible GMRES (FGMRES)…
We present a method for formulating gauge theories of chiral fermions in lattice field theory. The method makes use of a Wilson mass to remove doublers. Gauge invariance is then restored by modifying the theory in two ways: the magnitude of…
Different recently developed Krylov space methods for solving linear systems are studied and compared for the solution of the Dirac equation on the lattice. Stabilized Biconjugate Gradient (BiCGstab2) is shown to be a robust and efficient…
Lattice formulations of QCD with Wilson fermions and a chirally twisted quark mass matrix provide an attractive framework for non-perturbative numerical studies. Owing to reparameterization invariance, the limiting continuum theory is just…
Quantum simulators have the exciting prospect of giving access to real-time dynamics of lattice gauge theories, in particular in regimes that are difficult to compute on classical computers. Future progress towards scalable quantum…
Preconditioning is at the core of modern many-fermion Monte Carlo algorithms, such as Hybrid Monte Carlo, where the repeated solution of a linear problem involving an ill-conditioned matrix is needed. We report on a performance comparison…
We present a modified version of the PRESB preconditioner for two-by-two block system of linear equations with the coefficient matrix $$\textbf{A}=\left(\begin{array}{cc} F & -G^* G & F \end{array}\right),$$ where $F\in\mathbb{C}^{n\times…
The numerical and computational aspects of chiral fermions in lattice quantum chromodynamics are extremely demanding. In the overlap framework, the computation of the fermion propagator leads to a nested iteration where the matrix vector…
The construction of multigrid operators for disordered linear lattice operators, in particular the fermion matrix in lattice gauge theories, by means of algebraic multigrid and block LU decomposition is discussed. In this formalism, the…
Gauge fixing is an essential step in lattice QCD calculations, particularly for studying gauge-dependent observables. Traditional iterative algorithms are computationally expensive and often suffer from critical slowing down and scaling…
Within the framework of quenched lattice QCD and using O(a) improved Wilson fermions and non-perturbative renormalisation, a high statistics computation of low moments of the unpolarised nucleon structure functions is given. Particular…
We consider three improvements to seed methods for Hermitian linear systems with multiple right-hand sides: only the Krylov subspace for the first system is used for seeding subsequent right-hand sides, the first right-hand side is solved…
In this paper we present deflation and augmentation techniques that have been designed to accelerate the convergence of Krylov subspace methods for the solution of linear systems of equations. We review numerical approaches both for linear…
Preconditioned Krylov subspace (KSP) methods are widely used for solving large-scale sparse linear systems arising from numerical solutions of partial differential equations (PDEs). These linear systems are often nonsymmetric due to the…
We study equilibrium statistical mechanics of Fermion lattice systems which require a different treatment compared with spin lattice systems due to the non-commutativity of local algebras for disjoint regions. Our major result is the…
I review recent machine trends and algorithmic developments for dynamical lattice QCD simulations with the HMC algorithm for Wilson-type fermions. The topics include the trend toward multi-core processors and general purpose GPU (GPGPU)…
We test an algebraic algorithm based on the coordinate-space method, evaluating with high accuracy the critical mass for Wilson fermions in lattice QCD at two loops. We test the results by using different types of infrared regularization.