Related papers: Adaptive Multigrid Algorithm for Lattice QCD
The multigrid algorithm is an efficient numerical method for solving a variety of elliptic partial differential equations (PDEs). The method damps errors at progressively finer grid scales, resulting in faster convergence compared to…
We consider $4d$ compact lattice QED in the quenched approximation. First, we briefly summarize the spectrum of the staggered Dirac operator and its connection with random matrix theory. Afterwards we present results for the low-lying…
In this work, we propose a robust and easily implemented algebraic multigrid method as a stand-alone solver or a preconditioner in Krylov subspace methods for solving either symmetric and positive definite or saddle point linear systems of…
In lattice QCD it is possible, in principle, to determine the parameters in the effective chiral lagrangian (including weak interaction couplings) by performing numerical simulations in the $\epsilon$--regime, i.e. at quark masses where the…
Chiral dynamics makes definitive predictions for the electromagnetic polarizabilities of hadrons near the chiral limit; but, agreement with experiment is tenuous in some cases. We provide an overview of lattice QCD methods to compute the…
Efficient algorithms for the solution of partial differential equations on parallel computers are often based on domain decomposition methods. Schwarz preconditioners combined with standard Krylov space solvers are widely used in this…
Computing the trace of the inverse of large matrices is typically addressed through statistical methods. Deflating out the lowest eigenvectors or singular vectors of the matrix reduces the variance of the trace estimator. This work…
Lattice QCD using fermions whose Dirac operator obeys the Ginsparg-Wilson relation, is perhaps the best known formulation of QCD with a finite cutoff. It reproduces all the low energy QCD phenomenology associated with chiral symmetry at…
Multilevel techniques are efficient approaches for solving the large linear systems that arise from discretized partial differential equations and other problems. While geometric multigrid requires detailed knowledge about the underlying…
In this paper, we propose a class of adaptive multiresolution (also called adaptive sparse grid) ultra-weak discontinuous Galerkin (UWDG) methods for solving some nonlinear dispersive wave equations including the Korteweg-de Vries (KdV)…
We present numerical results for the 2-flavour Schwinger model with dynamical chiral lattice fermions. We insert an approximately chiral hypercube Dirac operator into the overlap formula to construct the overlap hypercube operator. This is…
We present a new multigrid scheme for solving the Poisson equation with Dirichlet boundary conditions on a Cartesian grid with irregular domain boundaries. This scheme was developed in the context of the Adaptive Mesh Refinement (AMR)…
In this paper, we introduce the overlap Dirac operator, which satisfies the Ginsparg-Wilson relation, to the matter sector of two-dimensional N=(2,2) lattice supersymmetric QCD (SQCD) with preserving one of the supercharges. It realizes the…
Efficient numerical solvers for sparse linear systems are crucial in science and engineering. One of the fastest methods for solving large-scale sparse linear systems is algebraic multigrid (AMG). The main challenge in the construction of…
It is possible to improve the chiral behaviour and the approach to the continuum limit of correlation functions in lattice QCD with standard and twisted Wilson fermions by taking arithmetic averages of correlators computed in theories…
We present simulation results for lattice QCD with light pions. For the quark fields we apply chirally symmetric lattice Dirac operators, in particular the overlap hypercube operator, along with the standard overlap operator for comparison.…
The multigrid-reduction-in-time (MGRIT) technique has proven to be successful in achieving higher run-time speedup by exploiting parallelism in time. The goal of this article is to develop and analyze a MGRIT algorithm, using FCF-relaxation…
We present evidence that multigrid works for wave equations in disordered systems, e.g. in the presence of gauge fields, no matter how strong the disorder, but one needs to introduce a "neural computations" point of view into large scale…
To numerically solve a generic elliptic equation on two-dimensional domains with rectangular Cartesian grids, we propose a cut-cell geometric multigrid method that features (1) general algorithmic steps that apply to two-dimensional…
Complete spectra of the staggered Dirac operator $\Dirac$ are determined in four-dimensional $SU(2)$ gauge fields with and without dynamical fermions. An attempt is made to relate the performance of multigrid and conjugate gradient…