---
title: 'GD–λAMG: Multigrid-Accelerated Davidson'
url: https://www.emergentmind.com/topics/multigrid-accelerated-davidson
type: topic
---

# GD–λAMG: Multigrid-Accelerated Davidson

Multigrid-accelerated Davidson (GD–λAMG) refers to a class of eigensolvers that combine the generalized Davidson method with a multigrid-accelerated inner iteration for efficiently computing extremal (typically smallest magnitude) eigenvalues and eigenvectors of large, sparse Hermitian matrices. This approach is particularly relevant in the spectral analysis of the Hermitian Wilson–Dirac operator in lattice quantum chromodynamics (QCD), where resolving the low-lying spectrum is crucial for various physical observables and for accelerating linear system solutions via deflation [2004.08146].

## 1. Problem Formulation and Context in Lattice QCD

The Hermitian Wilson–Dirac operator $Q$ on a four-dimensional lattice $\mathcal{L}$ of size $N_t\times N_s^3$ with quark mass $m_0$ is constructed as
\[
D\psi(x) = (m_0 + 4/a) \psi(x) - \frac{1}{2a} \sum_{\mu=0}^3 [ (I - \gamma_\mu)\otimes U_\mu(x)]\psi(x+\hat\mu) - \frac{1}{2a} \sum_{\mu=0}^3 [ (I + \gamma_\mu)\otimes U^\dagger_\mu(x-\hat\mu)] \psi(x-\hat\mu)
\]
with $\Gamma_5 = \gamma_0\gamma_1\gamma_2\gamma_3$. The Hermitian operator is defined as $Q := \Gamma_5 D$, so $Qx = \lambda x$ becomes the target eigenproblem with $x\in\mathbb{C}^N$, $\lambda\in\mathbb{R}$, specifically focusing on eigenpairs with small $|\lambda|$.

Small$-|\lambda|$ eigenpairs of $Q$ are essential in noise-reduction, all-to-all propagators, low-mode averaging, linear solve deflation, and in the study of spontaneous chiral symmetry breaking via the Banks–Casher relation. The spectral density near zero directly informs physical phenomena and the efficiency of many computational techniques [2004.08146].

## 2. Generalized Davidson Method: Core Algorithmic Components

The generalized Davidson method forms the outer iteration framework. An orthonormal basis $V_m = [v_1, \ldots, v_m]$ is maintained, spanning the current search space. At each step, harmonic Ritz pairs $(\theta, s)$ are found by the generalized eigenproblem
\[
(W_m^\dagger A V_m) s = \theta (W_m^\dagger V_m) s,\quad W_m = A V_m
\]
producing approximate eigenvalues $\theta$ and corresponding Ritz vectors $u = V_m s$.

The residual $r = Au - \theta u$ quantifies the approximation error. A new search direction $t$ is obtained by approximately solving the shifted equation (the correction equation)
\[
(A - \tau I)t = r
\]
with the shift $\tau$ typically chosen close to $\theta$. The new vector $t$ is orthogonalized against the current basis and any locked (converged) eigenvectors to expand the search space. This process is repeated until $k$ target eigenpairs are identified.

## 3. Multigrid-Accelerated DD–$\alpha$AMG Inner Iterations

At each Davidson step, the correction equation takes the form $(Q - \sigma I) y = r$ with shift $\sigma = \tau$, where $r$ is the residual. The algorithm employs left preconditioning with $\Gamma_5$, yielding $(D - \sigma \Gamma_5)t = \Gamma_5 r$, with $D = \Gamma_5 Q$.

A multigrid hierarchy is constructed using aggregation-based interpolation $P \in \mathbb{C}^{N\times N_c}$, built from $n_{tv}$ “test vectors” (approximations to the smallest$-|\lambda|$ modes). The restriction is $R = P^\dagger$, and the coarse operator is $Q_c = R Q P$. For the shifted system, $D_c - \sigma \Gamma_5^c$ with $D_c = P^\dagger D P$ and $\Gamma_5^c = P^\dagger \Gamma_5 P$ is used.

A two-level error propagator for $(Q-\sigma I)$ is constructed as
\[
E = (I - M(Q-\sigma I))^\nu \left[ I - P(Q_c-\sigma I)^{-1} R(Q-\sigma I) \right] (I - M(Q-\sigma I))^\mu
\]
with $M$ a smoother (e.g., SAP or GMRES) and $\mu,\nu$ specifying pre- and post-smoothing steps. The approximate inverse $M \approx (Q-\sigma I)^{-1}$ is realized by a single V-cycle: pre-smoothing, coarse correction, post-smoothing.

## 4. Dynamic Coarse Operator Updates and Local Coherence

“Local coherence” refers to the property that small$-|\lambda|$ eigenvectors exhibit similar structure on each aggregate. For a vector $v$, the coherence metric $lc(v) = \|P P^\dagger v\|/\|v\|$ is used; values near 1 indicate that the interpolation $P$ provides an accurate representation.

