Search

Recently, regime-switching option pricing based on fractional diffusion models has been used, which explains many significant empirical facts about financial markets better. There are many methods to solve the problem, but to the best of our knowledge, effective preconditioners for the second-order schemes have not been proposed. Thus, in this article, an implicit numerical scheme is developed for a regime-switching European option pricing problem under a multi-state tempered fractional model. The scheme is proven to be unconditionally stable and converges quadratically in space and linearly in time. Besides, the resulting linear system is solved using an iterative method, and a preconditioner is proposed to accelerate the rate of convergence. The preconditioner is constructed through circulant approximations to the Toeplitz blocks due to the coefficient matrix, which is is a block matrix with Toeplitz blocks. The spectral analysis of the preconditioned matrix is given, which demonstrates that the spectrum of the preconditioned matrix is clustered around 1. Numerical examples show the efficiency of the proposed method, and an empirical study is also provided.

To better simulate the prices of underlying assets and improve the accuracy of pricing financial derivatives, an increasing number of new models are being proposed. Among them, the Lévy process with jumps has received increasing attention because of its capacity to model sudden movements in asset prices. This paper explores the Hamilton–Jacobi–Bellman (HJB) equation with a fractional derivative and an integro-differential operator, which arise in the valuation of American options and stock loans based on the Lévy- α -stable process with jumps model. We design a fast solution strategy that includes the policy iteration method, Krylov subspace method, and banded preconditioner, aiming to solve this equation rapidly. To solve the resulting HJB equation, a finite difference method including an upwind scheme, shifted Grünwald approximation, and trapezoidal method is developed with stability and convergence analysis. Then, an algorithmic framework involving the policy iteration method and the Krylov subspace method is employed. To improve the performance of the above solver, a banded preconditioner is proposed with condition number analysis. Finally, two examples, sugar option pricing and stock loan valuation, are provided to illustrate the effectiveness of the considered model and the efficiency of the proposed preconditioned policy–Krylov subspace method. [ABSTRACT FROM AUTHOR]

Fractional derivatives and regime-switching models are widely used in various fields of finance because they can describe the nonlocal properties of the solutions and the changes in the market status, respectively. The regime-switching time-fractional diffusion equations that combine both advantages are also used in European option pricing; however, to our knowledge, American option pricing based on such models and their numerical methods is yet to be studied. Hence, a fast algorithm for solving the multi-state time-fractional linear complementary problem arising from the regime-switching time-fractional American option pricing models is developed in this paper. To construct the solution strategy, the original problem has been converted into a Hamilton–Jacobi–Bellman equation, and a nonlinear finite difference scheme has been proposed to discretize the problem with stability analysis. A policy-Krylov subspace method is developed to solve the nonlinear scheme. Further, to accelerate the convergence rate of the Krylov method, a tri-diagonal preconditioner is proposed with condition number analysis. Numerical experiments are presented to demonstrate the validity of the proposed nonlinear scheme and the efficiency of the proposed preconditioned policy-Krylov subspace method. [ABSTRACT FROM AUTHOR]

The anisotropic space-fractional diffusion equations in two dimensions are discretized by the Crank-Nicolson difference scheme with the weighted and shifted Grünwald formula, which is unconditionally stable and second-order convergence. The coefficient matrix of the discretized linear system possesses a two-level Toeplitz-like structure. Due to the anisotropy, the standard multigrid method converges slowly. By utilizing the GMRES method with a newly proposed tridiagonal preconditioner as a smoother, the convergence rate of the multigrid method can be accelerated significantly. The proposed tridiagonal preconditioner is shown to be invertible and a numerical experiment is given to demonstrate the efficiency of the proposed multigrid method with preconditioned smoother. [ABSTRACT FROM AUTHOR]

A second order alternating direction implicit scheme for time-dependent Riesz space distributed-order advection-dispersion equations is applied to higher dimensions with the Tensor-Train decomposition technique. The solutions are solved in compressed format, the Tensor-Train format, and the errors accumulated due to compressions are analyzed to ensure convergence. Problems with low-rank data are tested, the results illustrated a steeper growth in the ranks of the numerical solutions than that in related works. [ABSTRACT FROM AUTHOR]

