Abstract:This work develops Monte Carlo Euler adaptive time stepping methods for the weak approximation problem of jump diffusion driven stochastic differential equations. The main result is the derivation of a new expansion for the omputational error, with computable leading order term in a posteriori form, based on stochastic flows and discrete dual backward problems which extends the results in [STZ]. These expansions lead to efficient and accurate computation of error estimates. Adaptive algorithms for either stochastic time steps or quasi-deterministic time steps are described. Numerical examples show the performance of the proposed error approximation and of the described adaptive time-stepping methods.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to develop an effective adaptive time - step Monte Carlo Euler method for weakly approximating stochastic differential equations (SDEs) driven by jump - diffusions. Specifically, the paper focuses on how to accurately estimate the computational error and adjust the time - step based on these estimates to improve the accuracy and efficiency of the numerical solution. ### Specific description of the problem 1. **Model setting**: The paper considers a d - dimensional stochastic process \(X(t)=(X_1(t), \ldots, X_d(t))\) which satisfies the following stochastic differential equation (SDE): \[ X(t)=X(0)+\int_{0}^{t} a(s, X(s^{-})) \, ds+\sum_{\ell = 1}^{\ell_0}\int_{0}^{t} b_{\ell}(s, X(s^{-})) \, dW_{\ell}(s)+\int_{0}^{t}\int_{Z} c(s, X(s^{-}), z) \, p(ds, dz) \] where: - \(a(s, X(s^{-}))\) is the drift term. - \(b_{\ell}(s, X(s^{-}))\) is the diffusion term related to the Wiener process. - \(c(s, X(s^{-}), z)\) is the jump term related to the Poisson random measure. 2. **Objective**: For a given scalar function \(g: \mathbb{R}^d\rightarrow\mathbb{R}\), the objective is to construct an approximation of \(E[g(X(T))]\), that is, the expected value at time \(T\). To this end, the paper adopts the Monte Carlo Euler method for numerical simulation. 3. **Main challenges**: - **Error estimation**: It is necessary to accurately estimate the time - discretization error \(E_T\) and the statistical error \(E_S\), and ensure that the total computational error \(E_C\) satisfies the given tolerance \(\text{TOL}\). - **Adaptive time - step**: Through a posteriori error estimation, dynamically adjust the time - step to minimize the computational workload while ensuring that the error is within an acceptable range. ### Solution The main contributions of the paper are: - **New error expansion formula**: Derived a new expansion of the computational error, which contains computable leading terms, based on the stochastic flow and the discrete dual backward problem. - **Adaptive algorithms**: Proposed two adaptive time - step algorithms, namely "quasi - deterministic" time and "stochastic" time - step, which are respectively suitable for different types of stochastic processes. - **Efficient computation**: By introducing the dual function, avoided directly solving the complex forward evolution equation, thereby reducing the computational complexity. ### Conclusion The paper verifies the effectiveness and accuracy of the proposed adaptive Monte Carlo Euler method through theoretical analysis and numerical experiments. This method not only improves the computational efficiency of weakly approximating SDEs driven by jump - diffusions but also performs excellently in error control.

Adaptive Weak Approximation of Diffusions with Jumps

Variable-Step Euler–Maruyama Approximations of Regime-Switching Jump Diffusion Processes

Strong convergence of a class of adaptive numerical methods for SDEs with jumps

Lower error bounds and optimality of approximation for jump-diffusion SDEs with discontinuous drift

A higher order approximation method for jump-diffusion SDEs with discontinuous drift coefficient

On the convergence of adaptive approximations for stochastic differential equations

Numerical Solutions of Jump Diffusions with Markovian Switching

Adaptive Mesh Construction for the Numerical Solution of Stochastic Differential Equations with Markovian Switching

Higher-order adaptive methods for exit times of Itô diffusions

Randomized Milstein algorithm for approximation of solutions of jump-diffusion SDEs

Strong convergence of an adaptive time-stepping Milstein method for SDEs with monotone coefficients

Strong Approximation of Monotone Stochastic Partial Differential Equations Driven by Multiplicative Noise

High order weak approximation of Stochastic Differential Equations for bounded and measurable test functions

Weak approximation schemes for SDEs with super-linearly growing coefficients

On the Rate of Convergence of Weak Euler Approximation for Nondegenerate Itô Diffusion and Jump Processes

Weak convergence rates for temporal numerical approximations of stochastic wave equations with multiplicative noise

Tamed-adaptive Euler-Maruyama approximation for SDEs with superlinearly growing and piecewise continuous drift, superlinearly growing and locally Hölder continuous diffusion

Multilevel Monte Carlo method for jump-diffusion SDEs

An adaptive strong order 1 method for SDEs with discontinuous drift coefficient

On the $\varepsilon$-Euler-Maruyama scheme for time inhomogeneous jump-driven SDEs

Strongly asymptotically optimal schemes for the strong approximation of stochastic differential equations with respect to the supremum error