File Download
Links for fulltext
(May Require Subscription)
- Publisher Website: 10.1137/S1064827598310227
- Scopus: eid_2-s2.0-0033293884
- WOS: WOS:000084272500003
- Find via
Supplementary
- Citations:
- Appears in Collections:
Article: Galerkin Projection Methods for Solving Multiple Linear Systems
Title | Galerkin Projection Methods for Solving Multiple Linear Systems |
---|---|
Authors | |
Keywords | Multiple linear systems Krylov space Conjugate gradient method Galerkin projection |
Issue Date | 1999 |
Publisher | Society for Industrial and Applied Mathematics. The Journal's web site is located at http://www.siam.org/journals/sisc.php |
Citation | SIAM Journal on Scientific Computing, 1999, v. 21 n. 3, p. 836-850 How to Cite? |
Abstract | In this paper, we consider using conjugate gradient (CG) methods for solving multiple linear systems A^{(i)} x^{(i)} = b^{(i)} , for 1 ≤ i ≤ s, where the coefficient matrices A^{(i)} and the right-hand sides b( i) are different in general. In particular, we focus on the seed projection method which generates a Krylov subspace from a set of direction vectors obtained by solving one of the systems, called the seed system, by the CG method and then projects the residuals of other systems onto the generated Krylov subspace to get the approximate solutions. The whole process is repeated until all the systems are solved. Most papers in the literature [T. F. Chan and W. L. Wan, SIAM J. Sci. Comput., 18 (1997), pp. 16981721; B. Parlett Linear Algebra Appl., 29 (1980), pp. 323346; Y. Saad, Math. Comp., 48 (1987), pp. 651662; V. Simoncini and E. Gallopoulos, SIAM J. Sci. Comput., 16 (1995), pp. 917933; C. Smith, A. Peterson, and R. Mittra, IEEE Trans. Antennas and Propagation, 37 (1989), pp. 14901493] considered only the case where the coefficient matrices A( i) are the same but the right-hand sides are different. We extend and analyze the method to solve multiple linear systems with varying coefficient matrices and right-hand sides. A theoretical error bound is given for the approximation obtained from a projection process onto a Krylov subspace generated from solving a previous linear system. Finally, numerical results for multiple linear systems arising from image restorations and recursive least squares computations are reported to illustrate the effectiveness of the method. |
Persistent Identifier | http://hdl.handle.net/10722/42993 |
ISSN | 2019 Impact Factor: 1.976 2015 SCImago Journal Rankings: 2.166 |
ISI Accession Number ID |
DC Field | Value | Language |
---|---|---|
dc.contributor.author | Ng, MK | en_HK |
dc.contributor.author | Chan, TF | en_HK |
dc.date.accessioned | 2007-03-23T04:36:25Z | - |
dc.date.available | 2007-03-23T04:36:25Z | - |
dc.date.issued | 1999 | en_HK |
dc.identifier.citation | SIAM Journal on Scientific Computing, 1999, v. 21 n. 3, p. 836-850 | en_HK |
dc.identifier.issn | 1064-8275 | en_HK |
dc.identifier.uri | http://hdl.handle.net/10722/42993 | - |
dc.description.abstract | In this paper, we consider using conjugate gradient (CG) methods for solving multiple linear systems A<sup>(i)</sup> x<sup>(i)</sup> = b<sup>(i)</sup> , for 1 ≤ i ≤ s, where the coefficient matrices A<sup>(i)</sup> and the right-hand sides b( i) are different in general. In particular, we focus on the seed projection method which generates a Krylov subspace from a set of direction vectors obtained by solving one of the systems, called the seed system, by the CG method and then projects the residuals of other systems onto the generated Krylov subspace to get the approximate solutions. The whole process is repeated until all the systems are solved. Most papers in the literature [T. F. Chan and W. L. Wan, SIAM J. Sci. Comput., 18 (1997), pp. 16981721; B. Parlett Linear Algebra Appl., 29 (1980), pp. 323346; Y. Saad, Math. Comp., 48 (1987), pp. 651662; V. Simoncini and E. Gallopoulos, SIAM J. Sci. Comput., 16 (1995), pp. 917933; C. Smith, A. Peterson, and R. Mittra, IEEE Trans. Antennas and Propagation, 37 (1989), pp. 14901493] considered only the case where the coefficient matrices A( i) are the same but the right-hand sides are different. We extend and analyze the method to solve multiple linear systems with varying coefficient matrices and right-hand sides. A theoretical error bound is given for the approximation obtained from a projection process onto a Krylov subspace generated from solving a previous linear system. Finally, numerical results for multiple linear systems arising from image restorations and recursive least squares computations are reported to illustrate the effectiveness of the method. | en_HK |
dc.format.extent | 351104 bytes | - |
dc.format.extent | 26112 bytes | - |
dc.format.mimetype | application/pdf | - |
dc.format.mimetype | application/msword | - |
dc.language | eng | en_HK |
dc.publisher | Society for Industrial and Applied Mathematics. The Journal's web site is located at http://www.siam.org/journals/sisc.php | - |
dc.relation.ispartof | SIAM Journal on Scientific Computing | - |
dc.rights | © 1999 Society for Industrial and Applied Mathematics. First Published in SIAM Journal on Scientific Computing in volume 21, issue 3, published by the Society for Industrial and Applied Mathematics (SIAM). | - |
dc.subject | Multiple linear systems | en_HK |
dc.subject | Krylov space | en_HK |
dc.subject | Conjugate gradient method | en_HK |
dc.subject | Galerkin projection | en_HK |
dc.title | Galerkin Projection Methods for Solving Multiple Linear Systems | en_HK |
dc.type | Article | en_HK |
dc.description.nature | published_or_final_version | en_HK |
dc.identifier.doi | 10.1137/S1064827598310227 | - |
dc.identifier.scopus | eid_2-s2.0-0033293884 | - |
dc.identifier.hkuros | 52937 | - |
dc.identifier.volume | 21 | - |
dc.identifier.issue | 3 | - |
dc.identifier.spage | 836 | - |
dc.identifier.epage | 850 | - |
dc.identifier.isi | WOS:000084272500003 | - |