An efficient direct solver for the boundary concentrated FEM in 2D (Q1864158)

This paper analyzes the boundary concentrated finite element method which is well suited for the numerical approximation of elliptic boundary value problems with smooth coefficients and nonsmooth boundary conditions. This method is a variant of the \(hp\)-method. It uses very refined meshes in the neighbourhood of the boundary and coarse meshes away from the boundary. In the two-dimensional case, the authors obtain \(O(\log N)\) complexity estimates for the Cholesky factorization of the resulting stiffness matrix, where \(N\) denotes a size parameter of the problem. These theoretical estimates are nicely illustrated by various numerical experiments.