---
title: 'System Interpolation: Multi-Domain Approaches'
url: https://www.emergentmind.com/topics/system-interpolation
type: topic
---

# System Interpolation: Multi-Domain Approaches

System interpolation denotes a family of constructions in which interpolation is performed at the level of a system rather than merely at the level of a scalar function. In the cited literature, the term ranges from nonuniform sampled-data reconstruction with Chebyshev nodes and word-serial hardware architectures, to frequency-domain interpolation of transfer functions for model reduction and identification, to direct completion of missing LPV trajectories from data, to interpolation inequalities used to close nonlinear PDE estimates, to interpolation-derived integrable systems, to proof-theoretic effective interpolation, and to parameter-space interpolation between language-model checkpoints [1001.1185], [2506.03410], [2508.11332], [2409.18808], [1706.05155], [2601.02821], [2601.21414]. The unifying feature is that interpolation acts on a structured object—signals, trajectories, transfer functions, solution classes, proof predicates, or checkpoint parameters—and thereby changes the behavior of an entire system.

## 1. Conceptual range of the term

The literature uses “system interpolation” in several technically distinct but structurally related senses.

| Setting | Interpolated object | System-level role |
|---|---|---|
| Sampled-data and sensing | Signals, soundfields, virtual views | Reconstruction, hardware realization, remote sensing |
| Dynamical systems | Transfer functions, tangential data, spectral zeros | Model reduction, identification, passivity preservation |
| Behavioral/data-driven systems | Finite trajectories under scheduling | Missing-data completion without explicit model identification |
| PDE and weak-solution theory | Intermediate norms or solution regimes | Closing estimates or interpolating solution behavior |
| Integrable and logical systems | Determinant data, proof obligations | Generating dynamics or extracting separators/proofs |
| Parameter-space model control | Checkpoint parameters | Continuous control of reasoning intensity |

In sampled-data settings, interpolation typically reconstructs a signal or field from sparse or nonuniform observations. In control and model reduction, interpolation is usually frequency-domain and transfer-function based, with exact matching at selected frequencies or directions. In behavioral system theory, interpolation can mean completing a partially specified trajectory so that it belongs to the admissible behavior. In PDE analysis, interpolation often refers to inequalities between summability and smoothness spaces, or to a continuous parameter between conservative and dissipative solution concepts. In integrable systems and logic, interpolation is not primarily approximation; it is the device from which the system itself is derived or the proof-theoretic property is extracted [2206.12502], [2601.02821].

This diversity suggests that “system interpolation” is best understood as an umbrella term for constructions where interpolation is constrained by, or acts directly upon, an ambient system structure rather than an isolated function value table.

## 2. Sampled-data, acoustic, and viewpoint systems

A hardware-oriented meaning of system interpolation appears in the Chebyshev sampled-data setting of “Some Architectures for Chebyshev Interpolation” [1001.1185]. The signal is sampled over \(x\in[-1,1]\), but not on an equispaced grid; instead, the sample locations are the roots of \(T_{N+1}\),
\[
x_k = \cos\left(\frac{2k-1}{2n}\pi\right), \quad k=1,2,\dots,n.
\]
The reconstruction polynomial is written in a Chebyshev basis,
\[
P_N(x)=\sum_{i=0}^{N} c_i \bar{T}_i(x),
\]
with coefficients obtained by a DCT-like transform. The paper contrasts transform-domain and time-domain realizations, introduces Chebyshev Type Interpolation Functions \(\phi_i(x)\), and then proposes a word-serial systolic-array architecture for an 8-sample window. Its first stage computes \(\{c_i\}\), its second stage evaluates Chebyshev polynomials recursively using a reused multiplier-adder pair, and its output stage performs an FIR-like accumulation. The dependence-graph scheduling is summarized by
\[
T=[1\;\;0], \qquad s=[1\;\;1]^T.
\]
For \(N=8\), the reported latency is 16 cycles, with 100% hardware usage efficiency. The same paper also links interpolation to the sampling front end: for the bandlimited test signal \(\sin(4x)+0.5\sin(8x)\), error \(<1.1\%\) required 10 equispaced points but 8 Chebyshev points; for \(e^{-x}\sin(8x)\), error \(<4.1\%\) required 11 equispaced points but 8 Chebyshev points. Under an 8-bit flash/SAR comparison-count model, the resulting Flash-SAR hybrid consumed 1552 au per window versus 2560 au and 2816 au for the equispaced flash-only alternatives, corresponding to approximately 39% and 44% power savings.

