Hierarchical Schur complement preconditioner for the stochastic Galerkin finite element methods
Loading...
Links to Files
Collections
Author/Creator
Author/Creator ORCID
Date
2013-01-25
Type of Work
Department
Program
Citation of Original Publication
Sousedík, Bedřich; Ghanem, Roger G.; Phipps, Eric T.; Hierarchical Schur complement preconditioner for the stochastic Galerkin finite element methods; Numerical Linear Algebra with Applications, 21, 1, pages 136-151, 25 January, 2013; https://doi.org/10.1002/nla.1869
Rights
This item is likely protected under Title 17 of the U.S. Copyright Law. Unless on a Creative Commons license, for uses protected by Copyright Law, contact the copyright holder or the author.
This is the pre-peer reviewed version of the following article: Sousedík, Bedřich; Ghanem, Roger G.; Phipps, Eric T.; Hierarchical Schur complement preconditioner for the stochastic Galerkin finite element methods; Numerical Linear Algebra with Applications, 21, 1, pages 136-151, 25 January, 2013; https://doi.org/10.1002/nla.1869, which has been published in final form at https://onlinelibrary.wiley.com/doi/10.1002/nla.1869. This article may be used for noncommercial purposes in accordance with Wiley Terms and Conditions for Use of Self-Archived Versions
This is the pre-peer reviewed version of the following article: Sousedík, Bedřich; Ghanem, Roger G.; Phipps, Eric T.; Hierarchical Schur complement preconditioner for the stochastic Galerkin finite element methods; Numerical Linear Algebra with Applications, 21, 1, pages 136-151, 25 January, 2013; https://doi.org/10.1002/nla.1869, which has been published in final form at https://onlinelibrary.wiley.com/doi/10.1002/nla.1869. This article may be used for noncommercial purposes in accordance with Wiley Terms and Conditions for Use of Self-Archived Versions
Subjects
Abstract
Use of the stochastic Galerkin finite element methods leads to large systems of linear equations obtained by the discretization of tensor product solution spaces along their spatial and stochastic dimensions. These systems are typically solved iteratively by a Krylov subspace method. We propose a preconditioner, which takes an advantage of the recursive hierarchy in the structure of the global matrices. In particular, the matrices posses a recursive hierarchical two-by-two structure, with one of the submatrices block diagonal. Each of the diagonal blocks in this submatrix is closely related to the deterministic mean-value problem, and the action of its inverse is in the implementation approximated by inner loops of Krylov iterations. Thus, our hierarchical Schur complement preconditioner combines, on each level in the approximation of the hierarchical structure of the global matrix, the idea of Schur complement with loops for a number of mutually independent inner Krylov iterations, and several matrix–vector multiplications for the off-diagonal blocks. Neither the global matrix nor the matrix of the preconditioner need to be formed explicitly. The ingredients include only the number of stiffness matrices from the truncated Karhunen–Loève expansion and a good preconditioned for the mean-value deterministic problem. We provide a condition number bound for a model elliptic problem, and the performance of the method is illustrated by numerical experiments. Copyright © 2013 John Wiley & Sons, Ltd.