Boundary Concentrated FEM

A method with optimal (up to logarithmic terms) complexity for solving elliptic problems is proposed. The method relies on interior regularity but the solution may have globally weak regularity due to rough boundary data or geometries. Elliptic regularity results, high order approximation results, and an efficient preconditioner are presented.