Bootstrap algebraic multigrid for the 2D Wilson Dirac system

From MaRDI portal



Abstract: We develop an algebraic multigrid method for solving the non-Hermitian Wilson discretization of the 2-dimensional Dirac equation. The proposed approach uses a bootstrap setup algorithm based on a multigrid eigensolver. It computes test vectors which define the least squares interpolation operators by working mainly on coarse grids, leading to an efficient and integrated self learning process for defining algebraic multigrid interpolation. The algorithm is motivated by the gamma-symmetry of the Dirac equation, which carries over to the Wilson discretization. This discrete gamma-symmetry is used to reduce a general Petrov Galerkin bootstrap setup algorithm to a Galerkin method for the Hermitian and indefinite formulation of the Wilson matrix. Kaczmarz relaxation is used as the multigrid smoothing scheme in both the setup and solve phases of the resulting Galerkin algorithm. The overall method is applied to the odd-even reduced Wilson matrix, which also fulfills the discrete gamma-symmetry. Extensive numerical results are presented to motivate the design and demonstrate the effectiveness of the proposed approach.











This page was built for publication: Bootstrap algebraic multigrid for the 2D Wilson Dirac system

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2878959)