Abstract

We develop an algebraic multigrid method for solving the non-Hermitian Wilson discretization of the two-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_5$-symmetry of the Dirac equation, which carries over to the Wilson discretization. This discrete $\gamma_5$-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_5$-symmetry. ...

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call

Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.