arXiv paper· applied dynamical systems / data-driven operator theory· method paper

Functional dynamic mode decomposition: Learning infinite-dimensional systems from data

Stefan Klus, Eirini Ioannou

math.DS · cs.LG · stat.ML

Quick verdict

Functional DMD extends DMD to infinite-dimensional systems via functional data; elegant operator theory, but validation is synthetic-only, no error bounds.

Summary

The paper extends projected and exact dynamic mode decomposition (DMD) from matrices to infinite-dimensional Hilbert spaces, learning finite-rank evolution operators directly from functional snapshot data — observables, densities, or wavefunctions — rather than from vector-valued measurements. Training pairs of functions define functional data vectors whose Gram matrices yield projected functional DMD (a Galerkin projection, A = C_uu⁻¹C_uv), while the SVD and pseudoinverse of the associated operator give exact functional DMD (Ã = V U⁺). The authors prove the best-fit property (Theorem 2.9), verify the eigenfunction construction (Theorem 2.16), show that classical projected/exact DMD arise as special cases under grid discretization (Examples 2.6, 2.17), and relate the framework to generalized EDMD. Applications cover the 1-D heat equation (first four frequencies 1.00–4.02 recovered, the fifth degraded), random walks on graphons with metastable-set detection and graphon reconstruction up to a constant, overdamped Langevin dynamics on the Himmelblau potential (eigenvalues 0.851/0.695 versus EDMD 0.854/0.699), Koopman–von Neumann mechanics with unit-circle spectra, and the Kuramoto–Sivashinsky equation. Main limitations: all examples are synthetic, no convergence rates or error bounds are derived, validation is mostly visual, and the nonlinear PDE case remains exploratory.

Review complete

Assessment at a glance

applied dynamical systems / data-driven operator theory

3

risk flags identified

● high 0

● medium 2

● low 1

Reproducibility readinesshigh
Code & data availableyes
Review confidencehigh
AI text signal5/20 passages flagged
Deep research viableCode-reproducible

5 strengths · 5 weaknesses · 5 author questions

AI text signal = Binoculars zero-shot detector score over 20 sampled passages (Fast-DetectGPT co-reported) — absence does not prove human authorship, and a flag does not imply misconduct.

Levels reflect verification complexity — not a quality score. Evolve never rates your paper.

Key claims under review

  • Projected and exact DMD extend to separable Hilbert spaces, learning finite-rank operators directly from functional snapshot pairs via Gram matrices and operator pseudoinverses.
  • The projected functional DMD matrix is the minimizer of the training prediction error (Theorem 2.9), and projected/exact variants share identical eigenvalues because their matrices are similar (Section 2.3.4).
  • Classical projected and exact DMD algorithms are special cases of functional DMD obtained by discretizing the training functions on a grid (Examples 2.6, 2.17).
  • Functional DMD identifies invariant densities, metastable sets, and graphons (up to a multiplicative constant) from density-only data, without tracking individual trajectories.
  • The framework applies to Koopman, Perron–Frobenius, and Koopman–von Neumann operators for graphons, ODEs, and SDEs, with nonlinear infinite-dimensional extensions deferred to future work.

Strengths

  • Rigorous operator-theoretic development

    Lemma 2.7, Theorem 2.9 (best-fit minimization), Theorem 2.16 (verification that Algorithm 2.15 outputs eigenfunctions of à = V U⁺), and Lemma 2.20 (exact-DMD eigenfunctions project onto their projected counterparts) are all proven in-line; the similarity of A and à is shown explicitly in Section 2.3.4.

  • Unifying reduction to classical DMD

    Examples 2.6(ii) and 2.17 show that gridding the training functions turns C_uu, C_uv into UᵀU, UᵀV, recovering projected and exact DMD exactly — this clarifies the implicit infinite-dimensional meaning of the standard algorithms rather than merely re-deriving them.

  • Verifiable benchmarks with quantitative cross-checks

  • Density-only capability without particle tracking

  • Candid limitations and released code

Weaknesses

  • No convergence or perturbation theory

  • Validation is overwhelmingly visual

  • Synthetic examples only

  • Kuramoto–Sivashinsky case is underdeveloped

  • Thin, asymmetric baseline comparison

