Iterative solution of shifted positive-definite linear systems arising in a numerical method for the heat equation based on Laplace transformation and quadrature

DOI10.1017/S1446181112000107zbMATH Open1264.65046arXiv1111.5105OpenAlexW3105719630MaRDI QIDQ2911076FDOQ2911076

Authors: William Mclean, Vidar Thomée 

Publication date: 12 September 2012

Published in: The ANZIAM Journal (Search for Journal in Brave)

Abstract: In earlier work we have studied a method for discretization in time of a parabolic problem which consists in representing the exact solution as an integral in the complex plane and then applying a quadrature formula to this integral. In application to a spatially semidiscrete finite element version of the parabolic problem, at each quadrature point one then needs to solve a linear algebraic system having a positive definite matrix with a complex shift, and in this paper we study iterative methods for such systems. We first consider the basic and a preconditioned version of the Richardson algorithm, and then a conjugate gradient method as well as a preconditioned version thereof.

Full work available at URL: https://arxiv.org/abs/1111.5105

Recommendations

• scientific article; zbMATH DE number 219982
• scientific article; zbMATH DE number 592768
• scientific article; zbMATH DE number 1183050
• A polynomial based iterative method for linear parabolic equations
• Non-stationary iterative solution of positive definite linear systems

zbMATH Keywords

Laplace transformpreconditioningalgorithmfinite elementsheat equationsemidiscretizationconjugate gradient methodquadratureRichardson iteration

Mathematics Subject Classification ID

Preconditioners for iterative methods (65F08) Iterative numerical methods for linear systems (65F10) Laplace transform (44A10) Method of lines for initial value and initial-boundary value problems involving PDEs (65M20) Finite element, Rayleigh-Ritz and Galerkin methods for initial value and initial-boundary value problems involving PDEs (65M60) Numerical methods for integral transforms (65R10) Transform methods (e.g., integral transforms) applied to PDEs (35A22) Numerical solution of discretized equations for initial value and initial-boundary value problems involving PDEs (65M22)

Cites Work

• Algorithm 740: Fortran subroutines to compute improved incomplete Cholesky factorizations
• Algorithm 832
• Block preconditioning of real-valued iterative algorithms for complex linear systems
• A parallel method for time-discretization of parabolic problems based on contour integral representation and quadrature
• A parallel method for time discretization of parabolic equations based on Laplace transformation and quadrature
• Numerical solution via Laplace transforms of a fractional order evolution equation
• Solution of Sparse Indefinite Systems of Linear Equations
• A Contour Integral Method for the Black–Scholes and Heston Equations
• Recent computational developments in Krylov subspace methods for linear systems
• Time discretization of an evolution equation via Laplace transforms
• Maximum-norm error analysis of a numerical solution via Laplace transformation and quadrature of a fractional-order evolution equation
• Exponentially Convergent Algorithms for the Operator Exponential with Applications to Inhomogeneous Problems in Banach Spaces
• Time discretization via Laplace transformation of an integro-differential equation of parabolic type
• On conjugate gradient type methods and polynomial preconditioners for a class of complex non-Hermitian matrices
• Restarted full orthogonalization method for shifted linear systems
• The Solution of Parametrized Symmetric Linear Systems

Cited In (2)

• The approximation of parabolic equations involving fractional powers of elliptic operators
• Numerical approximation of fractional powers of elliptic operators

Uses Software

• PyAMG

This page was built for publication: Iterative solution of shifted positive-definite linear systems arising in a numerical method for the heat equation based on Laplace transformation and quadrature

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2911076)