ODE solvers using band-limited approximations

From MaRDI portal
Publication:348985

DOI10.1016/J.JCP.2014.02.001zbMATH Open1349.65692arXiv1208.3285OpenAlexW2076936268MaRDI QIDQ348985FDOQ348985

K. Sandberg, Gregory Beylkin

Publication date: 5 December 2016

Published in: Journal of Computational Physics (Search for Journal in Brave)

Abstract: We use generalized Gaussian quadratures for exponentials to develop a new ODE solver. Nodes and weights of these quadratures are computed for a given bandlimit c and user selected accuracy epsilon, so that they integrate functions eibx, for all |b|lec, with accuracy epsilon. Nodes of these quadratures do not concentrate excessively near the end points of an interval as those of the standard, polynomial-based Gaussian quadratures. Due to this property, the usual implicit Runge Kutta (IRK) collocation method may be used with a large number of nodes, as long as the method chosen for solving the nonlinear system of equations converges. We show that the resulting ODE solver is symplectic and demonstrate (numerically) that it is A-stable. We use this solver, dubbed Band-limited Collocation (BLC-IRK), in the problem of orbit determination. Since BLC-IRK minimizes the number of nodes needed to obtain the solution, in this problem we achieve speed close to that of explicit multistep methods.


Full work available at URL: https://arxiv.org/abs/1208.3285




Recommendations




Cites Work


Cited In (7)





This page was built for publication: ODE solvers using band-limited approximations

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