Related papers: Multigrid solver for clover fermions
Lattice regularization of chiral fermions is an important development of the theory of elementary particles. Nontheless, brute force computer simulations are very expensive, if not prohibitive. In this letter I exploit the non-interacting…
The improved Wilson quark action - the clover action - is constructed to have smaller discretization errors than the normal Wilson quark action. We test this in a quenched spectroscopy computation on 6 lattice ensembles with spacings from…
Unfitted finite element methods have emerged as a popular alternative to classical finite element methods for the solution of partial differential equations and allow modeling arbitrary geometries without the need for a boundary-conforming…
Clover improvement is the standard choice for lattice simulations of QCD, when the lattice artefacts coming from Wilson fermions have to be reduced. However, the clover improvement is not limited to QCD, but can be applied to a wider range…
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.…
We present preliminary results of light hadron spectroscopy using valence, tadpole-improved, Clover fermions on an ensemble of gauge configurations generated with 2 flavors of staggered fermions at a beta of 5.6. We compare the slope and…
We show that using the multi-splitting algorithm as a preconditioner for the domain wall Dirac linear operator, arising in lattice QCD, effectively reduces the inter-node communication cost, at the expense of performing more on-node…
Simulations at physical quark masses are affected by the critical slowing down of the solvers. Multigrid preconditioning has proved to deal effectively with this problem. Multigrid accelerated simulations at the physical value of the pion…
A modern Fortran implementation of three Dirac operators (Wilson, Brillouin, Susskind) in lattice QCD is presented, based on OpenMP shared-memory parallelization and SIMD pragmas. The main idea is to apply a Dirac operator to $N_v$ vectors…
The Adaptive Aggregation-based Domain Decomposition Multigrid method (arXiv:1303.1377) is extended for two degenerate flavors of twisted mass fermions. By fine-tuning the parameters we achieve a speed-up of the order of hundred times…
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…
A fast multigrid solver is presented for high-order accurate Stokes problems discretised by local discontinuous Galerkin (LDG) methods. The multigrid algorithm consists of a simple V-cycle, using an element-wise block Gauss-Seidel smoother.…
We develop an improved lattice action for heavy quarks based on Brillouin-type fermions, that have excellent energy-momentum dispersion relation. The leading discretization errors of $O(a)$ and $O(a^2)$ are eliminated at tree-level. We…
We describe a way to optimize the chiral behavior of Wilson-type lattice fermion actions by studying the low energy real eigenmodes of the Dirac operator. We find a candidate action, the clover action with fat links with a tuned clover…
We present a comparison of a number of iterative solvers of linear systems of equations for obtaining the fermion propagator in lattice QCD. In particular, we consider chirally invariant overlap and chirally improved Wilson (maximally)…
Considering Ginsparg-Wilson type fermions dynamically in lattice QCD simulations is a challenging task. The hope is to be able to approach smaller pion masses and to eventually reach physical situations. The price to pay is substantially…
The Riesz maps of the $L^2$ de Rham complex frequently arise as subproblems in the construction of fast preconditioners for more complicated problems. In this work we present multigrid solvers for high-order finite element discretizations…
Multigrid solvers face multiple challenges on parallel computers. Two fundamental ones read as follows: Multiplicative solvers issue coarse grid solves which exhibit low concurrency and many multigrid implementations suffer from an…
We consider two-flavor QCD in the lattice regularization with improved Wilson fermions. In this formulation chiral symmetry is explicitly broken at order a and hence the isovector axial currents require improvement as well as a finite…
We compare different conjugate gradient -- like matrix inversion methods (CG, BiCGstab1 and BiCGstab2) employing for this purpose the compact lattice quantum electrodynamics (QED) with Wilson fermions. The main goals of this investigation…