Many applications in contemporary science involve multiscale dynamics with
time and space scale separation of patterns of motion, with fewer slowly
evolving variables and much larger set of faster evolving variables. A direct
numerical simulation of the evolution of such dynamics is typically
computationally expensive, due to both the large number of fast variables and
necessity for a small discretization time step to resolve fast components of
dynamics.