Related papers: Multigrid solver for clover fermions
The goal of this primer is to provide a relatively short exposition of the basics of multigrid methods, simplified by focusing on fundamental concepts in a variational setting. This is done by way of a quadratic energy minimization…
We present a comparison of different multigrid approaches for the solution of systems arising from high-order continuous finite element discretizations of elliptic partial differential equations on complex geometries. We consider the…
We propose an adaptive multigrid preconditioning technology for solving linear systems arising from Discontinuous Petrov-Galerkin (DPG) discretizations. Unlike standard multigrid techniques, this preconditioner involves only trace spaces…
Overlap fermions are a powerful tool for investigating the chiral and topological structure of the vacuum and the thermal states of QCD. We study various chiral and topological aspects of the finite temperature phase transition of N_f=2…
A coarse grid correction (CGC) approach is proposed to enhance the efficiency of the matrix exponential and $\varphi$ matrix function evaluations. The approach is intended for iterative methods computing the matrix-vector products with…
We have implemented a variant of the force gradient integrator proposed by Kennedy et.al. and are using it in our production 2+1 flavor DWF simulations with pion masses of 180 MeV in (4.5fm)3 volumes. We find modest speed-ups (\sim 20%)…
Numerical climate- and weather-prediction requires the fast solution of the equations of fluid dynamics. Discontinuous Galerkin (DG) discretisations have several advantageous properties. They can be used for arbitrary domains and support a…
Solving discretized versions of the Dirac equation represents a large share of execution time in lattice Quantum Chromodynamics (QCD) simulations. Many high-performance computing (HPC) clusters use graphics processing units (GPUs) to offer…
We propose and analyze a Vanka-type multigrid solver for solving a sequence of complex-shifted Laplacian systems arising in diagonalization-based parallel-in-time algorithms for evolutionary equations. Under suitable assumption, local…
We present preliminary results from UKQCD simulations at light quark masses using two flavours of non-pertubatively improved Wilson fermions. We report on the performance of the standard HMC algorithm at these quark masses where m_pi/m_rho…
We expand the most general lattice Dirac operator D in a basis of simple operators. The Ginsparg-Wilson equation turns into a system of coupled quadratic equations for the expansion coefficients. Our expansion of D allows for a natural…
This paper analyzes the Krylov convergence rate of a Helmholtz problem preconditioned with Multigrid. The multigrid method is applied to the Helmholtz problem formulated on a complex contour and uses GMRES as a smoother substitute at each…
In this paper we discuss how the peculiar properties of twisted lattice QCD at maximal twist can be employed to set up a consistent computational scheme in which, despite the explicit breaking of chiral symmetry induced by the presence of…
Solving partial differential equations is crucial to analysing and predicting complex, large-scale physical systems but pushes conventional high-performance computers to their limits. Application specific photonic processors are an exciting…
We study exceptional modes of both the Wilson and the clover action in order to understand why quenched clover spectroscopy suffers so severely from exceptional configurations. We show that, in contrast to the case of the Wilson action, a…
We solve Poisson's equation using new multigrid algorithms that converge rapidly. The novel feature of the 2D and 3D algorithms are the use of extra diagonal grids in the multigrid hierarchy for a much richer and effective communication…
Classical iterative methods for tomographic reconstruction include the class of Algebraic Reconstruction Techniques (ART). Convergence of these stationary linear iterative methods is however notably slow. In this paper we propose the use of…
We show that it is possible to improve the chiral behaviour and the approach to the continuum limit of correlation functions in lattice QCD with Wilson fermions by taking arithmetic averages of correlators computed in theories regularized…
We compare the efficiency of four different algorithms to compute the overlap Dirac operator, both for the speed, i.e., time required to reach a desired numerical accuracy, and for the adaptability, i.e., the scaling of speed with the…
We present a simple and effective multigrid-based Poisson solver of second-order accuracy in both gravitational potential and forces in terms of the one, two and infinity norms. The method is especially suitable for numerical simulations…