Application of Randomized Quadrature Formulas to the Finite Element Method for Elliptic Equations


The implementation of the finite element method for linear elliptic equations requires to assemble the stiffness matrix and the load vector. In general, the entries of this matrix-vector system are not known explicitly but need to be approximated by quadrature rules. If the coefficient functions of the differential operator or the forcing term are irregular, then standard quadrature formulas, such as the barycentric quadrature rule, may not be reliable. In this paper we investigate the application of two randomized quadrature formulas to the finite element method for such elliptic boundary value problems with irregular coefficient functions. We give a detailed error analysis of these methods, discuss their implementation, and demonstrate their capabilities in several numerical experiments.

ArXiv Preprint
Raphael Kruse
Raphael Kruse

Prof. Dr. Raphael Kruse is the head of the working group “Numerik stochastischer Differentialgleichungen” at Martin-Luther-University Halle-Wittenberg. His research interests include numerical methods and stochastic analysis for stochastic evolution equations and Monte Carlo methods.