added csolve_sparse()

R uses cholmod for sparse factorization and solving which is fast, but Eigen's
SimplicialLLT Solver is magnitudes faster for what we need