Direct Solvers and Higher Order Methods for PDEs