---
title: Discrete Quasi-Trefftz Spaces
url: https://www.emergentmind.com/topics/discrete-quasi-trefftz-spaces
type: topic
---

# Discrete Quasi-Trefftz Spaces

Discrete quasi-Trefftz spaces are finite-dimensional, equation-dependent approximation spaces whose functions satisfy the governing partial differential equation exactly only in special cases and otherwise satisfy it approximately in a controlled local sense, typically through vanishing Taylor coefficients of the residual or through a relaxed weak local constraint. They were developed to extend classical Trefftz discretizations beyond the regime of linear homogeneous PDEs with piecewise-constant coefficients and explicitly available local solutions, while retaining the central Trefftz objective: replacing full polynomial spaces by much smaller problem-adapted spaces without sacrificing approximation order [1506.04521][2412.00806][2505.18480].

## 1. Definition and relation to classical Trefftz spaces

In the classical Helmholtz setting, a discrete Trefftz space is a finite-dimensional subspace of functions that satisfy the PDE pointwise on each mesh element; for a mesh \(\mathcal T_h\), the survey literature writes
\[
T(\mathcal T_h):=\{v\in H^1_{\rm pw}(\Omega):(-\Delta-k^2)v=0 \text{ pointwise in each }K\in\mathcal T_h\},
\]
and a discrete Trefftz space is any finite-dimensional subspace assembled from local spaces \(V_{p_K}(K)\subset T(\mathcal T_h)\) [1506.04521]. In that same survey, quasi-Trefftz is described as a typically larger space whose elements satisfy the PDE only approximately, for example up to higher-order terms in \(h_K\) [1506.04521].

Later formulations make that approximation notion precise. For a second-order scalar elliptic operator \(L v=-\nabla\!\cdot(\alpha\nabla v)+\beta\cdot\nabla v+\gamma v\), the unified framework defines on each element \(K\) the local quasi-Trefftz space
\[
V_h^{QT}(K):=\{v\in P^p(K):D^i(Lv)(x_K)=0\ \forall\,|i|\le p-2\}=\ker A_K,
\]
where \(x_K\) is a chosen center and \(A_K\) collects the scaled residual jets at that point [2412.00806]. The same structure appears in the general Taylor-based theory, where for an operator \(\mathcal L\) of order \(\gamma\) one defines
\[
\mathcal D_p(\Pi)=T_{p-\gamma}[\mathcal L\Pi],\qquad \mathbb{QT}_p(x_0):=\ker \mathcal D_p,
\]
so the local quasi-Trefftz space is the kernel of “Taylor truncation composed with the differential operator” [2505.18480].

This formulation places exact Trefftz and quasi-Trefftz spaces in a strict inclusion chain whenever exact polynomial solutions exist. For the elliptic framework,
\[
T_p(K)=\{v\in P^p(K):Lv=0\text{ exactly in }K\},\qquad
T_p(K)\subsetneq V_h^{QT}(K)\subset P^p(K),
\]
and \(\dim T_p(K)<\dim V_h^{QT}(K)\) in general [2412.00806]. A common misconception is that quasi-Trefftz spaces are merely ad hoc enrichments; the cited constructions instead characterize them as structured kernels of local residual operators.

## 2. Local algebraic structure, basis construction, and dimensions

The local algebraic core is a constrained polynomial space. In the general scalar theory, polynomial spaces are decomposed by total degree,
\[
\mathbb P_p=\bigoplus_{\ell=0}^p \widetilde{\mathbb P}_\ell,
\]
and the quasi-Trefftz operator \(\mathcal D_p:\mathbb P_p\to\mathbb P_{p-\gamma}\) is shown, under a minimal non-degeneracy assumption on the principal part of \(\mathcal L\), to be surjective [2505.18480]. This yields the dimension formula
\[
\dim \mathbb{QT}_p=\binom{p+d}{d}-\binom{p-\gamma+d}{d},
\]
which is the canonical dimension count for Taylor-based polynomial quasi-Trefftz spaces [2505.18480].

At the implementation level, several equivalent constructions recur across the literature. One may assemble the matrix of residual moments or truncated jets in a monomial basis and compute a nullspace basis by Gaussian elimination, SVD, or rank-revealing QR [2412.00806][2408.00392][2505.18480]. Alternatively, one may prescribe “Cauchy data” and recover the remaining coefficients by explicit recurrences ordered by homogeneous degree; this is the preferred construction in the space-time wave equation, the diffusion-advection-reaction equation, and several first-order or vector generalizations [2011.04617][2312.09919][2509.08936][2509.00193].

Representative local dimensions appearing in the literature are as follows.

