diff --git a/dune/solvers/common/algorithm.hh b/dune/solvers/common/algorithm.hh
index ecec1416d531c875c75db839e72f0f6d08d452cb..95e1d0b802196c55d9b903ff9b60216bc78274fc 100644
--- a/dune/solvers/common/algorithm.hh
+++ b/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);