37,303 research outputs found
On the ADI method for the Sylvester Equation and the optimal- points
The ADI iteration is closely related to the rational Krylov projection
methods for constructing low rank approximations to the solution of Sylvester
equation. In this paper we show that the ADI and rational Krylov approximations
are in fact equivalent when a special choice of shifts are employed in both
methods. We will call these shifts pseudo H2-optimal shifts. These shifts are
also optimal in the sense that for the Lyapunov equation, they yield a residual
which is orthogonal to the rational Krylov projection subspace. Via several
examples, we show that the pseudo H2-optimal shifts consistently yield nearly
optimal low rank approximations to the solutions of the Lyapunov equations
An ADI extrapolated Crank-Nicolson orthogonal spline collocation method for nonlinear reaction-diffusion systems: a computational study
An alternating direction implicit (ADI) orthogonal spline collocation (OSC)
method is described for the approximate solution of a class of nonlinear
reaction-diffusion systems. Its efficacy is demonstrated on the solution of
well-known examples of such systems, specifically the Brusselator, Gray-Scott,
Gierer-Meinhardt and Schnakenberg models, and comparisons are made with other
numerical techniques considered in the literature. The new ADI method is based
on an extrapolated Crank-Nicolson OSC method and is algebraically linear. It is
efficient, requiring at each time level only operations where
is the number of unknowns. Moreover,it is shown to produce
approximations which are of optimal global accuracy in various norms, and to
possess superconvergence properties
A Riemannian Trust Region Method for the Canonical Tensor Rank Approximation Problem
The canonical tensor rank approximation problem (TAP) consists of
approximating a real-valued tensor by one of low canonical rank, which is a
challenging non-linear, non-convex, constrained optimization problem, where the
constraint set forms a non-smooth semi-algebraic set. We introduce a Riemannian
Gauss-Newton method with trust region for solving small-scale, dense TAPs. The
novelty of our approach is threefold. First, we parametrize the constraint set
as the Cartesian product of Segre manifolds, hereby formulating the TAP as a
Riemannian optimization problem, and we argue why this parametrization is among
the theoretically best possible. Second, an original ST-HOSVD-based retraction
operator is proposed. Third, we introduce a hot restart mechanism that
efficiently detects when the optimization process is tending to an
ill-conditioned tensor rank decomposition and which often yields a quick escape
path from such spurious decompositions. Numerical experiments show improvements
of up to three orders of magnitude in terms of the expected time to compute a
successful solution over existing state-of-the-art methods
Application of Operator Splitting Methods in Finance
Financial derivatives pricing aims to find the fair value of a financial
contract on an underlying asset. Here we consider option pricing in the partial
differential equations framework. The contemporary models lead to
one-dimensional or multidimensional parabolic problems of the
convection-diffusion type and generalizations thereof. An overview of various
operator splitting methods is presented for the efficient numerical solution of
these problems.
Splitting schemes of the Alternating Direction Implicit (ADI) type are
discussed for multidimensional problems, e.g. given by stochastic volatility
(SV) models. For jump models Implicit-Explicit (IMEX) methods are considered
which efficiently treat the nonlocal jump operator. For American options an
easy-to-implement operator splitting method is described for the resulting
linear complementarity problems.
Numerical experiments are presented to illustrate the actual stability and
convergence of the splitting schemes. Here European and American put options
are considered under four asset price models: the classical Black-Scholes
model, the Merton jump-diffusion model, the Heston SV model, and the Bates SV
model with jumps
ADI splitting schemes for a fourth-order nonlinear partial differential equation from image processing
We present directional operator splitting schemes for the numerical solution of a fourth-order, nonlinear partial differential evolution equation which arises in image processing. This equation constitutes the H−1-gradient flow of the total variation and represents a prototype of higher-order equations of similar type which are popular in imaging for denoising, deblurring and inpainting problems. The efficient numerical solution of this equation is very challenging due to the stiffness of most numerical schemes. We show that the combination of directional splitting schemes with implicit time-stepping provides a stable and computationally cheap numerical realisation of the equation
An alternating direction implicit spectral method for solving two dimensional multi-term time fractional mixed diffusion and diffusion-wave equations
In this paper, we consider the initial boundary value problem of the two
dimensional multi-term time fractional mixed diffusion and diffusion-wave
equations. An alternating direction implicit (ADI) spectral method is developed
based on Legendre spectral approximation in space and finite difference
discretization in time. Numerical stability and convergence of the schemes are
proved, the optimal error is , where are the
polynomial degree, time step size and the regularity of the exact solution,
respectively. We also consider the non-smooth solution case by adding some
correction terms. Numerical experiments are presented to confirm our
theoretical analysis. These techniques can be used to model diffusion and
transport of viscoelastic non-Newtonian fluids
- …