We focus our attention on the family of General Linear Methods (GLMs), for the numerical solution of second order ordinary differential equations (ODEs). These are multivalue methods introduced in (D'Ambrosio et al., 2012) [3] with the aim to provide an unifying approach for the analysis of the properties of accuracy of numerical methods for second order ODEs. Our investigation is addressed to providing the building blocks useful to analyze the linear stability properties of GLMs for second order ODEs: thus, we present the extension of the classical notions of stability matrix, stability polynomial, stability and periodicity interval, A-stability and P-stability to the family of GLMs. Special attention will be given to the derivation of highly stable GLMs, whose stability properties depend on the stability polynomial of indirect Runge-Kutta-Nyström methods based on Gauss-Legendre collocation points, which are known to be P-stable. In this way, we are able to provide P-stable GLMs whose order of convergence is greater than that of the corresponding RKN method, without heightening the computational cost. We finally provide and discuss examples of P-stable irreducible GLMs satisfying the mentioned features. © 2013 Elsevier B.V. All rights reserved.

P-stable general Nystrom methods for y''= f(y(t))

D'Ambrosio, Raffaele;
2014-01-01

Abstract

We focus our attention on the family of General Linear Methods (GLMs), for the numerical solution of second order ordinary differential equations (ODEs). These are multivalue methods introduced in (D'Ambrosio et al., 2012) [3] with the aim to provide an unifying approach for the analysis of the properties of accuracy of numerical methods for second order ODEs. Our investigation is addressed to providing the building blocks useful to analyze the linear stability properties of GLMs for second order ODEs: thus, we present the extension of the classical notions of stability matrix, stability polynomial, stability and periodicity interval, A-stability and P-stability to the family of GLMs. Special attention will be given to the derivation of highly stable GLMs, whose stability properties depend on the stability polynomial of indirect Runge-Kutta-Nyström methods based on Gauss-Legendre collocation points, which are known to be P-stable. In this way, we are able to provide P-stable GLMs whose order of convergence is greater than that of the corresponding RKN method, without heightening the computational cost. We finally provide and discuss examples of P-stable irreducible GLMs satisfying the mentioned features. © 2013 Elsevier B.V. All rights reserved.
File in questo prodotto:
Non ci sono file associati a questo prodotto.
Pubblicazioni consigliate

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/11697/120002
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 13
  • ???jsp.display-item.citation.isi??? 5
social impact