A finite element method for surface PDEs: matrix properties