Adaptive Aggregation Based Domain Decomposition Multigrid for the Lattice Wilson Dirac Operator
arXiv:1303.1377
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.
Additional comparison to method of arXiv:1011.2775 and to mixed-precision odd-even preconditioned BiCGStab. Results of numerical experiments changed slightly due to more systematic use of odd-even preconditioning
References in corpus (6)
- Ab-initio Determination of Light Hadron Masses
- Adaptive multigrid algorithm for the lattice Wilson-Dirac operator
- Lattice QCD at the physical point: light quark masses
- Nucleon mass and sigma term from lattice QCD with two light fermion flavors
- Isospin splittings in the light baryon octet from lattice QCD and QED
- Adaptive Multigrid Algorithm for Lattice QCD
Cited by in corpus (7)
- Simulation of QCD with N_f=2+1 flavors of non-perturbatively improved Wilson fermions
- Lattice QCD with Domain Decomposition on Intel Xeon Phi Co-Processors
- Kenneth Wilson and lattice QCD
- An adaptive aggregation based domain decomposition multilevel method for the lattice wilson dirac operator: multilevel results
- Hierarchically deflated conjugate residual
- Adaptive algebraic multigrid on SIMD architectures
- Lattice QCD with 2+1 Flavors and Open Boundaries: First Results of the Baryon Spectrum