168 research outputs found
Random walk algorithm for the Dirichlet problem for parabolic integro-differential equation
We consider stochastic differential equations driven by a general Lévy 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. Mathematics Subject Classification (2010) 65C30 · 60H10 · 35R09 · 60H35 · 60J7
Numerical Methods for the Fractional Laplacian: a Finite Difference-quadrature Approach
The fractional Laplacian is a non-local operator which
depends on the parameter and recovers the usual Laplacian as . A numerical method for the fractional Laplacian is proposed, based on
the singular integral representation for the operator. The method combines
finite difference with numerical quadrature, to obtain a discrete convolution
operator with positive weights. The accuracy of the method is shown to be
. Convergence of the method is proven. The treatment of far
field boundary conditions using an asymptotic approximation to the integral is
used to obtain an accurate method. Numerical experiments on known exact
solutions validate the predicted convergence rates. Computational examples
include exponentially and algebraically decaying solution with varying
regularity. The generalization to nonlinear equations involving the operator is
discussed: the obstacle problem for the fractional Laplacian is computed.Comment: 29 pages, 9 figure
Monte Carlo method for parabolic equations involving fractional Laplacian
We apply the Monte Carlo method to solving the Dirichlet problem of linear
parabolic equations with fractional Laplacian. This method exploit- s the idea
of weak approximation of related stochastic differential equations driven by
the symmetric stable L\'evy process with jumps. We utilize the jump- adapted
scheme to approximate L\'evy process which gives exact exit time to the
boundary. When the solution has low regularity, we establish a numeri- cal
scheme by removing the small jumps of the L\'evy process and then show the
convergence order. When the solution has higher regularity, we build up a
higher-order numerical scheme by replacing small jumps with a simple process
and then display the higher convergence order. Finally, numerical experiments
including ten- and one hundred-dimensional cases are presented, which confirm
the theoretical estimates and show the numerical efficiency of the proposed
schemes for high dimensional parabolic equations.Comment: 30pages 2 figure
Superconvergence of a discontinuous Galerkin method for fractional diffusion and wave equations
We consider an initial-boundary value problem for
, that is, for a fractional
diffusion () or wave () equation. A numerical solution
is found by applying a piecewise-linear, discontinuous Galerkin method in time
combined with a piecewise-linear, conforming finite element method in space.
The time mesh is graded appropriately near , but the spatial mesh is
quasiuniform. Previously, we proved that the error, measured in the spatial
-norm, is of order , uniformly in , where
is the maximum time step, is the maximum diameter of the spatial finite
elements, and . Here,
we generalize a known result for the classical heat equation (i.e., the case
) by showing that at each time level the solution is
superconvergent with respect to : the error is of order
. Moreover, a simple postprocessing step
employing Lagrange interpolation yields a superconvergent approximation for any
. Numerical experiments indicate that our theoretical error bound is
pessimistic if . Ignoring logarithmic factors, we observe that the
error in the DG solution at , and after postprocessing at all , is of
order .Comment: 24 pages, 2 figure
Recommended from our members
Dirichlet Form Theory and its Applications
Theory of Dirichlet forms is one of the main achievements in modern probability theory. It provides a powerful connection between probabilistic and analytic potential theory. It is also an effective machinery for studying various stochastic models, especially those with non-smooth data, on fractal-like spaces or spaces of infinite dimensions. The Dirichlet form theory has numerous interactions with other areas of mathematics and sciences.
This workshop brought together top experts in Dirichlet form theory and related fields as well as promising young researchers, with the common theme of developing new foundational methods and their applications to specific areas of probability. It provided a unique opportunity for the interaction between the established scholars and young researchers
A Probabilistic Numerical Method for Fully Nonlinear Parabolic PDEs
We consider the probabilistic numerical scheme for fully nonlinear PDEs
suggested in \cite{cstv}, and show that it can be introduced naturally as a
combination of Monte Carlo and finite differences scheme without appealing to
the theory of backward stochastic differential equations. Our first main result
provides the convergence of the discrete-time approximation and derives a bound
on the discretization error in terms of the time step. An explicit
implementable scheme requires to approximate the conditional expectation
operators involved in the discretization. This induces a further Monte Carlo
error. Our second main result is to prove the convergence of the latter
approximation scheme, and to derive an upper bound on the approximation error.
Numerical experiments are performed for the approximation of the solution of
the mean curvature flow equation in dimensions two and three, and for two and
five-dimensional (plus time) fully-nonlinear Hamilton-Jacobi-Bellman equations
arising in the theory of portfolio optimization in financial mathematics
- …