Papers Submitted
Abstract:
There are two important classes of numerical methods for stochastic differential equations (SDEs): strong methods and weak methods. Strong methods construct numerical approximations to trajectories of the SDEs directly, using the Brownian motions driving the SDEs. Weak methods compute a numerical trajectory using a sequence of random variables independent of the Brownian motions. The convergence of weak methods is usually expressed indirectly in terms of the convergence of expected values of test functions of the trajectories. Here we present an alternative formulation of convergence for weak methods in terms of the well-known Prokhorov metric on spaces of random variables. For a general class of weak methods, we establish bounds on the rates of convergence in terms of the Prokhorov metric. In doing so, we revisit the original proofs of convergence for weak methods and show explicitly how the bounds on the error depend on the smoothness of the test functions. As an application of our result, we use the Strassen--Dudley theorem to show that the true solution to the system of SDEs and the approximation from the weak method can be embedded in the same probability space in such a way that values generated by the weak method converge there in a strong sense. We conclude with a review of the existing results for pathwise convergence of weak methods and the corresponding strong results available under embedding.