In active noise control, interpolation is used to synthesize virtual error microphones inside a region of interest from physical microphones outside it [2309.10605]. The proposed system is a simulated 3D free-field ANC setup with one primary source, two secondary sources, eight monitoring microphones, and two virtual ear microphones. The interpolated field \(\hat p(n,x,y,z)\) is produced by a PINN constrained by the homogeneous wave equation,
\[
\nabla^2 p - \frac{1}{c^2}\frac{\partial^2 p}{\partial t^2}=0,
\]
through a loss
\[
\mathfrak{L}=\mathfrak{L}_{\mathrm{data}}+\mathfrak{L}_{\mathrm{PDE}}.
\]
The network has one hidden layer with 16 neurons, \(\tanh\) activation, \(A=100\) collocation points, and \(5\times 10^5\) training epochs. Relative to a spherical-harmonic baseline, the PINN reduced interpolation error by approximately \(8\,\mathrm{dB}\) for \(r_s\in[0.2,0.4]\) m. When its interpolated virtual signals replaced physical in-ROI error microphones in FxLMS, the steady-state ear-level noise reduction improved by about \(13\,\mathrm{dB}\), and the residual field around the ear regions was approximately \(10\,\mathrm{dB}\) lower.

A related systems use of interpolation appears in live free-viewpoint video streaming [2112.10603]. There, view interpolation is not only an image-synthesis primitive but the mechanism that changes the architecture of the service. A synchronized 12-camera rig captures live video; a CNN-based midpoint interpolator, VINet, is recursively stacked so that \(n\) stages produce \(2^n-1\) views between adjacent cameras. The deployed system uses 4 stages, yielding 15 interpolated views between adjacent cameras and 177 total views across the 12-camera setup. These dense views are organized into multi-view cluster frames, encoded, segmented in time and viewpoint space, and streamed through HLS. The reported server-side average latencies were 12.85 ms for view interpolation, 1.27 ms for adaptive stitching, 5.62 ms for encoding, and 5.45 ms for scheduling, implying throughput of about 77 FPS. The architectural claim is that, once dense views are generated and packaged, server load becomes unrelated to the number of clients.

## 3. Frequency-domain interpolation for linear systems, identification, and passivity

In large-scale MIMO model reduction, system interpolation is explicitly frequency-domain and tangential [2506.03410]. For
\[
G(s)=C(sI-A)^{-1}B+D,
\]
the paper develops a left-tangential barycentric interpolant
\[
R(s)=M(s)^{-1}N(s)
\]
that satisfies
\[
U_i^*R(s_i)=U_i^*G(s_i)
\]
for any admissible weight matrices \(W_k\), while preserving the feedthrough \(D\). The free weights are not chosen by direct minimization of \(\|R-G\|\), but through a weighted \(H_2\) surrogate,
\[
\min_W \|M(R-G)\|_{H_2}^2=\min_W \|N-MG\|_{H_2}^2.
\]
The resulting quadratic problem has an explicit stationary condition
\[
\widetilde W\,\widetilde C\Theta \widetilde C^* = C\Theta \widetilde C^*,
\]
and the minimized surrogate value \(\gamma_r(\widetilde W_r)\) decreases monotonically with added interpolation rank:
\[
\gamma_{r+1}(\widetilde W_{r+1}) \le \gamma_r(\widetilde W_r).
\]
The framework is AAA-like in point addition, but uses low-rank tangential blocks, Gramian-based weight updates, and max-error, grid, or random point selection. On the 270-state ISS example, the max-error variant delivered \(L_\infty\)-error behavior close to balanced truncation and empirically produced stable reduced models on the tested stable systems.

