An hp-version discontinuous Galerkin method for integro-differential equations of parabolic type

We study the numerical solution of a class of parabolic integro-differential equations with weakly singular kernels. We use an $hp$-version discontinuous Galerkin (DG) method for the discretization in time. We derive optimal $hp$-version error estimates and show that exponential rates of convergence can be achieved for solutions with singular (temporal) behavior near $t=0$ caused by the weakly singular kernel. Moreover, we prove that by using nonuniformly refined time steps, optimal algebraic convergence rates can be achieved for the $h$-version DG method. We then combine the DG time-stepping method with a standard finite element discretization in space, and present an optimal error analysis of the resulting fully discrete scheme. Our theoretical results are numerically validated in a series of test problems

The obstacle problem for semilinear parabolic partial integro-differential equations

This paper presents a probabilistic interpretation for the weak Sobolev solution of the obstacle problem for semilinear parabolic partial integro-differential equations (PIDEs). The results of Leandre (1985) concerning the homeomorphic property for the solution of SDEs with jumps are used to construct random test functions for the variational equation for such PIDEs. This results in the natural connection with the associated Reflected Backward Stochastic Differential Equations with jumps (RBSDEs), namely Feynman Kac's formula for the solution of the PIDEs. Moreover it gives an application to the pricing and hedging of contingent claims with constraints in the wealth or portfolio processes in financial markets including jumps.Comment: 31 page

Regularity of the Optimal Stopping Problem for Jump Diffusions

The value function of an optimal stopping problem for jump diffusions is known to be a generalized solution of a variational inequality. Assuming that the diffusion component of the process is nondegenerate and a mild assumption on the singularity of the L\'{e}vy measure, this paper shows that the value function of this optimal stopping problem on an unbounded domain with finite/infinite variation jumps is in $W^{2,1}_{p, loc}$ with $p\in(1, \infty)$. As a consequence, the smooth-fit property holds.Comment: To Appear in the SIAM Journal on Control and Optimizatio

Variational Theory and Domain Decomposition for Nonlocal Problems

In this article we present the first results on domain decomposition methods for nonlocal operators. We present a nonlocal variational formulation for these operators and establish the well-posedness of associated boundary value problems, proving a nonlocal Poincar\'{e} inequality. To determine the conditioning of the discretized operator, we prove a spectral equivalence which leads to a mesh size independent upper bound for the condition number of the stiffness matrix. We then introduce a nonlocal two-domain variational formulation utilizing nonlocal transmission conditions, and prove equivalence with the single-domain formulation. A nonlocal Schur complement is introduced. We establish condition number bounds for the nonlocal stiffness and Schur complement matrices. Supporting numerical experiments demonstrating the conditioning of the nonlocal one- and two-domain problems are presented.Comment: Updated the technical part. In press in Applied Mathematics and Computatio

Optimal error estimates of a mixed finite element method for\ud parabolic integro-differential equations with non smooth initial data

In this article, a new mixed method is proposed and analyzed for parabolic integro-differential equations (PIDE) with nonsmooth initial data. Compared to mixed methods for PIDE, the present method does not bank on a reformulation using a resolvent operator. Based on energy arguments and without using parabolic type duality technique, optimal L2-error estimates are derived for semidiscrete approximations, when the initial data is in L2. Due to the presence of the integral term, it is, further, observed that estimate in dual of H(div)-space plays a role in our error analysis. Moreover, the proposed analysis follows the spirit of the proof technique used for deriving optimal error estimates of finite element approximations to PIDE with smooth data and therefore, it unifies both the theories, i.e., one for smooth data and other for nonsmooth data. Finally, the proposed analysis can be easily extended to other mixed method for PIDE with rough initial data and provides an improved result