| Setting | Local space | Dimension |
|---|---|---|
| Scalar operator of order \(\gamma\) in \(d\) dimensions | \(\mathbb{QT}_p\) | \(\binom{p+d}{d}-\binom{p-\gamma+d}{d}\) |
| Second-order scalar elliptic, \(d=2\) | \(V_h^{QT}(K)\) | \(2p+1\) |
| Space-time wave equation | \({}^p(K)\) | \(\binom{p+n}{n}+\binom{p-1+n}{n}\) |
| Second-order time-harmonic Maxwell, \(d=3\) | \(\mathbb Q\mathbb T_p\) | \(3p^2+10p-2\) |
| Embedded Trefftz on quadrilaterals | \(T(K)\) | \(2p+2\) |

These counts reflect the same structural principle: quasi-Trefftz spaces eliminate precisely the polynomial directions that generate low-order residuals. In two-dimensional second-order scalar problems, this produces the particularly important linear-in-\(p\) count \(2p+1\), in contrast with the quadratic growth of the full polynomial space [2412.00806][2312.09919].

## 3. Global discretization patterns

Most discrete quasi-Trefftz methods are formulated in discontinuous Galerkin form. In the unified framework one begins with the broken polynomial space
\[
V_h=\bigoplus_{K\in\mathcal T_h}P^p(K),
\]
then decomposes it as \(V_h=U_h\oplus T_h\), where \(T_h\) is the discontinuous global space whose elementwise restrictions lie in the local quasi-Trefftz spaces [2412.00806]. The method enforces two conditions simultaneously: local residual matching in the \(U_h\)-direction and the DG residual against all test functions in \(T_h\). The resulting \(2\times2\) block system simplifies in most quasi-Trefftz methods because the upper-right block vanishes, so one first solves local particular problems and then a reduced global DG problem on \(T_h\) [2412.00806].

For scalar elliptic problems, the global bilinear form is typically an interior-penalty formulation with upwind treatment of advection. The diffusion-advection-reaction quasi-Trefftz DG method uses broken trial and test spaces with a DG norm containing broken \(H^1\)-seminorms, jump penalization, and upwind terms, and proves consistency, discrete coercivity, and boundedness [2312.09919]. A closely related general elliptic formulation employs SIPG/upwind couplings for local spaces \(Q_k(T)\subset P_k(T)\) defined by vanishing residual jets at element centers [2408.00392].

The space-time wave equation provides a different global pattern. There, local quasi-Trefftz polynomial spaces on space-time elements are assembled into a DG space for the first-order acoustic variables \((w,\tau)\), and the choice of numerical fluxes makes the coupling causal, so no global space-time solve is required and one may march in time via element patches, even in parallel [2011.04617].

A distinct line is the embedded Trefftz DG method. On each anisotropic quadrilateral element \(K\), one defines a tensor-product space \(V(K)\), a smaller local test space \(L(K)\), and the relaxed local Trefftz condition
\[
(A_K v_h,q_h)_K=0\qquad\forall\,q_h\in L(K),
\]
with \(A_K(v)=-\varepsilon^2\Delta v+v\) [2606.03845]. The corresponding local embedded Trefftz space \(T(K)\) is not imposed by pointwise residual annihilation but by weak orthogonality to \(L(K)\); the global method then imposes the DG bilinear form only on the product space \(T_h\), again yielding a reduced global system [2606.03845].

Although quasi-Trefftz methods are predominantly discontinuous, a globally conforming variant now exists for Helmholtz. The Trefftz Continuous Galerkin method constructs a space \(V_{(N_e,N_n)}\subset H^1_\kappa(\Omega)\) from edge-based and node-based functions that are exact local Helmholtz solutions on each cell and continuous across the global Cartesian cut mesh, leading to a conforming Trefftz discretization rather than a DG coupling [2512.02145].

## 4. Approximation theory, stability, and convergence

The defining analytical feature of quasi-Trefftz spaces is that they preserve the approximation order of the surrounding full polynomial space. In the general scalar Taylor-based theory, if \(u\) is an exact local solution of \(\mathcal L u=0\) with \(u\in \mathcal C^{p+\gamma}\), then its Taylor polynomial \(T_pu\) belongs to \(\mathbb{QT}_p\), and the best-approximation error is \(O(h^{p+1})\) in \(L^\infty\) and \(O(h^p)\) in the gradient [2505.18480]. For local spaces \(Q_k(T)\) associated with a linear differential operator of order \(m\), one has
\[
\inf_{q\in Q_k(T)}\|u-q\|_{H^s(T)}\le C(d,m,k)\,h_T^{k+1-s}\|u\|_{H^{k+1}(T)},
\]
so the same \(h\)-rates as full \(P_k\)-approximation are recovered [2408.00392].

