We consider the application of symmetric Boundary Value Methods to linear autonomous Hamiltonian systems. The numerical approximation of the Hamiltonian function exhibits a superconvergence property, namely its order of convergence is p + 2 for a p order symmetric method. We exploit this result to define a natural projection procedure that slightly modifies the numerical solution so that, without changing the convergence properties of the numerical method, it provides orbits lying on the same quadratic manifold as the continuous ones. A numerical test is also reported.

A symplectic Lanczos-type algorithm to compute the eigenvalues of positive definite Hamiltonian matrices

AMODIO, Pierluigi
2003-01-01

Abstract

We consider the application of symmetric Boundary Value Methods to linear autonomous Hamiltonian systems. The numerical approximation of the Hamiltonian function exhibits a superconvergence property, namely its order of convergence is p + 2 for a p order symmetric method. We exploit this result to define a natural projection procedure that slightly modifies the numerical solution so that, without changing the convergence properties of the numerical method, it provides orbits lying on the same quadratic manifold as the continuous ones. A numerical test is also reported.
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: https://hdl.handle.net/11586/38502
 Attenzione

Attenzione! I dati visualizzati non sono stati sottoposti a validazione da parte dell'ateneo

Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 5
  • ???jsp.display-item.citation.isi??? 5
social impact