Abstract
In this paper, we study the block lower triangular Toeplitz-like with tri-diagonal blocks system which arises from the time-fractional partial differential equation. Existing fast numerical solver (e.g., fast approximate inversion method) cannot handle such linear system as the main diagonal blocks are different. The main contribution of this paper is to propose a fast direct method for solving this linear system, and to illustrate that the proposed method is much faster than the classical block forward substitution method for solving this linear system. Our idea is based on the divide-and-conquer strategy and together with the fast Fourier transforms for calculating Toeplitz matrix-vector multiplication. The complexity needs O(MNlog2M) arithmetic operations, where M is the number of blocks (the number of time steps) in the system and N is the size (number of spatial grid points) of each block. Numerical examples from the finite difference discretization of time-fractional partial differential equations are also given to demonstrate the efficiency of the proposed method.
Original language | English |
---|---|
Pages (from-to) | 203-211 |
Number of pages | 9 |
Journal | Journal of Computational Physics |
Volume | 303 |
DOIs | |
Publication status | Published - 2015 |
Scopus Subject Areas
- Numerical Analysis
- Modelling and Simulation
- Physics and Astronomy (miscellaneous)
- General Physics and Astronomy
- Computer Science Applications
- Computational Mathematics
- Applied Mathematics
User-Defined Keywords
- Block triangular Toeplitz-like matrix
- Direct methods
- Divide-and-conquer strategy
- Fast Fourier transform
- Fractional partial differential equations