Krylov integrators for Hamiltonian systems
Abstract
We consider Arnoldi like processes to obtain symplectic subspaces for Hamiltonian systems. Large systems are locally approximated by ones living in low dimensional subspaces; we especially consider Krylov subspaces and some extensions. This will be utilized in two ways: solve numerically local small dimensional systems or in a given numerical, e.g. exponential, integrator, use the subspace for approximations of necessary functions. In the former case one can expect an excellent energy preservation. For the latter this is so for linear systems. For some second order exponential integrators we consider these two approaches are shown to be equivalent. In numerical experiments with nonlinear Hamiltonian problems their behaviour seems promising.
- Publication:
-
arXiv e-prints
- Pub Date:
- December 2017
- DOI:
- 10.48550/arXiv.1712.04047
- arXiv:
- arXiv:1712.04047
- Bibcode:
- 2017arXiv171204047K
- Keywords:
-
- Mathematics - Numerical Analysis;
- 65P10;
- 37M15;
- 65F60;
- 65F10
- E-Print:
- 24 pages, 10 figures