Coarsest-level improvements in multigrid for lattice QCD on large-scale computers

05/18/2022
by   Jesus Espinoza-Valverde, et al.
0

Numerical simulations of quantum chromodynamics (QCD) on a lattice require the frequent solution of linear systems of equations with large, sparse and typically ill-conditioned matrices. Algebraic multigrid methods are meanwhile the standard for these difficult solves. Although the linear systems at the coarsest level of the multigrid hierarchy are much smaller than the ones at the finest level, they can be severely ill-conditioned, thus affecting the scalability of the whole solver. In this paper, we investigate different novel ways to enhance the coarsest-level solver and demonstrate their potential using DD-αAMG, one of the publicly available algebraic multigrid solvers for lattice QCD. We do this for two lattice discretizations, namely clover-improved Wilson and twisted mass. For both the combination of two of the investigated enhancements, deflation and polynomial preconditioning, yield significant improvements in the regime of small mass parameters. In the clover-improved Wilson case we observe a significantly improved insensitivity of the solver to conditioning, and for twisted mass we are able to get rid of a somewhat artificial increase of the twisted mass parameter on the coarsest level used so far to make the coarsest level solves converge more rapidly.

READ FULL TEXT
research
08/14/2017

Eigenvalue Solvers for Modeling Nuclear Reactors on Leadership Class Machines

Three complementary methods have been implemented in the code Denovo tha...
research
03/28/2022

Acceleration of multiple precision solver for ill-conditioned algebraic equations with lower precision eigensolver

There are some types of ill-conditioned algebraic equations that have di...
research
06/25/2021

Linear solvers for power grid optimization problems: a review of GPU-accelerated linear solvers

The linear equations that arise in interior methods for constrained opti...
research
11/09/2022

Solution of time-harmonic Maxwell's equations by a domain decomposition method based on PML transmission conditions

Numerical discretization of the large-scale Maxwell's equations leads to...
research
08/30/2023

Towards Robust Solvers for Nuclear Fusion Simulations Using JOREK: A Numerical Analysis Perspective

One of the most well-established codes for modeling non-linear Magnetohy...
research
02/14/2023

Are sketch-and-precondition least squares solvers numerically stable?

Sketch-and-precondition techniques are popular for solving large least s...

Please sign up or login with your details

Forgot password? Click here to reset