The unified DG framework turns this local property into global quasi-optimality. Discrete stability and Céa-type arguments yield
\[
\|u-u_h\|_{V_h}\lesssim \inf_{v_h\in V_h}\|u-v_h\|_{V_h}\lesssim \inf_{w_h\in T_h}\|u-w_h\|_{V_h},
\]
and for coercive second-order problems one obtains
\[
\|u-u_h\|_{DG}\lesssim h^p\|u\|_{H^{p+1}(\Omega)},\qquad
\|u-u_h\|_{L^2(\Omega)}\lesssim h^{p+1}\|u\|_{H^{p+1}(\Omega)}
\]
under the stated regularity assumptions [2412.00806]. The diffusion-advection-reaction analysis reaches the same conclusion: the quasi-Trefftz DG method is well-posed, consistent, stable, and high-order convergent, with local spaces smaller than full polynomial spaces of the same degree [2312.09919].

For the space-time wave equation with piecewise-smooth coefficients, the exact solution’s Taylor polynomial of degree \(p\) belongs to the local quasi-Trefftz space \({}^p(K)\), and the assembled DG method satisfies a quasi-optimal bound in the mesh-dependent energy norm with asymptotic order \(O(h^{p+1/2})\), together with \(h^{p+1}\)-type behavior at final time under duality arguments [2011.04617]. The first-order Helmholtz formulation has the analogous local statement: if \((p,v)\) is a sufficiently smooth exact solution, then \((\mathcal T_d p,\mathcal T_{d-1}v)\in \mathbb{QT}_d^F\), and the local approximation error scales like \(|x-x_0|^{d+1}\) for the fields and \(|x-x_0|^d\) for their gradients [2509.08936].

A misconception often attached to quasi-Trefftz spaces is that approximate local PDE satisfaction necessarily degrades stability. The available analyses point in the opposite direction: local residual annihilation is designed so that the global constants remain those of the full DG method, while the main effect is a reduction of the global dimension [2412.00806]. This suggests that quasi-Trefftz spaces should be understood less as weakened Trefftz spaces than as dimension-reduced polynomial spaces with PDE-adapted null directions.

## 5. Wave-propagation specializations

Wave problems have supplied the main impetus for quasi-Trefftz developments because the exact Trefftz philosophy is especially attractive there and especially fragile under variable coefficients. In homogeneous Helmholtz media, plane waves \(e^{ik d\cdot x}\), generalized harmonic polynomials, evanescent waves, fundamental solutions, and multipoles are exact local Trefftz functions [1506.04521]. As soon as the coefficients vary, “no closed-form family of exact local solutions is available,” which motivates replacing exact local solutions by high-order approximate ones [2508.09435].

One major specialization is the generalized plane wave construction. On each cell \(T\), one chooses a center \(x_T\) and basis functions of the form
\[
G_j(x)=\exp(i\phi_j(x)),\qquad
\phi_j(x)= i\kappa\,d_j\cdot(x-x_T)+\sum_{2\le |\alpha|\le p}\beta_\alpha^{(j)}(x-x_T)^\alpha,
\]
and determines the higher-order coefficients by imposing that the first \(q=p-2\) terms of the Taylor expansion of the PDE residual vanish [2508.09435]. The same paper explicitly relates GPW bases to polynomial quasi-Trefftz bases, making the phase-based and polynomial constructions part of a single framework.

Another specialization is the Trefftz Continuous Galerkin method for the two-dimensional Helmholtz equation. On each Cartesian rectangle \(K\), one solves a homogeneous Helmholtz–Dirichlet problem with nonzero trace on one chosen edge and zero trace on the others, obtaining single-edge modes
\[
\varphi_n(x,y)=
\sin(\kappa\nu_n x)\,
\frac{\sin(\kappa\sqrt{1-\nu_n^2}\,y)}
{\sin(\kappa\sqrt{1-\nu_n^2}\,h_2)},\qquad
\nu_n=\frac{n\pi}{\kappa h_1},
\]
which are propagative for \(\nu_n<1\) and evanescent for \(\nu_n>1\) [2512.02145]. These functions can also be written as linear combinations of four complex-direction plane waves, are transplanted to adjacent cells sharing an edge, and yield globally continuous edge and node spaces in \(H^1_\kappa(\Omega)\) [2512.02145]. The resulting method proves wavenumber-explicit best-approximation bounds, bounded coefficients, and spectral accuracy for analytic Helmholtz solutions, with exponential decay of the approximation error both at fixed frequency and along suitable high-frequency sequences [2512.02145].

