---
title: Local Taylor-based Polynomial Quasi-Trefftz Space
url: https://www.emergentmind.com/topics/local-taylor-based-polynomial-quasi-trefftz-space
type: topic
---

# Local Taylor-based Polynomial Quasi-Trefftz Space

Local Taylor-based polynomial quasi-Trefftz spaces are local, equation-dependent polynomial spaces whose elements are not exact local solutions of a partial differential equation, but local approximate solutions in the precise sense that a truncated Taylor expansion of the residual vanishes at a chosen point. In the scalar setting, the modern formulation identifies such a space as the kernel of the **quasi-Trefftz operator**, i.e. the composition of Taylor truncation with the differential operator. For scalar linear operators with smooth coefficients, this kernel-based viewpoint yields spaces that are substantially smaller than full polynomial spaces while retaining the same local approximation orders, together with fully explicit basis-construction procedures valid in arbitrary dimension and for operators of arbitrary order [2505.18480].

## 1. Definition and conceptual setting

Trefftz methods use local discrete spaces made of exact local solutions of the governing PDE. Local Taylor-based polynomial quasi-Trefftz spaces relax exactness: a polynomial is retained if, after applying the differential operator, its Taylor expansion up to the maximal admissible order vanishes at a fixed point, typically an element center or barycenter. For a scalar differential operator $\mathcal{L}$ of order $\gamma$, polynomial degree $p\ge \gamma$, and expansion point $\mathbf{x}_0$, the scalar local space is
\[
\mathbb{QT}_p
=
\left\{
\Pi \in \mathbb{P}_p \mid T_{p-\gamma}[\mathcal{L}\Pi]=0
\right\},
\]
equivalently $\mathbb{QT}_p=\ker(\mathcal{D}_p)$ for
\[
\mathcal{D}_p:\mathbb{P}_p\to \mathbb{P}_{p-\gamma},
\qquad
\Pi\mapsto T_{p-\gamma}[\mathcal{L}\Pi].
\]
The truncation order $p-\gamma$ is the natural one because an operator of order $\gamma$ applied to a degree-$p$ polynomial produces a polynomial of degree at most $p-\gamma$ [2505.18480].

A practically important inhomogeneous variant appears for general linear PDEs with smooth coefficients and source term. If $\mathcal{M}$ has order $m$ on an element $E$ with center $x^E$, the degree-$p$ local quasi-Trefftz space with source $f$ is
\[
{}_f^p(E)
:=
\left\{
v\in \mathbb{P}^p(E)
:
D^\beta(\mathcal{M}v-f)(x^E)=0
\quad
\forall\, |\beta|\le p-m
\right\}.
\]
This states that the PDE is satisfied to Taylor order $p-m$ at the chosen point. In this formulation, if $p<m$, every polynomial of degree at most $p$ is trivially quasi-Trefftz, so the construction becomes nontrivial only for $p\ge m$ [2408.00392].

The modifier “local” is essential. These spaces are defined pointwise through Taylor data and are therefore elementwise objects. They become usable in numerical discretizations only after a global coupling mechanism—most often discontinuous Galerkin, embedded Trefftz, or LDG coupling—is imposed across mesh interfaces [2412.00806].

## 2. Algebraic structure and explicit basis construction

The central structural fact is that the quasi-Trefftz operator is most naturally analyzed on the graded decomposition of polynomial spaces,
\[
\mathbb{P}_p=\bigoplus_{l=0}^{p}\widetilde{\mathbb{P}}_l,
\]
where $\widetilde{\mathbb{P}}_l$ denotes homogeneous polynomials of degree $l$. In this representation, the quasi-Trefftz operator has a block-triangular structure. Its leading block is determined by the principal part of the differential operator at the expansion point,
\[
\mathcal{L}_*
:=
\sum_{|\mathbf{j}|=\gamma} c_{\mathbf{j}}(\mathbf{x}_0)\,\partial^{\mathbf{j}},
\]
while all lower-order terms and coefficient variations contribute only to strictly lower blocks. This graded viewpoint is the basis of the systematic scalar theory and is the reason an explicit forward-substitution construction is possible [2505.18480].

