Working with sparse linear systems