This page is a copy of research/scientific_computing/former/prec_trf (Wed, 31 Aug 2022 15:00:53)
Preconditioners based on trigonometric transforms
Participants
- Kurt Otto
(coordinator), Dept. of Scientific Computing, Uppsala Univ.
- Sverker Holmgren, Dept. of Scientific Computing, Uppsala Univ.
- Elisabeth Larsson, Dept. of Scientific Computing, Uppsala Univ.
- Eva Mossberg, Dept. of Scientific Computing, Uppsala Univ.
Research
The state of the art for solving linear systems of equations arising from discretizations of PDEs is to employ some Krylov subspace method. In order to achieve an acceptable rate of convergence and, more importantly, a short total execution time, it is crucial to construct effective, parallelizable preconditioners.
We have designed preconditioners [1] based on the fast Fourier transform, which have been expediently used for second-order accurate discretizations of first-order systems of PDEs [4]. The solution procedures are highly parallelizable, and have already from the beginning been implemented on a variety of parallel computer architectures. Moreover, the convergence properties have been thoroughly analyzed [3,2,5]. The preconditioning technique has been generalized to several fast trigonometric transforms [6,#], which resulted in favorable convergence properties for high-order discretizations of a scalar first-order PDE [7] and for a second-order discretization of the Helmholtz equation [10]. It has also been successfully applied to discretizations of waveguide problems in underwater acoustics [11,13,14] and electromagnetics [12].
Publications
Refereed
-
Iterative solution methods and preconditioners for block-tridiagonal systems of equations
. In SIAM Journal on Matrix Analysis and Applications, volume 13, pp 863-886, 1992. (DOI
).
-
Semicirculant preconditioners for first-order partial differential equations
. In SIAM Journal on Scientific Computing, volume 15, pp 385-407, 1994. (DOI
).
-
Analysis of preconditioners for hyperbolic partial differential equations
. In SIAM Journal on Numerical Analysis, volume 33, pp 2131-2165, 1996. (DOI
).
-
Semicirculant solvers and boundary corrections for first-order partial differential equations
. In SIAM Journal on Scientific Computing, volume 17, pp 613-630, 1996. (DOI
).
-
Analysis of semi-Toeplitz preconditioners for first-order PDEs
. In SIAM Journal on Scientific Computing, volume 17, pp 47-64, 1996. (DOI
).
-
A framework for polynomial preconditioners based on fast transforms I: Theory
. In BIT Numerical Mathematics, volume 38, pp 544-559, 1998. (DOI
).
-
A framework for polynomial preconditioners based on fast transforms II: PDE applications
. In BIT Numerical Mathematics, volume 38, pp 721-736, 1998. (DOI
).
-
Object-Oriented Construction of Parallel PDE Solvers
. In Modern Software Tools for Scientific Computing, pp 203-226, Birkhäuser, Boston, MA, 1997.
-
Object-oriented software tools for the construction of preconditioners
. In Scientific Programming, volume 6, pp 285-295, 1997. (External link
).
-
Iterative solution of the Helmholtz equation by a second-order method
. In SIAM Journal on Matrix Analysis and Applications, volume 21, pp 209-229, 1999. (DOI
).
-
Iterative solution of the Helmholtz equation by a fourth-order method
. In suppl, volume 40:1 of Bollettino di Geofisica Teorica ed Applicata, pp 104-105, OGS, Trieste, Italy, 1999.
-
A domain decomposition method for the Helmholtz equation in a multilayer domain
. In SIAM Journal on Scientific Computing, volume 20, pp 1713-1731, 1999. (DOI
).
-
Helmholtz and parabolic equation solutions to a benchmark problem in ocean acoustics
. In Journal of the Acoustical Society of America, volume 113, pp 2446-2454, 2003. (DOI
).
-
Parallel solution of the Helmholtz equation in a multilayer domain
. In BIT Numerical Mathematics, volume 43, pp 387-411, 2003. (DOI
).
Doctoral theses
-
Construction and Analysis of Preconditioners for First-order PDE
. Ph.D. thesis, Comprehensive summaries of Uppsala dissertations from the Faculty of Science nr 439, Acta Universitatis Upsaliensis, Uppsala, 1993.
-
Domain Decomposition and Preconditioned Iterative Methods for the Helmholtz Equation
. Ph.D. thesis, Comprehensive Summaries of Uppsala Dissertations from the Faculty of Science and Technology nr 523, Acta Universitatis Upsaliensis, Uppsala, 2000.