diff --git a/dune/matrix-vector/triangularsolve.hh b/dune/matrix-vector/triangularsolve.hh
index 325a1ea70bf8c8d833fcd10099fc428f783527f1..baa2c52e7a13916e4db7a5ec42782901ff13ba72 100644
--- a/dune/matrix-vector/triangularsolve.hh
+++ b/dune/matrix-vector/triangularsolve.hh
@@ -16,6 +16,7 @@ namespace MatrixVector {
         auto cIt = it->begin();
         assert(cIt.index() == it.index());
         x[i] = b[i] / *cIt;
+        cIt++;
         for (; cIt != it->end(); ++cIt) {
           const size_t j = cIt.index();
           b[j] -= x[i] * *cIt;