Volume 396 - The 38th International Symposium on Lattice Field Theory (LATTICE2021) - Oral presentation
Implementation of Simultaneous Inversion of a Multi-shifted Dirac Matrix for Twisted-Mass Fermions within DD$\alpha$AMG
S. Yamamoto*, S. Bacchio and J. Finkenrath
Full text: pdf
Pre-published on: May 16, 2022
Published on:
Abstract
At physical light quark masses, efficient linear solvers are crucial for carrying out the millions of inversions of the Dirac matrix required for obtaining high statistics in quark correlation functions. Adaptive algebraic multi-grid methods have proven to be very efficient in such cases, exhibiting mild critical slowing down towards very light quark masses and outperforming traditional solver methods, such as the conjugate gradient method, at the physical point.
We will discuss our implementations of simultaneous inversion of a (degenerate) Dirac matrix for twisted-mass fermions for multiple right-hand-sides (rhs) with multi-shifts and block-Krylov solvers. The implementation is carried out within the community library DD$\alpha$AMG, which implements aggregation-based Domain Decomposition adaptive algebraic multi-grid methods. The block-Krylov solvers are provided via the Fast Accurate Block Linear krylOv Solver (Fabulous) library and can be used at coarser levels.
Our code inverts Dirac matrices with different twisted-mass terms and for multiple rhs simultaneously and is thus also suitable for components within a typical lattice QCD simulation workflow, such as the rational approximation. We show preliminary results on scalability and compare the performance of our implementation when using different Block-Krylov solver techniques.
DOI: https://doi.org/10.22323/1.396.0536
How to cite

Metadata are provided both in "article" format (very similar to INSPIRE) as this helps creating very compact bibliographies which can be beneficial to authors and readers, and in "proceeding" format which is more detailed and complete.

Open Access