The fractional diffusion equation is discretized by an implicit finite difference scheme with the shifted Grünwald formula, which is unconditionally stable. The coefficient matrix of the discretized linear system is equal to the sum of a scaled identity matrix and two diagonal-timesToeplitz matrices. Standard circulant preconditioners may not work for such Toeplitz-like linear systems. The main aim of this paper is to propose and develop approximate inverse preconditioners for such Toeplitz-like matrices. An approximate inverse preconditioner is constructed to approximate the inverses of weighted Toeplitz matrices by circulant matrices, and then combine them together rowby-row. Because of Toeplitz structure, both the discretized coefficient matrix and the preconditioner can be implemented very efficiently by using fast Fourier transforms. Theoretically, we show that the spectra of the resulting preconditioned matrices are clustered around one. Thus Krylov subspace methods with the proposed preconditioner converge very fast. Numerical examples are given to demonstrate the effectiveness of the proposed preconditioner and show that its performance is better than the other testing preconditioners.
Stochastic matrices play an important role in the study of probability theory
and statistics, and are often used in a variety of modeling problems in
economics, biology and operation research. Recently, the study
of tensors and
their applications
became a hot topic in numerical analysis and
optimization.
In this paper, we focus on studying stochastic tensors and, in
particular, we study the extreme points of a set of multi-stochastic tensors. Two
necessary and sufficient conditions for a multi-stochastic tensor to be an
extreme point are established.
These conditions characterize the “generators” of multi-stochastic tensors.
An algorithm to search the convex combination of extreme points for an arbitrary given
multi-stochastic tensor is developed. Based on our obtained results, some expression
properties for third-order and n-dimensional multi-stochastic tensors
(${n=3}$ and 4) are derived, and all extreme points of
3-dimensional and 4-dimensional triply-stochastic tensors can be produced in a simple
way. As an application, a new approach for the
partially filled square problem under the framework of multi-stochastic
tensors is given.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.