---
title: Hiptmair–Xu Preconditioners
url: https://www.emergentmind.com/topics/hiptmair-xu-preconditioners
type: topic
---

# Hiptmair–Xu Preconditioners

Hiptmair--Xu preconditioners are nodal auxiliary-space preconditioners for discrete \(H(\curl)\) and \(H(\div)\) operators. Their defining idea is to replace a difficult edge- or face-element solve by an additive combination of a simple smoother on the original space, one or more scalar or vector \(H^1\)-type auxiliary solves in nodal spaces, and exact-sequence corrections through discrete gradient or curl maps. In the positive Maxwell regime this produces coercive preconditioned operators, while in semidefinite curl-curl problems the same architecture resolves the large gradient kernel through auxiliary-space splitting rather than direct inversion of the singular operator [2306.13217][1708.05853].

## 1. Core auxiliary-space architecture

In the positive Maxwell setting on a bounded polyhedral domain \(\Omega\subset\mathbb R^3\), the discrete operator has the curl-curl-plus-mass form
\[
\langle \mathscr M(u),v\rangle
=
\int_\Omega \boldsymbol{curl}(u)\cdot \boldsymbol{curl}(v)+\gamma^2 u\cdot v\,dx,
\qquad u,v\in W(\Omega),
\]
with \(W(\Omega)\) the lowest-order Nédélec edge element space. The added positive mass term makes the operator coercive on \(W(\Omega)\), which is precisely the regime in which HX is classically formulated as a positive Maxwell preconditioner [2306.13217].

A standard HX formula recalled in later work is
\[
\tilde{\mathscr M}^{-1}
=
\operatorname{diag}(\mathscr M)^{-1}
+
G_\Omega \mathscr L^{-1} G_\Omega^*
+
\sum_{j=1,2,3} \Pi_\Omega^{e_j}\mathscr L^{-1}(\Pi_\Omega^{e_j})^* .
\]
Here \(\operatorname{diag}(\mathscr M)^{-1}\) is a Jacobi smoother on the edge space, \(G_\Omega:V(\Omega)\to W(\Omega)\) is the discrete gradient from the nodal \(P_1\) space, and each \(\Pi_\Omega^{e_j}\) is a nodal-to-edge transfer induced by edge interpolation. This realizes the characteristic HX decomposition into one edge-space smoothing term, one gradient correction, and three nodal lifting terms [2306.13217].

In the weighted Maxwell formulation with jump coefficients, the same three-part pattern is written as
\[
\mathbf B_h
=
\mathbf J_h^{-1}
+
\hat{\mathbf r}_h(-\mathbf \Delta_h)^{-1}\hat{\mathbf r}_h^\ast
+
\mathbf T_h(\mathbf A_h^\nabla)^{-1}\mathbf T_h^\ast .
\]
The factors are the Jacobi smoother \(\mathbf J_h^{-1}\), a vector nodal Laplacian solve transferred by the edge interpolation \(\hat{\mathbf r}_h\), and a gradient-space correction through the restriction of \(\mathbf A_h\) to \(\nabla Z_h(\Omega)\). The paper emphasizing jump coefficients states that implementation requires solving four scalar elliptic problems, which is one of the practical attractions of the method [1708.05853].

A closely related formula appears in recent time-harmonic practice for a nearby positive Maxwell block:
\[
H_{\text{HX}}
=
S^{-1}
+
P_{\text{curl}}(L+k^2M)^{-1}P_{\text{curl}}^T
+
(k^2)^{-1}G(-\Delta)^{-1}G^T .
\]
This variant is not a new HX theory; it is an implementation of the same auxiliary-space principle in a real-valued split formulation, with a smoother, a nodal/vector-Laplacian correction, and a scalar-gradient correction [2507.13066].

## 2. Regular decomposition, exact sequences, and robustness