Loewner rational interpolation provides a different realization-based paradigm for learning low-dimensional systems from frequency-response data [1910.00110]. In the SISO setting, the Loewner and shifted Loewner matrices are
\[
L_{ij}=\frac{H(\mu_i)-H(\gamma_j)}{\mu_i-\gamma_j}, \qquad
L^{(s)}_{ij}=\frac{\mu_i H(\mu_i)-\gamma_j H(\gamma_j)}{\mu_i-\gamma_j},
\]
leading to a descriptor interpolant with \(E=-L\) and \(A=-L^{(s)}\). Under relative complex Gaussian noise on the measurements, the paper proves that for each \(s\in\Omega\),
\[
|H_\sigma(s)-H(s)|\le C_s \sigma
\]
with probability at least
\[
1-4e^{-n/2},
\]
provided a small-noise condition involving \(\kappa_2(sE-A)\), \(\|F_E\|_2+\|F_A\|_2\), and \(\|B\|_\infty/\|B\|_2\) holds. The practical message is that interpolation-point selection, left/right partitioning, and realization scaling directly affect robustness.

A related but identification-specific use of barycentric interpolation is the rapid frequency-domain scheme of [2510.00525]. Here the transfer function is unknown and each \(G(j\omega_k)\) must be obtained by an expensive sinusoidal experiment. The identified interpolant is
\[
R(s)=M^{-1}(s)N(s),
\]
with
\[
M(s)=I+\sum_{k=0}^{\ell}\frac{w_k}{s-j\omega_k}, \qquad
N(s)=D+\sum_{k=0}^{\ell}\frac{w_k G(j\omega_k)}{s-j\omega_k},
\]
so that \(R(j\omega_k)=G(j\omega_k)\) when \(w_k\neq 0\). The paper exploits transient as well as steady-state data from each sinusoidal run and replaces the nonconvex output fit by the convex surrogate
\[
\min_{\widehat W}\sum_{k=1}^{m}\frac{1}{n_k}
\left\|
\begin{bmatrix}1 & \widehat W\end{bmatrix}
(\mathcal N u^{(k)}-\mathcal M y^{(k)})
\right\|_2^2.
\]
When \(\hat X_2\) is positive definite, the unconstrained optimum is
\[
\widehat W=-\hat X_0\hat X_2^{-1}.
\]
A Lyapunov-based LMI then enforces stability. On a 270-state ISS structural model, adaptive frequency selection outperformed gridded selection in \(\mathsf H^2\) error.

