Abstract
Algorithms are presented that enable the element matrices for the standard finite element space, consisting of continuous piecewise polynomials of degree $n$ on simplicial elements in $\mathbb{R}^d$, to be computed in optimal complexity $\mathcal{O}(n^{2d})$. The algorithms (i) take into account numerical quadrature; (ii) are applicable to nonlinear problems; and (iii) do not rely on precomputed arrays containing values of one-dimensional basis functions at quadrature points (although these can be used if desired). The elements are based on Bernstein polynomials and are the first to achieve optimal complexity for the standard finite element spaces on simplicial elements.
Original language | English |
---|---|
Pages (from-to) | 3087-3109 |
Number of pages | 23 |
Journal | SIAM Journal on Scientific Computing |
Volume | 33 |
Issue number | 6 |
Early online date | 1 Nov 2011 |
DOIs | |
Publication status | Published - 2011 |
Keywords
- finite elements
- optimal assembly
- polynomials