Preconditioning Stochastic Galerkin Saddle Point Systems
SIAM Journal on Matrix Analysis and Applications
SIAM Journal on Scientific Computing
SIAM Journal on Scientific Computing
Journal of Computational Physics
Hi-index | 0.01 |
The discretization of linear partial differential equations with random data by means of the stochastic Galerkin finite element method results in general in a large coupled linear system of equations. Using the stochastic diffusion equation as a model problem, we introduce and study a symmetric positive definite Kronecker product preconditioner for the Galerkin matrix. We compare the popular mean-based preconditioner with the proposed preconditioner which—in contrast to the mean-based construction—makes use of the entire information contained in the Galerkin matrix. We report on results of test problems, where the random diffusion coefficient is given in terms of a truncated Karhunen-Loève expansion or is a lognormal random field.