Much effort is put into the construction of general linear methods with the aim of achieving an excellent long-time behavior for the integration of Hamiltonian systems. In this article, a backward error analysis is presented, which permits to get sharp estimates for the parasitic solution components and for the error in the Hamiltonian. For carefully constructed methods (symmetric and zero growth parameters) the error in the parasitic components typically grows like $h^{p+4}exp(h^2Lt)$, where p is the order of the method, and L depends on the problem and on the coefficients of the method. This is confirmed by numerical experiments.

Long-term stability of multi-value methods for ordinary differential equations

D'AMBROSIO, RAFFAELE;
2014

Abstract

Much effort is put into the construction of general linear methods with the aim of achieving an excellent long-time behavior for the integration of Hamiltonian systems. In this article, a backward error analysis is presented, which permits to get sharp estimates for the parasitic solution components and for the error in the Hamiltonian. For carefully constructed methods (symmetric and zero growth parameters) the error in the parasitic components typically grows like $h^{p+4}exp(h^2Lt)$, where p is the order of the method, and L depends on the problem and on the coefficients of the method. This is confirmed by numerical experiments.
File in questo prodotto:
Non ci sono file associati a questo prodotto.

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: http://hdl.handle.net/11386/4243253
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus ND
  • ???jsp.display-item.citation.isi??? ND
social impact