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

by   Nicholas Paskal, et al.

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.


page 29

page 30


Optimal Control via Combined Inference and Numerical Optimization

Derivative based optimization methods are efficient at solving optimal c...

2nd-order Updates with 1st-order Complexity

It has long been a goal to efficiently compute and use second order info...

Higher Order Targeted Maximum Likelihood Estimation

Asymptotic efficiency of targeted maximum likelihood estimators (TMLE) o...

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 ...

The anatomy of Boris type solvers and large time-step plasma simulations

This work gives a Lie operator derivation of various Boris solvers via a...

Euclidean distance and maximum likelihood retractions by homotopy continuation

We define a new second-order retraction map for statistical models. We a...

Second-Order Occlusion-Aware Volumetric Radiance Caching

We present a second-order gradient analysis of light transport in partic...