Add MatrixIndexSet class

This is a drop in replacement for Dune::MatrixIndexSet from dune-istl. But it uses a sorted std::vector instead of a std::set internally for better performance. If you expect very dense rows, Dune::MatrixIndexSet may still perform better.

Using this in the operator assembler can improve performance measurably.