REVIEW 2 major objections 5 minor 1 cited by
Structure-Preserving Generalized Manifold Galerkin Reduction for Port-Hamiltonian Systems
T0 review · 2 major / 5 minor · reviewed 2026-08-04 · deepseek-v4-flash
Pith's one-line read A generalized manifold Galerkin reduction preserves the port-Hamiltonian form for any approximation map whose tangent space contains the port matrix and avoids the singular set of (J−R)^{-1}.
desk verdict Clean pH-preservation theorem built on GMG reduction; the generic non-degeneracy claim is unproved and the numerics are in-sample, but the core extension is real and worth refereeing. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing piece is the generalized manifold Galerkin (GMG) reduction map defined with structure matrix G = (J−R)^{−1}. For a given embedding Jacobian V = D_x ϕ, it produces the tangent projection W = (F_G(V))^T, where F_G(V) = (V^T G V)^{−1} V^T G. This projection becomes the 'reduction map' that turns the full-order residual into the reduced dynamics. The non-degeneracy condition V^T G V invertible and the port-inclusion condition span(B) ⊆ span(V) together force the reduced system to inherit the pH structure. The paper constructs embeddings of the form ϕ(ˇx) = B ˇx_1 + V η(ˇx_2) so that the port condition holds automatically, and instantiates η as linear and quadratic maps.
What would settle it
Compute det(V(ˇx(t))^T (J−R)^{−1} V(ˇx(t))) along the reduced trajectory for GMG-QM-ROM on a mass-spring-damper with a strong input; if it crosses zero at any time, the ROM is no longer a pH system (or the projection breaks), disproving the genericity claim.
Extended reading notes
Core claim
The core discovery is that the pH structure survives the GMG reduction under two explicitly checkable conditions on the approximation map's Jacobian V: span(B) ⊆ span(V) and V ∈ S_{(J−R)^{−1}}, i.e., V^T (J−R)^{−1} V invertible. Under these conditions the reduced system is exactly of port-Hamiltonian form with ˇJ = W^T J W, ˇR = W^T R W, ˇH = H∘ϕ, and ˇB = W^T B, where W = (F_{(J−R)^{−1}}(V))^T. The proof works by rewriting the GMG projection so that it factors into the desired pH structure, and the output equation is matched through the port-inclusion condition. This is the first such result allowing a completely general nonlinear embedding map.
Load-bearing premise
The construction requires that det(J−R) ≠ 0 and that, for every reduced state along the trajectory, the Jacobian of the nonlinear embedding satisfies det((D_xϕ)^T (J−R)^{−1} D_xϕ) ≠ 0; the paper assumes this non-degeneracy for all reduced states without giving a constructive guarantee that it holds away from the training data.
Editorial extensions
If this is right
- With the two conditions satisfied, reduced-order models of both linear and nonlinear pH systems remain provably port-Hamiltonian, preserving passivity and the dissipation inequality.
- The framework removes the need to restrict to linear embeddings or special separable ansatzes, opening the door to data-driven nonlinear embeddings (e.g., neural-network maps) as long as the tangent condition is enforced.
- The explicit reduced matrices (^J, ^R, ^B) and reduced Hamiltonian (^H) are given in closed form, so the resulting ROM can be simulated with standard port-Hamiltonian solvers.
- Combining the method with structure-preserving DEIM keeps the ROM cheap while retaining the pH structure, as the numerical examples illustrate.
- The quadratic embedding realization (GMG-QM-ROM) achieves lower state and output errors than the linear one and is close to the theoretical lower bound, indicating the approach can exploit curvature of the solution manifold.
Reading between the lines
- The paper asserts that the non-degeneracy condition is 'generically satisfied' but provides no proof; a measure-theoretic or probabilistic statement about the set of initial conditions/trajectories would turn this into a rigorous guarantee.
- The assumption det(J−R)≠0 excludes significant classes of pH systems, such as lossless systems with singular J or systems with dependent ports; extending the framework to differential-algebraic pH systems would broaden the reach.
- The quadratic embedding could be replaced by higher-degree polynomial or neural-network maps without any change to the structural argument, provided the tangent-space condition is enforced during training—this suggests a direct route to learning-based structure-preserving MOR.
- The numerical experiments are limited to mass-spring-damper systems; testing on fluid or electrical networks would check whether the claimed accuracy gains persist outside mechanical examples.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript proposes a structure-preserving model-order-reduction framework for port-Hamiltonian (pH) systems based on the generalized manifold Galerkin (GMG) reduction. For a (possibly nonlinear) embedding \phi:\mathbb{R}^r\to\mathbb{R}^N, the reduced model is defined by enforcing the residual to be orthogonal with respect to the GMG projector W(\hat x)=F_{(J-R)^{-1}}(D_{\hat x}\phi)^T. Theorem 3.1 gives sufficient conditions—span(B)⊆span(D_{\hat x}\phi) and D_{\hat x}\phi(\hat x)\in S_{(J-R)^{-1}} for every \hat x—under which the resulting ROM is again a pH system, with \hat J=W^TJW, \hat R=W^TRW, \hat H=H\circ\phi, and \hat B=W^TB. The authors instantiate the framework with a linear embedding (GMG-POD-ROM) and a quadratic embedding (GMG-QM-ROM), both built from snapshot data, and use structure-preserving DEIM for nonlinear Hamiltonians. Numerical comparisons on a linear and a nonlinear mass-spring-damper system report lower relative state and output errors than the SP1-POD-ROM and SP2-POD-ROM baselines.
Significance. If the key non-degeneracy hypothesis can be guaranteed or checked for the constructed embeddings, the theorem provides an elegant and quite general structure-preserving reduction framework. The algebraic derivation is clean and self-contained, the construction of the reduced Hamiltonian and the constant port matrix in the special cases is natural, and the numerical study includes state errors, output errors, and energy-balance errors. The main weakness is that the central sufficient condition of Theorem 3.1 is never verified for the data-driven embeddings, and the abstract's claim that the condition is 'generically satisfied' is not substantiated anywhere in the paper. Because the GMG map W in Eq. (9) is undefined if the condition fails, this gap is load-bearing for the claimed generality.
major comments (2)
- [Sec. 4.1, Thm. 3.1, Eqs. (12), (17), (18)-(22)] The non-degeneracy condition D_{\hat x}\phi(\hat x)\in S_{(J-R)^{-1}} is asserted but never proved for the proposed embeddings. After Eq. (12), Section 4.1 verifies only condition (i), namely span(B)⊆span(D_{\hat x}\phi); condition (ii) is not addressed. For the quadratic map (17), D_{\hat x}\phi is affine in \hat x_2, so det((D_{\hat x}\phi)^T(J-R)^{-1}D_{\hat x}\phi) is a polynomial in the reduced state. The constructions (18)-(22) impose no constraint on its zeros, and Algorithm 2 gives no check even for the linear map, where V=[B,\bar V] is built only from a POD of X-BB^\dagger X. If the determinant vanishes on the ROM trajectory, W in (9) is undefined and the ROM (10) is not a pH system. The abstract's 'generically satisfied' claim therefore needs either a rigorous proof or an algorithmic safeguard/check.
- [Sec. 5, Eqs. (26)-(27)] The numerical errors are computed on the same trajectories used to construct the approximation maps: the snapshot matrix X in (13) is used to build the POD bases in (18)-(19) and the quadratic lifting M in (22), and the same X is then used to evaluate e_{x,red} and e_y. Thus the reported 'lower relative reduction error' is a training/interpolation error, not a predictive error on unseen inputs or initial conditions. The comparison is still informative, but the claim should be qualified as an interpolation result and, ideally, supplemented by a hold-out test or cross-validation to support the advertised accuracy advantage.
minor comments (5)
- [Table 1, GMG-QM-ROM row] The displayed formula reads 'B\hat x_1 + V_1\hat x_1 + V_2M(\hat x_2\otimes\hat x_2)'; the second term should presumably be V_1\hat x_2.
- [Sec. 5.1] The sentence introducing SP2-POD-ROM contains a missing citation: 'the structure-preserving MOR method presented in ,'.
- [Eq. (28)] The formula for the lower bound appears to be missing the relative-error denominator; as printed, it is a product of two norms without a division by \sum_i\|x_i\|_2^2.
- [Sec. 2.1 and Thm. 3.1] The definition (1) assumes constant J,R,B, but the reduced system in Theorem 3.1 has state-dependent \hat J,\hat R,\hat B. A short remark reconciling this with the definition (or extending the definition to state-dependent matrices) would improve readability.
- [Sec. 2.1, abstract] The assumption det(J-R)\neq 0 is not part of the standard definition of pH systems and is imposed for the entire paper. The abstract and introduction currently claim applicability to 'pH systems' without this caveat; the scope restriction should be stated explicitly in the abstract and revisited in the conclusion.
Circularity Check
No significant circularity: the structure-preservation theorem is a self-contained algebraic derivation; the main caveats are an unproved 'generic' claim and in-sample numerical validation.
full rationale
The central result, Theorem 3.1, is derived directly from the GMG reduction definition (9) and the definitions of the reduced quantities W, J, R, H, and B. The proof uses only W^T D_ x phi = I and span(B) subset span(D_ x phi) to rewrite the projected vector field as (W^T J W - W^T R W) grad_ x H(phi(x)) plus W^T B u, and to match the output equation. This is an algebraic identity, not an equivalence with its own conclusion. Condition (ii), D_ x phi in S_{(J-R)^{-1}}, is exactly the invertibility required for W in (9) to exist, so the theorem is a conditional statement rather than a tautology. The cited works [3] and [4] are prior work by a co-author, but they only supply the GMG framework and an approximation lower bound; neither is used to justify pH preservation, and the proof above is self-contained. Two caveats are correctness/validation concerns rather than circularity: the abstract claims the non-degeneracy conditions are 'generically satisfied' without a proof in Sections 3-4, and the numerical basis V, V1, V2 is fitted to the same snapshot matrix X (13) used to evaluate errors (26)-(29), so the reported accuracy is in-sample rather than a predictive test on unseen inputs. Neither caveat makes the derivation equivalent to its inputs by construction.
Assumptions & free parameters
free parameters (5)
- quadratic lifting dimension r_n =
8 (linear example); input to Algorithm 3
- regularization parameter λ_reg =
10^{-3} (linear example); max(0.2*e_proj(r), 10^{-2.5}) (nonlinear example)
- DEIM tolerance ε_DEIM =
10^{-8}
- reduced order r =
5-20 (linear), 6-20 (nonlinear)
- POD bases V1/V2 and DEIM basis U =
POD of snapshot matrices X−BB†X and q(x_i)
assumptions (6)
- domain assumption J=-J^T, R=R^T⪰0, and det(J−R)≠0 for the pH system
- ad hoc to paper The non-degeneracy condition Dˇxϕ(ˇx) ∈ S_{(J−R)^{−1}} holds for all ˇx
- standard math GMG reduction framework from [4] (point projection property and tangent reduction map)
- domain assumption Structure-preserving DEIM approximation from [6] yields a valid reduced Hamiltonian
- domain assumption Gauss-Legendre time integration preserves pH structure
- domain assumption Snapshot matrix X is representative of the dynamics being reduced
Cite this review
Pith. "Pith review of Structure-Preserving Generalized Manifold Galerkin Reduction for Port-Hamiltonian Systems." pith.science (2026). https://pith.science/paper/3MISP6SX
@misc{pith2026260308656,
author = {Pith},
title = {Pith review of: Structure-Preserving Generalized Manifold Galerkin Reduction for Port-Hamiltonian Systems},
year = {2026},
howpublished = {\url{https://pith.science/paper/3MISP6SX}},
note = {Machine review of arXiv:2603.08656}
}
read the original abstract
This paper considers structure-preserving model order reduction (MOR) techniques for port-Hamiltonian (pH) systems, which are typically derived from energy-based modeling. To keep favorable properties of \pH systems such as passivity in a reduced order model (ROM), we use structure-preserving methods in the reduction process. Although projection-based structure-preserving MOR methods for nonlinear pH systems based on nonlinear approximation ansatzes have recently been proposed, existing approaches typically rely on specific structures of the approximation map and the underlying pH system. To address this limitation, we propose a \MOR framework based on generalized manifold Galerkin (GMG) reduction. The resulting framework can employ general nonlinear approximation maps while preserving the pH structure. We establish sufficient conditions for structure preservation, show that the associated non-degeneracy conditions are generically satisfied. We further present linear and quadratic approximation maps within the proposed framework. Numerical examples for a linear and a nonlinear mass-spring-damper system show that the proposed \MOR methods have lower relative reduction error compared to existing methods.
Figures
Figures from the paper (1 more)
Forward citations
Cited by 1 Pith paper
-
Structure-Preserving Discretization and Model Reduction for Energy-Based Models
A Petrov-Galerkin discretization framework preserves discrete dissipation inequalities for a general class of energy-based models, including circuits, Cahn-Hilliard, and doubly nonlinear diffusion.
Reference graph
Works this paper leans on
-
[4]
P . Buchfink, S. Glas, B. Haasdonk, and B. Unger. Model reduction on manifolds: A differential geometric framework. Physica D: Nonlinear Phenomena, 468, 2024.doi:10.1016/j.physd.2024.134299
arXiv 2024
-
[1]
J. Annoni. Modeling for Wind Farm Control. PhD thesis, University of Minnesota, 2016
2016
-
[2]
A. C. Antoulas, C. A. Beattie, and S. Gugercin. Interpolatory Model Reduction of Large-Scale Dynamical Systems, pages 3–58. Springer, 2010.doi:10.1007/978-1-4419-5757-3_1
-
[3]
P . Buchfink, S. Glas, and B. Haasdonk. Approximation bounds for model reduction on polynomially mapped manifolds. Comptes Rendus. Mathématique, 362:1881–1891, 12 2024.doi:10.5802/crmath.632
-
[5]
F . Califano, R. Rashad, F . P . Schuller, and S. Stramigioli. Geometric and energy-aware decomposition of the Navier-Stokes equations: A port-Hamiltonian approach. Physics of Fluids, 33, 4 2021. doi:10.1063/5. 0048359
doi:10.1063/5 2021
-
[6]
Chaturantabut, C
S. Chaturantabut, C. A. Beattie, and S. Gugercin. Structure-preserving model reduction for nonlinear port-Hamiltonian systems. SIAM Journal on Scientific Computing, 38:B837–B865, 2016. doi:10.1137/ 15M1055085
2016
-
[7]
S. Chaturantabut and D. C. Sorensen. Nonlinear model reduction via discrete empirical interpolation. SIAM Journal on Scientific Computing, 32:2737–2764, 2010.doi:10.1137/090766498
-
[8]
Duindam, A
V . Duindam, A. Macchelli, S. Stramigioli, and H. Bruyninckx. Modeling and control of complex physical systems the port-Hamiltonian approach. Springer Science & Business Media, 2009. doi:10.1007/ 978-3-642-03196-0
2009
Show all 31 references
-
[9]
Falaize and T
A. Falaize and T . Hélie. Passive guaranteed simulation of analog audio circuits: A port-Hamiltonian approach. Applied Sciences (Switzerland), 6, 9 2016.doi:10.3390/app6100273
2016 doi
-
[10]
Y. Geng, L. Ju, B. Kramer, and Z. Wang. Data-driven reduced-order models for port-Hamiltonian systems with operator inference. Computer Methods in Applied Mechanics and Engineering, 442:118042, 7 2025. doi:10.1016/j.cma.2025.118042. 13 Structure-preserving model reduction on m...
2025
-
[11]
Greif and K
C. Greif and K. Urban. Decay of the Kolmogorov N-width for wave problems. Applied Mathematics Letters, 96:216–222, 2019.doi:10.1016/j.aml.2019.05.013
2019 doi
-
[12]
Gugercin, R
S. Gugercin, R. V . Polyuga, C. Beattie, and A. J. van der Schaft. Structure-preserving tangential interpola- tion for model reduction of port-Hamiltonian systems. Automatica, 48:1963–1974, 2012. doi:10.1016/j. automatica.2012.05.052
1963 doi
-
[13]
T . C. Ionescu and A. Astolfi. Families of moment matching based, structure preserving approximations for linear port-Hamiltonian systems. Automatica, 49:2424–2434, 2013.doi:10.1016/j.automatica.2013.05.006
2013 doi
-
[14]
T . C. Ionescu and A. Astolfi. Moment matching for nonlinear port-Hamiltonian and gradient systems. In IFAC-Proceedings, volume 46, pages 395–399, 2013.doi:10.3182/20130904-3-RO-2023.00060
2013
-
[15]
Jacob and H
B. Jacob and H. Zwart. Linear Port-Hamiltonian Systems on Infinite-dimensional Spaces. Springer Science & Business Media, 2012.doi:10.1007/978-3-0348-0399-1
2012 doi
-
[16]
Kawano and J
Y. Kawano and J. M. Scherpen. Structure preserving truncation of nonlinear port-Hamiltonian systems. IEEE Transactions on Automatic Control, 63:4286–4293, 2018.doi:10.1109/TAC.2018.2811787
2018
-
[17]
Kotyczka and L
P . Kotyczka and L. Lefevre. Discrete-time port-Hamiltonian systems based on Gauss-Legendre collocation. IFAC-PapersOnLine, 51(3):125–130, 2018
2018
-
[18]
Mehrmann and B
V . Mehrmann and B. Unger. Control of port-Hamiltonian differential-algebraic systems and applications.Acta Numerica, 32:395–515, 2023
2023
-
[19]
B. Moore. Principal component analysis in linear systems: Controllability observability and model reduction. IEEE Transactions on Automatic Control, 26:17–32, 1981.doi:10.1109/TAC.1981.1102568
1981
-
[20]
Morandin, J
R. Morandin, J. Nicodemus, and B. Unger. Port-Hamiltonian dynamic mode decomposition. SIAM Journal on Scientific Computing, 45:A1690–A1710, 2023.doi:10.1137/22M149329X
2023 doi
-
[21]
Ohlberger and S
M. Ohlberger and S. Rave. Reduced basis methods: Success, limitations and future challenges. In Proceedings of ALGORITMY, pages 1–12, 2016
2016
-
[22]
Pagliantini and F
C. Pagliantini and F . Vismara. Gradient-preserving hyper-reduction of nonlinear dynamical systems via discrete empirical interpolation. SIAM Journal on Scientific Computing, 45:A2725–A2754, 10 2023. doi: 10.1137/22M1503890
2023 doi
-
[23]
R. V . Polyuga and A. J. van der Schaft. Effort- and flow-constraint reduction methods for structure preserving model reduction of port-Hamiltonian systems. Systems and Control Letters, 61:412–421, 2012. doi:10.1016/ j.sysconle.2011.12.008
2012
-
[24]
Rashad, F
R. Rashad, F . Califano, A. J. van der Schaft, and S. Stramigioli. Twenty years of distributed port-Hamiltonian systems: A literature review, 12 2020.doi:10.1093/imamci/dnaa018
2020 doi
-
[25]
Rettberg, J
J. Rettberg, J. Kneifl, J. Herb, P . Buchfink, J. Fehr, and B. Haasdonk. Data-driven identification of latent port-Hamiltonian systems. arXiv preprint arXiv:2408.08185, 8 2024
2024 arXiv
-
[26]
K. Sato. Riemannian optimal model reduction of linear port-Hamiltonian systems. Automatica, 93:428–434, 7 2018.doi:10.1016/j.automatica.2018.03.051
2018 doi
-
[27]
P . Schulze. Structure-preserving model reduction for port-Hamiltonian systems based on separable nonlinear approximation ansatzes. Frontiers in Applied Mathematics and Statistics, 9:1160250, 2023. doi:10.3389/ fams.2023.1160250
2023
-
[28]
Schwerdtner and M
P . Schwerdtner and M. Voigt. SOBMOR: Structured optimization-based model order reduction.SIAM Journal on Scientific Computing, 45:A502–A529, 2023.doi:10.1137/20M1380235
2023 doi
-
[29]
Sirovich
L. Sirovich. Turbulence and the dynamics of coherent structures. ii. Symmetries and transformations.Quarterly of Applied mathematics, 45:573–582, 1987.doi:10.1090/qam/910463
1987 doi
-
[30]
A. J. van der Schaft and D. Jeltsema. Port-Hamiltonian systems theory: An introductory overview. Foundations and Trends in Systems and Control, 1:173–378, 2014.doi:10.1561/2600000002
2014 doi
-
[31]
T . Wolf, B. Lohmann, R. Eid, and P . Kotyczka. Passivity and structure preserving order reduction of linear port-Hamiltonian systems using Krylov subspaces. European Journal of Control, 16:401–406, 2010. doi: 10.3166/ejc.16.401-406. 14
2010 doi
Reviewed August 4, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.