The lack of continuity and the role of infinite and infinitesimal in numerical methods for ODEs: The case of symplecticity
The lack of continuity and the role of infinite and infinitesimal in numerical methods for ODEs: The case of symplecticity
复制标题
DOI:
10.1016/j.amc.2011.03.022
复制
发表时间:
2010-10
期刊:
影响因子:
--
通讯作者:
L. Brugnano;F. Iavernaro;D. Trigiante
中科院分区:
文献类型:
--
作者:
L. Brugnano;F. Iavernaro;D. Trigiante
When numerically integrating canonical Hamiltonian systems, the long-term conservation of some of its invariants, for example the Hamiltonian function itself, assumes a central role. The classical approach to this problem has led to the definition of symplectic methods, among which we mention Gauss–Legendre collocation formulae. Indeed, in the continuous setting, energy conservation is derived from symplecticity via an infinite number of infinitesimal contact transformations. However, this infinite process cannot be directly transferred to the discrete setting. By following a different approach, in this paper we describe a sequence of methods, sharing the same essential spectrum (and, then, the same essential properties), which are energy preserving starting from a certain element of the sequence on, i.e., after a finite number of steps.