English

On the complexity of solving initial value problems

Numerical Analysis 2017-01-18 v1 Computational Complexity

Abstract

In this paper we prove that computing the solution of an initial-value problem y˙=p(y)\dot{y}=p(y) with initial condition y(t0)=y0Rdy(t_0)=y_0\in\R^d at time t0+Tt_0+T with precision eμe^{-\mu} where pp is a vector of polynomials can be done in time polynomial in the value of TT, μ\mu and Y=supt0uT\infnormy(u)Y=\sup_{t_0\leqslant u\leqslant T}\infnorm{y(u)}. Contrary to existing results, our algorithm works for any vector of polynomials pp over any bounded or unbounded domain and has a guaranteed complexity and precision. In particular we do not assume pp to be fixed, nor the solution to lie in a compact domain, nor we assume that pp has a Lipschitz constant.

Keywords

Cite

@article{arxiv.1202.4407,
  title  = {On the complexity of solving initial value problems},
  author = {Olivier Bournez and Daniel S. Graça and Amaury Pouly},
  journal= {arXiv preprint arXiv:1202.4407},
  year   = {2017}
}

Comments

8 pages (two columns per page), submitted to ISSAC'12 conference

R2 v1 2026-06-21T20:22:21.801Z