Iterative Krylov Subspace Methods for Linear Port‐Hamiltonian Systems

Authors

Stefan Maier, Nicole Marheineke

Abstract

              In this work, we present a structure‐preserving Krylov subspace iteration scheme for solving the equation systems that arise from the Gauss integration of linear energy‐conserving and dissipative differential systems (e.g., Poisson systems, gradient systems, and port‐Hamiltonian systems). Exploiting the relation between Gauss integrators and diagonal Padé approximations, the ‐Arnoldi process yields iterates that are not only energy‐preserving up to convergence but also on each iteration level. We extend the approach to cover also energy dissipation. The use of the ‐Arnoldi approximation for Gauss integration enhances the computational efficiency as it allows a termination of the iteration without loss of energy‐associated structures as soon as the desired accuracy of the numerical integrator is reached. We investigate the performance in splitting schemes for linear port‐Hamiltonian systems.

Citation

• Journal: Proceedings in Applied Mathematics and Mechanics
• Year: 2026
• Volume: 26
• Issue: 2
• Pages:
• Publisher: Wiley
• DOI: 10.1002/pamm.70139

BibTeX

@article{Maier_2026,
  title={{Iterative Krylov Subspace Methods for Linear Port‐Hamiltonian Systems}},
  volume={26},
  ISSN={1617-7061},
  DOI={10.1002/pamm.70139},
  number={2},
  journal={Proceedings in Applied Mathematics and Mechanics},
  publisher={Wiley},
  author={Maier, Stefan and Marheineke, Nicole},
  year={2026}
}

Download the bib file

References

• van der Schaft A, Jeltsema D (2014) Port-Hamiltonian Systems Theory: An Introductory Overvie – 10.1561/9781601987877
• McLachlan RI (1995) On the Numerical Integration of Ordinary Differential Equations by Symmetric Composition Methods. SIAM J Sci Comput 16(1):151–168. https://doi.org/10.1137/091601 – 10.1137/0916010
• Liesen J, Strakos Z (2012) Krylov Subspace Method – 10.1093/acprof:oso/9780199655410.001.0001
• Sogabe T (2022) Krylov Subspace Methods for Linear Systems. Springer Nature Singapor – 10.1007/978-981-19-8532-4
• Frommer A., Progress in Industrial Mathematics at ECMI 2023 (2026)
• Mönch M, Marheineke N (2025) Commutator-based operator splitting for linear port-Hamiltonian systems. Applied Numerical Mathematics 210:25–38. https://doi.org/10.1016/j.apnum.2024.12.00 – 10.1016/j.apnum.2024.12.007
• Hairer E., Geometric Numerical Integration: Structure Preserving Algorithms for Ordinary Differential Equations (2006)
• Maier S, Marheineke N, Frommer A (2025) Energy-preserving iteration schemes for Gauss collocation integrators. Linear Algebra and its Applications. https://doi.org/10.1016/j.laa.2025.09.00 – 10.1016/j.laa.2025.09.005
• McLachlan RI, Quispel GRW (2002) Splitting methods. Acta Numerica 11:341–434. https://doi.org/10.1017/s096249290200005 – 10.1017/s0962492902000053
• Sexton JC, Weingarten DH (1992) Hamiltonian evolution for the hybrid Monte Carlo algorithm. Nuclear Physics B 380(3):665–677. https://doi.org/10.1016/0550-3213(92)90263- – 10.1016/0550-3213(92)90263-b
• Baker GA, Graves-Morris P (1996) Padé Approximants Second Editio – 10.1017/cbo9780511530074
• Gautschi W., Numerical Analysis: An Introduction (1997)
• Higham NJ (2008) Functions of Matrice – 10.1137/1.9780898717778
• Lancaster P., The Theory of Matrices: With Applications (1985)
• Hairer E., Solving Ordinary Differential Equations II. Stiff and Differential‐Algebraic Problems (1996)
• Frommer A, Simoncini V (2008) Matrix Functions. Mathematics in Industry 275–30 – 10.1007/978-3-540-78841-6_13
• Lopez L, Simoncini V (2006) Analysis of Projection Methods for Rational Function Approximation to the Matrix Exponential. SIAM J Numer Anal 44(2):613–635. https://doi.org/10.1137/0506259 – 10.1137/05062590
• Gugercin S, Polyuga RV, Beattie C, van der Schaft A (2012) Structure-preserving tangential interpolation for model reduction of port-Hamiltonian systems. Automatica 48(9):1963–1974. https://doi.org/10.1016/j.automatica.2012.05.05 – 10.1016/j.automatica.2012.05.052
• Al-Mohy AH, Higham NJ (2011) Computing the Action of the Matrix Exponential, with an Application to Exponential Integrators. SIAM J Sci Comput 33(2):488–511. https://doi.org/10.1137/10078886 – 10.1137/100788860