Interpolation also underlies scalar positive-real synthesis. The modified Riccati method of [1811.11348] solves rational analytic interpolation with positivity and derivative constraints by reusing the covariance extension equation in the form
\[
P=\Gamma(P-Phh'P)\Gamma' + (u+U\sigma+U\Gamma Ph)(u+U\sigma+U\Gamma Ph)'.
\]
The data enter through
\[
D=(W+\tfrac12I)^{-1}(W-\tfrac12I), \qquad
\begin{bmatrix}u & U\end{bmatrix}=
\begin{bmatrix}0 & I_n\end{bmatrix}V^{-1}DV,
\]
and the interpolant is reconstructed from
\[
a=(I-U)(\Gamma Ph+\sigma)-u, \qquad
b=(I+U)(\Gamma Ph+\sigma)+u.
\]
For each spectral-zero polynomial \(\sigma\) and feasible generalized Pick data, the solution is unique; moreover,
\[
\operatorname{rank}P=\deg f.
\]

Passivity-preserving interpolation is treated differently in [2308.03500]. There the transfer function \(Z(s)\) is square and passive, and the key object is the para-Hermitian function
\[
\Phi(s)=Z^{\mathsf T}(-s)+Z(s).
\]
Interpolation points and directions are obtained from a deflating subspace of the corresponding system pencil \(S(s)\). The reduced model satisfies the bitangential conditions
\[
Z(\lambda_j)Wr_j=\hat Z(\lambda_j)Wr_j, \qquad
r_j^{\mathsf T}W^{\mathsf T}Z(-\lambda_j)=
r_j^{\mathsf T}W^{\mathsf T}\hat Z(-\lambda_j).
\]
The paper then parameterizes a shifted family
\[
Z_\xi(s)=Z\!\left(s-\frac{\xi}{2}\right)-\frac{\xi}{2}I_m
\]
and shows that the corresponding reduced model has a stronger passivity certificate,
\[
(\hat X_\xi,\hat{\mathcal M})\succeq \xi\,\mathrm{diag}(\hat X_\xi,I_m),
\]
which yields a normalized port-Hamiltonian realization with passivity-radius lower bound \(\xi/2\).

## 4. Structured nonlinear and parametric system interpolation

For polynomial and quadratic-bilinear systems, a single linear transfer function is inadequate, so interpolation is formulated on multivariate generalized transfer functions. In [1904.11891], polynomial parametric systems
\[
E(p)\dot x=A(p)x+\sum_{\ell=2}^{d}H_\ell(p)x^{\otimes \ell}
+\sum_{\eta=1}^{d}N_\eta(p)\big(u\otimes x^{\otimes \eta}\big)+B(p)u
\]
are associated with
\[
F_L(s_1,p)=C(p)\Phi(s_1,p)B(p),
\]
\[
F_H^{(\ell)}(s_1,\dots,s_{\ell+1},p)
=C(p)\Phi(s_{\ell+1},p)H_\ell(p)
\big(\Phi(s_\ell,p)B(p)\otimes\cdots\otimes \Phi(s_1,p)B(p)\big),
\]
and analogous \(F_N^{(\eta)}\). Projection spaces are then constructed so that the reduced model satisfies two-sided tangential Hermite interpolation, not only in the frequency variables \(s_j\) but also in the parameter \(p\). A Loewner-inspired realization is built from
\[
L=-W^TEV, \qquad L_s=-W^TAV,
\]
followed by SVD compression and effective projections \(V_{\mathrm{eff}},W_{\mathrm{eff}}\). The numerical examples show that directly reducing the original cubic Chafee–Infante system is more accurate than first rewriting it in quadratic-bilinear form, with up to 3 orders of magnitude better accuracy; a parametric Chafee–Infante ROM of order 5 reproduces dynamics well across \(p\in[0.25,2]\); and CUR approximation can reduce online cost for nonlinear evaluation.

A structurally parallel program is developed for quadratic-bilinear systems with internal differential structure [2304.14292]. The paper uses generic operators
\[
K(s),\quad B(s),\quad C(s),\quad N(s),\quad H(s_1,s_2)
\]
to cover unstructured first-order QB systems, structured second-order systems, and time-delay QB systems. It defines structured symmetric subsystem transfer functions and structured generalized transfer functions, for example
\[
G_{\mathrm{sym},1}(s_1)=C(s_1)K(s_1)^{-1}B(s_1),
\]
and derives exact interpolation conditions by requiring the projection bases to contain recursively defined subsystem responses such as \(K(\sigma)^{-1}B(\sigma)\), bilinear images \(K(\cdot)^{-1}N(\cdot)(I_m\otimes V)\), and quadratic images \(K(\cdot)^{-1}H(\cdot,\cdot)(V\otimes V)\). Structure preservation is built in:
\[
K(s)=\sum_{j=1}^{n_K} h_j(s)K_j
\quad \Rightarrow \quad
\widehat K(s)=\sum_{j=1}^{n_K} h_j(s)\widehat K_j.
\]
On a quadratic time-delayed reaction-diffusion model, oversampled two-sided structured interpolation outperformed POD by about four orders of magnitude in time-domain accuracy and about six orders of magnitude on the first symmetric transfer function. On the Toda lattice molecular-dynamics example, preserving the internal block structure through split congruence transformations was important for stable reduced simulations.

## 5. Behavioral and direct data-driven trajectory interpolation

A distinct use of system interpolation appears in the behavioral LPV setting of [2508.11332]. The systems considered are discrete-time LPV systems with shifted-affine dependence,
\[
\sum_{i=0}^{n_\ell} r_i(q^i p)\,q^i w = 0, \qquad
r_i(q^i p)=r_{i,0}+\sum_{j=1}^{n_p} r_{i,j}\,q^i p_j,
\]
and the goal is to recover missing values in a partially specified finite trajectory without identifying an explicit parametric model. Given a data dictionary
\[
D=(\breve w_{[1,N_d]},\breve p_{[1,N_d]})\in \mathcal B_{[1,N_d]},
\]
the key representation theorem states that every admissible length-\(L\) trajectory is characterized by a Hankel equation involving both \(\breve w\) and the lifted product signal \(\breve w^{\breve p}\), under the generalized persistence of excitation condition
\[
\operatorname{rank}
\begin{pmatrix}
\mathcal H_L(\breve w_{[1,N_d]})\\
\mathcal H_L(\breve w^{\breve p}_{[1,N_d]})
\end{pmatrix}
= n + (m+n_p n_w)L.
\]

Interpolation is then formulated as exact missing-data completion for a given scheduling trajectory \(\tilde p_{[1,L]}\) and a set \(I_g\) of known entries. Existence is characterized by a rank/column-space compatibility condition, and uniqueness by
\[
\operatorname{rank}\big(\mathcal H_L(w_{[1,N_d]})|_{I_g}\,\mathcal N_{\tilde p}\big)
=
\operatorname{rank}\big(\mathcal H_L(w_{[1,N_d]})\,\mathcal N_{\tilde p}\big)
= n+mL,
\]
where \(\mathcal N_{\tilde p}\) spans the kernel of
\[
N_{\tilde p}:=\mathcal H_L(w^{\breve p}_{[1,N_d]})-
\tilde{\mathcal P}\mathcal H_L(w_{[1,N_d]}).
\]
Algorithmically, one solves a linear system for \(g\),
\[
\begin{bmatrix}
\mathcal H_L(w_{[1,N_d]})|_{I_g}\\
\mathcal H_L(w^p_{[1,N_d]})-\tilde{\mathcal P}\mathcal H_L(w_{[1,N_d]})
\end{bmatrix} g
=
\begin{bmatrix}
\operatorname{vec}(w_{I_g})\\
0
\end{bmatrix},
\]
and then reconstructs
\[
\operatorname{vec}(\tilde w_{[1,L]})=\mathcal H_L(w_{[1,N_d]})g.
\]

The mass-spring-damper example illustrates the behavioral interpretation. With \(m=d=1\), \(\kappa_0=10\), \(\kappa_1=9\), \(\tau=0.1\), the paper identifies \(n_w=2\), \(n_p=1\), \(n=2\), \(L=2\), and \(m=1\). For horizon \(L=30\), at least \(N_d=121\) data points were needed. With \(K=35\) known samples, Conditions 1–3 held and the interpolated trajectory matched the true trajectory up to numerical precision; with \(K=10\), existence held but uniqueness failed, yielding infinitely many admissible interpolants. This formulation treats simulation as a special case of interpolation and extends naturally to approximation when exact compatibility fails.

## 6. Interpolation in PDE regularity and solution families

In PDE analysis, system interpolation often means interpolation between function spaces rather than sample values. For the stationary 3D Navier–Stokes system in a bounded smooth domain,
\[
-\nu\Delta v +(v\cdot \nabla)v+\nabla p=f,\qquad \nabla\cdot v=0,\qquad v|_{\partial\Omega}=0,
\]
the paper [2409.18808] uses the inequality
\[
|u|^{(l)}_{\Omega}\le
C\bigl(|u|^{(2+\alpha)}_{\Omega}\bigr)^{\omega}
\bigl(\|u\|_{L^q(\Omega)}\bigr)^{1-\omega},
\qquad
\omega=\frac{lq+3}{q(2+\alpha)+3},
\]
to interpolate between \(L^q(\Omega)\) and \(C^{2+\alpha}(\Omega)\). Taking \(q=6\), obtained from \(W_2^1(\Omega)\hookrightarrow L^6(\Omega)\), the paper estimates the convection term \((v\cdot\nabla)v\) in \(C^\alpha(\Omega)\) and derives
\[
|(v\cdot \nabla)v|^{(\alpha)}_\Omega
\le C\bigl(|v|^{(2+\alpha)}_\Omega\bigr)^{a_1}
\bigl(\|v\|_{L^6(\Omega)}\bigr)^{a_2},
\qquad
a_1=\frac{6(2+\alpha)+6}{6(2+\alpha)+3}<1.
\]
Young’s inequality then yields
\[
|(v\cdot \nabla)v|^{(\alpha)}_\Omega
\le \varepsilon |v|^{(2+\alpha)}_\Omega
+ C_\varepsilon \bigl(\|v\|_{L^6(\Omega)}\bigr)^B,
\]
with
\[
B = 2-\frac{6}{6(2+\alpha)+3},
\]
and the final a priori bound
\[
|v|^{(2+\alpha)}_{\Omega}+|\nabla p|^{(\alpha)}_{\Omega}
\le C|f|^{(\alpha)}_{\Omega}
+ C\bigl(\|v\|_{L^6(\Omega)}\bigr)^B.
\]
Here interpolation is the bridge between weak Sobolev control and Schauder regularity for a nonlinear system.

A different PDE meaning appears in the two-component Camassa–Holm system [1402.1060]. There, interpolation is between conservative and dissipative global weak solutions at wave breaking. In Lagrangian coordinates, the energy variable is split into \(h\) and an effective part \(\bar h\). At each breaking time \(\tau_n(\xi)\),
\[
\bar h(\tau_n(\xi),\xi)=(1-\alpha)\bar h(\tau_n(\xi)-0,\xi),\qquad
h(\tau_n(\xi),\xi)=h(\tau_n(\xi)-0,\xi),
\]
with \(\alpha\in[0,1]\). Thus
\[
\alpha=0 \Rightarrow \text{conservative},\qquad
\alpha=1 \Rightarrow \text{fully dissipative}.
\]
In Eulerian variables, the singular part of the energy measure is multiplied by \(1-\alpha\) across breaking:
\[
\mu_{\mathrm{ac}}=\mu^-_{\mathrm{ac}},\qquad
\mu_{\mathrm{s}}=(1-\alpha)\mu^-_{\mathrm{s}}.
\]
All \(\alpha\)-dissipative solutions solve the same 2CH system weakly; what changes is the fraction of concentrated energy retained after collision. This is a literal continuous interpolation of solution behavior at the system level.

## 7. Interpolation-generated dynamics, proof systems, and parameter-space cognition

In integrable-systems theory, interpolation can be the source from which the system itself is derived. The elliptic interpolation construction of [1706.05155] starts from
\[
\psi(z)=\prod_{i=1}^{2N}\frac{\Gamma_{p,q}(u_i/z)}{\Gamma_{p,q}(k/u_i z)},
\qquad \prod_{i=1}^{2N}u_i=k^N,
\]
and asks for elliptic functions \(P(z),Q(z)\) such that
\[
\psi(q^{-s})=\frac{P(q^{-s})}{Q(q^{-s})},\qquad s=0,1,\dots,N.
\]
Casorati determinants of \((P,\psi Q)^T\) factor into residual spectral polynomials \(F(z)\) and \(G(z)\), whose zeros become dynamical variables. After gauge normalization, the contiguous relations \(L_2,L_3\) form a Lax pair, and compatibility yields the elliptic Garnier system. For \(N=3\), the construction reduces to the elliptic Painlevé equation. The multiple interpolation problem of [2206.12502] plays a closely related role for discrete integrable systems. There a determinant solution of
\[
Y_1(x_s)f_1(x_s)+\cdots+Y_m(x_s)f_m(x_s)=0,
\qquad s=1,\dots,|n|+m-1,
\]
produces \(\tau\)-functions satisfying Hirota’s discrete KP equations and, after the interpolation-specific constraint is imposed, a non-autonomous multidimensional Toda system. In the \(m=2\) reduction, the resulting equations become the non-autonomous discrete-time Toda lattice.

Proof complexity introduces a logical meaning of effective system interpolation [2601.02821]. The paper studies the uniform effective disjunction property and the uniform effective interpolation property for \(EF\) and stronger normal proof systems \(S\). If \(EF\) has the uniform effective disjunction property, then every sufficiently strong \(S\) corresponding to a theory \(T\) in the language of \(V_1^1\) also has it. If \(EF\) has the uniform effective interpolation property, then so does \(S\). A major consequence is that, under this hypothesis, every disjoint \(NE\)-pair has a separator in \(E\), and therefore
\[
NE\cap coNE = E.
\]
The paper also proves that if \(A_1,A_2\) is an \(NE\)-pair with \(A_1\cup A_2=\mathbb N\), then an exponential-time algorithm can, for each input \(n\), find some \(i\in\{1,2\}\) such that \(n\in A_i\). Here “interpolation” is not approximation but effective proof-theoretic extraction.

A final, very different use appears in checkpoint interpolation for large language models [2601.21414]. The paper defines a parameter-space interpolation between an Instruct checkpoint \(\Theta^{(\mathrm I)}\) and a Thinking checkpoint \(\Theta^{(\mathrm T)}\),
\[
\Theta^{(\mathrm M)}=\lambda \Theta^{(\mathrm T)}+(1-\lambda)\Theta^{(\mathrm I)},
\]
and then makes the coefficient query-dependent,
\[
\Theta^{(\mathrm M)}(q)=\lambda(q)\Theta^{(\mathrm T)}+(1-\lambda(q))\Theta^{(\mathrm I)}.
\]
Its pilot study reports a convex, monotonic Pareto frontier between accuracy and token usage, with representation continuity supported by a correlation \(r=0.974\) between the first PCA component of the \(\langle\text{think}\rangle\) representation and \(\lambda\). The DAMI framework estimates \(\lambda(q)\) either by preference learning over accuracy/cost tradeoffs or by a zero-shot confidence discrepancy signal,
\[
S_{\mathrm{final}}=S_{\mathrm{ambi}}+S_{\mathrm{dis}},
\qquad
\lambda(q)=\sigma\left(\frac{S_{\mathrm{final}}(q)-\mu}{\tau}\right).
\]
On five mathematical reasoning benchmarks, the paper reports that DAMI-Pref improves accuracy over the Thinking model by 3.4 points with 29% fewer tokens, while DAMI-Conf improves accuracy by 2.5 points with 40% fewer tokens. This suggests a parameter-space notion of system interpolation in which the interpolated object is not a transfer function or trajectory but the model’s cognitive configuration itself.

Across these domains, system interpolation is less a single technique than a recurring structural idea: interpolation becomes “systemic” when it is governed by architecture, dynamics, admissibility, proof theory, or parameter geometry, and when its output modifies how the system as a whole behaves.

Source: https://www.emergentmind.com/topics/system-interpolation