Efficient time stepping for numerical integration using reinforcement
  learning
  Many problems in science and engineering require an efficient numerical
approximation of integrals or solutions to differential equations. For systems
with rapidly changing dynamics, an equidistant discretization is often
inadvisable as it either results in prohibitively large errors or computational
effort. To this end, adaptive schemes, such as solvers based on Runge--Kutta
pairs, have been developed which adapt the step size based on local error
estimations at each step. While the classical schemes apply very generally and
are highly efficient on regular systems, they can behave sub-optimal when an
inefficient step rejection mechanism is triggered by structurally complex
systems such as chaotic systems. To overcome these issues, we propose a method
to tailor numerical schemes to the problem class at hand. This is achieved by
combining simple, classical quadrature rules or ODE solvers with data-driven
time-stepping controllers. Compared with learning solution operators to ODEs
directly, it generalises better to unseen initial data as our approach employs
classical numerical schemes as base methods. At the same time it can make use
of identified structures of a problem class and, therefore, outperforms
state-of-the-art adaptive schemes. Several examples demonstrate superior
efficiency. Source code is available at
https://github.com/lueckem/quadrature-ML.
