dc.contributor.author | Haslinger, Jaroslav | |
dc.contributor.author | Kozubek, Tomáš | |
dc.contributor.author | Kučera, Radek | |
dc.contributor.author | Peichl, Gunther | |
dc.date.accessioned | 2007-12-17T13:19:01Z | |
dc.date.available | 2007-12-17T13:19:01Z | |
dc.date.issued | 2007 | |
dc.identifier.citation | Numerical Linear Algebra with Applications. 2007, vol. 14, issues 9, p. 713-739. | en |
dc.identifier.issn | 1070-5325 | |
dc.identifier.issn | 1099-1506 | |
dc.identifier.uri | http://hdl.handle.net/10084/64507 | |
dc.language.iso | en | en |
dc.publisher | Wiley | en |
dc.relation.ispartofseries | Numerical Linear Algebra with Applications | en |
dc.relation.uri | http://dx.doi.org/10.1002/nla.550 | en |
dc.subject | saddle-point system | en |
dc.subject | fictitious domain method | en |
dc.subject | Schur complement | en |
dc.subject | orthogonal projectors | en |
dc.subject | BiCGSTAB algorithm | en |
dc.subject | multigrid | en |
dc.title | Projected Schur complement method for solving non-symmetric systems arising from a smooth fictitious domain approach | en |
dc.type | article | en |
dc.identifier.location | Není ve fondu ÚK | en |
dc.description.abstract-en | This paper deals with a fast method for solving large-scale algebraic saddle-point systems arising from fictitious domain formulations of elliptic boundary value problems. A new variant of the fictitious domain approach is analyzed. Boundary conditions are enforced by control variables introduced on an auxiliary boundary located outside the original domain. This approach has a significantly higher convergence rate; however, the algebraic systems resulting from finite element discretizations are typically non-symmetric. The presented method is based on the Schur complement reduction. If the stiffness matrix is singular, the reduced system can be formulated again as another saddle-point problem. Its modification by orthogonal projectors leads to an equation that can be efficiently solved using a projected Krylov subspace method for non-symmetric operators. For this purpose, the projected variant of the BiCGSTAB algorithm is derived from the non-projected one. The behavior of the method is illustrated by examples, in which the BiCGSTAB iterations are accelerated by a multigrid strategy. | en |
dc.identifier.doi | 10.1002/nla.550 | |
dc.identifier.wos | 000250976900003 | |