High dimensional conservative spatial distributed-order fractional diffusion equation is discretized by midpoint quadrature rule, Crank–Nicolson method, and a finite volume approximation, with alternating direction implicit scheme. The resulting scheme is shown to be consistent and unconditionally stable, hence convergent with order 3 − α , where α is the maximum of the involving fractional orders. Moreover, if the initial condition and source term possess Tensor-Train format (TT-format) with low TT-ranks, the scheme can be solved in TT-format, such that higher dimensional cases can be considered. Perturbation analysis ensures that the accumulated errors due to data recompression do not affect the overall convergence order. Numerical examples with low TT-ranks initial conditions and source terms, and with dimensions up to 20 are tested. [ABSTRACT FROM AUTHOR]

In this article, a fast‐iterative method and a fast‐direct method is proposed for solving one‐dimensional and two‐dimensional tempered fractional diffusion equations with constant coefficients. The proposed iterative method is accelerated by circulant preconditioning which is shown to converge superlinearly while the proposed direct method is based on circulant and skew‐circulant representation for Toeplitz matrix inversion. In one‐dimensional case, the operation cost of the proposed methods are both shown to be O ( N log N ) with O ( N ) memory requirement in each time step, where N is the number of spatial nodes. With the alternating direction implicit method, it is proven that the proposed fast solution algorithms can be extended to handle two‐dimensional tempered fractional diffusion equations with O ( N 2 log N ) operation cost and O ( N 2 ) memory requirement in each time step, where the number of spatial nodes in x ‐direction and y ‐direction both equal to N . Numerical examples are provided to illustrate the effectiveness and efficiency of the proposed methods. [ABSTRACT FROM AUTHOR]

In this paper, fast numerical methods for solving space-fractional diffusion equations are studied in two stages. Firstly, a fast direct solver for an implicit finite difference scheme proposed by Haoet al.[A fourth-order approximation of fractional derivatives with its applications, J. Comput. Phys. 281 (2015), pp. 787–805], which is fourth-order accurate in space and second-order accurate in time, is developed based on a circulant-and-skew-circulant (CS) representation of Toeplitz matrix inversion. Secondly, boundary value method with spatial discretization of Haoet al.[A fourth-order approximation of fractional derivatives with its applications, J. Comput. Phys. 281 (2015), pp. 787–805] is adopted to produce a numerical solution with higher order accuracy in time. Particularly, a method with fourth-order accuracy in both space and time can be achieved. GMRES method is employed for solving the discretized linear system with two preconditioners. Based on the CS representation of Toeplitz matrix inversion, the two preconditioners can be applied efficiently, and the convergence rate of the preconditioned GMRES method is proven to be fast. Numerical examples are given to support the theoretical analysis. [ABSTRACT FROM PUBLISHER]

In this paper, finite difference schemes for differential equations with both temporal and spatial fractional derivatives are studied. When the order of the time fractional derivative is in, an alternating direction implicit (ADI) scheme with second-order accuracy in both space and time is constructed. For equations with time fractional derivatives of order lying in, a scheme is derived and solved by the generalized minimal residual method. We also propose a preconditioner to improve the efficiency for the implementation of the scheme in this situation. [ABSTRACT FROM PUBLISHER]

In recent years, considerable literature has proposed the more general class of exponential Lévy processes as the underlying model for prices of financial quantities, which thus better explain many important empirical facts of financial markets. Finite moment log stable, Carr–Geman–Madan–Yor and KoBoL models are chosen from those above-mentioned models as the dynamics of underlying equity prices in this paper. With such models pricing barrier options, one kind of financial derivatives is transformed to solve specific fractional partial differential equations (FPDEs). This study focuses on numerically solving these FPDEs via the fully implicit scheme, with the shifted Grünwald approximation. The circulant preconditioned generalized minimal residual method which converges very fast with theoretical proof is incorporated for solving resultant linear systems. Numerical examples are given to demonstrate the effectiveness of the proposed preconditioner and show the accuracy of our method compared with that done by the Fourier cosine expansion method as a benchmark. [ABSTRACT FROM PUBLISHER]

