Related papers: Multigrid solver for clover fermions
We compare the quenched hadron spectrum on blocked and unblocked lattices for the Wilson quark action, the clover action and the tadpole-improved clover action. The latter gives a spectrum markedly closer to the original one, even though…
In Hybrid Monte Carlo simulations for full QCD, the gauge fields evolve smoothly as a function of Molecular Dynamics time. Here we investigate improved methods of estimating the trial or starting solutions for the Dirac matrix inversion as…
We calculate the Sheikholeslami-Wohlert coefficient of the O(a) improvement term for Wilson fermions in any representation of the gauge group SU(N) perturbatively at the one-loop level. The result applies to QCD with adjoint quarks and to…
We propose a new two-grid approach based on Bellman-Kalaba quasilinearization and Axelsson-Xu finite element two-grid method for the solution of singularly perturbed reaction-diffusion equations. The algorithms involve solving one…
In the light-front formulation of field theory, it is possible to write down a chirally invariant mass term. It thus appears as if one could solve the species doubling problem on a light-front quantized transverse lattice in a chirally…
The past decade has witnessed a dramatic acceleration of lattice quantum chromodynamics calculations in nuclear and particle physics. This has been due to both significant progress in accelerating the iterative linear solvers using…
We discuss the implementation of a Sheikholeslami-Wohlert term for simulations of lattice QCD with dynamical Wilson fermions as required by Symanzik's improvement program. We show that for the Hybrid Monte Carlo or Kramers equation…
We present a detailed design of a (P)HMC simulation algorithm for N_f = 2 + 1 + 1 maximally twisted Wilson quark flavours. The algorithm retains even/odd and mass-shift preconditionings combined with multiple Molecular Dynamics time scales…
We report calculations of the parameter B_K appearing in the Delta S=2 neutral kaon mixing matrix element, whose uncertainty limits the power of unitarity triangle constraints for testing the standard model or looking for new physics. We…
We report on calculations of the charmonium and bottomonium spectrum in lattice QCD. We use ensembles of gauge fields with three flavors of sea quarks, simulated with the asqtad improved action for staggered fermions. For the heavy quarks…
It is shown how various ideas that are well established for the solution of Poisson's equation using plane wave and multigrid methods can be combined with wavelet concepts. The combination of wavelet concepts and multigrid techniques turns…
The combination of a non-overlapping Schwarz preconditioner and the Hybrid Monte Carlo (HMC) algorithm is shown to yield an efficient simulation algorithm for two-flavour lattice QCD with Wilson quarks. Extensive tests are performed, on…
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…
A Maxwell solver derived from finite element method with \mathcal{O}(N) computing cost is developed to improve the numerical dispersion properties in relativistic particle-in-cell (PIC) simulations. The correction of the dispersion relation…
Representing massless Dirac fermions on a spatial lattice poses a potential challenge known as the Fermion Doubling problem. Addition of a quadratic term to the Dirac Hamiltonian circumvents this problem. We show that the modified…
Computing quark propagators with overlap fermions requires the solution of a shifted unitary linear system. Jagels and Reichel have shown that for such systems it is possible to construct a minimal residual algorithm by short recurrences.…
It is shown the central field Dirac equation can be simplified through the use of real conjugate spinors to substitute for the upper and lower components of the bi-spinor eigensolutions. This substitution reduces the Dirac equation for the…
We introduce a class of efficient multiple right-hand side multigrid algorithm for domain wall fermions. The simultaneous solution for a modest number of right hand sides concurrently allows for a significant reduction in the time spent…
The Helmholtz equation arises when modeling wave propagation in the frequency domain. The equation is discretized as an indefinite linear system, which is difficult to solve at high wave numbers. In many applications, the solution of the…
Iterative solvers preconditioned with algebraic multigrid have been devised as an optimal technology to speed up the response of large sparse linear systems. In this work, this technique was implemented in the framework of the dual…