Sparse matrix computation