Full Review Report

$1.99 · one-time

Unlock the complete preliminary analysis Evolve produced for this paper.

Report preview

Scope & Claims

The paper (arXiv 2609.29159; math.DS/cs.LG/stat.ML) proposes *functional DMD*: a reformulation of projected and exact dynamic mode decomposition in which the learned object is a finite-rank operator on a separable Hilbert space H (L², Sobolev, or RKHS) rather than an m×m matrix built from vector snapshots. Central claims: (i) projected functional DMD is a Galerkin projection with representation matrix A = C_uu⁻¹C_uv (Section 2.2); (ii) exact functional DMD uses the SVD/pseudoinverse of the 'functional data operator' U, giving à = V U⁺ = Σ⁻¹ΘᵀC_uv ΘΣ⁻¹ (Section 2.3); (iii) classical DMD algorithms are special cases obtained by gridding the training functions (Examples 2.6, 2.17); (iv) the framework covers Koopman, Perron–Frobenius, and Koopman–von Neumann operators for graphons, ODEs, and SDEs. Claims are appropriately scoped to linear dynamics; nonlinear extensions are deferred.

The remaining sections cover Methodology Assessment, Statistical Evidence, Results & Reproducibility, Limitations & Risks, and Verdict & Recommendations.

  • All 5 strengths & 5 weaknesses, analyzed in detail
  • 3 risk flags with severity, section & reviewer questions
  • Methods & statistical evidence assessment
  • 5 questions a peer reviewer would ask
  • Full six-dimension evaluation report & recommendations

One-time purchase · this paper only

One-time purchase, this paper only. AI-generated preliminary analysis — not peer review.

Frequently asked questions

What is functional dynamic mode decomposition (functional DMD)?

A variant of DMD that works directly on functions — observables, probability densities, or wavefunctions — living in an infinite-dimensional Hilbert space. Instead of building matrices from vector snapshots, it learns a finite-rank evolution operator from Gram matrices of the functional training data, so systems like PDEs need not be turned into finite-dimensional problems before analysis.

How does functional DMD relate to standard DMD?

When the training functions are represented by their values on a spatial grid, the functional Gram matrices become the usual products UᵀU and UᵀV, and projected/exact functional DMD reduce exactly to classical projected and exact DMD. The paper proves that classical DMD algorithms are special cases of the functional versions.

What can functional DMD be used for?

Estimating spectra and forecasting for heat-type PDEs, detecting metastable sets and reconstructing graphons from random-walk density data, finding invariant densities and timescales of Langevin (SDE) dynamics from density estimates alone — without tracking individual particles — and analyzing Koopman–von Neumann mechanics, including potential quantum-simulation applications.

Is there code available for functional DMD?

Yes. The authors state the supporting code is available in the open-source d3s repository (github.com/sklus/d3s). All examples in the paper are synthetic benchmarks (heat equation, graphon random walks, Himmelblau potential, Koopman–von Neumann mechanics, Kuramoto–Sivashinsky) that are reproducible on a standard CPU.

About this review

This is an AI-generated preliminary review by Evolve.Pub, produced from the openly available text of the paper. It is not peer review and does not represent the authors' views. The canonical record of the paper is at the publisher link above.

  • Figures are available only as captions in the extracted text; visual claims (eigenfunction accuracy, graphon reconstruction quality, KS eigenvalue plots) could not be independently verified.
  • Some quantitative values exist only inside figures (e.g., eigenvalue spreads in Figures 4a, 7b, 8b) and were taken from captions and surrounding text.
  • Several numeric fractions lost their fraction bars in extraction (e.g., inverse temperature and lag times in Examples 3.6/3.10); they were interpreted from context and not relied upon for judgments.
  • Full text including references appears complete, supporting high confidence despite the figure limitation.

Go deeper on this paper

Wrote a paper? Stress-test it before reviewers do.

Upload your manuscript for the same structured scrutiny — reasoning, statistics, evidence, figures, and reproducibility.

Scan your paper
Functional dynamic mode decomposition: Learning infinite-dimensional systems from data | Evolve.Pub