English

ODE Solvers Using Bandlimited Approximations

Numerical Analysis 2013-11-21 v2 Classical Analysis and ODEs

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 cc and user selected accuracy ϵ\epsilon, so that they integrate functions eibxe^{ibx}, for all bc|b|\le c, 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.

Keywords

Cite

@article{arxiv.1208.3285,
  title  = {ODE Solvers Using Bandlimited Approximations},
  author = {Gregory Beylkin and Kristian Sandberg},
  journal= {arXiv preprint arXiv:1208.3285},
  year   = {2013}
}

Comments

25 pages, 2 figures

R2 v1 2026-06-21T21:51:19.780Z