3,505 research outputs found
The automatic solution of partial differential equations using a global spectral method
A spectral method for solving linear partial differential equations (PDEs)
with variable coefficients and general boundary conditions defined on
rectangular domains is described, based on separable representations of partial
differential operators and the one-dimensional ultraspherical spectral method.
If a partial differential operator is of splitting rank , such as the
operator associated with Poisson or Helmholtz, the corresponding PDE is solved
via a generalized Sylvester matrix equation, and a bivariate polynomial
approximation of the solution of degree is computed in
operations. Partial differential operators of
splitting rank are solved via a linear system involving a block-banded
matrix in operations. Numerical
examples demonstrate the applicability of our 2D spectral method to a broad
class of PDEs, which includes elliptic and dispersive time-evolution equations.
The resulting PDE solver is written in MATLAB and is publicly available as part
of CHEBFUN. It can resolve solutions requiring over a million degrees of
freedom in under seconds. An experimental implementation in the Julia
language can currently perform the same solve in seconds.Comment: 22 page
Non-negativity preserving numerical algorithms for stochastic differential equations
Construction of splitting-step methods and properties of related
non-negativity and boundary preserving numerical algorithms for solving
stochastic differential equations (SDEs) of Ito-type are discussed. We present
convergence proofs for a newly designed splitting-step algorithm and simulation
studies for numerous numerical examples ranging from stochastic dynamics
occurring in asset pricing theory in mathematical finance (SDEs of CIR and CEV
models) to measure-valued diffusion and superBrownian motion (SPDEs) as met in
biology and physics.Comment: 23 pages, 7 figures. Figures 6.2 and 6.3 in low resolution due to
upload size restrictions. Original resolution at
http://gisc.uc3m.es/~moro/profesional.htm
A Parallel Algorithm for solving BSDEs - Application to the pricing and hedging of American options
We present a parallel algorithm for solving backward stochastic differential
equations (BSDEs in short) which are very useful theoretic tools to deal with
many financial problems ranging from option pricing option to risk management.
Our algorithm based on Gobet and Labart (2010) exploits the link between BSDEs
and non linear partial differential equations (PDEs in short) and hence enables
to solve high dimensional non linear PDEs. In this work, we apply it to the
pricing and hedging of American options in high dimensional local volatility
models, which remains very computationally demanding. We have tested our
algorithm up to dimension 10 on a cluster of 512 CPUs and we obtained linear
speedups which proves the scalability of our implementationComment: 25 page
- …