Explicit cost bounds of stochastic Galerkin approximations for parameterized PDEs with random coefficients

From MaRDI portal
Publication:2007265

DOI10.1016/J.CAMWA.2015.12.005zbMATH Open1443.65200arXiv1507.05545OpenAlexW1464641790MaRDI QIDQ2007265FDOQ2007265

Clayton G. Webster, Guannan Zhang, Nick C. Dexter

Publication date: 12 October 2020

Published in: Computers & Mathematics with Applications (Search for Journal in Brave)

Abstract: This work analyzes the overall computational complexity of the stochastic Galerkin finite element method (SGFEM) for approximating the solution of parameterized elliptic partial differential equations with both affine and non-affine random coefficients. To compute the fully discrete solution, such approaches employ a Galerkin projection in both the deterministic and stochastic domains, produced here by a combination of finite elements and a global orthogonal basis, defined on an isotopic total degree index set, respectively. To account for the sparsity of the resulting system, we present a rigorous cost analysis that considers the total number of coupled finite element systems that must be simultaneously solved in the SGFEM. However, to maintain sparsity as the coefficient becomes increasingly nonlinear in the parameterization, it is necessary to also approximate the coefficient by an additional orthogonal expansion. In this case we prove a rigorous complexity estimate for the number of floating point operations (FLOPs) required per matrix-vector multiplication of the coupled system. Based on such complexity estimates we also develop explicit cost bounds in terms of FLOPs to solve the stochastic Galerkin (SG) systems to a prescribed tolerance, which are used to compare with the minimal complexity estimates of a stochastic collocation finite element method (SCFEM), shown in our previous work [16]. Finally, computational evidence complements the theoretical estimates and supports our conclusion that, in the case that the coefficient is affine, the coupled SG system can be solved more efficiently than the decoupled SC systems. However, as the coefficient becomes more nonlinear, it becomes prohibitively expensive to obtain an approximation with the SGFEM.


Full work available at URL: https://arxiv.org/abs/1507.05545




Recommendations




Cites Work


Cited In (2)

Uses Software





This page was built for publication: Explicit cost bounds of stochastic Galerkin approximations for parameterized PDEs with random coefficients

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2007265)