This paper deals with the numerical solution of the Lax system L' = [B(L), L], L(0) = L-0 (*), where L-0 is a constant symmetric matrix, B(-) maps symmetric matrices into skew-symmetric matrices, and [B(L), L] is the commutator of B(L) and L. Here two different procedures, based on the approach recently proposed by Calvo, Iserles and Zanna (the MGLRK methods), are suggested. Such an approach is a computational form for the Flaschka formulation of (*). Our numerical procedures consist in solving (*) by a Runge-Kutta method, then, a single step of a Gauss-Legendre Runge-Kutta (GLRK) method may be applied to the Flaschka formulation of (*). In the first procedure we compute the approximation of the Lax system by a continuous explicit RK method, instead, the second procedure computes the approximation of the Lax system by a GLRK method (the same method used for the Flaschka system). The computational costs have been derived and compared with the ones of the MGLRK methods. Finally, several numerical tests and computational comparisons will be shown.

Numerical Procedures Based on Runge-Kutta Methods for Solving Isospectral Flows / Lopez, L.; Politi, T.. - In: APPLIED NUMERICAL MATHEMATICS. - ISSN 0168-9274. - STAMPA. - 25:4(1997), pp. 443-459. [10.1016/S0168-9274(97)00051-2]

Numerical Procedures Based on Runge-Kutta Methods for Solving Isospectral Flows

Politi, T.
1997-01-01

Abstract

This paper deals with the numerical solution of the Lax system L' = [B(L), L], L(0) = L-0 (*), where L-0 is a constant symmetric matrix, B(-) maps symmetric matrices into skew-symmetric matrices, and [B(L), L] is the commutator of B(L) and L. Here two different procedures, based on the approach recently proposed by Calvo, Iserles and Zanna (the MGLRK methods), are suggested. Such an approach is a computational form for the Flaschka formulation of (*). Our numerical procedures consist in solving (*) by a Runge-Kutta method, then, a single step of a Gauss-Legendre Runge-Kutta (GLRK) method may be applied to the Flaschka formulation of (*). In the first procedure we compute the approximation of the Lax system by a continuous explicit RK method, instead, the second procedure computes the approximation of the Lax system by a GLRK method (the same method used for the Flaschka system). The computational costs have been derived and compared with the ones of the MGLRK methods. Finally, several numerical tests and computational comparisons will be shown.
1997
Numerical Procedures Based on Runge-Kutta Methods for Solving Isospectral Flows / Lopez, L.; Politi, T.. - In: APPLIED NUMERICAL MATHEMATICS. - ISSN 0168-9274. - STAMPA. - 25:4(1997), pp. 443-459. [10.1016/S0168-9274(97)00051-2]
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/11589/10195
Citazioni
  • Scopus 7
  • ???jsp.display-item.citation.isi??? 9
social impact