科研进展
一阶/二阶时间发展方程预条件时空离散矩阵的一致谱分析(周涛)
发布时间:2022-12-16 |来源:

  Solving evolutionary equations in a parallel-in-time manner is an attractive topic. The iterative algorithm based on the block $\alpha$-circulant preconditioning technique has shown promising advantages, especially for hyperbolic problems. By fast Fourier transform for factorizing the involved circulant matrices, the preconditioned iteration can be computed efficiently via the so-called diagonalization technique, which yields a direct parallel implementation across all time levels. In recent years, considerable efforts have been devoted to exploring the spectral property of the iteration matrix arising from the used time-integrator, which leads to many case-by-case studies. Denoting by $\mathcal{K} $ and $\mathcal{P}_\alpha$ the all-at-once matrix of the evolutionary PDEs and the corresponding block $\alpha$-circulant preconditioner, we will present a systematic spectral analysis for the matrix $\mathcal{P}_\alpha^{-1}\mathcal{K}$ for both the first-order and second-order evolutionary problems. For the first-order problems our analysis works for all stable single-step time-integrators, while for the second-order problems our analysis works for a large class of symmetric two-step methods which could be arbitrarily high-order. Illustrative numerical experiments are presented to complement our theory. 

   

  Publication: 

  SIAM Journal on Matrix Analysis and Applications, Volume 43, Issue 3, pp: 1331-1353. 

    

  Author: 

  Shulin Wu 

  School of Mathematics and Statistics, Northeast Normal University, Changchun, 130024, China 

  Tao Zhou 

  LSEC, Institute of Computational Mathematics and Scientific/Engineering Computing, AMSS, Chinese Academy of Sciences, Beijing, 100190, China  

  Email: tzhou@lsec.cc.ac.cn 

  Zhi Zhou 

  Department of Applied Mathematics, The Hong Kong Polytechnic University, Kowloon, Hong Kong, China 


附件下载:

    联系我们
    参考
    相关文章