The analytical foundation of HX is a discrete regular decomposition. In the edge-element setting one seeks, for \( \mathbf v_h\in V_h(\Omega)\), a splitting of the form
\[
\mathbf v_h=\nabla p_h+\mathbf r_h\mathbf w_h+\mathbf R_h ,
\]
with \(p_h\) in a nodal space, \(\mathbf w_h\) in a vector nodal space, and \(\mathbf R_h\) a remainder controlled by smoothing. Later auxiliary-space analyses identify this decomposition as the exact mechanism matching the three blocks of the preconditioner [1708.05850][1708.05853].

For heterogeneous Maxwell problems with strongly discontinuous coefficients, the relevant norms are coefficient-weighted:
\[
\|\mathbf v\|_{H^\ast(\mathrm{curl};\Omega)}
:=
\Big(\sum_{r=1}^{N_0}\alpha_r\|\nabla\times \mathbf v\|^2_{0,\Omega_r}
+
\beta_r\|\mathbf v\|^2_{0,\Omega_r}\Big)^{1/2}.
\]
The key difficulty is that classical unweighted decompositions do not control the effect of large jumps in \(\alpha\) and \(\beta\). The weighted theory therefore introduces upper intersection sets
\[
\mathcal S_k
=
\bigcup_{\alpha_l\ge \alpha_k,\ l\ne k}(\partial\Omega_k\cap\partial\Omega_l)
\ \cup\
(\partial\Omega_k\cap\partial\Omega),
\]
together with the geometric notion of thorny vertices, which are the obstruction to full-space weighted stability [1708.05853].

The no-thorny case yields the sharpest HX results. If the coefficients belong to the class \(\mathbb P_M(\tau)\), every \(\mathbf v_h\in V_h(\Omega)\) admits
\[
\mathbf v_h=\nabla p_h+\mathbf r_h\mathbf w_h+\mathbf R_h
\]
with
\[
\|p_h\|_{H^1_\beta(\Omega)}
\ \stlm\
\rho^{m_0}(h)\,\|\mathbf v_h\|_{H^\ast(\mathrm{curl};\Omega)},
\]
and
\[
\|\mathbf w_h\|_{H^1_\ast(\Omega)}+h^{-1}\|\mathbf R_h\|_{L^2_\beta(\Omega)}
\ \stlm\
\rho^{m_0}(h)\,\|\mathbf v_h\|_{H^\ast(\mathrm{curl};\Omega)}.
\]
This leads to
\[
\mathrm{cond}(\mathbf B_h^{-1}\mathbf A_h)\ \stlm\ \rho^{m_0}(h),
\]
and in the best geometric case,
\[
\mathrm{cond}(\mathbf B_h^{-1}\mathbf A_h)\ \stlm\ 1.
\]
The hidden constants depend on geometry, mesh regularity, and the class parameter \(\tau\), but not on the jump magnitudes themselves [1708.05853].

The companion decomposition paper develops the local ingredients behind these estimates. It proves regular Helmholtz-type decompositions for lowest-order Nédélec functions on ordinary polyhedra, on unions of polyhedra meeting at an edge or at a vertex, and under constraints prescribing zero tangential degrees of freedom on selected faces and edges. The characteristic estimates are of the form
\[
\|\mathbf w_h\|_{1,G}+h^{-1}\|{\bf R}_h\|_{0,G}
\lesssim
\log(1/h)\|\mathbf v_h\|_{H(\curl;G)},
\]
with the logarithmic factor removable in several favorable cases. This provides the interface-compatible splitting toolkit used later in the weighted HX convergence theory [1708.05850].

## 3. Semidefinite curl-curl operators and saddle-point systems

HX is not restricted to strictly positive operators. An auxiliary-space theory for semidefinite linear systems formulates the preconditioner abstractly as
\[
B=\Pi\,\widetilde B\,\Pi^t,
\]
with \(\Pi:\widetilde V\to V\) surjective and \(\widetilde B\) simpler than \(B\). For semi-SPD \(A\), the theory replaces ordinary condition-number formulas by range-restricted identities involving an infimum over \(\phi\in\mathcal N(A)\). This makes the kernel treatment explicit rather than implicit [2509.07179].

