Add the sparseRange() utility function

This helper allows to use `for`-loops with
structured bindings for sparse ranges `R`
whose iterators `it` provide an
`it->index()` method, like this

```cpp
for(auto&& [A_i, i] : sparseRange(R))
  doSomethingWithValueAndIndex(A_i, i);
```