Runge-Kutta type methods for orthogonal integration

From MaRDI portal





This paper is concerned with the numerical solution of matrix differential equations that possess orthogonal lows by means of onestep methods that preserve the orthogonality along the numerical solution. Such differential systems have the form \(Y'=F(Y)Y\), where \(F\) is a matrix function that maps orthogonal matrices to skew-symmetric matrices.NEWLINENEWLINENEWLINEThe author shows that the linearly-implicit modification of the implicit midpoint rule given by \(Y_{n+1}=Y_n+(h/2)F(Y_n)\) \((Y_n+Y_{n+1})\) preserves orthogonality and has order one, while the implicit midpoint rule does not preserve this property. A second-order linearly-implicit method which preserves orthogonality is also given.NEWLINENEWLINENEWLINEFinally, some numerical experiments are presented to show the effects of the lack of orthogonality of some numerical methods in the numerical solution of two selected nonlinear matrix problems.




Cited in
(21)








This page was built for publication: Runge-Kutta type methods for orthogonal integration

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5961742)