In that framework, the \(H(\curl)\) HX preconditioner for the SPD case \(\epsilon>0\) is written as
\[
B_{\curl}
=
D_{\curl}^{-1}
+
\Pi_h^{\curl} A_{\grad^3}^{-1} (\Pi_h^{\curl})^t
+
\epsilon^{-1}\grad A_{\grad}^{-1}\grad^t,
\]
while in the semidefinite case \(\epsilon=0\) the gradient correction is dropped:
\[
B_{\curl}
=
D_{\curl}^{-1}
+
\Pi_h^{\curl} A_{\grad^3}^{-1} (\Pi_h^{\curl})^t .
\]
The paper proves
\[
\kappa(B_{\curl}A_{\curl})\lesssim 1
\]
and, analogously for \(H(\div)\),
\[
\kappa(B_{\div}A_{\div})\lesssim 1,
\]
thereby recovering uniform optimality in both SPD and semidefinite regimes through the same auxiliary-space variational identities [2509.07179].

A distinct semidefinite application arises in the mixed finite element discretization of time-harmonic Maxwell saddle-point systems at vanishing wave number. For
\[
\mathcal A
=
\begin{pmatrix}
A & B^T\\
B & 0
\end{pmatrix},
\]
with \(A\) the singular discrete curl-curl matrix, \(C\) spanning \(\ker(A)\), and \(L=BC\), the exact inverse formula yields
\[
p=L^{-1}C^T f,
\]
and reduces the saddle-point problem to
\[
Au=(I-B^TL^{-1}C^T)f,\qquad Bu=g.
\]
The crucial point is that the right-hand side lies in \(\mathcal R(A)\), so the singular curl-curl equation is consistent. The paper then states that one may compute an arbitrary particular solution \(u_0\) by CG with the Hiptmair-Xu preconditioner and recover the full mixed solution by explicit kernel corrections. This is not a new HX construction for the full saddle-point matrix; it is an exact algebraic reduction showing when HX is the natural solver for the singular curl-curl block [1610.03196].

## 4. Substructured and interface HX preconditioners

Classical HX is a volume preconditioner, but it can be transferred to a nonoverlapping substructuring setting. For positive Maxwell problems, a central exact identity for the scalar Schur complement is
\[
(R^*T_LR)^{-1}=\mathscr B\mathscr L^{-1}\mathscr B^*,
\]
where \(T_L^{-1}=BL^{-1}B^*\) is the block-diagonal local Schur inverse and \(R\) is the interface assembly operator. This identity says that tracing the inverse of the global conforming scalar volume operator to the skeleton gives exactly the inverse of the assembled scalar Schur complement, not merely a spectrally equivalent approximation [2306.13217].

Once this identity is available, the volume HX decomposition can be pushed term-by-term to the interface. The resulting skeleton analogue is
\[
\underline{\mathscr B}\,\tilde{\mathscr M}^{-1}\,\underline{\mathscr B}^*
=
\underline{\mathscr B}\,\operatorname{diag}(\mathscr M)^{-1}\,\underline{\mathscr B}^*
+
G_\Sigma(R^*T_LR)^{-1}G_\Sigma^*
+
\sum_{j=1,2,3}(\Pi_\Sigma^{e_j})(R^*T_LR)^{-1}(\Pi_\Sigma^{e_j})^* .
\]
Replacing the exact scalar Schur inverse by a Neumann–Neumann preconditioner
\[
Q_{nn}
=
D_{\Sigma}^{-1}R^{*}DT_{L}^{-1}DRD_{\Sigma}^{-1}
\]
gives the practical substructured HX preconditioner
\[
Q_{hx}
=
\underline{\mathscr B}\operatorname{diag}(\mathscr M)^{-1}\underline{\mathscr B}^{*}
+
G_{\Sigma}Q_{nn}G_{\Sigma}^*
+
\sum_{j=1,2,3}(\Pi_\Sigma^{e_j})Q_{nn}(\Pi_\Sigma^{e_j})^* .
\]
The architecture is recognizably HX, but every nodal inverse has been replaced by a scalar interface Schur-complement preconditioner [2306.13217].

