ترغب بنشر مسار تعليمي؟ اضغط هنا

Random walk algorithm for the Dirichlet problem for parabolic integro-differential equation

74   0   0.0 ( 0 )
 نشر من قبل Michael Tretyakov
 تاريخ النشر 2020
  مجال البحث الهندسة المعلوماتية
والبحث باللغة English




اسأل ChatGPT حول البحث

We consider stochastic differential equations driven by a general Levy processes (SDEs) with infinite activity and the related, via the Feynman-Kac formula, Dirichlet problem for parabolic integro-differential equation (PIDE). We approximate the solution of PIDE using a numerical method for the SDEs. The method is based on three ingredients: (i) we approximate small jumps by a diffusion; (ii) we use restricted jump-adaptive time-stepping; and (iii) between the jumps we exploit a weak Euler approximation. We prove weak convergence of the considered algorithm and present an in-depth analysis of how its error and computational cost depend on the jump activity level. Results of some numerical experiments, including pricing of barrier basket currency options, are presented.



قيم البحث

اقرأ أيضاً

This paper studies the solvability of a class of Dirichlet problem associated with non-linear integro-differential operator. The main ingredient is the probabilistic construction of continuous supersolution via the identification of the continuity se t of the exit time operators under Skorohod topology.
A simple-to-implement weak-sense numerical method to approximate reflected stochastic differential equations (RSDEs) is proposed and analysed. It is proved that the method has the first order of weak convergence. Together with the Monte Carlo techniq ue, it can be used to numerically solve linear parabolic and elliptic PDEs with Robin boundary condition. One of the key results of this paper is the use of the proposed method for computing ergodic limits, i.e. expectations with respect to the invariant law of RSDEs, both inside a domain in $mathbb{R}^{d}$ and on its boundary. This allows to efficiently sample from distributions with compact support. Both time-averaging and ensemble-averaging estimators are considered and analysed. A number of extensions are considered including a second-order weak approximation, the case of arbitrary oblique direction of reflection, and a new adaptive weak scheme to solve a Poisson PDE with Neumann boundary condition. The presented theoretical results are supported by several numerical experiments.
A convexification-based numerical method for a Coefficient Inverse Problem for a parabolic PDE is presented. The key element of this method is the presence of the so-called Carleman Weight Function in the numerical scheme. Convergence analysis ensure s the global convergence of this method, as opposed to the local convergence of the conventional least squares minimization techniques. Numerical results demonstrate a good performance.
In this paper, several two-grid finite element algorithms for solving parabolic integro-differential equations (PIDEs) with nonlinear memory are presented. Analysis of these algorithms is given assuming a fully implicit time discretization. It is sho wn that these algorithms are as stable as the standard fully discrete finite element algorithm, and can achieve the same accuracy as the standard algorithm if the coarse grid size $H$ and the fine grid size $h$ satisfy $H=O(h^{frac{r-1}{r}})$. Especially for PIDEs with nonlinear memory defined by a lower order nonlinear operator, our two-grid algorithm can save significant storage and computing time. Numerical experiments are given to confirm the theoretical results.
201 - Jing Sun , Weihua Deng , Daxin Nie 2021
We make the split of the integral fractional Laplacian as $(-Delta)^s u=(-Delta)(-Delta)^{s-1}u$, where $sin(0,frac{1}{2})cup(frac{1}{2},1)$. Based on this splitting, we respectively discretize the one- and two-dimensional integral fractional Laplaci an with the inhomogeneous Dirichlet boundary condition and give the corresponding truncation errors with the help of the interpolation estimate. Moreover, the suitable corrections are proposed to guarantee the convergence in solving the inhomogeneous fractional Dirichlet problem and an $mathcal{O}(h^{1+alpha-2s})$ convergence rate is obtained when the solution $uin C^{1,alpha}(bar{Omega}^{delta}_{n})$, where $n$ is the dimension of the space, $alphain(max(0,2s-1),1]$, $delta$ is a fixed positive constant, and $h$ denotes mesh size. Finally, the performed numerical experiments confirm the theoretical results.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا