New numerical integrators based on solvability and splitting
A new family of integration schemes for the linear matrix differential equation \(Y'(t) = A(t) Y(t)\), \(Y(0) = I\), is presented. It is based on representing the approximate solutions as a product of upper and lower triangular matrices. These triangular matrices are expressed in terms of multivariate integrals. A brief theoretical analysis of the algorithm is given. In particular it turns out that the approximations evolve in SL\((n)\) if the trace of the coefficient matrix \(A(t)\) vanishes. Thus it is possible to interpret the method as a new type of volume-preserving algorithm. Moreover the author discusses suitable ways to choose the numerical integration rules for the above mentioned multivariate integrals.
- High order optimized geometric integrators for linear differential equations
- scientific article; zbMATH DE number 3136073 (Why is no real title available?)
- scientific article; zbMATH DE number 1745051 (Why is no real title available?)
- scientific article; zbMATH DE number 1795895 (Why is no real title available?)
- scientific article; zbMATH DE number 3200675 (Why is no real title available?)
- Improved high order integrators based on the Magnus expansion
- Magnus and Fer expansions for matrix differential equations: the convergence problem
- Nineteen Dubious Ways to Compute the Exponential of a Matrix, Twenty-Five Years Later
- On the exponential solution of differential equations for a linear operator
- On the global error of discretization methods for highly-oscillatory ordinary differential equations
- On the solution of linear differential equations in Lie groups
- Solving Ordinary Differential Equations I
- Volume-preserving algorithms for source-free dynamical systems
This page was built for publication: New numerical integrators based on solvability and splitting
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q885931)