For scalar linear equations, the crucial solvability input is that, for every $n\ge 0$, the map
\[
\mathcal{L}_*:\widetilde{\mathbb{P}}_{n+\gamma}\to \widetilde{\mathbb{P}}_{n}
\]
is surjective for every nontrivial linear constant-coefficient homogeneous differential operator. The construction then proceeds degree by degree. The homogeneous components below degree $\gamma$ are free. For each higher degree, one solves
\[
\mathcal{L}_*(\Pi_{n+\gamma})
=
-\big[\text{known lower-degree contributions}\big],
\]
so that the residual vanishes at the required Taylor order. The outcome is a fully explicit local basis with minimal computational cost for scalar operators with smooth coefficients [2505.18480].

Earlier constructive formulations already exhibited the same recursive pattern, although without the later fully general kernel formalism. In the general smooth-coefficient setting for a PDE of order $m$, one may impose a non-degeneracy condition such as
\[
\alpha_{e_1\cdot m}(x^E)\neq 0,
\]
choose as free “Cauchy data” the coefficients with first index $\alpha_1<m$, and recover the remaining coefficients recursively from the Taylor constraints. Each admissible choice of Cauchy data determines a unique quasi-Trefftz polynomial, and a basis is obtained by systematically selecting a basis for those data [2408.00392]. A closely related recursive construction appears for diffusion-advection-reaction operators with piecewise-smooth coefficients, where the polynomial coefficients are computed from local Taylor expansions of the PDE coefficients by a simple completion algorithm [2312.09919].

## 3. Dimension, economy, and approximation properties

The main appeal of local Taylor-based polynomial quasi-Trefftz spaces is that they preserve the approximation order of full polynomial spaces while drastically reducing local dimension. In the scalar kernel formulation,
\[
\dim(\mathbb{QT}_p)=\dim(\mathbb{P}_p)-\dim(\mathbb{P}_{p-\gamma}),
\]
so the dimension is reduced by exactly the number of independent Taylor residual constraints [2505.18480].

For smooth exact solutions, the Taylor polynomial itself belongs to the local quasi-Trefftz space. In the scalar homogeneous case, if $u$ solves $\mathcal{L}u=0$, then $T_p[u]\in \mathbb{QT}_p$, and the best local quasi-Trefftz approximation reproduces Taylor accuracy:
\[
|u(x)-u_a(x)|\le C|x-x_0|^{p+1},
\qquad
|\nabla u(x)-u_a(x)|\le C|x-x_0|^{p}.
\]
In the inhomogeneous smooth-coefficient formulation, the degree-$p$ Taylor polynomial of the exact solution at $x^E$ belongs to ${}_f^p(E)$, and the local best-approximation estimate matches that of the full polynomial space in $C^q$ seminorms [2505.18480; 2408.00392].

The following representative formulas illustrate the resulting economies.

| Setting | Local defining condition | Dimension result |
|---|---|---|
| Scalar operator of order $\gamma$ | $\mathbb{QT}_p=\{\Pi\in\mathbb{P}_p:T_{p-\gamma}[\mathcal{L}\Pi]=0\}$ | $\dim(\mathbb{P}_p)-\dim(\mathbb{P}_{p-\gamma})$ [2505.18480] |
| Second-order elliptic scalar PDE | ${}_0^p(E)=\{v\in\mathbb{P}^p(E):D^\beta(\mathcal{M}v)(x^E)=0,\ |\beta|\le p-2\}$ | $2p+1$ in 2D, $(p+1)^2$ in 3D [2408.00392] |
| Space-time heat equation | $\mathbb{Q}^p(K)=\{v\in\mathbb{P}^p(K):D^{(\alpha,i_t)}\mathcal{H}v(x_K,t_K)=0,\ |\alpha|+i_t\le p-2\}$ | $\binom{p+d}{d}+\binom{p-1+d}{d}$ [2411.14819] |
| Second-order Maxwell system | $T_{p-2}[\nabla\times\nabla\times\mathbf{\Pi}-\epsilon\mathbf{\Pi}]=0$ and $T_{p-1}[\nabla\cdot(\epsilon\mathbf{\Pi})]=0$ | $3p^2+10p-2$ [2509.00193] |

These formulas make precise the statement that quasi-Trefftz spaces are “much smaller” than full polynomial spaces. In the heat equation case, for example, the local dimension scales as $p^d$, whereas full space-time polynomials scale as $p^{d+1}$ [2411.14819]. In elliptic diffusion-advection-reaction problems, the same $h^p$ convergence as full polynomial DG can be obtained with these reduced spaces [2408.00392].

