An efficient jet marcher for computing the quasipotential for 2D SDEs

09/08/2021
by   Nicholas Paskal, et al.
0

We present a new algorithm, the efficient jet marching method (EJM), for computing the quasipotential and its gradient for two-dimensional SDEs. The quasipotential is a potential-like function for nongradient SDEs that gives asymptotic estimates for the invariant probability measure, expected escape times from basins of attractors, and maximum likelihood escape paths. The quasipotential is a solution to an optimal control problem with an anisotropic cost function which can be solved for numerically via Dijkstra-like label-setting methods. Previous Dijkstra-like quasipotential solvers have displayed in general 1st order accuracy in the mesh spacing. However, by utilizing higher order interpolations of the quasipotential as well as more accurate approximations of the minimum action paths (MAPs), EJM achieves second-order accuracy for the quasipotential and nearly second-order for its gradient. Moreover, by using targeted search neighborhoods for the fastest characteristics following the ideas of Mirebeau, EJM also enjoys a reduction in computation time. This highly accurate solver enables us to compute the prefactor for the WKB approximation for the invariant probability measure and the Bouchet-Reygner sharp estimate for the expected escape time for the Maier-Stein SDE.

READ FULL TEXT

page 29

page 30

research
09/23/2021

Optimal Control via Combined Inference and Numerical Optimization

Derivative based optimization methods are efficient at solving optimal c...
research
05/24/2021

2nd-order Updates with 1st-order Complexity

It has long been a goal to efficiently compute and use second order info...
research
01/15/2021

Higher Order Targeted Maximum Likelihood Estimation

Asymptotic efficiency of targeted maximum likelihood estimators (TMLE) o...
research
08/12/2022

A second order accurate numerical method for the Poisson-Nernst-Planck system in the energetic variational formulation

A second order accurate (in time) numerical scheme is proposed and analy...
research
02/18/2022

Analysis and approximations of an optimal control problem for the Allen-Cahn equation

The scope of this paper is the analysis and approximation of an optimal ...
research
12/13/2022

An adaptive spectral method for oscillatory second-order linear ODEs with frequency-independent cost

We introduce an efficient numerical method for second order linear ODEs ...
research
05/02/2018

Computing tensor Z-eigenvectors with dynamical systems

We present a new framework for computing Z-eigenvectors of general tenso...

Please sign up or login with your details

Forgot password? Click here to reset