On Finite difference schemes for partial integro-differential equations of L\'evy type

TL;DR

This paper introduces a finite difference scheme for Le9vy-type integro-differential equations, achieving improved error estimates by treating the operator as second order on the entire domain without restrictions on the Le9vy measure.

Contribution

It presents a novel finite difference approximation that handles the entire Le9vy operator as second order, removing previous limitations and improving error bounds.

Findings

Achieves error estimates of order (h + τ^k) with k in {1/2, 1}.

Eliminates additional errors present in previous methods.

No conditions required on the Le9vy measure.

Abstract

In this article we introduce a finite difference approximation for integro-differential operators of L\'evy type. We approximate solutions of integro-differential equations, where the second order operator is allowed to degenerate. In the existing literature, the L\'evy operator is treated as a zero/first order operator outside of a centered ball of radius $\delta$, leading to error estimates of order $\xi (\delta)+N(\delta)(h+\sqrt{\tau})$, where $h$ and $\tau$ are the spatial and temporal discretization parameters respectively. In these estimates $\xi (\delta) \downarrow 0$, but $N(\delta )\uparrow \infty$ as $\delta \downarrow 0$. In contrast, we treat the integro-differential operator as a second order operator on the whole unit ball. By this method we obtain error estimates of order $(h+\tau^k)$ for $k\in \{1/2,1\}$, eliminating the additional errors and the blowing up constants. Moreover, we do not pose any conditions on the L\'evy measure.