The corresponding estimate is
\[
\operatorname{cond}(Q_{hx}\,\underline R^*T_M\underline R)
\le
\operatorname{cond}(Q_{nn}R^*T_LR)\cdot
\operatorname{cond}(\tilde{\mathscr M}^{-1}\mathscr M).
\]
This shows that the interface Maxwell preconditioner inherits the HX quality factor from the volume problem and multiplies it by the scalar Schur-complement quality factor. The same paper explicitly notes that because the scalar ingredient is only one-level Neumann–Neumann, the method is not expected to be fully scalable in the strong sense; adding a coarse space is identified as the natural improvement [2306.13217].

## 5. Generalizations to nonstandard complexes and meshes

The HX paradigm has been extended well beyond lowest-order simplicial volume discretizations. For lowest-order mixed virtual element methods on polytopal meshes, the discrete complex is exact,
\[
dV_h^k=\ker_d(V_h^{k+1}),
\]
the cochain projection commutes,
\[
\Pi_h^{k+1}d=d\Pi_h^k,
\]
and one has a VEM regular decomposition
\[
v_h=\tilde v_h+\Pi_h^k\psi_h+dp_h
\]
with
\[
\|h^{-1}\tilde v_h\|+\|\psi_h\|_1\lesssim \|dv_h\|,
\qquad
\|p_h\|_{V^{k-1}}\lesssim \|v_h\|_{V^k}.
\]
This leads, in the 3D edge case, to the direct HX analogue
\[
B=(S^1)^{-1}+\Pi_h^1A_1^{-1}(\Pi_h^1)^*+\nabla A_2^{-1}\nabla^*,
\]
and to an \(H(\div)\) facet preconditioner with a recursive \(\curl(\cdot)\curl^*\) structure. The spectral condition number is proved to be bounded independently of \(h\), while robustness with respect to high aspect ratio is reported as a numerical observation rather than a theorem [2404.12823].

On triangulated surfaces, the geometry changes the auxiliary space itself. A central observation is that any vector field tangential to a triangulated surface \(\mathcal M_h\) cannot be continuous. The surface HX construction therefore uses the ambient nodal space \([H_h(\nabla_h)]^3\) together with a modified interpolation \(\bar\pi_h^{\nabla\cdot}\) for the \(H(\div)\) branch. The additive surface preconditioner is
\[
B_h^{\mathrm d}
=
S_h^{\mathrm d}
+
\bar\pi_h^{\mathrm d}\big[(A_h^\nabla)^{-1}\big]^3(\bar\pi_h^{\mathrm d})'
+
c^{-1}\,\mathrm d_h^-(A_h^{\mathrm d^-})^{-1}(\mathrm d_h^-)'
\]
and satisfies
\[
(A_h^{\mathrm d})^{-1}\lesssim B_h^{\mathrm d}\lesssim (A_h^{\mathrm d})^{-1},
\qquad
\kappa(B_h^{\mathrm d}A_h^{\mathrm d})\lesssim 1.
\]
For 2D surfaces the method uses four scalar surface Laplacian inverses: three from the ambient nodal-vector block and one from the exact-sequence term [2107.07978].

A further extension replaces the fixed-dimensional de Rham complex by a mixed-dimensional one on networks of manifolds. There the mixed-dimensional regular decomposition takes the form
\[
\mathfrak q_h=\Pi_h^k \mathfrak a_h+\mathfrak b_h+\mathfrak d \mathfrak f_h,
\]
with
\[
\|\mathfrak a_h\|_{H^1\mathfrak L^k}
+
\|h^{-1}\mathfrak b_h\|_{L^2\mathfrak L^k}
+
\|\mathfrak f_h\|_{H\mathfrak L^{k-1}}
\lesssim
\|\mathfrak q_h\|_{H\mathfrak L^k}.
\]
The corresponding mixed-dimensional HX preconditioner is
\[
\mathfrak B^k
=
(\mathfrak S^k)^{-1}
+
\Pi_h^k (\mathfrak A^k_{\mathrm{reg}})^{-1}(\Pi_h^k)^*
+
\mathfrak d (\mathfrak A^{k-1})^{-1}\mathfrak d^*,
\]
with
\[
\kappa(\mathfrak B^k\mathfrak A^k)\lesssim 1.
\]
This suggests that the HX mechanism is fundamentally a de Rham-complex construction rather than a feature specific to standard 3D \(H(\curl)\) discretizations [1910.04704].

