I am currently using SimplicialLDLT from Eigen (which I believe uses the Cholesky decomposition) to solve my large over-constrained matrix (30,000 by 1,000). It is taking aroun