## 4. Role in discontinuous Galerkin and Trefftz-like discretizations

Local Taylor-based polynomial quasi-Trefftz spaces are mainly used as building blocks inside discontinuous formulations. For elliptic diffusion-advection-reaction problems with smooth coefficients, polynomial quasi-Trefftz DG restricts a classical SIPG-plus-upwind bilinear form to the quasi-Trefftz trial and test spaces. In the non-homogeneous case, a local quasi-Trefftz particular solution is constructed elementwise and the homogeneous correction is solved for globally. The resulting scheme is stable, achieves high-order convergence, and its stated practical advantage over standard DG is higher accuracy for comparable numbers of degrees of freedom [2408.00392].

For piecewise-smooth coefficients, a closely related quasi-Trefftz DG method for the homogeneous diffusion-advection-reaction equation proves that the local quasi-Trefftz space has smaller dimension than the full polynomial space of the same degree while yielding the same optimal convergence rates. The method is shown to be well-posed, consistent, stable, and high-order convergent, with numerical experiments in two dimensions exhibiting the expected approximation and convergence behavior [2312.09919].

In space-time wave discretization, the quasi-Trefftz condition is imposed on polynomials whose Taylor coefficients satisfy the variable-coefficient wave operator up to order $p-2$ at the element center:
\[
{}^p\mathbb{U}(K)
=
\left\{
f\in\mathbb{P}^p(K):
D^i(\square_{\rho,G}f)(\mathbf{x}_K,t_K)=0,\ \forall |i|\le p-2
\right\}.
\]
The basis is generated recursively from Cauchy data $f(\cdot,t_K)$ and $\partial_t f(\cdot,t_K)$, and the DG scheme built from the associated first-order variables is proved stable and high-order convergent for piecewise-smooth coefficients [2011.04617].

A broader abstract interpretation is supplied by the unified framework for Trefftz-like methods. There the discrete space is decomposed as
\[
V_h=\mathbb{U}_h\oplus \mathbb{T}_h,
\]
with the quasi-Trefftz space appearing as the globally coupled component satisfying local Taylor conditions. A notable analytical point is that the resulting error estimate is formulated in terms of the best approximation over the full discrete space $V_h$, not only over the quasi-Trefftz subspace. This framework covers quasi-Trefftz methods for second-order scalar elliptic PDEs and a scalar reaction-advection problem and clarifies their relation to embedded Trefftz DG and other local/global decompositions [2412.00806].

For parabolic problems, an inf-sup stable space-time LDG method for the heat equation admits quasi-Trefftz polynomial spaces as local trial spaces. The weaker inf-sup theory applies because quasi-Trefftz spaces generally do not satisfy $\partial_t V(K)\subseteq V(K)$, but the method still yields optimal algebraic rates in the natural energy-type norm with substantially reduced local dimension [2411.14819].

## 5. Extensions to wave problems, first-order systems, and Maxwell equations

The polynomial quasi-Trefftz idea first matured in wave-oriented contexts where exact polynomial Trefftz spaces are unavailable or too restrictive. For the 3D convected Helmholtz equation, a polynomial quasi-Trefftz function is defined by requiring the PDE residual to satisfy
\[
cH(x)=O(|x-x_C|^q),
\]
which is equivalent to vanishing of the Taylor coefficients of the residual up to degree $q-1$. The corresponding coefficient system is linear and triangular by degree, the dimension of the space is $(q+2)^2$, and the basis is generated by initializing one free monomial coefficient for multi-indices with $j_1\in\{0,1\}$ and solving recursively for the remainder. The same work reports that the polynomial basis avoids the ill-conditioning typical of wave-like quasi-Trefftz bases: in one test case, at order $n=8$, the condition number was $1.14\times 10^{17}$ for the plane-wave basis and $7\times 10^3$ for the polynomial quasi-Trefftz basis [2201.12993].

For variable-coefficient wave propagation, generalized plane waves (GPWs) provide a nonlinear quasi-Trefftz alternative. A GPW basis function has the form $e^{P(\mathbf{x})}$ and is required to satisfy a Taylor-truncated residual condition of the form
\[
T_{p-\gamma}\{\mathcal{L}[e^P]\}=0.
\]
The polynomial quasi-Trefftz space is the linear analog of this construction: it uses the ansatz $\varphi\in P_p$ and the linear condition $T_{p-\gamma}(\mathcal{L}\varphi)=0$. The two constructions are therefore closely related, but polynomial quasi-Trefftz spaces retain a finite-dimensional linear-algebraic structure [2508.09435].