## 6. Practical use, indefinite Maxwell problems, and abstract operator-preconditioning

In large-scale time-harmonic Maxwell computations, HX is often used indirectly rather than as a preconditioner for the original indefinite and non-Hermitian operator. In one recent parallel study, the complex system \(A=C-M-iB\) is rewritten in split real-imaginary form and preconditioned by the block-diagonal positive surrogate
\[
P=
\begin{pmatrix}
C+M+B & 0\\
0 & C+M+B
\end{pmatrix}.
\]
Each positive block is then approximately inverted by hypre’s AMS, which the paper identifies as HX in that software stack. The outer solver is right-preconditioned FGMRES because the inner block inverses are computed inexactly [2507.13066].

That study reports two HX variants. In the first, each positive block solve uses CG preconditioned by HX with maximum \(20\) inner iterations and relative residual \(10^{-2}\). In the second, AMS is used directly as the approximate inner solver. For \(k=1\), the outer FGMRES iteration counts for the inner-CG variant are \(28,26,26,26,26\) across refinements from \(2\times316\) to \(2\times4073\) unknowns, which the paper interprets as essentially mesh-independent. For \(k=10\), the outer counts remain roughly mesh-independent, around \(154\)–\(174\), but the method deteriorates with increasing wavenumber and even more strongly with increasing physical domain size. The authors state that the number of FGMRES iterations is at least proportional to \(k^{0.86}\), and on an enlarged-domain test the iteration count is multiplied by \(28\) relative to the small-domain case [2507.13066].

These observations clarify a common misconception. HX is tailored to positive Maxwell equations, not to the original indefinite time-harmonic scattering operator. In that setting it is therefore being used as inner technology for a nearby positive block problem, not as a direct preconditioner for \(C-M-iB\). The same paper nevertheless finds HX and Block Low-Rank preconditioners to be the most promising iterative approaches among those tested, while also showing that BLR can be substantially faster on a large-domain benchmark [2507.13066].

At a more abstract level, operator-preconditioning theory explains why HX often tolerates inexact auxiliary solves. A Petrov–Galerkin extension of operator preconditioning considers a perturbed main operator \(A_{h,\nu}\) and a perturbed preconditioner \(P_{h,\mu}\), with the key estimate
\[
\kappa_S(\mathbf P_\mu \mathbf A_\nu)
\le
K_\star
\left(\frac{1+\mu}{1-\mu}\right)
\left(\frac{1+\nu}{1-\nu}\right).
\]
The main implication is that discretization accuracy and preconditioning accuracy can be decoupled: one may require \(\nu\) small for PDE accuracy while allowing \(\mu=O(1)\) for a cheap or compressed auxiliary preconditioner. This is not an HX construction theorem, but it is directly relevant to practical HX implementations that replace exact nodal inverses by AMG cycles, incomplete factorizations, or other approximations [2011.05028].

Taken together, these developments portray Hiptmair--Xu preconditioners as a stable auxiliary-space paradigm whose essential ingredients are exact-sequence structure, a regular decomposition, and nodal \(H^1\)-type auxiliary solvers. The same pattern appears in positive Maxwell problems, semidefinite curl-curl systems, saddle-point reductions, Schur-complement substructuring, virtual elements, surface complexes, and mixed-dimensional geometries. What changes from one setting to another is not the core mechanism, but the decomposition theory and transfer operators needed to make that mechanism stable and robust [2306.13217][2509.07179].

Source: https://www.emergentmind.com/topics/hiptmair-xu-preconditioners