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 条
[11]  
Benzi M(1988)An optimal circulant preconditioner for Toeplitz systems SIAM J. Sci. Statist. Comput. 9 766-771
[12]  
Breiten T(2014)Fourth order accurate scheme for the space fractional diffusion equations SIAM J. Numer. Anal. 52 1418-1438
[13]  
Simoncini V(2013)Superlinearly convergent algorithms for the two-dimensional space–time Caputo–Riesz fractional diffusion equation Appl. Numer. Math. 70 22-41
[14]  
Stoll M(2018)Boundary problems for the fractional and tempered fractional operators Multiscale Model. Simul. 16 125-149
[15]  
Carreras BA(2016)Fast tensor product solvers for optimization problems with fractional differential equations as constraints Appl. Math. Comput. 273 604-623
[16]  
Lynch VE(2016)Spectral analysis and structure preserving preconditioners for fractional diffusion equations J. Comput. Phys. 307 262-279
[17]  
Zaslavsky GM(2010)Hierarchical singular value decomposition of tensors SIAM J. Matrix Anal. Appl. 31 2029-2054
[18]  
Chan R(2018)Tempered fractional diffusion equations for princing multi-asset options under CGMYe process Comput. Math. Appl. 76 1500-1514
[19]  
Jin X(2009)A new scheme for the tensor representation J. Fourier Anal. Appl. 15 706-722
[20]  
Chan R(2015)A fourth-order approximation of fractional derivatives with its applications J. Comput. Phys. 281 787-805