Related papers: Adaptive Multigrid Algorithm for the QCD Dirac-Wil…
MGARD (MultiGrid Adaptive Reduction of Data) is an algorithm for compressing and refactoring scientific data, based on the theory of multigrid methods. The core algorithm is built around stable multilevel decompositions of conforming…
NOTE: this is a shortened version of the abstract of the paper. Multigrid methods for propagators in gauge fields are investigated. Gauge fields are incorporated in algorithms in a covariant way. This avoids the necessity for gauge fixing…
An algorithm for generating optimal nonuniform grids for solving the two-body Schr\"odinger equation is developed and implemented. The shape of the grid is optimized to accurately reproduce the low-energy part of the spectrum of the…
New exact upper and lower bounds are derived on the spectrum of the square of the hermitian Wilson Dirac operator. It is hoped that the derivations and the results will be of help in the search for ways to reduce the cost of simulations…
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 consider a Wilson-Dirac operator with improved chiral properties. We show that, for arbitrarily rough gauge fields, it satisfies the index theorem if we identify the zero modes with the small real eigenvalues of the fermion operator and…
The aim of this paper is to design an efficient multigrid method for constrained convex optimization problems arising from discretization of some underlying infinite dimensional problems. Due to problem dependency of this approach, we only…
We present a convergent and scalable multigrid solver for high-frequency Helmholtz equations. Standard multigrid methods do not converge for high-frequency Helmholtz problems, and a common cure is adding a complex shift and using the…
We prove quadratic estimates for complex perturbations of Dirac-type operators, and thereby show that such operators have a bounded functional calculus. As an application we show that spectral projections of the Hodge--Dirac operator on…
In this article we consider two-grid finite element methods for solving semilinear interface problems in d space dimensions, for d=2 or d=3. We first describe in some detail the target problem class with discontinuous diffusion…
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…
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…
A new non-Abelian gauge transformation for two-forms is introduced. Construction is based on a fixed map from the spacetime to the loop space which attachs a closed loop to each point of the spacetime. It is argued that this set-up is…
In this article we investigate and solve exactly the modified Dirac oscillator in curved spacetime with spin and pseudospin symmetries through an algebraic approach. By focusing on the radial part of this problem, we use the Schr\"odinger…
We consider supersymmetric quantum mechanical systems in arbitrary dimensions on curved spaces with nontrivial gauge fields. The square of the Dirac operator serves as Hamiltonian. We derive a relation between the number of supercharges…
An idealized multigrid algorithm for the computation of propagators of staggered fermions is investigated. Exemplified in four-dimensional $SU(2)$ gauge fields, it is shown that the idealized algorithm preserves criticality under…
The intertwining technique has been widely used to study the Schr\"odinger equation and to generate new Hamiltonians with known spectra. This technique can be adapted to find the bound states of certain Dirac Hamiltonians. In this paper the…
At physical light quark masses, efficient linear solvers are crucial for carrying out the millions of inversions of the Dirac matrix required for obtaining high statistics in quark correlation functions. Adaptive algebraic multi-grid…
We report on an ongoing project to parametrize the Fixed-Point Dirac operator for massless quarks, using a very general construction which has arbitrarily many fermion offsets and gauge paths, the complete Clifford algebra and satisfies all…
Two new relaxation schemes are proposed for the smoothing step in the geometric multigrid solution of PDEs on 2D and 3D stretched structured grids. The new schemes are characterized by efficient line relaxation on branched sets of lines of…