Preconditioners for fractional diffusion equations based on the spectral symbol
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS(2022)
摘要
It is well known that the discretization of fractional diffusion equations with fractional derivatives alpha is an element of (1, 2), using the so-called weighted and shifted Grunwald formula, leads to linear systems whose coefficient matrices show a Toeplitz-like structure. More precisely, in the case of variable coefficients, the related matrix sequences belong to the so-called generalized locally Toeplitz class. Conversely, when the given FDE has constant coefficients, using a suitable discretization, we encounter a Toeplitz structure associated to a nonnegative function f(alpha), called the spectral symbol, having a unique zero at zero of real positive order between one and two. For the fast solution of such systems by preconditioned Krylov methods, several preconditioning techniques have been proposed in both the one- and two-dimensional cases. In this article we propose a new preconditioner denoted by p(p alpha) which belongs to the tau-algebra and it is based on the spectral symbol f(alpha). Comparing with some of the previously proposed preconditioners, we show that although the low band structure preserving preconditioners are more effective in the one-dimensional case, the new preconditioner performs better in the more challenging multi-dimensional setting.
更多查看译文
关键词
fractional differential equations,fractional order zero,GMRES,multi-level Toeplitz matrix,sine transform based preconditioner
AI 理解论文
溯源树
样例
![](https://originalfileserver.aminer.cn/sys/aminer/pubs/mrt_preview.jpeg)
生成溯源树,研究论文发展脉络
Chat Paper
正在生成论文摘要