Solution of Time-Dependent PDE Through Rapid Estimation of Block Gaussian Quadrature Nodes
Document Type
Article
Publication Date
3-1-2015
Department
Mathematics
School
Mathematics and Natural Sciences
Abstract
The stiffness of systems of ODEs that arise from spatial discretization of PDEs causes difficulties for both explicit and implicit time-stepping methods. Krylov Subspace Spectral (KSS) methods present a balance between the efficiency of explicit methods and the stability of implicit methods by computing each Fourier coefficient from an individualized approximation of the solution operator of the PDE. While KSS methods are explicit methods that exhibit a high order of accuracy and stability similar to that of implicit methods, their efficiency needs to be improved. A previous asymptotic study of block Lanczos iteration yielded estimates of extremal block Gaussian quadrature nodes for each Fourier component and led to an improvement in efficiency. In this paper, a more detailed asymptotic study is performed in order to rapidly estimate all nodes, thus drastically reducing computational expense without sacrificing accuracy. Numerical results demonstrate that the new node estimation scheme does in fact accomplish these aims.
Publication Title
Linear Algebra and its Applications
Volume
468
First Page
233
Last Page
259
Recommended Citation
Palchak, E. M.,
Cibotarica, A.,
Lambers, J. V.
(2015). Solution of Time-Dependent PDE Through Rapid Estimation of Block Gaussian Quadrature Nodes. Linear Algebra and its Applications, 468, 233-259.
Available at: https://aquila.usm.edu/fac_pubs/15277