526 research outputs found
Fast and Efficient Numerical Methods for an Extended Black-Scholes Model
An efficient linear solver plays an important role while solving partial
differential equations (PDEs) and partial integro-differential equations
(PIDEs) type mathematical models. In most cases, the efficiency depends on the
stability and accuracy of the numerical scheme considered. In this article we
consider a PIDE that arises in option pricing theory (financial problems) as
well as in various scientific modeling and deal with two different topics. In
the first part of the article, we study several iterative techniques
(preconditioned) for the PIDE model. A wavelet basis and a Fourier sine basis
have been used to design various preconditioners to improve the convergence
criteria of iterative solvers. We implement a multigrid (MG) iterative method.
In fact, we approximate the problem using a finite difference scheme, then
implement a few preconditioned Krylov subspace methods as well as a MG method
to speed up the computation. Then, in the second part in this study, we analyze
the stability and the accuracy of two different one step schemes to approximate
the model.Comment: 29 pages; 10 figure
Wavelet Galerkin method for fractional elliptic differential equations
Under the guidance of the general theory developed for classical partial
differential equations (PDEs), we investigate the Riesz bases of wavelets in
the spaces where fractional PDEs usually work, and their applications in
numerically solving fractional elliptic differential equations (FEDEs). The
technique issues are solved and the detailed algorithm descriptions are
provided. Compared with the ordinary Galerkin methods, the wavelet Galerkin
method we propose for FEDEs has the striking benefit of efficiency, since the
condition numbers of the corresponding stiffness matrixes are small and
uniformly bounded; and the Toeplitz structure of the matrix still can be used
to reduce cost. Numerical results and comparison with the ordinary Galerkin
methods are presented to demonstrate the advantages of the wavelet Galerkin
method we provide.Comment: 20 pages, 0 figure
Mini-Workshop: Adaptive Methods for Control Problems Constrained by Time-Dependent PDEs
Optimization problems constrained by time-dependent PDEs (Partial Differential Equations) are challenging from a computational point of view: even in the simplest case, one needs to solve a system of PDEs coupled globally in time and space for the unknown solutions (the state, the costate and the control of the system). Typical and practically relevant examples are the control of nonlinear heat equations as they appear in laser hardening or the thermic control of flow problems (Boussinesq equations). Specifically for PDEs with a long time horizon, conventional time-stepping methods require an enormous storage of the respective other variables. In contrast, adaptive methods aim at distributing the available degrees of freedom in an a-posteriori-fashion to capture singularities and are, therefore, most promising
An odyssey into local refinement and multilevel preconditioning III: Implementation and numerical experiments
In this paper, we examine a number of additive and multiplicative multilevel iterative methods and preconditioners in the setting of two-dimensional local mesh refinement. While standard multilevel methods are effective for uniform refinement-based discretizations of elliptic equations, they tend to be less effective for algebraic systems, which arise from discretizations on locally refined meshes, losing their optimal behavior in both storage and computational complexity. Our primary focus here is on Bramble, Pasciak, and Xu (BPX)-style additive and multiplicative multilevel preconditioners, and on various stabilizations of the additive and multiplicative hierarchical basis (HB) method, and their use in the local mesh refinement setting. In parts I and II of this trilogy, it was shown that both BPX and wavelet stabilizations of HB have uniformly bounded condition numbers on several classes of locally refined two- and three-dimensional meshes based on fairly standard (and easily implementable) red and red-green mesh refinement algorithms. In this third part of the trilogy, we describe in detail the implementation of these types of algorithms, including detailed discussions of the data structures and traversal algorithms we employ for obtaining optimal storage and computational complexity in our implementations. We show how each of the algorithms can be implemented using standard data types, available in languages such as C and FORTRAN, so that the resulting algorithms have optimal (linear) storage requirements, and so that the resulting multilevel method or preconditioner can be applied with optimal (linear) computational costs. We have successfully used these data structure ideas for both MATLAB and C implementations using the FEtk, an open source finite element software package. We finish the paper with a sequence of numerical experiments illustrating the effectiveness of a number of BPX and stabilized HB variants for several examples requiring local refinement
- …