English

Low-Rank Solvers for Energy-Conserving Hamiltonian Boundary Value Methods

Numerical Analysis 2026-05-18 v2 Numerical Analysis

Abstract

We study energy-conserving Hamiltonian Boundary Value Methods (HBVMs) for Hamiltonian systems, which arise in applications where long-term preservation of energy and symplecticity is essential. HBVMs are multi-stage schemes whose stage equations reformulate as matrix equations with a low-rank right-hand side. For linear systems, we exploit this structure directly via Krylov projection solvers. For nonlinear systems, we leverage it within simplified Newton iterations and as a preconditioner in a Newton--Krylov framework, combined with adaptive time-stepping for robust convergence. Numerical experiments on semi-discretized wave equations demonstrate the efficiency and robustness of the proposed approach.

Keywords

Cite

@article{arxiv.2511.21597,
  title  = {Low-Rank Solvers for Energy-Conserving Hamiltonian Boundary Value Methods},
  author = {Fabio Durastante and Mariarosa Mazza},
  journal= {arXiv preprint arXiv:2511.21597},
  year   = {2026}
}
R2 v1 2026-07-01T07:56:37.091Z