Tensor-Train Format Solution with Preconditioned Iterative Method for High Dimensional Time-Dependent Space-Fractional Diffusion Equations with Error Analysis

被引：0
作者：
Lot-Kei Chou
Siu-Long Lei
机构：
[1] University of Macau,Department of Mathematics
[2] Avenida da Universidade,undefined
来源：
Journal of Scientific Computing | 2019年 / 80卷
关键词：
High dimensional fractional diffusion equation; Tensor-Train decomposition; Krylov subspace method; Preconditioner; 65N22; 65F10; 26A33; 41A63;
D O I：
暂无
中图分类号：
学科分类号：
摘要：
In this paper, a first order implicit finite difference scheme with Krylov subspace linear system solver is employed to solving time-dependent space-fractional diffusion equations in high dimensions where the initial condition and source term are in tensor-train (TT) format with low TT-ranks. In the time-marching process, TT-format of the solution is maintained and the increment of TT-ranks due to addition is moderated by rounding. The error introduced by rounding is shown to be consistent with the first order finite difference scheme. On the other hand, the linear systems involved in the solution process are shown to possess Toeplitz-like structure so that the complexity and required memory for Krylov subspace solver can be optimized. Further reduction in complexity is made by utilizing a circulant preconditioner which accelerates the convergence rate of Krylov subspace method dramatically. Numerical examples for problems up to 20 dimensions are presented.
引用
收藏
页码：1731 / 1763
页数：32
相关论文
共 140 条
[31]  
Li B(2006)Finite difference approximations for two-sided space-fractional partial differential equations Appl. Numer. Math. 56 80-90
[32]  
Tian W(2006)Finite difference methods for two-dimensional fractional dispersion equation J. Comput. Phys. 211 249-261
[33]  
Zhang P(2017)Spectral analysis and multigrid preconditioners for two-dimensional space-fractional diffusion equations J. Comput. Phys. 350 992-1011
[34]  
Dolgov S(2013)A banded preconditioner for the two-sided, nonlinear space-fractional diffusion equation Comput. Math. Appl. 66 659-667
[35]  
Pearson JW(2011)Tensor-train decomposition SIAM J. Sci. Comput. 33 2295-2317
[36]  
Savostyanov DV(2013)Constructive representation of functions in low-rank tensor formats Constr. Approx. 37 1-18
[37]  
Stoll M(2009)Breaking the curse of dimensionality, or how to use SVD in many dimensions SIAM J. Sci. Comput. 31 3744-3759
[38]  
Donatelli M(2014)Preconditioning techniques for diagonal-times-Toeplitz matrices in fractional diffusion equations SIAM J. Sci. Comput. 36 A2698-2719
[39]  
Mazza M(2012)Multigrid method for fractional diffusion equations J. Comput. Phys. 231 693-703
[40]  
Serra-Capizzano S(2014)Circulant and skew-circulant splitting iteration for fractional advection-diffusion equations Int. J. Comput. Math. 91 2232-2242
12345678910