Adaptive Aggregation Based Domain Decomposition Multigrid for the Lattice Wilson Dirac Operator
read the original abstract
In lattice QCD computations a substantial amount of work is spent in solving discretized versions of the Dirac equation. Conventional Krylov solvers show critical slowing down for large system sizes and physically interesting parameter regions. We present a domain decomposition adaptive algebraic multigrid method used as a precondtioner to solve the "clover improved" Wilson discretization of the Dirac equation. This approach combines and improves two approaches, namely domain decomposition and adaptive algebraic multigrid, that have been used seperately in lattice QCD before. We show in extensive numerical test conducted with a parallel production code implementation that considerable speed-up over conventional Krylov subspace methods, domain decomposition methods and other hierarchical approaches for realistic system sizes can be achieved.
This paper has not been read by Pith yet.
Forward citations
Cited by 9 Pith papers
-
Complete lattice QCD calculation of $K^{-}\to \ell^{-}\bar{\nu}_{\ell}\ell^{'+}\ell^{'-}$ form factors
First complete lattice QCD determination of the four structure-dependent form factors for K- → ℓ- ν̄ℓ ℓ'+ ℓ'- decays at physical quark masses with controlled statistical and systematic errors.
-
Charmonium radiative transitions to dileptons from lattice QCD: The case of $h_c \to \eta_c \ell^+\ell^-$ and $\chi_{c1} \to J/\psi\,\ell^+\ell^-$
First fully dynamical lattice QCD yields Γ(h_c → η_c e⁺e⁻) = 5.45(19) keV (3σ above BESIII) and Γ(χ_c1 → J/ψ e⁺e⁻) = 2.869(90) keV, with continuum-extrapolated results and q² distributions.
-
Scalar and Tensor Form Factors for $\Lambda \rightarrow p\ell \bar{\nu}_\ell$ from Lattice QCD
Lattice QCD yields the scalar and tensor form factors for Λ→pℓν̄ℓ as functions of q², providing a model-independent input to constrain non-standard charged-current interactions via the predicted R^{μe} ratio compared ...
-
Mellin Moments of Pion and Kaon Unpolarized PDFs from Nonlocal Operators in Lattice QCD
First lattice QCD results for Mellin moments of pion and kaon unpolarized PDFs from nonlocal operator matrix elements on a 32^3×64 twisted-mass ensemble at m_π=260 MeV, extracted at NNLO and μ=2 GeV.
-
Mellin Moments of the Unpolarized Gluon PDF in the Proton from Nonlocal Operators in Lattice QCD
Lattice QCD extracts the ratio of the third to first Mellin moment of the gluon PDF at 2 GeV from nonlocal operators on an Nf=2+1+1 ensemble.
-
A Framework for Quantum Simulations of Energy-Loss and Hadronization in Non-Abelian Gauge Theories: SU(2) Lattice Gauge Theory in 1+1D
A quantum simulation framework is developed and demonstrated for energy loss and hadronization of a heavy quark in 1+1D SU(2) lattice gauge theory on 18 qubits of IBM hardware, with results matching classical simulations.
-
Pion and Kaon PDFs from Lattice QCD via Large Momentum Effective Theory and Short-Distance Factorization
Lattice QCD calculation of pion and kaon unpolarized quark PDFs on a 32^3×64 ensemble with 260 MeV pion mass, using LaMET and SDF matching at boosts up to 2.07 GeV.
-
Performance of Low Mode Averaging on Twisted-Mass Fermion Ensembles at the physical pion mass point
Numerical tests show low-mode averaging reduces noise in light meson and baryon observables on physical-pion twisted-mass ensembles, yielding renormalized chiral condensate 269.5(4.5) MeV.
-
Variance reduction strategies for lattice QCD
Variance reduction schemes based on decompositions of quark propagators have proven useful for precision lattice QCD observables and may help reduce the computational cost of reaching large volumes.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.