AMF-type W-methods for parabolic problems with mixed derivatives
Severiano González-Pinto, Ernst Hairer, Domingo Hernández-Abreu
and Soledad Pérez-Rodríguez,
Abstract. The time integration of differential equations obtained by the
space discretization via Finite Differences of evolution parabolic PDEs with mixed
derivatives in the elliptic operator is considered (MOL approach).
W-methods (Rosenbrock-type methods) are combined with
the Approximate Matrix Factorization technique (AMF), which is considered in
alternating direction implicit (ADI) sense.
The focus of the paper is a stability analysis, which is based on a scalar test
equation that is relevant for the class of problems when periodic
or homogeneous Dirichlet boundary conditions are considered.
Unconditional stability, independent of the number of space dimensions $m$,
is proved for a variety of AMF-type W-methods.
Numerical experiments with linear parabolic problems in dimension m=3 and m=4,
as well as with the Heston problem from financial option pricing are presented.
Key Words. Parabolic PDEs, mixed derivatives, time integration, W-methods,
Approximate Matrix Factorization, Alternating Direction Implicit schemes,
unconditional stability, Heston model.