Abstract
In this paper, we propose a new type of preconditioners for solving finite section Wiener-Hopf integral equations (αI + Aτ)xτ = g by the preconditioned conjugate gradient algorithm. We show that for an integer u > 1, the operator αI + Aτ> can be decomposed into a sum of operators αI + Pτ(u,v) for 0 ≤ v < u. Here Pτ(u,v) are gwvcirculant matrices. For u - 1, our preconditioners are defined as ( 1 u)∑v(αI+Pτ(u,v))-1. Thus the way the preconditioners are constructed is very similar to the approach used in the additive Schwarz method for elliptic problems. As for the convergence rate, we prove that the spectra of the resulting preconditioned operators ( 1 u)∑v(αI+Pτ(u,v))-1][αI+Aτ are clustered around 1 and thus the algorithm converges sufficiently fast. Finally, we discretize the resulting preconditioned equations by rectangular rule. Numerical results show that our methods converges faster than those preconditioned by using circulant integral operators.
Original language | English |
---|---|
Pages (from-to) | 77-96 |
Number of pages | 20 |
Journal | Applied Mathematics and Computation |
Volume | 72 |
Issue number | 1 |
DOIs | |
Publication status | Published - 15 Sept 1995 |
Externally published | Yes |
Funding
This research was supported in part by HKRGC grantno. CUHK 178/93E.