2013/03/06 by Andreas Frommer, Karsten Kahl, Frommer, Andreas +7 · 8 citations
Engineering · Computer Science · Mathematics · #Advanced Numerical Methods in Computational Mathematics #Matrix Theory and Algorithms #Numerical methods for differential equations
paper · pdf · doi:10.48550/arxiv.1303.1377
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.