From f4f369271964d227171e6ef045e6d78a3f986e0c Mon Sep 17 00:00:00 2001
From: Max Kahnt <max.kahnt@fu-berlin.de>
Date: Tue, 26 Jul 2016 13:39:45 +0200
Subject: [PATCH] Fix: Starting residual computation with wrong iterate.
---
dune/matrix-vector/triangularsolve.hh | 1 +
1 file changed, 1 insertion(+)
diff --git a/dune/matrix-vector/triangularsolve.hh b/dune/matrix-vector/triangularsolve.hh
index 325a1ea..baa2c52 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;
--
GitLab