pith. sign in

arxiv: 1704.06339 · v2 · pith:55T2YS2Knew · submitted 2017-04-20 · 🧮 math.NA · cs.NA

A Monte Carlo approach to computing stiffness matrices arising in polynomial chaos approximations

classification 🧮 math.NA cs.NA
keywords matricescarlomonteequationsmethodpolynomialapproachapproximations
0
0 comments X
read the original abstract

We use a Monte Carlo method to assemble finite element matrices for polynomial Chaos approximations of elliptic equations with random coefficients. In this approach, all required expectations are approximated by a Monte Carlo method. The resulting methodology requires dealing with sparse block-diagonal matrices instead of block-full matrices. This leads to the solution of a coupled system of elliptic equations where the coupling is given by a Kronecker product matrix involving polynomial evaluation matrices. This generalizes the Classical Monte Carlo approximation and Collocation method for approximating functionals of solutions of these equations.

This paper has not been read by Pith yet.

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.