Preconditioning of spectral methods via Birkhoff interpolation

被引:0
作者
Conor McCoid
Manfred R. Trummer
机构
[1] Simon Fraser University,Department of Mathematics
来源
Numerical Algorithms | 2018年 / 79卷
关键词
Two-point boundary value problems; Spectral methods; Preconditioning; Birkhoff interpolation; 65N35; 65L10;
D O I
暂无
中图分类号
学科分类号
摘要
High-order differentiation matrices as calculated in spectral collocation methods usually include a large round-off error and have a large condition number (Baltensperger and Berrut Computers and Mathematics with Applications 37(1), 41–48 1999; Baltensperger and Trummer SIAM J. Sci. Comput. 24(5), 1465–1487 2003; Costa and Don Appl. Numer. Math. 33(1), 151–159 2000). Wang et al. (Wang et al. SIAM J. Sci. Comput. 36(3), A907–A929 2014) present a method to precondition these matrices using Birkhoff interpolation. We generalize this method for all orders and boundary conditions and allowing arbitrary rows of the system matrix to be replaced by the boundary conditions. The preconditioner is an exact inverse of the highest-order differentiation matrix in the equation; thus, its product with that matrix can be replaced by the identity matrix. We show the benefits of the method for high-order differential equations. These include improved condition number and, more importantly, higher accuracy of solutions compared to other methods.
引用
收藏
页码:555 / 573
页数:18
相关论文
共 21 条
  • [1] Baltensperger R(1999)The errors in calculating the pseudospectral differentiation matrices for C̆ebys̆ev-Gauss-Lobatto points Computers and Mathematics with Applications 37 41-48
  • [2] Berrut JP(2003)Spectral differencing with a twist SIAM J. Sci. Comput. 24 1465-1487
  • [3] Baltensperger R(1906)General mean value and remainder theorems with applications to mechanical differentiation and quadrature Trans. Am. Math. Soc. 7 107-136
  • [4] Trummer MR(2000)On the computation of high order pseudospectral derivatives Appl. Numer. Math. 33 151-159
  • [5] Birkhoff GD(2016)Rectangular spectral collocation IMA J. Numer. Anal. 36 108-132
  • [6] Costa B(2016)On well-conditioned spectral collocation and spectral methods by the integral reformulation SIAM J. Sci Comput. 38 A3247-A3263
  • [7] Don WS(1992)The pseudospectral method for third-order differential equations SIAM J. Numer. Anal. 29 1626-1647
  • [8] Driscoll TA(2013)A fast and well-conditioned spectral method SIAM Rev. 55 462-489
  • [9] Hale N(1966)On Hermite-Birkhoff interpolation J. Math. Anal. Appl. 16 538-543
  • [10] Du K(1994)Efficient spectral-galerkin method i. direct solvers of second- and fourth-order equations using legendre polynomials SIAM J. Sci. Comput. 15 1489-1505