Spezialize sparseRangeFor for MultiTypeBlockVector not *Matrix

dune/solvers/common/algorithm.hh

@@ -246,7 +246,7 @@ constexpr auto hybridSize(const T& t)
  * \brief Hybrid for loop over sparse range
  */
 template<class... T, class F>
-void sparseRangeFor(Dune::MultiTypeBlockMatrix<T...>& range, F&& f)
+void sparseRangeFor(const Dune::MultiTypeBlockVector<T...>& range, F&& f)
 {
   integralRangeFor<std::size_t>(Indices::_0, hybridSize(range), [&](auto&& i) {
       f(range[i], i);