REVIEW 4 major objections 4 minor 1 cited by
This paper claims that the memory behaviors of modern state space models can be rebuilt as explicit, interpretable orthogonal-polynomial mechanisms.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
T0 review · deepseek-v4-flash
2026-08-02 21:02 UTC pith:DVF3VRGT
load-bearing objection A competent, cleanly written paper that adds five interpretable HiPPO extensions, but the abstract's claim about realizing modern SSM capabilities is too broad given the scalar time-warp limitation of Salience HiPPO and the thin synthetic evidence. the 4 major comments →
HiPPO Zoo: Explicit Memory Mechanisms for Interpretable State Space Models
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central claim is that capabilities typically associated with modern state space models can be realized through explicit, interpretable polynomial memory structures. HiPPO represents the past by projecting a signal onto orthogonal polynomials under a measure, with coefficients updated by a structured linear ODE. The paper shows five modifications: Volterra HiPPO puts nonlinear interactions into polynomial kernels in the readout; Salience HiPPO reshapes the history measure via a scalar time warp so memory is allocated to informative inputs; Associative Memory HiPPO adds a separate OP-based key-value bank with continuous addresses; Multiscale HiPPO embeds a continuum of timescales as a scal
What carries the argument
The load-bearing object is the HiPPO history representation: a signal's recent past is encoded as coefficients s(t) of an orthogonal polynomial basis defined by a measure ω(τ), updated online via a linear ODE ˙s = A s + b f. The paper's five extensions modify this substrate—the readout (Volterra kernels), the measure via a scalar time warp (Salience), the state space augmented with a separate OP key-value bank (Associative Memory), the state itself parameterized over inverse timescale (Multiscale), and the metric on histories induced by a forecasting map (Forecasting). The orthogonal polynomial basis is what keeps each mechanism explicit: coefficients correspond to named functionals of the p
Load-bearing premise
The load-bearing premise is that a single scalar salience signal multiplying the entire HiPPO dynamics can express the full range of selective memory that modern SSMs achieve with high-dimensional, input-dependent state updates; if real selectivity requires richer, coordinate-dependent gating, the Salience mechanism cannot capture it.
What would settle it
Give the model a task where two input channels must be remembered with different priorities at the same time step (e.g., attend to channel A while ignoring channel B at the same timestamp). Since Salience HiPPO scales all channels by the same scalar g(t), it cannot represent such channel-specific selectivity; a model with per-channel input-dependent gates would succeed. If the paper's mechanism fails this test while a gated SSM succeeds, the claim that scalar time-warping realizes modern selective memory is falsified.
If this is right
- If the claims hold, selective and associative memory behaviors in SSMs can be studied and controlled directly through the history measure and the OP basis, instead of being inferred from learned state transitions.
- Salience HiPPO provides an interpretable recipe for input-dependent memory: a learned scalar warps time, which is equivalent to adapting the memory measure; this could enable principled design of attention-like gating in streaming settings.
- Multiscale HiPPO shows a single linear system can represent a continuum of timescales with a fixed coefficient budget, so a model need not commit to a horizon a priori.
- Forecasting HiPPO shows that the training objective itself selects memory geometry: different horizons induce different predictive memories, which can be visualized as leading eigenfunctions of the induced metric.
- The Volterra HiPPO readout shows that nonlinear system identification can be done with interpretable kernels, with faster convergence than an MLP readout in the paper's experiment.
Where Pith is reading between the lines
- Inference: The scalar time-warping equivalence in Salience HiPPO suggests a precise diagnostic for how 'selective' any SSM actually is: if a model's selectivity can be mimicked by a single warped time axis, it is not using richer gating; future work could test this by fitting a Salience HiPPO to a trained selective SSM's input-output map.
- Inference: The OP associative memory's reproducing-kernel interference structure is a concrete model of crosstalk between stored associations; this could be used to predict memory failures in tasks with many overlapping write addresses, a testable quantitative prediction.
- Inference: The Multiscale construction's decoupling via spectral decomposition of the scale-coupling matrix implies that its expressivity is equivalent to a bank of fixed-timescale HiPPO systems; so its value is queryability and parameter sharing, not new dynamics—a claim the paper itself nearly states, and that could be probed by comparing to a tuned bank with the same total dimension.
- Inference: The paper leaves open whether these explicit mechanisms scale to real benchmarks; a natural extension would be to plug Volterra or Salience mechanisms into a modern SSM architecture and measure whether interpretability survives while performance holds.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper revisits the HiPPO framework and proposes five explicit memory mechanisms: Volterra HiPPO (polynomial nonlinear readouts), Salience HiPPO (scalar time-warping of the history measure), Associative Memory HiPPO (OP-based key-value memory with minimum-norm writes), Multiscale HiPPO (polynomial embedding of inverse timescale), and Forecasting HiPPO (predictive memory geometry Q=T^T W T). Each extension is paired with a synthetic experiment and interpretability diagnostics. The central claim is that capabilities usually implicit in learned SSMs—adaptive memory allocation, associative recall, multiscale coverage, and forecast-induced memory geometry—can be realized explicitly in polynomial memory structures.
Significance. If the claims hold, the paper provides a useful, transparent toolbox for interpreting memory in state space models, and the derivations are mostly elementary and checkable. The Volterra substitution, the time-warp identity, the associative-memory write/read construction, and the Q=T^T W T forecasting geometry are clean and clearly explained. I also credit the authors for explicitly disclaiming expressivity gains for Multiscale HiPPO (App. D.4). However, the breadth of the central claim exceeds what the five proof-of-concept experiments establish, especially for Salience HiPPO, where the mechanism is a scalar time reparameterization and is only tested on a task requiring temporal, not feature-selective, prioritization. The contribution is still valuable as an explicit-mechanism framework, but the claim should be calibrated to the demonstrated scope.
major comments (4)
- [Abstract, §3.2, Eq. (2)] The claim that Salience HiPPO realizes adaptive memory allocation "typically associated with modern SSMs" is broader than the mechanism. Eq. (2) multiplies the entire state dynamics by a single scalar g(t), and Eqs. (3)–(4) show this is equivalent to a standard HiPPO system on a reparameterized time axis. Consequently, all memory prioritization is temporal: there is no mechanism to gate different channels or input features differently at the same timestep, unlike the per-channel/input-dependent updates in Mamba-type models. The selective-copying task (App. C.2) interleaves informative and uninformative tokens at different times, so scalar time warping suffices; it never tests a task where salience differs across features within a single timestep. I recommend either adding a per-channel or coordinate-dependent salience extension and a task that requires within-timestep selectivity, or exp
- [§3.2, Eqs. (3)–(4), App. C.2] The time-warping equivalence treats g as an exogenous function of time. In the experiment, g is produced by a learned network that conditions on the current HiPPO state, so g is state-dependent. Then φ(t)=∫g(s)ds depends on the trajectory itself, and Eq. (4) is not a model-level equivalence to a fixed HiPPO system; it is a post hoc reparameterization of each realized trajectory. This does not invalidate the adaptive-memory idea, but the paper should clarify that the "equivalent standard HiPPO" statement is only exact for exogenous g, or reformulate the discussion as a state-dependent time warp.
- [App. C.3, Eq. (C.3)] The paper states that the associative-memory write solves a constrained minimum-norm problem that achieves the desired value at the key address, so that m_j(x_key)=y_t[j]. However, the reported update Δ = α(y−ŷ)k with α = g_write/(||k||^2+ε) does not exactly satisfy that constraint for ε>0. The exact solution of the stated constrained problem is Δ = g_write(y−ŷ)k/||k||^2. As written, the update is approximate. If ε is a numerical regularizer, the text should say so; otherwise the "closed-form minimum-norm write and minimal interference" claim is stronger than what is implemented.
- [§3.2 / App. C.2, discretization] Eq. (2) is a continuous-time equation, but the experiments are discrete-time. For a time-varying g(t), the exact zero-order-hold discretization is not the usual constant-g ZOH update; for piecewise-constant g within a step it is s_{t+1}=exp(Δt g_t A)s_t + (g_t A)^{-1}(exp(Δt g_t A)-I)b f_t. The appendix does not specify how g_t enters the discrete update actually used for training and evaluation, nor whether an Euler or exact scheme was used. Since the time-warp interpretation and the induced-measure visualizations depend on continuous-time integration, this implementation detail should be documented and the exact discrete update stated.
minor comments (4)
- [§3.2, Eq. (4)] In Eq. (4), the input should be written as f(φ^{-1}(t1)) rather than f(t1) unless a new symbol for the warped-time input is introduced.
- [§3.1] The sentence defining the polynomial sum says "i0, . . . , ik" but the sum indexes i1,...,ik; this is a typo.
- [Table 1] The table entries "bnbk" should be typeset as products b_n b_k (and similarly for the Leg-S entries), and the sign convention for Leg-T should be checked against the standard HiPPO definitions.
- [App. C.4] The multiscale base timescale τ0 is introduced but its numerical value in the experiments is never stated. Since the query mapping g(L)=τ0/L depends on it, the value should be reported.
Circularity Check
No significant circularity: each zoo extension is constructed from explicit definitions and fitted parameters are used for interpretability, not as independent predictions.
full rationale
The paper's derivation chain is self-contained and does not reduce to its inputs. Volterra HiPPO defines a polynomial readout of the HiPPO state and then fits the Volterra coefficients to a synthetic system; the fitted coefficients are visualized, not used as an independent prediction. Salience HiPPO's time-warp equivalence (Eqs. 2–4) is an exact change-of-variables identity, and the salience signal is learned, with the induced history measure computed post hoc from the trained g(t). Associative Memory HiPPO constructs an explicit OP key-value bank with a closed-form minimum-norm write update; solving the synthetic recall task is the intended behavior of that construction, not a separate derived result. Multiscale HiPPO substitutes a polynomial expansion into the HiPPO dynamics to derive coupled linear dynamics, and Appendix D.4 explicitly acknowledges that this is no more expressive than a bank of HiPPO systems—an honest statement that prevents overclaiming. Forecasting HiPPO fits a reduced-rank map T and then defines Q=T^T T as the induced history metric; this is an interpretability map of the fitted object, not a second, independently predicted result. The paper does not relabel fitted quantities as predictions, and there are no load-bearing self-citations: the cited HiPPO and SSM works are external and not authored by the present paper's authors. All experiments fit their own parameters and report the resulting behavior, which is consistent with the stated goal of illustration and interpretability rather than independent validation. Therefore, no circular step is identifiable under the specified criteria.
Axiom & Free-Parameter Ledger
free parameters (5)
- Leg-S timescale factor sigma (alpha=0.375) =
0.375
- Wray-Green system constants a,m,k,tau_max,alpha =
a=2, m=0.3, k=0.08, tau_max=50, alpha=0.004 (scaled 4e-3)
- Multiscale HiPPO base timescale tau0 =
not explicitly stated numerically
- Forecasting HiPPO horizons H =
H=4 and H=32
- Rank d of reduced-rank regression =
d not reported in main text or appendix
axioms (5)
- domain assumption Volterra series converge / truncating to first few kernels is adequate (Sec. 3.1).
- domain assumption The HiPPO state truncation f(t-tau) ≈ Σ s_n(t)P_n(tau) is accurate enough to support the Volterra readout (Eq. 3.1, C.1).
- domain assumption The reproducing-kernel / minimum-norm write update is a sufficient model of associative binding (Eq. C.3).
- standard math Jordan/matrix-exponential machinery for G=exp(J) is numerically well-conditioned (App. D.3).
- domain assumption The RBF-mixture signal model and streaming statistics give a valid estimate of the optimal linear predictor (C.5).
invented entities (1)
-
'predictive memory' and Q=T^T T as 'objective-induced history geometry'
no independent evidence
read the original abstract
Representing the past in a compressed, efficient, and informative manner is a central problem for systems trained on sequential data. The HiPPO framework, originally proposed by Gu & Dao et al., provides a principled approach to sequential compression by projecting signals onto orthogonal polynomial (OP) bases via structured linear ordinary differential equations. Subsequent works have embedded these dynamics in state space models (SSMs), where HiPPO structure serves as an initialization. Nonlinear successors of these SSM methods such as Mamba are state-of-the-art for many tasks with long-range dependencies, but the mechanisms by which they represent and prioritize history remain largely implicit. In this work, we revisit the HiPPO framework with the goal of making these mechanisms explicit. We show how polynomial representations of history can be extended to support capabilities of modern SSMs such as adaptive memory allocation and associative memory, while retaining direct interpretability in the OP basis. We introduce a unified framework comprising five such extensions, which we collectively refer to as a "HiPPO zoo." Each extension exposes a specific modeling capability through an explicit, interpretable modification of the HiPPO framework. The resulting models adapt their memory online and train in streaming settings with efficient updates. We illustrate the behaviors and modeling advantages of these extensions through a range of synthetic sequence modeling tasks, demonstrating that capabilities typically associated with modern SSMs can be realized through explicit, interpretable polynomial memory structures.
Figures
Forward citations
Cited by 1 Pith paper
-
SHiPPO: Recurrent Memory with Transported Polynomial Projections
SHiPPO transports HiPPO coefficient memories via right actions into Sylvester dynamics, and diagnostics show this recovers order-sensitive memory changes that high-rank writes cannot.
Reference graph
Works this paper leans on
-
[1]
We refer the reader to the survey of Totik and the book of Chihara for more information on orthogonal polynomials (Totik, 2005; Chihara, 1978). Orthonormal Legendre polynomials on [0,1] .Let Pn(x) denote the standard Legendre polynomials on [−1,1] , defined by the recurrence (n+ 1)P n+1(x) = (2n+ 1)xP n(x)−nP n−1(x), P 0(x) = 1, P1(x) =x, and orthogonal w...
2005
-
[6]
The solution follows from a Lagrange multiplier and yields the closed-form update Ct+1[j] =C t[j] +α t yt[j]−ˆyt[j] kt, α t = gwrite ∥kt∥2 2 +ε
and a full write (gwrite = 1). The solution follows from a Lagrange multiplier and yields the closed-form update Ct+1[j] =C t[j] +α t yt[j]−ˆyt[j] kt, α t = gwrite ∥kt∥2 2 +ε . Because the basis is orthonormal, ∥kt∥2 2 =⟨k t, kt⟩ is well-conditioned. Interpreted in function space, Eq. (C.3) adds a localized bump proportional to the truncated reproducing k...
2000
-
[7]
Leg-T corresponds to the truncated Legendre HiPPO system defined on a finite window, while Leg-S corresponds to the exponentially weighted (sliding) Legendre HiPPO system
The nonlinear benchmark introduced by Wray & Green (1994) is a purely second-order V olterra system with a separable 11 HiPPO Zoo ω(τ)P n(τ)b n Ank Leg-T U[0,1]L n(τ) √2n+ 1 bnbk(−1)n−k+1, n > k bnbk, n=k 0, n < k Leg-S exp(−τ)L n(e−τ ) √2n+ 1 0, n < k −bnbk n+ 1 2n+ 1 , n=k −bnbk, n > k Table 1.Orthogonal Polynomial HiPPO systems used i...
1994
-
[1978]
ISBN 0-677-04150-0. Dao, T. and Gu, A. Transformers are ssms: Generalized models and efficient algorithms through structured state space duality.arXiv preprint arXiv:2405.21060,
-
[2022]
Flow- state: Sampling rate invariant time series forecasting
Graf, L., Ortner, T., Wo´Ls ¸niak, S., Pantazi, A., et al. Flow- state: Sampling rate invariant time series forecasting. arXiv preprint arXiv:2508.05287,
-
[2023]
Orthogonal polynomials.arXiv preprint math/0512424,
Totik, V . Orthogonal polynomials.arXiv preprint math/0512424,
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.