The extension from scalar equations to systems is nontrivial. For the second-order time-harmonic Maxwell equation with variable coefficient $\epsilon(\mathbf{x})$, the local Taylor-based polynomial quasi-Trefftz space is defined by two simultaneous conditions:
\[
T_{p-2}\left[\nabla\times\nabla\times\mathbf{\Pi}-\epsilon\mathbf{\Pi}\right]=0,
\qquad
T_{p-1}\left[\nabla\cdot(\epsilon\mathbf{\Pi})\right]=0.
\]
The second, quasi-divergence condition is essential because of the structure of the curl-curl operator. A Helmholtz decomposition of homogeneous polynomial vector fields yields both an explicit dimension formula,
\[
\dim \mathbb{Q}\mathbb{T}_p = 3p^2+10p-2,
\]
and a recursive construction procedure for the local basis. This is identified as the first study of local Taylor-based polynomial quasi-Trefftz spaces for a system of PDEs [2509.00193].

A complementary first-order development appears for a first-order formulation of the Helmholtz equation in 2D. There the local space consists of triples $(p,v_x,v_y)\in\mathbb{P}^d\times(\mathbb{P}^{d-1})^2$ satisfying Taylor-truncated first-order residual constraints, and the local dimension is $2d+1$. Two explicit recursive constructions are given: a coupled system algorithm and a decoupled pressure algorithm based on the associated scalar Helmholtz equation [2509.08936].

## 6. Distinctions, limitations, and related constructions

Several distinctions are fundamental. First, quasi-Trefftz does not mean Trefftz. In a classical Trefftz space, the local residual vanishes identically in the element; in a local Taylor-based polynomial quasi-Trefftz space, only a finite Taylor jet of the residual vanishes at a point. The two coincide only in special constant-coefficient situations where exact polynomial solutions exist or where the truncated condition is sufficient to reconstruct an exact local polynomial solution [2505.18480].

Second, the local condition is point-centered rather than globally conforming. By itself, a quasi-Trefftz space does not enforce continuity or interface transmission. This is why quasi-Trefftz spaces are predominantly paired with DG-type couplings, embedded Trefftz constructions, or other weak interface formulations [2412.00806].

Third, the family is generally not nested in the polynomial degree: for diffusion-advection-reaction quasi-Trefftz spaces one has, in general, $\mathcal{Q}^p(T)\not\subseteq \mathcal{Q}^{p+1}(T)$. This distinguishes them from standard polynomial hierarchies and affects basis design, $p$-adaptivity, and conditioning strategies [2312.09919].

A further limitation is exposed by the Helmholtz equation itself. For nonzero wavenumber, no nonvanishing piecewise-polynomial exact Trefftz function exists because the Helmholtz operator contains both second- and zero-order terms. This is one of the original motivations for quasi-Trefftz constructions and for non-polynomial exact Trefftz spaces such as plane waves or generalized harmonic polynomials [1506.04521].

Finally, not every local PDE-adapted polynomial space should be identified with a Taylor-based quasi-Trefftz space. In the embedded Trefftz DG method for the Oseen problem, the analysis relies on an explicitly constructed local complement space
\[
X(T)=T_w(T)\oplus L(T),
\]
where $L(T)$ is polynomial, local, and designed so that the Oseen operator is stably invertible on it. The paper explicitly states that for Oseen this construction is polynomial and geometric, not strictly Taylor-based, even though it is closely related to earlier Taylor-based quasi-Trefftz ideas in the constant-coefficient Stokes case [2606.13229]. This distinction delineates the quasi-Trefftz concept from the broader class of local PDE-tailored polynomial decompositions.

In current usage, the term therefore denotes a precise algebraic and local notion: a polynomial kernel space defined by Taylor-truncated residual annihilation. Its significance lies in combining PDE adaptation, explicit local construction, reduced dimension, and high-order approximation for variable-coefficient problems where exact local Trefftz spaces are unavailable or impractical [2505.18480].

Source: https://www.emergentmind.com/topics/local-taylor-based-polynomial-quasi-trefftz-space