Quasi-Trefftz ideas have also been extended to first-order and vector systems. For a first-order formulation of the Helmholtz equation, the local space
\[
\mathbb{QT}^F_d\subset \mathbb P^d\times(\mathbb P^{d-1})^2
\]
is characterized equivalently by the scalar condition \(\mathcal T_{d-2}(\mathcal L(p))=0\) together with \(v=\frac1{\omega\rho}\nabla p\), and has dimension \(2d+1\) [2509.08936]. For the second-order time-harmonic Maxwell equation with variable coefficients, the vector-valued space \(\mathbb Q\mathbb T_p\) is defined through coupled curl-curl and divergence truncation constraints, its dimension is \(3p^2+10p-2\), and its constructive theory relies on a Helmholtz decomposition of homogeneous vector polynomial fields into solenoidal, irrotational, and harmonic components [2509.00193].

## 6. Computational aspects, numerical behavior, and active directions

Computationally, discrete quasi-Trefftz spaces trade larger local basis-construction work for smaller global systems. For scalar polynomial spaces, the construction is entirely local and can be done from monomial matrices, nullspace extraction, or block-triangular recursions [2408.00392][2505.18480]. For the first-order Helmholtz formulation, explicit recurrences cost \(\mathcal O(d^4)\) operations per function and \(\mathcal O(d^5)\) to build the full \(2d+1\) basis, whereas SVD-based algebraic approaches cost \(\mathcal O(d^6)\) flops [2509.08936]. In the space-time wave equation, the local basis is built by a single pass of dimension-\(O(p^n)\) linear recurrences and its conditioning remains moderate in the reported experiments [2011.04617].

Assembly may also benefit from the PDE-adapted basis. The Helmholtz survey emphasizes that Trefftz formulations often involve only face or boundary integrals and that plane-wave integrals can be evaluated in closed form on straight segments and polygonal faces [1506.04521]. The TCG Helmholtz construction sharpens this point: each local basis function expands into four complex-direction plane waves on each adjacent cell, and since the required cell and boundary integrals have closed-form expressions, all entries of the stiffness, mass, and boundary-term matrices can be computed exactly without quadrature for polygonal cells [2512.02145].

Conditioning remains a central issue. The Helmholtz survey records that plane-wave bases become nearly linearly dependent when \(kh_K\) is small or directions cluster, and reviews remedies such as limiting the number of plane waves per element, local orthonormalization by Gram–Schmidt or SVD, oversampling or least-squares formulations, hybrid bases, and directional adaptivity [1506.04521]. Quasi-Trefftz spaces inherit part of this challenge, but several constructions report favorable coefficient control. In the TCG Helmholtz method, edge and node expansion coefficients remain uniformly bounded in \(\ell^2\) as \(N_e\to\infty\) at fixed \(N_n\), grow at most polynomially in \(N_n\) at fixed frequency, and remain modest numerically in the high-frequency regime when \(N_e,N_n\sim \kappa h\) [2512.02145].

The numerical record across applications is consistent. The diffusion-advection-reaction method shows \(H^1\)-errors of order \(O(h^p)\), \(L^2\)-errors of order \(O(h^{p+1})\), and exponential \(p\)-convergence for smooth tests while using a local dimension \(2p+1\) in two dimensions rather than \((p+1)(p+2)/2\) [2312.09919]. The embedded Trefftz DG method on anisotropic quadrilateral meshes reproduces the convergence rates of the full DG method while reducing the size of the global system, with reported tests showing roughly half as many global unknowns or up to \(50\%\) fewer global unknowns for comparable accuracy [2606.03845]. The TCG Helmholtz method reports exponential convergence for analytic targets and robust high-frequency behavior, including \(8\)–\(10\) digits on complex polygonal domains with only a few DOFs per wavelength in least-squares Petrov-Galerkin solves [2512.02145].

Current research directions follow directly from these patterns. The 2015 Helmholtz survey identifies \(k\)-robust discretizations, variable media, complex geometries, adaptive directional refinement, and improved preconditioning as central open problems for Trefftz-type methods [1506.04521]. Subsequent quasi-Trefftz work has already moved into general scalar operators, space-time hyperbolic problems, first-order systems, Maxwell equations, generalized plane waves, anisotropic meshes, and globally conforming Helmholtz constructions [2505.18480][2509.00193][2606.03845][2512.02145]. A plausible implication is that “discrete quasi-Trefftz spaces” now designate not a single method class, but a unifying local approximation principle: encode the PDE directly into the finite-dimensional space by annihilating low-order residual structure, then exploit that space in DG, embedded, least-squares, or conforming global couplings.

Source: https://www.emergentmind.com/topics/discrete-quasi-trefftz-spaces