DeepAI AI Chat
Log In Sign Up

A general framework for the rigorous computation of invariant densities and the coarse-fine strategy

by   Stefano Galatolo, et al.

In this paper we present a general, axiomatical framework for the rigorous approximation of invariant densities and other important statistical features of dynamics. We approximate the system trough a finite element reduction, by composing the associated transfer operator with a suitable finite dimensional projection (a discretization scheme) as in the well-known Ulam method. We introduce a general framework based on a list of properties (of the system and of the projection) that need to be verified so that we can take advantage of a so-called “coarse-fine” strategy. This strategy is a novel method in which we exploit information coming from a coarser approximation of the system to get useful information on a finer approximation, speeding up the computation. This coarse-fine strategy allows a precise estimation of invariant densities and also allows to estimate rigorously the speed of mixing of the system by the speed of mixing of a coarse approximation of it, which can easily be estimated by the computer. The estimates obtained here are rigourous, i.e., they come with exact error bounds that are guaranteed to hold and take into account both the discretiazation and the approximations induced by finite-precision arithmetic. We apply this framework to several discretization schemes and examples of invariant density computation from previous works, obtaining a remarkable reduction in computation time. We have implemented the numerical methods described here in the Julia programming language, and released our implementation publicly as a Julia package.


page 1

page 2

page 3

page 4


DNN Approximation of Nonlinear Finite Element Equations

We investigate the potential of applying (D)NN ((deep) neural networks) ...

Projection error-based guaranteed L2 error bounds for finite element approximations of Laplace eigenfunctions

For conforming finite element approximations of the Laplacian eigenfunct...

An efficient two level approach for simulating Bose-Einstein condensates

In this work, we consider the numerical computation of ground states and...

Error bounds for Lanczos-based matrix function approximation

We analyze the Lanczos method for matrix function approximation (Lanczos...

Another Approximation of the First-Passage Time Densities for the Ratcliff Diffusion Decision Model

We present a novel method for approximating the probability density func...

Advanced Modeling of Rectangular Waveguide Devices with Smooth Profiles by Hierarchical Model Reduction

We present a new method for the analysis of smoothly varying tapers, tra...