An implicit second-order finite difference scheme, which is unconditionally stable, is employed to discretize fractional advection–diffusion equations with constant coefficients. The resulting systems are full, unsymmetric, and possess Toeplitz structure. Circulant and skew-circulant splitting iteration is employed for solving the Toeplitz system. The method is proved to be convergent unconditionally to the solution of the linear system. Numerical examples show that the convergence rate of the method is fast. [ABSTRACT FROM PUBLISHER]

We consider a discrete-time orthogonal spline collocation scheme for solving Schrödinger equation with wave operator. The scheme is proposed recently by Wang et al. (J Comput Appl Math 235 (2011), 1993-2005) and is showed to have high-order convergence rate when a parameter θ in the scheme is not less than \documentclass{article} \usepackage{amsmath,amsfonts, amssymb}\pagestyle{empty}\begin{document}$\frac{1}{4}$\end{document}. In this article, we show that the result can be extended to include \documentclass{article} \usepackage{amsmath,amsfonts, amssymb}\pagestyle{empty}\begin{document}$\theta\in(0,\frac{1}{4})$\end{document} under an assumption. Numerical example is given to justify the theoretical result. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013 [ABSTRACT FROM AUTHOR]

We consider the solution of differential equations with multidelays by using boundary value methods (BVMs). These methods require the solution of some nonsymmetric, large and sparse linear systems. The GMRES method with the Strang-type block-circulant preconditioner is proposed to solve these linear systems. If an Ak1,k2-stable BVM is used, we show that our preconditioner is invertible and the spectrum of the preconditioned matrix is clustered. It follows that when the GMRES method is applied to solving the preconditioned systems, the method would converge fast. Numerical results are given to show the effectiveness of our methods. [Copyright &y& Elsevier]

In this paper, we study nonsymmetric and highly nondiagonally dominant linear systems that arise from discretizations of constant-coefficient first-order partial differential equations (PDEs). We apply the generalized minimal residual method [Y. Saad, Iterative Methods for Sparse Linear Systems, PWS Publishing Company, Boston] for solving the system with a preconditioner based on the fast sine transform. An analytic formula for the eigenvalues of the preconditioned matrix is derived and it is shown that the eigenvalues are clustered around 1 except some outliers. The outlier eigenvalues are bounded and well separated from the origin when the size of system increases. In numerical experiments, we compare our preconditioner with the semi-Toeplitz preconditioner proposed in [SIAM J. Sci. Comput. 17 (1996) 47]. We refer to [J. Numer. Linear Algebra Appl. 1 (1992) 77, Numer. Math. J. Chinese Univ. 2 (1993) 116, BIT 32 (1992) 650, Linear Algebra Appl. 293 (1999) 85] for the early works on preconditioning techniques for PDEs. [Copyright &y& Elsevier]

In recent years, fractional partial differential equation (FPDE) has been widely applied in options pricing problems, which better explains many important empirical facts of financial markets. However, the vast majority of the literatures focus on pricing the single asset option under the FPDE framework. In this paper, a two-dimensional FPDE governing the valuation of rainbow options is established, where two underlying assets are assumed to follow independent exponential Lévy processes, and its boundary conditions are determined by solving one-dimensional FPDEs. A second-order accurate finite difference scheme is proposed to discretize the two-dimensional FPDE. Given the block Toeplitz with Toeplitz block structure of the coefficient matrix, a fast preconditioned Krylov subspace method is developed for solving the resulting linear system with O (N log N) computational complexity per iteration, where N is the matrix size. The proposed preconditioner accelerates the convergence of the iterative method with theoretical analysis. Numerical examples are given to demonstrate the accuracy and efficiency of our proposed numerical methods. [ABSTRACT FROM AUTHOR]