If the multigrid interpolation $P$ is constructed from the smallest$-|\lambda|$ eigenmodes, as the Davidson method advances and focuses on eigenvalues further from zero, $lc$ drops, resulting in reduced efficiency of the inner multigrid solve. GD–$\lambda$AMG addresses this by dynamically updating $P$ and $Q_c$. When new eigenvectors $x_j$ (corresponding to converged eigenvalues $\lambda_j$) become available, the $n_{tv}$ vectors with eigenvalues closest to the current target $\tau$ are selected to rebuild $P$—specifically,
\[
P = \text{block-diag}[\,x_j\,] \text{ aggregated over sublattices},\quad Q_c = R Q P,\, R = P^\dagger.
\]
This ensures that local coherence remains high throughout the Davidson iterations, maintaining the efficacy of the multigrid acceleration [2004.08146].

## 5. Detailed Algorithmic Workflow (GD–λAMG)

The overall algorithm proceeds as follows:

1. Initialize $V = [v_1, ..., v_{n_{tv}}]$ with random/test vectors; build $P$, $Q_c$; set $X = \varnothing$, $\Lambda = \varnothing$.
2. For each $m = n_{tv}+1, ...$:
    - Orthogonalize new guess $t$ to $V$ and $X$; normalize to obtain $v_m$; augment $V$.
    - Compute harmonic Ritz pairs $(\theta_i, s_i)$ from $(QV)^\dagger(QV)s = \theta (QV)^\dagger V s$; pick $\theta$ closest to 0 not yet locked.
    - Set $u = V s$, $r = Q u - \theta u$.
    - If $\|r\| \leq \varepsilon_\text{outer}$, lock $(\theta, u)$: $\Lambda \leftarrow \Lambda \cup \{\theta\}$, $X \leftarrow [X, u]$; select $n_{tv}$ locked vectors nearest to next $\tau$; rebuild $P$, $Q_c$; restart subspace (size $\leq m_{min}$).
    - Solve preconditioned correction by one V-cycle of DD–$\alpha$AMG: $t \leftarrow \text{MG\_solve}((D-\theta \Gamma_5)t = \Gamma_5 r)$.
3. Repeat until $k$ eigenpairs are obtained.

This process tightly couples the outer eigensolver with the inner multigrid-accelerated correction step, with frequent adaptation of the coarse-level components.

## 6. Performance Characteristics and Comparative Benchmarks

Performance tests on large lattices, e.g., $64 \times 32^3$, show that GD–$\lambda$AMG achieves substantial speedups over conventional eigensolvers. Key findings include:

- $\Gamma_5$-preconditioning yields an approximate $2\times$ reduction in solve time (e.g., 83 core-h to 41 core-h for 100 eigenpairs).
- DD–$\alpha$AMG preconditioning reduces inner iteration counts per eigenvalue from $300$–$1200$ (for pure GMRESR) to $30$ (for one V-cycle), corresponding to a $\sim7\times$ speed-up.
- Without interpolation/coarse-grid updates, cost scales as $k^{1.8}$ for $k$ eigenpairs; with updates, cost is near-linear in $k$, halving total time at $k=300$.
- Benchmarking against PARPACK (with optimal polynomial filter) and PRIMME (with the same AMG), GD–$\lambda$AMG is up to $10\times$ faster than PARPACK and $1.2$–$1.5\times$ faster than PRIMME (speedup increases with lattice size and $k$).
- For $k$ between $50$ and $1000$, GD–$\lambda$AMG maintains near-linear scaling, with the performance advantage over PRIMME increasing from $1.5\times$ to over $3\times$ [2004.08146].

## 7. Spectral Analysis and Physical Relevance

Computed low-lying eigenvalues enable construction of the spectral density $\rho(\lambda)$, a key observable in lattice QCD. For a $64 \times 32^3$ lattice with lattice spacing $a=0.071~\mathrm{fm}$ and $N_f=2$, the density from the smallest 1000 eigenvalues of $Q$ is compared to the Banks–Casher relation, which in the continuum (for $\lambda>m$, where $m$ is the bare quark mass, $m\approx8.22~\mathrm{MeV}$) predicts
\[
\rho(\lambda) \simeq \frac{2\lambda\Sigma}{\pi\sqrt{\lambda^2-m^2}}
\]
with $\Sigma$ the chiral condensate. The numerical results display nonzero density below $\lambda=m$ and a smooth interpolation near $m$ rather than the divergence expected in the infinite-volume, continuum limit. This quantifies the impact of lattice discretization and volume artifacts in practical QCD simulations [2004.08146].

In summary, GD–$\lambda$AMG couples harmonic-Davidson iteration with dynamically updated DD–$\alpha$AMG inner solves, exploiting local coherence and preconditioning to yield up to order-of-magnitude acceleration and near-linear scaling in both number of eigenpairs and lattice size, while preserving physically accurate spectral information critical for modern lattice QCD.

Source: https://www.emergentmind.com/topics/multigrid-accelerated-davidson