New shape functions for triangular \(p\)-FEM using integrated Jacobi polynomials (Q2494375)

The authors report a remarkable result concerning the computational cost in evaluating the element stiffness matrix and element mass matrix in the \(p\)-version of finite element method (FEM) applied to a linear second-order elliptic boundary value problem. In case of piecewise constant coefficient equations they construct a two dimensional basis functions such that the matrix vector multiplication and the generation of stiffness and mass matrices can be done in a number of arithmetic operations which equals the order of the number of elements in the finite element space.