Sparse matrix performance improvement -- NB: API unstable
Add benchmarks as example