Sparse matrix wrapper for Eigen
Since we ship Eigen anyway, it would be great if we had a sparse matrix wrapper for their linear solvers (thanks @he66coqe for suggesting this!). This would also yield access to sparse direct solvers without the need to install PETSc. Clearly, the Eigen implementations are not well-suited for (massively) parallel computations.