Circulant preconditioners for convolution-like integral equations with higher-order quadrature rules
Electronic transactions on numerical analysis, Tome 5 (1997), pp. 18-28
In this paper, we consider solving matrix systems arising from the discretization of convolution-like integral equations by preconditioned conjugate gradient (PCG) methods. Circulant integral operators as preconditioners have been proposed and studied. However, the discretization of these circulant preconditioned equations by employing higher-order quadratures leads to matrix systems that cannot be solved efficiently by using fast Fourier transforms (FFTs). The aim of this paper is to propose "inverted" circulant preconditioners for convolution-like integral equations. The discretization of these preconditioned integral equations by higher-order quadratures leads to matrix systems that involve only Toeplitz, circulant and diagonal matrix-vector multiplications, and hence can be computed efficiently by FFTs in each iteration. Numerical examples are given to illustrate the fast convergence of the method and the improvement of the accuracy of the computed solutions with using higher-order quadratures.
Classification :
45E10, 45L10, 65R20, 65J10
Keywords: integral equations, displacement kernel, quadratures, circulant matrices, Toeplitz matrices, fast Fourier transforms, signal processing
Keywords: integral equations, displacement kernel, quadratures, circulant matrices, Toeplitz matrices, fast Fourier transforms, signal processing
@article{ETNA_1997__5__a4,
author = {Ng, Michael K.},
title = {Circulant preconditioners for convolution-like integral equations with higher-order quadrature rules},
journal = {Electronic transactions on numerical analysis},
pages = {18--28},
year = {1997},
volume = {5},
zbl = {0910.65112},
language = {en},
url = {http://geodesic.mathdoc.fr/item/ETNA_1997__5__a4/}
}
TY - JOUR AU - Ng, Michael K. TI - Circulant preconditioners for convolution-like integral equations with higher-order quadrature rules JO - Electronic transactions on numerical analysis PY - 1997 SP - 18 EP - 28 VL - 5 UR - http://geodesic.mathdoc.fr/item/ETNA_1997__5__a4/ LA - en ID - ETNA_1997__5__a4 ER -
Ng, Michael K. Circulant preconditioners for convolution-like integral equations with higher-order quadrature rules. Electronic transactions on numerical analysis, Tome 5 (1997), pp. 18-28. http://geodesic.mathdoc.fr/item/ETNA_1997__5__a4/