Abstract
Given a stochastic differential equation (SDE) inRnwhose solution is constrained to lie in some manifoldM⊂Rn, we identify a class of numerical schemes for the SDE whose iterates remain close toMto high order. These schemes approximate a geometrically invariant scheme, which gives perfect solutions for any SDE that is diffeomorphic ton-dimensional Brownian motion. Unlike projection-based methods, they may be implemented without explicit knowledge ofM. They can even be implemented if the solution merely remains close toM, without being exactly confined to it. Our approach does not require simulating any iterated Itô integrals beyond those needed to implement the Euler–Maryuama scheme. We prove that the schemes converge under a standard set of assumptions, and illustrate their geometric advantages in a variety of numerical contexts, including Monte Carlo simulation of the Riemannian Langevin equation.
Sign-in/Register to access full text options 
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a callMore From: Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.
