From 19660295784f28080e0eca3b0e003e36e71360e1 Mon Sep 17 00:00:00 2001
From: Oliver Sander <oliver.sander@tu-dresden.de>
Date: Sun, 22 Nov 2020 11:31:38 +0100
Subject: [PATCH] Fix bug in method lowerTriangularSolve
This makes triangularsolvetest.cc pass.
---
dune/matrix-vector/triangularsolve.hh | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/dune/matrix-vector/triangularsolve.hh b/dune/matrix-vector/triangularsolve.hh
index 4d7c42e..5d21f49 100644
--- a/dune/matrix-vector/triangularsolve.hh
+++ b/dune/matrix-vector/triangularsolve.hh
@@ -37,7 +37,7 @@ namespace MatrixVector {
// Note: We could drop the check for ignore nodes here bcs. b[j] will
// be ignored anyway due to the check above.
if (ignore == nullptr or (*ignore)[j].none())
- b[j] -= x[i] * *cIt;
+ cIt->mmv(x[i], b[j]);
}
}
} else {
--
GitLab